Free Access
Issue
A&A
Volume 626, June 2019
Article Number L4
Number of page(s) 10
Section Letters to the Editor
DOI https://doi.org/10.1051/0004-6361/201935500
Published online 17 June 2019

© ESO 2019

1. Introduction

A stellar occultation occurs when a solar system object passes in front of a star from the point of view of an observer on Earth, causing a temporary drop in the observed flux of the star. This technique allows the determination of sizes and shapes with kilometre precision and to obtain characteristics of the object, such as its albedo, the presence of an atmosphere, rings, jets, or other structures around the body (Sicardy et al. 2011, 2016; Braga-Ribas et al. 2013, 2014; Dias-Oliveira et al. 2015; Benedetti-Rossi et al. 2016; Ortiz et al. 2015, 2017; Leiva et al. 2017; Bérard et al. 2017), or even the detection of topographic features (Dias-Oliveira et al. 2017).

Stellar occultations can also provide very accurate astrometric measurements of the occulted body, with uncertainties that can be as low as 5–10 km or even smaller for some objects (Leiva et al. 2017; Desmars et al. 2019). Compared with other methods in the context of the Galilean moons, classical CCD astrometry enables us to obtain positions with uncertainties in the 300–450 km level (Kiseleva et al. 2008), and relative positions between two close satellites achieve uncertainties of 100 km (Peng et al. 2012). Positions obtained using mutual phenomena have uncertainties at a level of 15–60 km (Saquet et al. 2018; Arlot et al. 2014; Dias-Oliveira et al. 2013; Emelyanov 2009), but they only occur at every equinox of the host planet (every six years in the case of Jupiter). The technique of mutual approximations also provides positions with uncertainties between 15–60 km, and this does not depend on the equinox of Jupiter (see Morgado et al. 2016, 2019 for details). Stellar occultation is then the only ground-based technique that can furnish astrometric measurements that are comparable with space probes, which usually have uncertainties smaller than 5 km (Tajeddine et al. 2015). In addition, the positions and sizes that can be obtained with stellar occultations are independent of reflectance models, which may have systematic errors due to variations on the satellite surface (Lindegren 1977).

In the context of the Galilean satellites, stellar occultations can provide shapes and sizes with uncertainties that are comparable with those of space probe images, but that are not affected by albedo variations on the satellite surface or by limb fitting, which is highly affected by the solar phase angle. From an astrometric point of view, these events can provide the best ground-based astrometry of these moons, with uncertainties smaller than 1 mas. This is at least one order of magnitude better than other methods.

Between 2019 and 2020, Jupiter will be in a very dense star region, the Galactic centre will be its background. This will only occur again in 2031. The probability of a stellar occultation by the Jovian moons increases dramatically (Gomes-Júnior et al. 2016). This is a great opportunity to observe stellar occultations by the Galilean moons, determine their positions, improve their ephemerides, and measure their shapes independently of probes. Accurate orbits help to prepare space missions targeting the Jovian system (Dirkx et al. 2016, 2017), and can help in the study of tides (Lainey et al. 2009; Lainey 2016). In the near future, we will have the ESA JUICE1 and the NASA Europa Clipper2 missions, which are scheduled for launch in the next decade (2020s). Moreover, only two stellar occultations by Ganymede have been observed in the past: one in 1972 (Carlson et al. 1973) and another in 2016 (D’Aversa et al. 2017). No occultations by Io, Europa, or Callisto are reported. Here we present the results obtained by the first observation of a stellar occultation by Europa, which occurred on 2017 March 31. These results serve as a proof of concept for the 2019–2021 campaign that is being organised.

The prediction method and the observational campaign are outlined in Sect. 2, while the data analysis and the occultation results are given in Sect. 3. The upcoming events for the period 2019–2021 can be found in Sect. 4. Our final remarks are contained in Sect. 5.

2. Prediction and observational campaign

The passage of Jupiter in a region that has the Galactic plane as background, as seen from Earth in the period between 2019–2020, creates a number of opportunities for observing stellar occultations by its satellites because this region has a high density of stars (Gomes-Júnior et al. 2016). This passage motivates the search for events involving the Galilean satellites because the satellites need stars with a magnitude of G = 11.5 or lower, so that the magnitude drop (Δmag) can be higher than 0.005 mag. This drop is otherwise very hard to observe with current equipment and techniques.

Following similar procedures as described in Assafin et al. (2010, 2012), we predicted that Europa would occult a star (mag. G = 9.5) on 2017 March 31 at 06:44 UTC. The shadow path would be crossing South America with a velocity of 17.78 km s−1. The prediction was made using the Europa ephemerides jup3103 with DE435 furnished by the Jet Propulsion Laboratory (JPL), and the stellar position was obtained from the Gaia DR1 catalogue (Clementini et al. 2016) because the occultation occurred before the launch of Gaia DR2 in April 2018 (Lindegren et al. 2018). The star was then identified in Gaia DR2. Its position (right ascension α, declination –δ, and errors) was updated using its proper motion and parallax for the occultation time (2017 March 17 06:44:00 UTC), and its G mag,

α star = 13 h 12 15 . 5430 ± 0.16 mas , δ star = 05 ° 56 48 . 7526 ± 0.12 mas , G mag = 9.5065 . $$ \begin{aligned}&\alpha _{\rm star} = 13^\mathrm{h}12\prime 15.\!\!^{\prime \prime }5430 \pm 0.16\,\mathrm{mas}, \nonumber \\&\delta _{\rm star} = -05^{\circ }56\prime 48.\!\!^{\prime \prime }7526 \pm 0.12\,\mathrm{mas}, \nonumber \\&G\,\mathrm{mag} = 9.5065. \end{aligned} $$(1)

This occultation was favourable for several potential observers in several countries in South America (Fig. 1). We then organised a campaign to observe this event, and a total of nine stations tried to observe the occultation. They are listed in Table 1. However, due to bad weather conditions, only four were able to obtain data (OPD/B&C, OPD/ROB, FOZ, and OBSPA)4. The observation made at Itajubá with the 40 cm telescope was discarded because the signal-to-noise ratio (S/N) was very low, while the other three (OPD/B&C, FOZ, and OBSPA) obtained positive detection.

thumbnail Fig. 1.

Occultation by Europa on 2017 March 31. The centre of the body at the time of closest approach between the shadow and the geocentre (CA), at 06:44:36 UTC, is represented by the large red dot. Blue lines represent the size limit of Europa (1561.2 ± 3.6 km). The small red dots represent the centre of the body, each separated by 1 min from the reference instant CA. The blue dots are the observers with positive detection, and the yellow dots are observers that were unable to obtain the light curve due to bad weather conditions. The offset between our result and the prediction (using the JPL jup310 ephemeris) was −0.90 ± 3.1 km for fc and −12.3 ± 2.0 km for gc, see Sect. 3. The black arrow at the bottom represents the direction of motion of the shadow.

Table 1.

Observational stations, technical details, and circumstance.

Figure 1 presents the post-occultation map. It shows the shadow path of Europa through South America, and the positions of the observation stations. The blue lines represent the size of Europa obtained after the fit (see Sect. 3).

3. Data analysis and results

Images from all data sets were corrected for bias and flat-field using standard IRAF5 procedures. A stacking image procedure6 developed by us in PYTHON (Astropy Collaboration 2013), following similar procedures as Zhu et al. (2018), was used in FOZ images in order to obtain an S/N that was high enough to emphasise the magnitude drop. The low time-resolution at OPD/ROB (cycle ≈ 7 s) did not allow the use of this technique, and it was not possible to determine the magnitude drop.

Differential aperture photometry was applied with the PRAIA package (Assafin et al. 2011). The calibrator and target apertures were automatically chosen to maximise the S/N for each frame. Io (J1) was set as calibrator in all sets of images because there were no other star in the field of view with a good S/N that could be used for calibration. The light fluxes of Europa and of the occulted star were measured together in the same aperture. The normalised light curves are displayed in Fig. 2, where the black lines are the observations and the red lines are the modelled fitted light curves.

thumbnail Fig. 2.

Three normalised light curves obtained with positive detection. They are shifted vertically for better viewing. Observational circumstances are given in Table 1. The difference in the depth of the curves is due to the use of different filters.

The ingress (ti) and egress (te) times were determined with a standard χ2 procedure between the observational light curve and the model. The model considers a sharp-edge occultation model convolved with Fresnel diffraction, stellar diameter (projected at the body distance, 0.633 km in our case Bourges et al. 2017), CCD bandwidth, and finite integration time (more details in Braga-Ribas et al. 2013). Table 2 contains the ingress and egress UTC times and the errors for each station with positive detection. Their respective uncertainties are given in seconds and kilometres (calculated using the event velocity of 17.78 km s−1). More information about the light-curve fit procedure is presented in Appendix B.

Table 2.

Occultation ingress (ti) and egress (te) times.

Each of the ingress and egress times were converted into a star position (f, g) regarding the body centre, with f and g being measured positively toward local celestial east and celestial north, respectively. Each pair of positions, from the same site, is a chord, and each position is a point of which we can fit the five parameters that defines an ellipse: (i and ii) the ellipse centre (fc, gc), (iii) apparent semi-major axis (a′), (iv) the apparent oblateness (ϵ′ = (a′−b′)/a′, where b′ is the apparent semi-minor axis), and (v) the position angle of the pole Pp of b′. This position angle was fixed as 24.2534°, derived from the pole position reported by Archinal et al. (2018). In Fig. 3 we show the three chords obtained from the occultation and the fitted ellipse. The ellipse parameters are

f c = 0.9 ± 3.1 km , g c = 12.3 ± 2.0 km , a = 1562.0 ± 3.6 km , ϵ = 0.0010 ± 0.0028 , P p = 24 . 2534 ° . $$ \begin{aligned} f_{\rm c}&= -0.9 \pm 3.1\,\mathrm{km}, \nonumber \\ {g}_{\rm c}&= -12.3 \pm 2.0\,\mathrm{km}, \nonumber \\ a\prime&= 1562.0 \pm 3.6\,\mathrm{km}, \nonumber \\ \epsilon \prime&= 0.0010 \pm 0.0028, \nonumber \\ P_{\rm p}&= 24.2534^{\circ }. \end{aligned} $$(2)

thumbnail Fig. 3.

Best elliptical limb fit (in black) using the three occultation chords. Blue, cyan, and green lines are the chords from the positive detection obtained at OPD/B&C, FOZ, and OBSPA, respectively. The small red segments at the end of each segment indicate the 1σ uncertainty on each chord extremity, derived from the time uncertainties provided in Table 2. The position of the centre of the fitted ellipse is showed as an X, with values of fc = −0.9 ± 3.1 km and gc = −12.3 ± 2.0 km, and represents the offset with respect to the JPL geocentric ephemeris. The fitted parameters of the ellipse are listed in Table 2.

It is important to highlight that the centre position (fc, gc) and the apparent size and shape of the ellipse (a′,ϵ′) are correlated. This correlation decreases with the number of obtained chords, and moreover, with the distribution of these chords in the shadow. Chords in both hemispheres are required to achieve the best results. In our case, the higher correlation was 0.81 between gc and ϵ′, which is expected because we only had chords in the southern hemisphere of the shadow.

We can compare our results with the values determined by Galileo mission images (Nimmo et al. 2007). Our apparent semi-major axis (a′ = 1562.0 ± 3.6 km) was between the axis a (1562.6 km) and b (1560.3 km) of the fitted ellipsoid. Our apparent semi-minor axis (b′ = 1560.4 ± 5.7 km) is nearly the c axis of the ellipsoid (1559.5 km). The expected oblateness was between 0.0019 and 0.0005, in comparison with 0.0010 ± 0.0028 that we obtained. The equivalent radius ( R eq = a 1 ϵ = 1561.2 ± 3.6 $ R_{\mathrm{eq}}=a{\prime}\sqrt{1 - \epsilon{\prime}} = 1561.2 \pm 3.6 $ km) agrees with the nominal equivalent radius of 1560.8 ± 0.3 km of Europa. According to Kay et al. (2019), the topographic features of Europa are at a level of some hundreds of metres, therefore it is not possible to infer any topography from our observations.

With the limb fitting, we also obtained a geocentric position for the Europa centre on 2017 March 31 at 06:44:00 UTC of

α J 2 = 13 h 12 15 . 548372 ± 0.96 mas , δ J 2 = 05 ° 56 48 . 687034 ± 0.62 mas . $$ \begin{aligned} \alpha _{J2}&= 13^\mathrm{h}12\prime 15.\!\!^{\prime \prime }548372 \pm 0.96\,\mathrm{mas}, \nonumber \\ \delta _{J2}&= -05^{\circ }56\prime 48.\!\!^{\prime \prime }687034 \pm 0.62\,\mathrm{mas}. \end{aligned} $$(3)

This position has an offset with respect to the JPL jup310 ephemerides of +0.28 mas and +3.81 mas for Δα cos δ and Δδ, respectively, while the offset with respect to the IMCCE NOE-5-2010-GAL (Lainey et al. 2009) is −2.59 mas and −6.01 mas, respectively, both in the sense of observation minus ephemeris. More details can be found in Appendix C.

4. Future events

Observing very many occultations by the Galilean satellites is required for better understanding their 3D shape (oblatness and pole positions). This in turn enables obtaining a highly accurate absolute position (with the help of Gaia DR2 catalogue, which furnishes the star positions) that are to be used in dynamical models and with which high-quality ephemerides can be obtained.

We conducted a search for occultation events that would occur between 2017 and 2021, using the GDR2 catalogue and the JPL jup310 and DE435 ephemerides. In Table 3 we summarise the occultation predictions between 2019 and 2021. The occultation presented here is highlighted. We also list some parameters that are necessary to produce occultation maps such as the closest approach instant (UTC) of the prediction; C/A, the apparent geocentric distance between the satellite and the star at the moment of the geocentric closest approach, in arcseconds; P/A, the satellite position angle regarding the star to be occulted at C/A, in degrees (zero at north of the star, increasing clockwise); magnitude G of the occulted GDR2 star; and finally, the expected magnitude drop (Δmag) for each event. For more information about the definition and use of these stellar occultation geometric elements, see Assafin et al. (2010, 2012). We present in Table 3 three events by Io, four by Europa, one by Ganymede, and three by Callisto with Δmag higher than 0.005 mag. The prediction maps are shown in the Appendix D.

Table 3.

Predicted stellar occultations by the Galilean moons between 2019 and 2021.

We draw the attention to the occultation by Io of a mag. G 5.8 star on 2012 April 2. The magnitude drop will be about 0.740 mag. This event can be observed in the South of North America, Central America, and in the North of South America (Fig. 4).

thumbnail Fig. 4.

Prediction of the occultation by Io on 02 April 2021 10:24 UTC, using JPL DE435 and jup310 ephemeris. The blue lines represent the expected size limit of Europa. The blue dashed lines are an uncertainty of 20 mas in Europa’s position. The red dots represents the centre of the body for a given time, each separated by 1 min. The arrow represents the sense of the shadow’s motion. The dark grey and the light grey area represent the night and the day part of the globe, respectively.

This type of observation is challenging because of the low magnitude drop. We suggest the use of red filters (or a narrow-band methane filter, centred on 889 nm with a width of 15 nm) in these observations. Jupiter should be placed outside the field of view of the camera when possible, which will reduce noise effects caused by scattered light from Jupiter. It is very important to test the equipment configuration (gain, exposure time, etc.) some days before the occultation in order to ensure the best S/N possible.

We also suggest the use of fast cameras with high-cadence images and minimum readout time. Stacking images techniques, such as the one applied here (detailed in Appendix A), can be used to increase the S/N of the images. Telescopes as small as 20 cm (8″) can provide useful data. Adaptive optics is also a possibility, with adequate exposure time. For a telescope with a large aperture and a sensitive camera, this could result in the direct detection of a thin atmosphere around a Galilean satellite. Observations with resolutions higher than the Fresnel scale (about 0.44 km at Jupiter) can result in the direct measurement of refraction spikes by a thin atmosphere such as the spikes observed in occultations by Pluto (Dias-Oliveira et al. 2015; Sicardy et al. 2016; Meza et al. 2019).

5. Discussion

Stellar occultation is a ground-based technique capable of determining positions, sizes, and shapes with uncertainties of some kilometres. These values are comparable with data obtained with space probes. In this work we took advantage of the beginning of the Jupiter passage through a very dense stellar region, with the Galactic centre as background, and organised an observational campaign to obtain data from a predicted stellar occultation by the Galilean moon Europa. Its ephemeris indicates that it would pass in front of a mag. G = 9.5 GDR2 star on 2017 March 31. Out of the nine stations in eight sites across South America, three obtained a positive detection of the event.

The fitted ellipse for Europa gives an area equivalent radius of 1561.2 ± 3.6 km. Its absolute position uncertainty is 0.80 mas (2.55 km), representing an offset with respect to the JPL jup310 ephemeris of 3.82 mas and of 6.54 mas with respect to the NOE-5-2010-GAL ephemeris. This position for Europa at the time of the occultation took advantage of the high-precision position of Gaia DR2 stars and the fact that stellar occultations render accurate body-star relative positions.

This is the first reported observation of a stellar occultation by Europa. The favourable configuration of Jupiter, which has the Galactic plane as background, increases the chances of observing other bright star occultations by its main satellites, as discussed in Sect. 4. This configuration will only occur again in 2031. For the future events predicted in the 2019 and 2021, campaigns will be organised in due time.

We encourage the observational community to observe the predicted future stellar occultations. Even the amateur community with telescopes as small as 20 cm (8″) can help provide useful data. A stellar occultation provides an apparent size and shape at a specific moment. Combined with the observations of the upcoming events, a 3D shape of these moons can be acquired with kilometre precision and complement space mission data. This information can aid in the study of planetary formation, evolution, and the influence of tides raised by the moons’ primary in their orbit. This highlights the importance of these events.

Finally, all this information can also be used together with dynamical models to ensure highly accurate orbits for these moons. These orbits can be helpful for future space probes aimed at the Jovian system, such as JUICE and the Europa Clipper mission.


4

Alias, defined in Table 1.

6

More details in Appendix A.

Acknowledgments

We thank our anonymous referee for helpful comments. This study was financed by the Coordenação de Aperfeiçoamento de Pessoal de Nível Superior – Brasil (CAPES) – Finance Code 001. Part of this research is suported by INCT do e-Universo, Brazil (CNPQ grants 465376/2014-2). Based in part on observations made at the Laboratório Nacional de Astrofísica (LNA), Itajubá-MG, Brazil. BM thanks the CAPES/Cofecub-394/2016-05 grant. G.B.R. acknowledges the support of the CAPES and FAPERJ/PAPDRJ (E26/203.173/2016) grants. ARGJ thanks FAPESP proc. 2018/11239-8. MA thanks CNPq (Grants 427700/2018-3, 310683/2017-3 and 473002/2013-2) and FAPERJ (Grant E-26/111.488/2013). VL’s research was supported by an appointment to the NASA Postdoctoral Program at the NASA Jet Propulsion Laboratory, California Institute of Technology, administered by Universities Space Research Association under contract with NASA. RVM acknowledges the grants CNPq-306885/2013, CAPES/Cofecub-2506/2015, FAPERJ/PAPDRJ-45/2013 and FAPERJ/CNE/05-2015. JIBC acknowledges CNPq grant 308150/2016-3. FBR acknowledges CNPq support, proc. 309578/2017-5. RS and OCW acknowledge Fapesp proc. 2016/24561-0, CNPq proc. 312813/2013-9 and 305737/2015-5.

References

  1. Anderson, J. D., Schubert, G., Jacobson, R. A., et al. 1998, Science, 281, 2019 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  2. Archinal, B. A., Acton, C. H., A’Hearn, M. F., et al. 2018, Celestial Mech. Dyn. Astron., 130, 22 [Google Scholar]
  3. Arlot, J. E., Emelyanov, N., Varfolomeev, M. I., et al. 2014, A&A, 572, A120 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  4. Assafin, M., Camargo, J. I. B., Vieira Martins, R., et al. 2010, A&A, 515, A32 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  5. Assafin, M., Vieira Martins, R., Camargo, J. I. B., et al. 2011, Gaia Follow-up Network for the Solar System Objects: Gaia FUN-SSO Workshop Proc., 85 [Google Scholar]
  6. Assafin, M., Camargo, J. I. B., Vieira Martins, R., et al. 2012, A&A, 541, A142 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  7. Astropy Collaboration (Robitaille, T. P., et al.) 2013, A&A, 558, A33 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  8. Benedetti-Rossi, G., Sicardy, B., Buie, M. W., et al. 2016, AJ, 152, 156 [NASA ADS] [CrossRef] [Google Scholar]
  9. Bérard, D., Sicardy, B., Camargo, J. I. B., et al. 2017, AJ, 154, 144 [NASA ADS] [CrossRef] [Google Scholar]
  10. Bourges, L., Mella, G., Lafrasse, S., et al. 2017, VizieR Online Data Catalog: II/346 [Google Scholar]
  11. Braga-Ribas, F., Sicardy, B., Ortiz, J. L., et al. 2013, ApJ, 773, 26 [NASA ADS] [CrossRef] [Google Scholar]
  12. Braga-Ribas, F., Sicardy, B., Ortiz, J. L., et al. 2014, Nature, 508, 72 [NASA ADS] [CrossRef] [Google Scholar]
  13. Carlson, R. W., Bhattacharyya, J. C., Smith, B. A., et al. 1973, Science, 182, 53 [NASA ADS] [CrossRef] [Google Scholar]
  14. Chow, G. C. 1960, Econometrica, 28, 591 [CrossRef] [Google Scholar]
  15. Clementini, G., Ripepi, V., Leccia, S., et al. 2016, A&A, 595, A133 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  16. D’Aversa, E., Oliva, F., Sindoni, G., et al. 2017, Eur. Planet. Sci. Congr., 850 [Google Scholar]
  17. Desmars, J., Meza, E., Sicardy, B., et al. 2019, A&A, 625, A43 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  18. Dias-Oliveira, A., Vieira-Martins, R., Assafin, M., et al. 2013, MNRAS, 432, 225 [NASA ADS] [CrossRef] [Google Scholar]
  19. Dias-Oliveira, A., Sicardy, B., Lellouch, E., et al. 2015, ApJ, 811, 53 [NASA ADS] [CrossRef] [Google Scholar]
  20. Dias-Oliveira, A., Sicardy, B., Ortiz, J. L., et al. 2017, AJ, 154, 22 [NASA ADS] [CrossRef] [Google Scholar]
  21. Dirkx, D., Lainey, V., Gurvits, L. I., & Visser, P. N. A. M. 2016, Planet. Space Sci., 134, 82 [NASA ADS] [CrossRef] [Google Scholar]
  22. Dirkx, D., Gurvits, L. I., Lainey, V., et al. 2017, Planet. Space Sci., 147, 14 [NASA ADS] [CrossRef] [Google Scholar]
  23. Emelyanov, N. V. 2009, MNRAS, 394, 1037 [NASA ADS] [CrossRef] [Google Scholar]
  24. Gomes-Júnior, A. R., Assafin, M., Beauvalet, L., et al. 2016, MNRAS, 462, 1351 [NASA ADS] [CrossRef] [Google Scholar]
  25. Jia, X., Kivelson, M. G., Khurana, K. K., & Kurth, W. S. 2018, Nat. Astron., 2, 459 [NASA ADS] [CrossRef] [Google Scholar]
  26. Kay, J., Schenk, P., & Prockter, L. 2019, Lunar Planet. Sci. Conf., 2975 [Google Scholar]
  27. Kiseleva, T. P., Kiselev, A. A., Kalinichenko, O. A., Vasilyeva, N. A., & Khovricheva, M. L. 2008, Sol. Syst. Res., 42, 414 [NASA ADS] [CrossRef] [Google Scholar]
  28. Lainey, V. 2016, Celestial Mech. Dyn. Astron., 126, 145 [NASA ADS] [CrossRef] [Google Scholar]
  29. Lainey, V., Arlot, J.-E., Karatekin, Ö., & van Hoolst, T. 2009, Nature, 459, 957 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  30. Leiva, R., Sicardy, B., Camargo, J. I. B., et al. 2017, AJ, 154, 159 [NASA ADS] [CrossRef] [Google Scholar]
  31. Lindegren, L. 1977, A&A, 57, 55 [NASA ADS] [Google Scholar]
  32. Lindegren, L., Hernández, J., Bombrun, A., et al. 2018, A&A, 616, A2 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  33. Meza, E., Sicardy, B., Assafin, M., et al. 2019, A&A, 625, A42 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  34. Morgado, B., Assafin, M., Vieira-Martins, R., et al. 2016, MNRAS, 460, 4086 [NASA ADS] [CrossRef] [Google Scholar]
  35. Morgado, B., Vieira-Martins, R., Assafin, M., et al. 2019, MNRAS, 482, 5190 [NASA ADS] [CrossRef] [Google Scholar]
  36. Nimmo, F., Thomas, P. C., Pappalardo, R. T., & Moore, W. B. 2007, Icarus, 191, 183 [NASA ADS] [CrossRef] [Google Scholar]
  37. Ortiz, J. L., Duffard, R., Pinilla-Alonso, N., et al. 2015, A&A, 576, A18 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  38. Ortiz, J. L., Santos-Sanz, P., Sicardy, B., et al. 2017, Nature, 550, 219 [NASA ADS] [CrossRef] [Google Scholar]
  39. Peng, Q. Y., He, H. F., Lainey, V., & Vienne, A. 2012, MNRAS, 419, 1977 [NASA ADS] [CrossRef] [Google Scholar]
  40. Roth, L., Saur, J., Retherford, K. D., et al. 2014, Science, 343, 171 [NASA ADS] [CrossRef] [Google Scholar]
  41. Saquet, E., Emelyanov, N., Robert, V., et al. 2018, MNRAS, 474, 4730 [NASA ADS] [CrossRef] [Google Scholar]
  42. Seber, G. A. F., & Lee, A. J. 2003, Linear Regression Analysis (Wiley) [CrossRef] [Google Scholar]
  43. Sicardy, B., Ortiz, J. L., Assafin, M., et al. 2011, Nature, 478, 493 [NASA ADS] [CrossRef] [Google Scholar]
  44. Sicardy, B., Talbot, J., Meza, E., et al. 2016, ApJ, 819, L38 [NASA ADS] [CrossRef] [Google Scholar]
  45. Sparks, W. B., Hand, K. P., McGrath, M. A., et al. 2016, ApJ, 829, 121 [NASA ADS] [CrossRef] [Google Scholar]
  46. Sparks, W. B., Richter, M., deWitt, C., et al. 2019, ApJ, 871, L5 [NASA ADS] [CrossRef] [Google Scholar]
  47. Tajeddine, R., Lainey, V., Cooper, N. J., & Murray, C. D. 2015, A&A, 575, A73 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  48. Zhu, B., Zhao, Y., Ge, J., & Ma, J. 2018, MNRAS, 481, 1469 [NASA ADS] [CrossRef] [Google Scholar]

Appendix A: Stacking image procedure

thumbnail Fig. A.1.

Normalised light flux from FOZ images with (black line) and without (light gray line) the stacking image procedure. The red line is the modelled light curve. Each stacked image contains the data of 20 single images stacked together.

The stacking image procedure was developed with the aim to attenuate the image noise, thereby increasing the image S/N when is not possible to just integrate over more time. This is the case when the Galilean moons are observed without a filter, as the brightness of Jupiter or even the satellite itself would quickly saturate the CCD image before we achieve an adequate S/N.

The first step in this procedure is to chose a well-sampled object in the images that is to be used as reference to align the images for the stacking. This object should not have a significant motion (in the sky plane) between the first and last image. In our case, in the absence of stars, we used Europa to align the images. Each individual image had an exposure of 0.05 seconds and the same cycle, and the motion of Europa was 5.39 mas s−1 in the sky plane, corresponding to a motion of 0.0074 pixels s−1 on the CCD.

We measured the object centroid (x, y) with a 2D circular symmetric Gaussian fit over pixels within one full-width at half-maximum (FWHM ∝ seeing) from the centre. This was done using the PRAIA package (Assafin et al. 2011). The uncertainty in the centroid measurement was about 36 mas (0.05 pixels). The alignment consists of vertical and horizontal shifts for each image (Δxy) relative to a chosen reference image; in our case, the first image.

We performed a stack procedure over each of the 20 images, which resulted in a temporal resolution of 1.0 s. Figure A.1 contains the light curve we obtained using all individual frames (light gray line), the light curve using the stacked images (black line), and the modelled light curve after the fit (red line). Before our procedure, the normalised light curve rms was 0.073, while after the stacking, we were able to achieve 0.016. It is important to highlight that the expected drop would be about 0.021 in flux. Without the stacking image procedure, it would be impossible to determine the ingress and/or the egress times.

Appendix B: Light-curve fitting procedure

thumbnail Fig. B.1.

Normalised light curve obtained at OPD/B&C with a zoom-in at the moment of the occultation. The black dots show the data, while the red lines represent the sharp-edge model for each observation. The dotted blue line shows the trapezoidal model.

The most common geological model of Europa is an ice surface on top of a global ocean (Anderson et al. 1998). The presence of plumes of water vapour (Roth et al. 2014; Jia et al. 2018; Sparks et al. 2016, 2019) is a strong evidence of this model.

The OPD/B&C light curve presents an apparent smooth flux decay on the ingress and egress times during the passage of the body in front of the star, and this slope would be compatible with the common picture of the Europa geology and water vapour evidence. These features must be analysed with caution, however, because the measurements show some level of correlated noise. To objectively assess the significance of this slope, we tested two scenarios: one simple (sharp-edge model) scenario in which we assumed that there are no slopes during the ingress and egress times, and one more complex scenario (trapezoidal model) with a smooth ingress and egress. It is also straightforward to show that the sharp-edge model is a particular case of a trapezoidal model, allowing us to analyse the decrease in residuals when the number of fitted parameters is increased.

thumbnail Fig. B.2.

Same as Fig. B.1 for the light curve obtained at FOZ.

thumbnail Fig. B.3.

Same as Fig. B.1 for the light curve obtained at OBSPA.

Because the more complex model has fewer degrees of freedom, it should at least fit the observables as well as the more straightforward case. Having set the characteristics of the problem, we used the Fisher-Snedecor F-test (Chow 1960; Seber & Lee 2003) to estimate the probability that the trapezoidal model fits the data better than a more simple sharp-edge model. Under the null hypothesis, the trapezoidal case does not produce a better fit than the standard case, or it could be argued that the slope coefficients are similar. The test follows an F-distribution, and the null hypothesis is rejected when the calculated F value is higher than the expected critical value (or significance).

As an example, we discuss the light curve obtained by OPD/B&C. The normalised chi-square statistics obtained from the sharp edge model was 1.06, whereas the value obtained with the trapezoidal model was 1.04. These models were applied to the data within an interval of five minutes centred at the central instant for each station. Using the F distribution, we can predict that within a 95% confidence level, the computed chi-square value for the plume model should lie around 0.74. In this case, increasing the number of model parameters leads only to a better fit of the noise because we found out that the probability of the trapezoidal model to explain the data better is about 17%. Our measurements therefore do not support any further information about the detection of a thin atmosphere (or a plume). The FOZ and OBSPA light curves have lower cadences and therefore a lower time resolution during ingress and egress of the event. The probabilities for these sites are even lower, about 1% and 3% respectively, as expected.

Figures B.1, B.2, and B.3 contain the individual light curves for each observational site (OPD/B&C, FOZ, and OBSPA, respectively) with a zoom (bottom graphs) of two minutes centred on the ingress and egress times. The x-axes are the time in minutes relative to 06:44:00 UTC, and the y-axes are the normalised light flux. The red lines are the sharp-edge model for each observation, and the dotted blue lines are the trapezoidal model.

Appendix C: Comparison between ephemeris

Table C.1.

Offset between the geocentric position of Europa obtained from the stellar occultation on 2017 March 31 at 06:44 UTC and different ephemerides in the sense observation minus ephemeris.

The absolute position of Europa obtained here has an uncertainty smaller than 1 mas (∼3 km). This level of accuracy is of the same order of magnitude as those obtained from space mission images (Tajeddine et al. 2015). However, our position is completely independent of any variation in the satellite albedo and is not affected by phase angle effects. Our geocentric position for Europa is given in Eq. (C.1),

α J 2 = 13 h 12 15 . 548372 ± 0.96 mas , δ J 2 = 05 ° 56 48 . 687034 ± 0.62 mas . $$ \begin{aligned} \alpha _{J2}&= 13^\mathrm{h}12\prime 15.\!\!^{\prime \prime }548372 \pm 0.96\,\mathrm{mas}, \nonumber \\ \delta _{J2}&= -05^{\circ }56\prime 48.\!\!^{\prime \prime }687034 \pm 0.62\,\mathrm{mas}. \end{aligned} $$(C.1)

Here we compare this position (αJ2, δJ2) with different geocentric ephemeris (αephem, δephem) using Eq. (C.2),

Δ α = α J 2 α ephem , Δ δ = δ J 2 δ ephem , s Δ α 2 cos 2 ( δ J 2 + δ ephem 2 ) + Δ δ 2 . $$ \begin{aligned} \Delta \alpha&= \alpha _{J2} - \alpha _{\rm ephem}, \nonumber \\ \Delta \delta&= \delta _{J2} - \delta _{\rm ephem}, \nonumber \\ s \simeq&\sqrt{\Delta \alpha ^2\cos ^2 \left(\dfrac{\delta _{J2} + \delta _{\rm ephem}}{2}\right) + \Delta \delta ^2} . \end{aligned} $$(C.2)

In Table C.1 we present the offset between the Europa position as derived based on this stellar occultation and the ephemeris. We compare our result with 12 ephemerides combinations. We use 6 planetary ephemerides (from the JPL, DE438s, DE436s, DE435, and DE430; and from the IMCCE, INPOP17a and INPOP13c) and two satellite ephemerides (from the JPL, jup310; and from the IMCCE, NOE-5-2010-GAL).

Appendix D: Future event maps

In this appendix we present the stellar occultation maps for the predicted events that will occur between 2019 and 2021 (see Table 3 for more details). We indicate in each figure the satellite that will occult the star, and the date and time of the occultation in UTC. The blue lines represent the expected size limit of the occulting satellite. The red dots are the centre of the body for a given time, each separated by 1 min. The large red dot represents the centre of the body at the time of closest approach (CA) indicated in the figure label. The arrow represents the direction of motion of the shadow. The dark grey and light grey areas represent the night and day part of the globe, respectively.

thumbnail Fig. D.1.

Occultation of a mag. G 10.9 star by Europa on 2019 May 5, 20:32 UTC. The predicted relative velocity of the event is 23.1 km s−1.

thumbnail Fig. D.2.

Occultation of a mag. G 9.1 star by Europa on 2019 June 4, 02:26 UTC. The predicted relative velocity of the event is 27.3 km s−1.

thumbnail Fig. D.3.

Occultation of a mag. G 10.2 star by Callisto on 2019 June 5, 23:12 UTC. The predicted relative velocity of the event is 21.0 km s−1.

thumbnail Fig. D.4.

Occultation of a mag. G 11.0 star by Io on 2019 September 9, 02:33 UTC. The predicted relative velocity of the event is 18.1 km s−1.

thumbnail Fig. D.5.

Occultation of a mag. G 11.3 star by Io on 2019 September 21, 13:08 UTC. The predicted relative velocity of the event is 26.5 km s−1.

thumbnail Fig. D.6.

Occultation of a mag. G 10.9 star by Callisto on 2020 June 20, 16:03 UTC. The predicted relative velocity of the event is 18.9 km s−1.

thumbnail Fig. D.7.

Occultation of a mag. G 11.3 star by Europa on 2020 June 22, 02:09 UTC. The predicted relative velocity of the event is 26.1 km s−1.

thumbnail Fig. D.8.

Occultation of a mag. G 5.8 star by Io on 2021 April 2, 10:24 UTC. The predicted relative velocity of the event is 16.5 km s−1.

thumbnail Fig. D.9.

Occultation of a mag. G 11.1 star by Ganymede on 2021 April 25, 07:55 UTC. The predicted relative velocity of the event is 27.7 km s−1.

thumbnail Fig. D.10.

Occultation of a mag. G 10.4 star by Callisto on 2021 May 4, 23:01 UTC. The predicted relative velocity of the event is 16.3 km s−1.

All Tables

Table 1.

Observational stations, technical details, and circumstance.

Table 2.

Occultation ingress (ti) and egress (te) times.

Table 3.

Predicted stellar occultations by the Galilean moons between 2019 and 2021.

Table C.1.

Offset between the geocentric position of Europa obtained from the stellar occultation on 2017 March 31 at 06:44 UTC and different ephemerides in the sense observation minus ephemeris.

All Figures

thumbnail Fig. 1.

Occultation by Europa on 2017 March 31. The centre of the body at the time of closest approach between the shadow and the geocentre (CA), at 06:44:36 UTC, is represented by the large red dot. Blue lines represent the size limit of Europa (1561.2 ± 3.6 km). The small red dots represent the centre of the body, each separated by 1 min from the reference instant CA. The blue dots are the observers with positive detection, and the yellow dots are observers that were unable to obtain the light curve due to bad weather conditions. The offset between our result and the prediction (using the JPL jup310 ephemeris) was −0.90 ± 3.1 km for fc and −12.3 ± 2.0 km for gc, see Sect. 3. The black arrow at the bottom represents the direction of motion of the shadow.

In the text
thumbnail Fig. 2.

Three normalised light curves obtained with positive detection. They are shifted vertically for better viewing. Observational circumstances are given in Table 1. The difference in the depth of the curves is due to the use of different filters.

In the text
thumbnail Fig. 3.

Best elliptical limb fit (in black) using the three occultation chords. Blue, cyan, and green lines are the chords from the positive detection obtained at OPD/B&C, FOZ, and OBSPA, respectively. The small red segments at the end of each segment indicate the 1σ uncertainty on each chord extremity, derived from the time uncertainties provided in Table 2. The position of the centre of the fitted ellipse is showed as an X, with values of fc = −0.9 ± 3.1 km and gc = −12.3 ± 2.0 km, and represents the offset with respect to the JPL geocentric ephemeris. The fitted parameters of the ellipse are listed in Table 2.

In the text
thumbnail Fig. 4.

Prediction of the occultation by Io on 02 April 2021 10:24 UTC, using JPL DE435 and jup310 ephemeris. The blue lines represent the expected size limit of Europa. The blue dashed lines are an uncertainty of 20 mas in Europa’s position. The red dots represents the centre of the body for a given time, each separated by 1 min. The arrow represents the sense of the shadow’s motion. The dark grey and the light grey area represent the night and the day part of the globe, respectively.

In the text
thumbnail Fig. A.1.

Normalised light flux from FOZ images with (black line) and without (light gray line) the stacking image procedure. The red line is the modelled light curve. Each stacked image contains the data of 20 single images stacked together.

In the text
thumbnail Fig. B.1.

Normalised light curve obtained at OPD/B&C with a zoom-in at the moment of the occultation. The black dots show the data, while the red lines represent the sharp-edge model for each observation. The dotted blue line shows the trapezoidal model.

In the text
thumbnail Fig. B.2.

Same as Fig. B.1 for the light curve obtained at FOZ.

In the text
thumbnail Fig. B.3.

Same as Fig. B.1 for the light curve obtained at OBSPA.

In the text
thumbnail Fig. D.1.

Occultation of a mag. G 10.9 star by Europa on 2019 May 5, 20:32 UTC. The predicted relative velocity of the event is 23.1 km s−1.

In the text
thumbnail Fig. D.2.

Occultation of a mag. G 9.1 star by Europa on 2019 June 4, 02:26 UTC. The predicted relative velocity of the event is 27.3 km s−1.

In the text
thumbnail Fig. D.3.

Occultation of a mag. G 10.2 star by Callisto on 2019 June 5, 23:12 UTC. The predicted relative velocity of the event is 21.0 km s−1.

In the text
thumbnail Fig. D.4.

Occultation of a mag. G 11.0 star by Io on 2019 September 9, 02:33 UTC. The predicted relative velocity of the event is 18.1 km s−1.

In the text
thumbnail Fig. D.5.

Occultation of a mag. G 11.3 star by Io on 2019 September 21, 13:08 UTC. The predicted relative velocity of the event is 26.5 km s−1.

In the text
thumbnail Fig. D.6.

Occultation of a mag. G 10.9 star by Callisto on 2020 June 20, 16:03 UTC. The predicted relative velocity of the event is 18.9 km s−1.

In the text
thumbnail Fig. D.7.

Occultation of a mag. G 11.3 star by Europa on 2020 June 22, 02:09 UTC. The predicted relative velocity of the event is 26.1 km s−1.

In the text
thumbnail Fig. D.8.

Occultation of a mag. G 5.8 star by Io on 2021 April 2, 10:24 UTC. The predicted relative velocity of the event is 16.5 km s−1.

In the text
thumbnail Fig. D.9.

Occultation of a mag. G 11.1 star by Ganymede on 2021 April 25, 07:55 UTC. The predicted relative velocity of the event is 27.7 km s−1.

In the text
thumbnail Fig. D.10.

Occultation of a mag. G 10.4 star by Callisto on 2021 May 4, 23:01 UTC. The predicted relative velocity of the event is 16.3 km s−1.

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.