Get 20M+ Full-Text Papers For Less Than $1.50/day. Start a 14-Day Trial for You or Your Team.

Learn More →

The 2014-2015 Brazilian Mutual Phenomena campaign for the Jovian satellites and improved results for the 2009 events

The 2014-2015 Brazilian Mutual Phenomena campaign for the Jovian satellites and improved results... The 2014-2015 Brazilian Mutual Phenomena campaign for the Jovian satellites and improved results for the 2009 events. a,b,c,< a,b,c b,c a,b,d e,f B. Morgado , R. Vieira-Martins , M. Assafin , A. Dias-Oliveira , D. I. Machado , J. I. a,b g h h i,a,b a,b B. Camargo , M. Malacarne , R. Sfair , O. C. Winter , F. Braga-Ribas , G. Benedetti-Rossi , h h h a,b,c,h g L. A. Boldrin , B. C. B. Camargo , H. S. Gaspar , A. R. Gomes-Júnior , J. O. Miranda , T. h f de Santana and L. L. Trabuco Observatório Nacional/MCTIC, R. General José Cristino 77, Rio de Janeiro, RJ 20.921-400, Brazil Laboratório Interinstitucional de e-Astronomia - LIneA and INCT do e-Universo, Rua Gal. José Cristino 77, Rio de Janeiro, RJ 20921-400, Brazil Observatório do Valongo/UFRJ, Ladeira Pedro Antonio 43, Rio de Janeiro, RJ 20080-090, Brazil Escola SESC de Ensino Médio, Avenida Ayrton Senna, 5677, Rio de Janeiro - RJ, 22775-004, Brazil Universidade Estadual do Oeste do Paraná (Unioeste), Avenida Tarquínio Joslin dos Santos 1300, Foz do Iguaçu, PR 85870-650, Brazil Polo Astronômico Casimiro Montenegro Filho/FPTI-BR, Avenida Tancredo Neves 6731, Foz do Iguaçu, PR 85867-900, Brazil Universidade Federal do Espírito Santo, Av. Fernando Ferrari 514, Vitória, ES 29075-910, Brasil UNESP - São Paulo State University, Grupo de Dinâmica Orbital e Planetologia, CEP 12516-410, Guaratinguetá, SP 12516-410, Brazil Federal University of Technology - Paraná (UTFPR/DAFIS), Av. Sete de Setembro, 3165, Curitiba, PR 80230-901, Brazil A R T I C L E I N F O A B S T R A C T Keywords: Progress in astrometry and orbital modelling of planetary moons in the last decade enabled better Methods: Data analysis determinations of their orbits. These studies need accurate positions spread over extended periods. Astrometry We present the results of the 2014-2015 Brazilian campaign for 40 mutual events from 47 observed Planets and satellites: individual: Io, light curves by the Galilean satellites plus one eclipse of Amalthea by Ganymede. We also reanalysed Europa, Ganymede, Callisto, Amalthea. and updated results for 25 mutual events observed in the 2009 campaign. All telescopes were equipped with narrow-band filters centred at 889 nm with a width of 15 nm to eliminate the scattered light from Jupiter. The albedos’ ratio was determined using images before and after each event. We simulated images of moons, umbra, and penumbra in the sky plane, and integrated their fluxes to compute albedos, simulate light curves and fit them to the observed ones using a chi-square fitting procedure. For that, we used the complete version of the Oren-Nayer reflectance model. The relative satellite positions mean uncertainty was 11.2 mas (í35 km) and 10.1 mas (í31 km) for the 2014-2015 and 2009 campaigns respectively. The simulated and observed ASCII light curve files are freely available in electronic form at the Natural Satellites DataBase (NSDB). The 40/25 mutual events from our 2014-2015/2009 campaigns represent a significant contribu- tion of 17%/15% in comparison with the PHEMU campaigns lead by the IMCCE. Besides that, our th result for the eclipse of Amalthea is only the 4 such measurement ever published after the three ones observed by the 2014-2015 international PHEMU campaign. Our results are suitable for new orbital/ephemeris determinations for the Galilean moons and Amalthea. 2009, 2017). 1. Introduction The uncertainty of the positions obtained from mutual Mutual phenomena between natural satellites – occul- phenomena is usually smaller than the ones obtained by other tations and eclipses – have been successfully used to im- methods. For instance, classical CCD astrometry achieves prove the orbital studies of these moons. For the Galilean uncertainties around 100 mas (í 300 km) (Kiseleva et al., satellites, they have been systematically observed since 1976 2008). For satellite-pair distances, the uncertainties are at (Aksnes and Franklin, 1976). These phenomena occur as the the 30 mas level (í 90 km) (Peng et al., 2012). Mutual ap- Earth and the Sun cross the orbital plane of the satellites. For proximations, based in the same geometrical configuration Jupiter, they happen every six years. of mutual occultations, achieve uncertainties at the 10 mas The photometry of these events offers a reliable source of level (í 30 km) (Morgado et al., 2016, 2019). very precise relative positions between two satellites. They In this paper, we present results for 47 light curves, 31 often achieve uncertainties bellow 5 mas (í 15 km) (Emelyanov, occultations and 16 eclipses, representing 40 mutual events 2009; Dias-Oliveira et al., 2013; Arlot et al., 2014a; Saquet between the Galilean moons observed by three stations in et al., 2018). These relative positions can constrain the or- Brazil, during the 2014-2015 mutual phenomena campaign. bital studies of these moons and give us hints about their We also present one event, an eclipse involving the inner structure and formation processes (Lainey et al., 2004b,a, satellite Amalthea (J5). We also used our improved methods to re-analyse 25 light curves, 13 occultations and 12 eclipses, Based in part on observations made at the Laboratório Nacional de Astrofísica (LNA), Itajubá-MG, Brazil. of 25 mutual phenomena observed by our group during the Corresponding author 2009 mutual phenomena campaign. We compared the new Morgado.fis@gmail.com (B. Morgado) results with those by Dias-Oliveira et al. (2013), Arlot et al. ORCID(s): 0000-0003-0088-1808 (B. Morgado) B. Morgado, R. Vieira-Martins, M. Assafin et al.: Preprint submitted to Elsevier Page 1 of 11 arXiv:1909.05238v1 [astro-ph.EP] 11 Sep 2019 Mutual Phenomena between the Jovian moons (2014a) and Morgado et al. (2016). In Section 2 we detail the observational campaigns. In Section 3 we present the photometry used to produce the ob- served light curves and describe the new, improved light- curve fitting procedures developed and used in this work. Section 4.1 contains new results from the re-analysis of 25 mutual events observed in 2009 and the comparison with the older results. In Section 4.2, we present the results for the 47 light curves involving 40 mutual events observed from Brazil during the 2014-2015 campaign. In Section 4.3, we present the result for the eclipse involving Amalthea. Our conclusions are set on Section 5. Figure 1: Image of Jupiter, Io (J1), Europa (J2) and Callisto (J4) obtained with the 0.6 m diameter Zeiss telescope of the 2. Mutual phenomena campaign details Observatório Pico dos Dias, equipped with a methane filter Every six years, during Jupiter equinox, we can observe on 2014 November 19. The planet and the satellites present mutual occultations and eclipses between Jupiter’s regular about the same brightness due to the use of the narrow-band satellites. The results presented here come from the col- filter, centred at  = 889 nm with 15 nm width. laboration between five Brazilian institutes. The prediction of these events was provided by the Institut de Mécanique Céleste et de Calcul des Ephémérides (IMCCE) (Arlot et al., this filter because in this spectral region, Jupiter’s albedo 2014a,b). drops to 0.05 due to the absorption in the upper atmosphere The 2009 mutual phenomena campaign was the first large (Karkoschka, 1994, 1998). Figure 1 shows an example of an attempt of the kind carried out in Brazil for the Galilean image obtained with the 0.6 m Zeiss telescope from OPD. moons. Observations and instruments are described in de- This filter has been successfully used in the 2009 mutual tail in Dias-Oliveira et al. (2013). We re-analysed 25 mutual phenomena campaign (Dias-Oliveira et al., 2013) and the events encompassing 25 light curves, 13 from occultations mutual approximation campaigns started in 2014 (Morgado and 12 from eclipses, and obtained new results for this cam- et al., 2016, 2019). paign. Discrepancies between the results obtained by Dias- Oliveira et al. (2013) and by Arlot et al. (2014a) motivated 3. Light curve analysis this re-analysis, see more details in Section 4.1. The last Brazilian mutual phenomena campaign of 2014- In mutual phenomena, one can determine relative posi- 2015 obtained data from three telescopes spread on the South tions between the satellites through the analysis of the events’ and South-East of Brazil, with apertures ranging between light curves. In our procedure, we simulate theoretical light 28 and 60 cm. We obtained 47 light curves, 31 for occulta- curves and use them to fit the observed ones. tions and 16 for eclipses, from 40 events observed between The parameters of interest are: (i) the impact parameter November 2014 and June 2015. In Table 1 we present the (s ), the smallest apparent angular distance in the sky plane stations, observers and instrumental details of each station. between both satellite’s centres in the case of occultations It also contains the number of light curves obtained by each or between the eclipsed satellite centre and the centre of the observatory. Moreover, we added the Minor Planet Center eclipsing shadow in the sky plane for eclipses, both cases in (MPC) observatory code of the station (XXX for the station a topocentric frame; (ii) the central instant (t ), the instant of without a code). time that this smallest distance occurs; and (iii) the apparent In Table 2, we list the observational details of each ob- relative velocity (v ) between both satellites in the sky plane. served event. It contains the date of the event and the satel- In the supplementary material, we also provide the inter- lites’ pairs in the form "S oS " for occultation and "S eS " 1 2 1 2 satellite tangential coordinates (X, Y ). For occultations, these for eclipses, where 1 stands for Io, 2 for Europa, 3 for Ganymede coordinates between both satellitesâĂŹ centres are in a topocen- and 4 for Callisto. We furnish the sites involved in each ob- tric frame. For mutual eclipses, these coordinates are in servation (using the alias defined in Table 1). For each event, a topocentric frame and the mean difference between the we give the solar phase angle (i) and the zenith distance (z), eclipsed satellite centre and the centre of the eclipsing satel- both in degrees. In the last column, we list the instrumental liteâĂŹs shadow in the sky plane. albedos’ ratio of the involved satellites (and its uncertainty), determined by using images before and after the event. This 3.1. Obtaining the observed light-curves albedos’ ratio is only needed for occultations. Firstly all images were corrected by Bias and Flat-Field In all observations we used a narrow band filter cen- using standard procedures with the Image Reduction and tred at 889 nm with a width of 15 nm. This bandpass is in Analysis Facility (IRAF) (Butcher and Stevens, 1981). Then, the methane absorption region of the spectrum. We chose we determined the light flux of the targets in the images by differential aperture photometry using the PRAIA package Website: http://nsdb.imcce.fr/multisat/nsszph515he.htm B. Morgado, R. Vieira-Martins, M. Assafin et al.: Preprint submitted to Elsevier Page 2 of 11 Mutual Phenomena between the Jovian moons Table 1 2014-2015 mutual phenomena campaign observer list. City/ Country Longitude Observers Telescope N of positive Station alias Latitude CCD detections MPC code Height Itajubá/ MG, Brazil 45 34’ 57" W B. Morgado 60 cm 20 OPD 22 32’ 07" S H. S. Gaspar Andor/IKon-L 874 1.864 km R. Sfair B. C. B. Camargo T. Santana L. A. Boldrin M. E. Huaman G. Benedetti-Rossi A. R. Gomes-Júnior Foz do Iguaçu/ PR, Brazil 54 35’ 37" W D. I. Machado 28 cm 22 FOZ 25 26’ 05" S L. L. Trabuco SBIG/ST-7X-ME X57 0.184 km Vitória/ ES, Brazil 40 19’ 00" W M. Malacarne 35 cm 5 GOA 20 17’ 52" S J. O. Miranda SBIG/ST-8X-ME XXX 0.026 km (Assafin et al., 2011). During an occultation, both satellites the time and the dynamical parameters using equation (1) are measured together in the same aperture, and a third satel- (Assafin et al., 2009). lite is used as calibrator. In the case of eclipses, the eclipsed satellite is measured alone in the aperture and the eclipsing 2 2 2 s.t/ = s + v .t * t / (1) satellite (or any other) is used as calibrator. The light curve 0 0 is then normalised by a polynomial fit so that the flux ratio We also need a reflectance model to take into account the outside the flux drop gets equal to 1.0, and the flux drop can phase effect and how the surface of the satellite will reflect be adequately evaluated. the sunlight. For eclipses, we further need some information During the photometry of mutual events, one must take about the Sun, such as its radius and a model to consider the care with the possibility of a parasitic flux, as pointed out Sun’s limb darkening. by Emel’yanov (2017) and Arlot et al. (2017). The origin We could do simulations with triaxial bodies with vary- of this flux is likely to be the background (mostly Jupiter’s ing albedo, but not for practical purposes, due to photometry scattered light) or the CCD detector. In our case, we attenu- limitations. Thus, satellites are considered as spheres with a ate this parasitic flux with the Methane filter and a rigorous known appropriate radius, and the albedo is uniform along calibration process. Tests showed that the parasitic flux in the surface. The relative velocity between both satellites is our images is one order of magnitude below the noise of our constant during the mutual event (a few minutes only). For observations. notation, the occulting/eclipsing is denoted Sat and Sat is 1 2 the occulted/eclipsed one. 3.2. Simulating light curves In the simulations of occultations and eclipses, we used The procedures utilised here follow the same principles the same geometric relations described in detail in Dias-Oliveira outlined in Assafin et al. (2009) and Dias-Oliveira et al. (2013). et al. (2013). However, improvements were made in almost every step, as explained in sections 3.2.1, 3.2.2 and 3.3. 3.2.1. Occultation The parameters needed in the modelling of a mutual oc- The first step in simulating a light curve of a mutual oc- cultation or eclipse can be separated into two complemen- cultation is the production of a 2D satellite apparent pro- tary types. The first refers to physical characteristics: sizes file, simulating how the body reflects the sunlight as seen by and shapes of each satellite, and the satellites’ albedo for oc- an observer on Earth. As pointed out by Vasundhara et al. cultations. The second type relates to dynamics, and are the (2017), it is essential to use a realistic intensity distribution parameters of interest: s , t and v . 0 0 0 for the satellite. Physical parameters such as radius and shape are known However, this approach demands previous knowledge about from space probes’ data. Albedos are determined from aux- the satellite surface (albedo maps) that can change with time iliary observations made before and after events with the or even for different effective wavelengths, it is important to same instrument setup. The apparent separation in the sky highlight that these maps are not know for the wavelength plane between both satellites can be written as a function of of our observations ( = 889 nm). The same applies to the B. Morgado, R. Vieira-Martins, M. Assafin et al.: Preprint submitted to Elsevier Page 3 of 11 Mutual Phenomena between the Jovian moons Hapke scattering law (Hapke, 1981; Hapke and Wells, 1981; curve was then normalised by the sums of individual fluxes Hapke, 1984, 1986, 2002, 2008, 2012) used by Emelyanov (F + F ). Figure 2 shows a simulated light curve for the 1 2 (2009), Arlot et al. (2014a) and Saquet et al. (2018), which event when Europa occulted Io on February 22, 2015. The requires unknown parameters in the wavelength band of our black dots indicate seven instants for which the respective observations. 2D profiles are displayed on each corresponding box. For We successfully solved these issues by adopting a gener- this event the albedo ratio was 0.960, as determined before alisation of Lambert scattering law given by Oren and Nayar the event. (1994). The Oren-Nayer model takes into account the direc- 3.2.2. Eclipse tion of radiance and the roughness of the surface in a natural Following Dias-Oliveira et al. (2013), we developed a way, so that the reflectance depends only of the albedo and numerical 2D mask that incorporates the two regions of the in one more parameter that very smoothly tunes a wide range eclipsing satellite’s shadow, the umbra and the penumbra. of surface roughness, and most importantly, regardless of the This mask was then applied to the 2D profile of the eclipsed wavelength. This model realistically reproduces the illumi- satellite, considering the separation between both satellites nation of an object in modern computer graphic scenes for as seen from the heliocentre. Once again, the spatial resolu- movies and for the full Moon. (Oren and Nayar, 1994). In tion was set as 1 mas (í 3 km). Dias-Oliveira et al. (2013) a simplified version of the model For the penumbra region, the fraction of sunlight that was used. Here, we implemented the complete version in reaches the eclipsed satellite was computed by using a nu- Oren and Nayar (1994), taking into account the direct illu- merical method. The solar limb darkening was taken into mination and all inter-reflection components of the radiance. account by using Hestroffer and Magnan (1998) empirical The albedo ratio between the satellites is determined in- law, with parameters set for the 889 nm spectral region. dependently, by using observations right before and after The light flux was numerically integrated for a given in- the mutual event with the same instrument setup. In Dias- stant by using the profile after the mask was applied (F.t/). Oliveira et al. (2013), analytic expressions involving the ter- Then the light curve was normalised using the Light flux of minator were used to take solar phase angle effects into ac- the eclipsed satellite (F ). Figure 3 is a simulated light curve count in the evaluation of the flux measurements of albedo for the event when Europa eclipsed Ganymede on April 12, observations. Here, following a more rigorous approach, 2015. The black dots represent seven instants for which the we also simulated the 2D profiles of the satellites for these respective 2D profiles are displayed. observations, to better determine the effective area and re- flectance of the satellites to compute more accurate albedos. 3.3. Fitting procedure In fact, first we measured the light fluxes between both We took a somewhat different approach from that in Dias- satellites (F , F ) separately using images right before or af- 1 2 Oliveira et al. (2013). Here, the parameters of interest were ter the occultation. In the other hand, we compute the 2D determined by the minimisation of the Chi-square test, Eq. profile of each satellites for the given instants and obtained (3), where the simulated light curve is compared with the the simulated light flux for each satellite (F , F ), this S1 S2 observed one. simulations already take into consideration the size of the satellite and the phase angle. The ratio of albedo (A _A ) 1 2 can be determined using Equation (2). We considered each .LC * LC / obs model satellite’s albedo as uniform, due to the lack of information .N * P/ = (3) about albedo variations for the wavelength of our observa- tions. N is the number of observations used in the process, and P is the number of parameters fitted (P = 3).  is the light curve’s standard deviation outside the event (the noise). The F A F 1 1 s1 = : : (2) parameters (s ,t and v ) for which the chi-square is mini- 0 0 0 F A F 2 2 s2 2 mum ( =  ) were set as the solution for the obser- min vations. The normalised  is expected to be around 1 for Similar to Dias-Oliveira et al. (2013) we discretise the good fittings. satellite into a 2D profile with elements of a given spatial Initially, a large range of parameter values was tested, resolution. However, unlike Dias-Oliveira et al. (2013), we then the ranges were narrowed as the iterative process pro- created 2D satellite profiles with much better spatial reso- ceeded. For computing speed, the spatial resolution was de- lutions, 1 mas (í 3 km), avoiding eventual round off errors graded in the first steps and then is set to the nominal value in the simulated flux counts. The profiles were positioned of 1 mas as we approached the  minimum. for each instant of time (t) using the separation between the The uncertainty of each parameter (1 error bar) was de- satellites, overlapping the occulted satellite with the occult- termined by changing that parameter from its nominal solu- ing one when necessary, obtaining a new combined 2D pro- 2 2 2 tion value, so that  changes from  to  + 1. This file of both satellites. The light flux (F.t/) was numerically min min procedure is repeated for each parameter: t , s and v . integrated over the combined 2D profile for a given instant 0 0 0 The procedures described in Section 3 were developed (t). Then, this is repeated for every instant (t) of the event as a PYTHON software that analyses and fits observed light to produce the simulated light curve. This simulated light B. Morgado, R. Vieira-Martins, M. Assafin et al.: Preprint submitted to Elsevier Page 4 of 11 Mutual Phenomena between the Jovian moons Figure 2: Occultation simulation. Io was occulted by Europa on February 22, 2015, at 02:07 UTC as seen at OPD. The bottom panel shows the simulated light curve and seven instants (black dots) were highlighted. The 2D profiles for each instant are displayed in the top panel. The albedo ratio was 0.960, as determined before the event. The profile resolution is 1 mas (í 3.2 km). Figure 3: Eclipse simulation. Ganymede was eclipsed by Europa on April 12, 2015 at 01:46 UTC as seen at FOZ. The bottom panel shows the simulated light curve and seven instants (black dots) were highlighted. The 2D profiles for each instant are displayed in the top panel. The profile resolution was 1 mas (í 3.6 km). curves. This software uses functions from NUMPY, ASTROPY 4. Results of the mutual phenomena (Astropy Collaboration et al., 2013), SCIPY and MATPLOTLIB campaigns libraries. 4.1. Re-analysis of the mutual phenomena between the Galilean moons - 2009 From April to October of 2009, 25 light curves for 25 mutual events between the Galilean moons, 13 occultations and 12 eclipses, were observed by our group with the 60 cm Zeiss telescope of the Observatório Pico dos Dias (OPD) us- ing the methane filter. These events were analysed by Dias- B. Morgado, R. Vieira-Martins, M. Assafin et al.: Preprint submitted to Elsevier Page 5 of 11 Mutual Phenomena between the Jovian moons Oliveira et al. (2013) and by Arlot et al. (2014a). However, the results obtained by both presented a significant difference (higher than 3). More recently, Morgado et al. (2016) analysed 5 mutual approximations that were auxiliary observations designed for albedo determinations made immediately before and af- ter mutual occultations from the 2009 campaign. The mu- tual approximation technique is based on the same geometri- cal parameters that describe a mutual occultation (Morgado et al., 2016, 2019). From the measurements of the separation between both satellites, this technique allows for determin- ing the time of maximum apparent approach between these satellites, which can be directly compared with the central instant from the occultation. The results obtained with this independent technique agree within the errors with the re- sults obtained by Arlot et al. (2014a). The discrepancy of the results between Dias-Oliveira et al. Figure 4: Light curve of the event where Io was occulted by Europa on February 22 2015, observed with the 0.60 m tele- (2013) and Arlot et al. (2014a), the agreement between the scope at OPD. The measured normalised flux is denoted by later with Morgado et al. (2016) for 5 events, and the sig- black dots and the fitted model represented by the red line. In nificant improvements in our light curve fitting procedure the bottom panel, the red dots are the residuals of the fitting motivated us to re-analyse our 2009 campaign data. in the sense observation minus fit. Here, we present our updated results for the 2009 events in Table 3. Table 3 is organised as follows: the date of the event and the satellites’ pairs in the form "S oS " for occul- 1 2 (2014a) and Morgado et al. (2016) within 1. We have a rms tation and "S eS " for eclipses, where 1 stands for Io, 2 for 1 2 of 9.9 mas (í 30 km) and 14.8 mas (í 45 km) in comparison Europa, 3 for Ganymede and 4 for Callisto. We furnish the 3 4 with the JPL’s and the IMCCE’s ephemeris. sites involved in each observation (using the alias defined in Table 1). For each event, we give the obtained central in- 4.2. Mutual phenomena between the Galilean stant (t ) and its uncertainty (t ) in seconds of time, and the 0 0 moons - 2014-2015 ephemeris offset (t ) in mas; the impact parameter (s ), its 0 0 Here we present the results concerning the latest cam- uncertainty (s ) and its offset (s ) in mas; and the rela- 0 0 paign. We obtained new 47 light curves, 31 occultations tive velocity (v ), its uncertainty (v ) and its offset (v ) in 0 0 0 and 16 eclipses, from 40 events observed by 3 stations in mas per second. All times are UTC. In the last columns, we the South and South-East of Brazil. have the rms between the observed light fluxes and the fitted An example is the event where Europa occulted Io on ones, the number of images utilised (N) and the normalised February 22 2015. The observed light curve is illustrated of our fit. in Figure 4. In the upper panel, the black dots are the light The corresponding inter-satellite tangential coordinates flux observed and the red line the model fitted. The bottom (X and Y ) in the sense occulting/eclipsing satellite minus oc- panel contains the residuals in the sense observation minus culted/eclipsed satellite for the central instant can be found in model. For this event, the offset for the central instant was the supplementary material, such form is the same presented +6.5 mas (í 21 km) and for the impact parameter -0.8 mas by Emelyanov and Gilbert (2006); Emelyanov (2009); Ar- (í 3 km). The offsets regard to the JPL’s jup310 and the lot et al. (2014a); Saquet et al. (2018). The plots of the DE435 ephemeris. re-fitted light curves are available as online material in the A second example is the case when Europa eclipsed Ganymede supplementary material. The simulated and observed ASCII on April 12 2015. The observed light curve is illustrated in light curve files are freely available in electronic form at the Figure 5. For this event, the offset for the central instant was NSDB . +13.6 mas (í 49 km), and for the impact parameter +9.6 The re-analysis resulted in a mean uncertainty of 15.3 mas (í 35 km). mas (í 46 km) for the impact parameter and 4.9 mas (í 15 The multiple coverage observational strategy reduced the km) for the central instant. In Table 4, we compare the up- number of events lost by overcast weather or instrumental is- dated results with the ones from Arlot et al. (2014a), Dias- sues. An example was the Io occultation by Ganymede on Oliveira et al. (2013) and Morgado et al. (2016). The error March 03 2015, OPD and FOZ observed this event. Also, of each parameter normalises the differences. If the value is two other stations in the USA observed this event, one in less than one, both results agree within 1. At the bottom, Arnold (AAC) and another in Scottsdale (SCO). These ob- we have the mean difference and the standard deviation for servations were made in the context of the international mu- each parameter. The JPL ephemeris utilised was jup310 and DE435. The updated results now agree with those by Arlot et al. The IMCCE ephemeris utilised was NOE-5-2010-GAL and DE435. Website: http://nsdb.imcce.fr/nsdb/home.html B. Morgado, R. Vieira-Martins, M. Assafin et al.: Preprint submitted to Elsevier Page 6 of 11 Mutual Phenomena between the Jovian moons the satellites’ pairs in the form "S oS " for occultation and 1 2 "S eS " for eclipses, where 1 stands for Io, 2 for Europa, 3 1 2 for Ganymede and 4 for Callisto. We furnish the sites in- volved in each observation (using the alias defined in Ta- ble 1). For each event, we give the obtained central instant (t ) and its uncertainty (t ) in seconds of time, and the 0 0 ephemeris offset (t ) in mas; the impact parameter (s ), 0 0 its uncertainty (s ) and the offset (s ) in mas; and the rel- 0 0 ative velocity (v ), its uncertainty (v ) and the offset (v ) 0 0 0 in mas per second. All times are UTC. In the last columns, we have the rms between the observational curves and the fitted ones, the number of images utilised (N) and the nor- malised  of our fit. The corresponding inter-satellite tangential coordinates (X and Y ) in the sense occulting/eclipsing satellite minus oc- culted/eclipsed satellite for the central instant can be found in Figure 5: Light curve of the event when Ganymede was the supplementary material, such formalism is the same pre- eclipsed by Europa on April 12 2015, observed with the 0.28 sented by Emelyanov and Gilbert (2006); Emelyanov (2009); m telescope at FOZ. Arlot et al. (2014a); Saquet et al. (2018). The plots of the fit- ted light curves are available as online material in the supple- mentary material. The simulated and observed light curve ASCII files are freely available in electronic form at the NSDB. The mean uncertainty of our results is 14.8 mas (í 45 km) for the impact parameter and 7.5 mas (í 23 km) for the central instant. The rms relative to JPL ephemeris was 9.2 mas (í 28 km) and 13.5 mas (í 40 km) relative to IMCCE ephemeris. From the 2014-2015 events, 10 were also analysed using a different procedure and published by Saquet et al. (2018). In average, the comparison between this procedure and ours agrees within 1. 4.3. Amalthea eclipsed by Ganymede - 02 March One particular event in our 2014-2015 observational cam- paign was the eclipse of Amalthea by Ganymede. The as- trometry of this inner satellite is not easy to be done due to its proximity to Jupiter (major semi-axis equal to 2.54 Jupiter’s radius). Often, coronagraphy techniques are needed to sep- arate this object from Jupiter’s scattered light (Kulyk et al., 2002; Veiga and Vieira Martins, 2005; Robert et al., 2017). The positional uncertainty of classical astrometry for this Figure 6: Light curves of the event where Io was occulted by satellite is in the 120 mas level (í 360 km). Ganymede on March 03 2015, observed at AAC, SCO, OPD The observation of mutual eclipses involving Galilean and FOZ, respectively. The AAC and SCO light curves were reduced and analysed by Saquet et al. (2018). moons and inner satellite was strongly advocated by Vachier et al. (2002). The first registration of this kind of event was given by Christou et al. (2010), regarding three eclipses of Amalthea observed during the 2009 mutual phenomena cam- tual phenomena campaign PHEMU15, (Saquet et al., 2018; paign. More recently, Saquet et al. (2016) also analysed three Emel’yanov, 2017). Both light curves are available at the more eclipses of Amalthea and the first observation of an NSDB. In the Figure 6 we compare our light curves (OPD eclipse of Thebe during the 2014-2015 campaign. and FOZ) with the ones analysed by Saquet et al. (2018) Here we present the results of one eclipse of Amalthea (AAC and SCO). The central instant obtained by the obser- by Ganymede observed on March, 2 of 2015 at the 1.6 m vations agrees within 2. Notice that all curves present sim- Perkin-Elmer telescope of the Observatório Pico dos Dias ilar features and we highlight the small residual in our light (OPD, MPC code: 874). This observation was made using curves. the IKon-L CCD camera with the narrow Methane filter . The results for these events are presented in Table 5. Table 5 is organised as follows: the date of the event and 5 Centred at 889 nm with a width of 15 nm. B. Morgado, R. Vieira-Martins, M. Assafin et al.: Preprint submitted to Elsevier Page 7 of 11 Mutual Phenomena between the Jovian moons Figure 7: Digital coronagraphy on an image containing part of Jupiter and its inner satellite Amalthea, as obtained on 02 March, 2015, with the 1.6 m Perkin-Elmer telescope. The left panel (i) displays the original image. The central panel (ii) shows the clean bright object profile (Jupiter) obtained. The right panel (iii) displays the final coronagraphed image (see text). After correcting by Bias and Flat-Field using the same procedure described above, we applied a digital coronagra- phy technique to reduce the influence of Jupiter brightness in the images, this coronagraphy was done using the PRAIA package (Assafin et al., 2008, 2009). Briefly, the procedure is as follows. The centroid of the bright object is iteratively determined. Concentric rings with radius R are formed for each image pixel at a distance R to the centroid. Quartile statistics of weighted fluxes inside each ring are performed, and for each image pixel an average count is assigned. The result is an improved profile with cleaner pixel counts that better represent the bright object. The profile is then sub- tracted from the original image, resulting in the final coron- agraphed (science) image, see Figure 7. Aperture photometry was done using the PRAIA package, where the size of the aperture was manually determined to Figure 8: Light curve of the event when Amalthea was eclipsed maximise the signal to noise ratio. The light curve simula- by Ganymede on March 02 2015, observed with the 1.6 m tion and fitting procedure were the same described in Sec- telescope at OPD. tions 3.2.2 and 3.3. Notice that Amalthea’s triaxial shape is 125  73  64 km (uncertainty of 2 km in each axis; Thomas et al. (1998)) and its rotation phase during the event was un- 5. Conclusions known. However, without any loss of precision, in our sim- ulations, we considered Amalthea as an equivalent sphere We presented in this paper the results for 40 mutual events with a radius equal to 83.5 km. Because of the time resolu- from the observation and analysis of 47 light curves, 31 oc- tion of the observations (8 seconds), the spherical shape was cultations and 16 eclipses, obtained during the 2014-2015 indistinguishable from the elliptical one. mutual phenomena campaign between the Galilean satellites. The light curve of this event is illustrated in Fig. 8. We The observations were made at three stations in the South obtained a central instant with an uncertainty of 19.4 mas and South-East of Brazil, using telescopes with diameters (í 58.2 km) and an impact parameter with an uncertainty of ranging between 28 and 60 cm. We also obtained updated 76.3 mas (í 228 km). This corresponds to a mean uncer- results from the re-analysis of 25 mutual events, 13 occul- tainty of 47.8 mas (í 143 km). The result of this event is tations and 12 eclipses, observed in Brazil by our group in displayed in the last line of Table 5. The positions obtained 2009 with a 60 cm aperture telescope. In all observations, by Christou et al. (2010) had mean uncertainty of 82 mas (í we used a narrow band methane filter centred at 889 nm with 246 km) and the ones obtained by Saquet et al. (2016) had a a width of 15 nm, that eliminates Jupiter’s scattered light. mean uncertainty of 45 mas (í 135 km). We used the Oren-Nayer model (Oren and Nayar, 1994) to characterise the reflectance of the surface of the satellites. It is a generalisation of Lambert’s scattering law. The main advantage of this model is that it does not require previous B. Morgado, R. Vieira-Martins, M. Assafin et al.: Preprint submitted to Elsevier Page 8 of 11 Mutual Phenomena between the Jovian moons knowledge about the satellite surface, and delivers excellent sion Europa Clipper , scheduled to be launched in the next results. The albedo ratio was instrumentally obtained by us- decade (2020s). ing satellite observations before and after the mutual occulta- tions, and the same light curve simulation routines that take Acknowledgements into account solar phase angle and surface reflectance. In We thank our anonymous referee and N. Emelyanov for our procedure, the simulated light curves fitted to the ob- helpful comments. This study was financed by the Coorde- served ones had a normalised chi-square very close to 1.0, nação de Aperfeiçoamento de Pessoal de Nível Superior - indicating good agreement of our model to the data. Brasil (CAPES) - Finance Code 001. Part of this research The updated results for 2009 now agree within 1 with is suported by INCT do e-Universo, Brazil (CNPQ grants those from Arlot et al. (2014a) and from Morgado et al. (2016). 465376/2014-2). Based in part on observations made at the The internal mean uncertainty was 10.1 mas (í 31 km). For Laboratório Nacional de Astrofísica (LNA), Itajubá-MG, Brazil. the 2014-2015 campaign, the mean uncertainty was 11.2 mas BM thanks the CAPES/Cofecub-394/2016-05 grant. RVM (í 35 km). There is no significant difference between the acknowledges the grants: CNPq-304544/2017-5, 401903/2016- quality of the data for both campaigns. Our result is compa- 8, Capes/Cofecub-2506/2015, Faperj: PAPDRJ-45/2013 and rable with other independent observations during the same E-26/203.026/2015. MA thanks CNPq (Grants 427700/2018- campaign, such as the 23 light curves observed by Vasund- 3, 310683/2017-3 and 473002/2013-2) and FAPERJ (Grant hara et al. (2017) and the 21 light curves observed by Zhang E-26/111.488/2013). JIBC acknowledges CNPq grants 308489/2013- et al. (2019), which report uncertainties in the 10 and 20 mas 6 and 308150/2016-3. RS and OCW acknowledges Fapesp level (30 and 60 km), respectively. Saquet et al. (2018) pub- proc. 2016/24561-0 and 2011/08171-3, CNPq proc. 312813/2013- lished the light curves of the international campaign organ- 9 and 305737/2015-5. FBR acknowledges CNPq support, ised by the IMCCE. It reunites 609 light curves with a mean proc. 309578/2017-5. GBR thanks to the support of the internal uncertainty of 24 mas (í 75 km) . From these 609 CAPES and FAPERJ/PAPDRJ (E26/203.173/2016) grants. light curves, 10 were observed by our group, the parameters ARGJ thanks FAPESP proc. 2018/11239-8. This collabora- obtained by Saquet et al. (2018) agree, on average, within 1 tion, as part of the Encelade working group, has been sup- with the results presented here. ported by the International Space Sciences Institute (ISSI) Compared to the 236 mutual events covered by the 2014- in Bern, Switzerland. 2015 international PHEMU campaign lead by the IMCCE, Observatoire de Paris (Saquet et al., 2018), the 40 ones cov- ered by our 2014-2015 campaign (with only 5 events in com- References mon) represent a significant contribution of about 17%. No- Aksnes, K., Franklin, F.A., 1976. Mutual phenomena of the Galilean satel- tice that this campaign was favourable for the north hemi- lites in 1973. III. Final results from 91 light curves. Astronomical Journal sphere, enhancing the weight of our southern results due to 81, 464–481. doi:10.1086/111908. parallax effects. Equally, our new results for the 25 events Arlot, J.E., Cooper, N., Emelyanov, N., Lainey, V., Meunier, L.E., Mur- ray, C., Oberst, J., Pascu, D., Pasewaldt, A., Robert, V., Tajeddine, R., represent about 15% of the 172 mutual events covered by Willner, K., 2017. Natural satellites astrometric data from either space the international PHEMU campaign of 2009 (Arlot et al., probes and ground-based observatories produced by the European con- 2014a). In a similar way, our result for the eclipse by Amalthea sortium “ESPaCE”. Notes Scientifiques et Techniques de l’Institut de th is only the 4 such measurement ever published for the 2014- Mecanique Celeste 105. 2015 campaign, after the 3 ones observed by Saquet et al. Arlot, J.E., Emelyanov, N., 2019. Natural satellites mutual phenomena ob- (2016), representing a significant contribution to the orbit of servations: Achievements and future. Planetary and Space Science 169, 70–77. doi:10.1016/j.pss.2019.02.004. this inner satellite of Jupiter. All the data are freely available Arlot, J.E., Emelyanov, N., Varfolomeev, M.I., Amossé, A., Arena, C., As- to anyone at NSDB for further research and orbital fitting. safin, M., Barbieri, L., Bolzoni, S., Bragas-Ribas, F., Camargo, J.I.B., All these results can be used to improve the orbit and Casarramona, F., Casas, R., Christou, A., Colas, F., Collard, A., Combe, ephemeris of the Galilean satellites (plus Amalthea) taking S., Constantinescu, M., Dangl, G., De Cat, P., Degenhardt, S., Delcroix, M., Dias-Oliveira, A., Dourneau, G., Douvris, A., Druon, C., Ellington, into account the tidal forces, as pointed out by Lainey et al. C.K., Estraviz, G., Farissier, P., Farmakopoulos, A., Garlitz, J., Gault, (2009). D., George, T., Gorda, S.Y., Grismore, J., Guo, D.F., Herald, D., Ida, The next mutual phenomena events for the Galilean satel- M., Ishida, M., Ivanov, A.V., Klemt, B., Koshkin, N., Le Campion, J.F., lites will occur in 2021 and will favour the southern hemi- Liakos, A., Liao, S.L., Li, S.N., Loader, B., Lopresti, C., Lo Savio, E., sphere, due to Jupiter’s declination. An observational cam- Marchini, A., Marino, G., Masi, G., Massallé, A., Maulella, R., McFar- land, J., Miyashita, K., Napoli, C., Noyelles, B., Pauwels, T., Pavlov, paign such as this one will be organised in due time, the H., Peng, Q.Y., Perelló, C., Priban, V., Prost, J., Razemon, S., Rous- prediction of these events are already in the IMCCE web- selle, J.P., Rovira, J., Ruisi, R., Ruocco, N., Salvaggio, F., Sbarufatti, site (Arlot and Emelyanov, 2019). These campaigns can G., Shakun, L., Scheck, A., Sciuto, C., da Silva Neto, D.N., Sinyaeva, increase the accuracy and precision of ephemeris and can be N.V., Sofia, A., Sonka, A., Talbot, J., Tang, Z.H., Tejfel, V.G., Thuillot, helpful to space missions aimed at the Jovian system. For ex- W., Tigani, K., Timerson, B., Tontodonati, E., Tsamis, V., Unwin, M., Venable, R., Vieira-Martins, R., Vilar, J., Vingerhoets, P., Watanabe, H., ample, we have the ESA mission JUICE and NASA’s mis- Yin, H.X., Yu, Y., Zambelli, R., 2014a. The PHEMU09 catalogue and The standard deviation after fitting the light curves. astrometric results of the observations of the mutual occultations and Website: http://nsdb.imcce.fr/multisat/nssephme.htm Website: https://www.nasa.gov/europa/. Website: http://sci.esa.int/juice/. B. Morgado, R. Vieira-Martins, M. Assafin et al.: Preprint submitted to Elsevier Page 9 of 11 Mutual Phenomena between the Jovian moons eclipses of the Galilean satellites of Jupiter made in 2009. Astronomy particle phase function hockey stick relation. Icarus 221, 1079–1083. and Astrophysics 572, A120. doi:10.1051/0004-6361/201423854. doi:10.1016/j.icarus.2012.10.022. Arlot, J.E., Saquet, E., Robert, V., Lainey, V., 2014b. The Phemu 2015 Hapke, B., Wells, E., 1981. Bidirectional reflectance spectroscopy. 2. Ex- campaign of observations of the mutual events of the Galilean satellites periments and observations. Journal of Geophysical Research 86, 3055– of Jupiter, in: European Planetary Science Congress, pp. EPSC2014–57. 3060. doi:10.1029/JB086iB04p03055. Assafin, M., Campos, R.P., Vieira Martins, R., da Silva Neto, D.N., Ca- Hestroffer, D., Magnan, C., 1998. Wavelength dependency of the Solar margo, J.I.B., Andrei, A.H., 2008. Instrumental and digital coronagra- limb darkening. Astronomy and Astrophysics 333, 338–342. phy for the observation of the Uranus satellites’ upcoming mutual events. Karkoschka, E., 1994. Spectrophotometry of the Jovian Planets and Titan Planetary and Space Science 56, 1882–1887. doi:10.1016/j.pss.2007. at 300- to 1000-nm Wavelength: The Methane Spectrum. Icarus 111, 05.030. 174–192. doi:10.1006/icar.1994.1139. Assafin, M., Vieira-Martins, R., Braga-Ribas, F., Camargo, J.I.B., da Silva Karkoschka, E., 1998. Methane, Ammonia, and Temperature Measure- Neto, D.N., Andrei, A.H., 2009. Observations and Analysis of Mutual ments of the Jovian Planets and Titan from CCD-Spectrophotometry. Events between the Uranus Main Satellites. Astronomical Journal 137, Icarus 133, 134–146. doi:10.1006/icar.1998.5913. 4046–4053. doi:10.1088/0004-6256/137/4/4046. Kiseleva, T.P., Kiselev, A.A., Kalinichenko, O.A., Vasilyeva, N.A., Assafin, M., Vieira Martins, R., Camargo, J.I.B., Andrei, A.H., Da Silva Khovricheva, M.L., 2008. Results of astrometric observations of Neto, D.N., Braga-Ribas, F., 2011. PRAIA - Platform for Reduction Jupiter’s Galilean satellites at the Pulkovo Observatory from 1986 of Astronomical Images Automatically, in: Gaia follow-up network for to 2005. Solar System Research 42, 414–433. doi:10.1134/ the solar system objects : Gaia FUN-SSO workshop proceedings, pp. S0038094608050055. 85–88. Kulyk, I., Jockers, K., Karpov, N., Sergeev, A., 2002. Astrometric CCD Astropy Collaboration, Robitaille, T.P., Tollerud, E.J., Greenfield, P., observations of the inner Jovian satellites in 1999-2000. Astronomy and Droettboom, M., Bray, E., Aldcroft, T., Davis, M., Ginsburg, A., Price- Astrophysics 383, 724–728. doi:10.1051/0004-6361:20011770. Whelan, A.M., Kerzendorf, W.E., Conley, A., Crighton, N., Barbary, Lainey, V., Arlot, J.E., Karatekin, Ö., van Hoolst, T., 2009. Strong tidal K., Muna, D., Ferguson, H., Grollier, F., Parikh, M.M., Nair, P.H., Un- dissipation in Io and Jupiter from astrometric observations. Nature 459, ther, H.M., Deil, C., Woillez, J., Conseil, S., Kramer, R., Turner, J.E.H., 957–959. doi:10.1038/nature08108. Singer, L., Fox, R., Weaver, B.A., Zabalza, V., Edwards, Z.I., Azalee Lainey, V., Arlot, J.E., Vienne, A., 2004a. New accurate ephemerides for Bostroem, K., Burke, D.J., Casey, A.R., Crawford, S.M., Dencheva, N., the Galilean satellites of Jupiter. II. Fitting the observations. Astronomy Ely, J., Jenness, T., Labrie, K., Lim, P.L., Pierfederici, F., Pontzen, A., and Astrophysics 427, 371–376. doi:10.1051/0004-6361:20041271. Ptak, A., Refsdal, B., Servillat, M., Streicher, O., 2013. Astropy: A Lainey, V., Duriez, L., Vienne, A., 2004b. New accurate ephemerides for community Python package for astronomy. Astronomy and Astrophysics the Galilean satellites of Jupiter. I. Numerical integration of elaborated 558, A33. doi:10.1051/0004-6361/201322068, arXiv:1307.6212. equations of motion. Astronomy and Astrophysics 420, 1171–1183. Butcher, H., Stevens, R., 1981. Image Reduction and Analysis Facility doi:10.1051/0004-6361:20034565. Development. Kitt Peak National Observatory Newsletter 16, 6. Lainey, V., Jacobson, R.A., Tajeddine, R., Cooper, N.J., Murray, C., Robert, Christou, A.A., Lewis, F., Roche, P., Hidas, M.G., Brown, T.M., 2010. Ob- V., Tobie, G., Guillot, T., Mathis, S., Remus, F., Desmars, J., Arlot, J.E., servational detection of eclipses of J5 Amalthea by the Galilean satel- De Cuyper, J.P., Dehant, V., Pascu, D., Thuillot, W., Le Poncin-Lafitte, lites. Astronomy and Astrophysics 522, A6. doi:10.1051/0004-6361/ C., Zahn, J.P., 2017. New constraints on Saturn’s interior from Cassini 201014822, arXiv:1104.0042. astrometric data. Icarus 281, 286–296. doi:10.1016/j.icarus.2016.07. Dias-Oliveira, A., Vieira-Martins, R., Assafin, M., Camargo, J.I.B., Braga- 014, arXiv:1510.05870. Ribas, F., da Silva Neto, D.N., Gaspar, H.S., Pires dos Santos, P.M., Morgado, B., Assafin, M., Vieira-Martins, R., Camargo, J.I.B., Dias- Domingos, R.C., Boldrin, L.A.G., Izidoro, A., Carvalho, J.P.S., Sfair, Oliveira, A., Gomes-Júnior, A.R., 2016. Astrometry of mutual approx- R., Sampaio, J.C., Winter, O.C., 2013. Analysis of 25 mutual eclipses imations between natural satellites. Application to the Galilean moons. and occultations between the Galilean satellites observed from Brazil in Monthly Notices of the Royal Astronomical Society 460, 4086–4097. 2009. Monthly Notices of the Royal Astronomical Society 432, 225– doi:10.1093/mnras/stw1244, arXiv:1605.06573. 242. doi:10.1093/mnras/stt447. Morgado, B., Vieira-Martins, R., Assafin, M., Machado, D.I., Camargo, Emelyanov, N.V., 2009. Mutual occultations and eclipses of the Galilean J.I.B., Sfair, R., Malacarne, M., Braga-Ribas, F., Robert, V., Bassallo, satellites of Jupiter in 2002-2003: final astrometric results. Monthly T., Benedetti-Rossi, G., Boldrin, L.A., Borderes-Motta, G., Camargo, Notices of the Royal Astronomical Society 394, 1037–1044. doi:10. B.C.B., Crispim, A., Dias-Oliveira, A., Gomes-Júnior, A.R., Lainey, V., 1111/j.1365-2966.2009.14398.x. Miranda, J.O., Moura, T.S., Ribeiro, F.K., Santana, T., Santos-Filho, Emel’yanov, N.V., 2017. Current problems of dynamics of moons of planets S., Trabuco, L.L., Winter, O.C., Yamashita, T.A.R., 2019. APPROX and binary asteroids based on observations. Solar System Research 51, - mutual approximations between the Galilean moons: the 2016-2018 20–37. doi:10.1134/S0038094617010014. observational campaign. Monthly Notices of the Royal Astronomical Emelyanov, N.V., Gilbert, R., 2006. Astrometric results of observations Society 482, 5190–5200. doi:10.1093/mnras/sty3040, arXiv:1811.02913. of mutual occultations and eclipses of the Galilean satellites of Jupiter Oren, M., Nayar, S.K., 1994. Generalization of lambert’s reflectance model, in 2003. Astronomy and Astrophysics 453, 1141–1149. doi:10.1051/ in: Proceedings of the 21st Annual Conference on Computer Graph- 0004-6361:20064810. ics and Interactive Techniques, ACM, New York, NY, USA. pp. 239– Hapke, B., 1981. Bidirectional reflectance spectroscopy. 1. Theory. Journal 246. URL: http://doi.acm.org/10.1145/192161.192213, doi:10.1145/ of Geophysical Research 86, 4571–4586. 192161.192213. Hapke, B., 1984. Bidirectional reflectance spectroscopy 3. Correction for Peng, Q.Y., He, H.F., Lainey, V., Vienne, A., 2012. Precise CCD positions macroscopic roughness. Icarus 59, 41–59. doi:10.1016/0019-1035(84) of Galilean satellite-pairs. Monthly Notices of the Royal Astronomical 90054-X. Society 419, 1977–1982. doi:10.1111/j.1365-2966.2011.19852.x. Hapke, B., 1986. Bidirectional reflectance spectroscopy 4. The extinction Robert, V., Saquet, E., Colas, F., Arlot, J.E., 2017. CCD astrometric obser- coefficient and the opposition effect. Icarus 67, 264–280. doi:10.1016/ vations of Amalthea and Thebe in the Gaia era. Monthly Notices of the 0019-1035(86)90108-9. Royal Astronomical Society 467, 694–698. doi:10.1093/mnras/stx123. Hapke, B., 2002. Bidirectional Reflectance Spectroscopy. 5. The Coherent Saquet, E., Emelyanov, N., Colas, F., Arlot, J.E., Robert, V., Christophe, B., Backscatter Opposition Effect and Anisotropic Scattering. Icarus 157, Dechambre, O., 2016. Eclipses of the inner satellites of Jupiter observed 523–534. doi:10.1006/icar.2002.6853. in 2015. Astronomy and Astrophysics 591, A42. doi:10.1051/0004-6361/ Hapke, B., 2008. Bidirectional reflectance spectroscopy. 6. Effects of poros- 201628246, arXiv:1605.06935. ity. Icarus 195, 918–926. doi:10.1016/j.icarus.2008.01.003. Saquet, E., Emelyanov, N., Robert, V., Arlot, J.E., Anbazhagan, P., Bail- Hapke, B., 2012. Bidirectional reflectance spectroscopy 7. The single lié, K., Bardecker, J., Berezhnoy, A.A., Bretton, M., Campos, F., Ca- B. Morgado, R. Vieira-Martins, M. Assafin et al.: Preprint submitted to Elsevier Page 10 of 11 Mutual Phenomena between the Jovian moons pannoli, L., Carry, B., Castet, M., Charbonnier, Y., Chernikov, M.M., Christou, A., Colas, F., Coliac, J.F., Dangl, G., Dechambre, O., Del- croix, M., Dias-Oliveira, A., Drillaud, C., Duchemin, Y., Dunford, R., Dupouy, P., Ellington, C., Fabre, P., Filippov, V.A., Finnegan, J., Foglia, Table 2 S., Font, D., Gaillard, B., Galli, G., Garlitz, J., Gasmi, A., Gaspar, H.S., Mutual events and observation conditions. Gault, D., Gazeas, K., George, T., Gorda, S.Y., Gorshanov, D.L., Gual- doni, C., Guhl, K., Halir, K., Hanna, W., Henry, X., Herald, D., Houdin, Date Event Obs. i z Ratio of o o G., Ito, Y., Izmailov, I.S., Jacobsen, J., Jones, A., Kamoun, S., Karda- yy-mm-dd S xS ( ) ( ) Albedo 1 2 sis, E., Karimov, A.M., Khovritchev, M.Y., Kulikova, A.M., Laborde, 14-11-02 4o1 OPD 10.60 73.49 3.56 , 0.03 J., Lainey, V., Lavayssiere, M., Le Guen, P., Leroy, A., Loader, B., 14-11-19 4o2 OPD 10.69 43.44 3.59 , 0.02 Lopez, O.C., Lyashenko, A.Y., Lyssenko, P.G., Machado, D.I., Maig- 14-12-20 2o1 FOZ 8.63 50.72 0.97 , 0.03 urova, N., Manek, J., Marchini, A., Midavaine, T., Montier, J., Mor- 14-12-21 4e1 FOZ 8.52 77.19 – gado, B.E., Naumov, K.N., Nedelcu, A., Newman, J., Ohlert, J.M., Ok- 14-12-21 3o1 FOZ 8.52 65.78 1.65 , 0.19 sanen, A., Pavlov, H., Petrescu, E., Pomazan, A., Popescu, M., Pratt, A., Raskhozhev, V.N., Resch, J.M., Robilliard, D., Roschina, E., Rothen- 14-12-24 2e3 FOZ 8.13 40.85 – berg, E., Rottenborn, M., Rusov, S.A., Saby, F., Saya, L.F., Selvaku- 15-01-21 2e1 FOZ 3.45 46.09 – mar, G., Signoret, F., Slesarenko, V.Y., Sokov, E.N., Soldateschi, J., 15-02-02 3o2 FOZ 0.95 67.15 1.53 , 0.09 Sonka, A., Soulie, G., Talbot, J., Tejfel, V.G., Thuillot, W., Timerson, 15-02-22 2o1 OPD 3.20 41.50 0.99 , 0.04 B., Toma, R., Torsellini, S., Trabuco, L.L., Traverse, P., Tsamis, V., Un- 15-02-22 2e1 OPD 3.20 43.44 – win, M., Abbeel, F.V.D., Vand enbruaene, H., Vasundhara, R., Velikod- 15-03-01 2o1 FOZ 4.57 57.20 1.01 , 0.06 sky, Y.I., Vienne, A., Vilar, J., Vugnon, J.M., Wuensche, N., Zeleny, P., 15-03-01 2e1 FOZ 4.58 65.89 – 2018. The PHEMU15 catalogue and astrometric results of the Jupiter’s 15-03-03 3o1 OPD 4.94 62.37 1.56 , 0.03 Galilean satellite mutual occultation and eclipse observations made in 15-03-03 3o1 FOZ 4.94 58.30 1.56 , 0.03 2014-2015. Monthly Notices of the Royal Astronomical Society 474, 15-03-06 1e2 OPD 5.46 43.08 – 4730–4739. doi:10.1093/mnras/stx2957. Thomas, P.C., Burns, J.A., Rossier, L., Simonelli, D., Veverka, J., Chap- 15-03-09 3e2 OPD 6.14 44.04 – man, C.R., Klaasen, K., Johnson, T.V., Belton, M.J.S., Galileo Solid 15-03-13 1e2 FOZ 6.79 48.89 – State Imaging Team, 1998. The Small Inner Satellites of Jupiter. Icarus 15-03-13 1e3 FOZ 6.66 60.76 – 135, 360–371. doi:10.1006/icar.1998.5976. 15-03-16 4o2 FOZ 7.11 48.84 3.58 , 0.03 Vachier, F., Arlot, J.E., Thuillot, W., 2002. Mutual phenomena involving J5 15-03-17 3e2 FOZ 7.27 58.36 – Amalthea in 2002-2003. Astronomy and Astrophysics 394, L19–L21. 15-03-18 2e1 GOA 7.54 41.46 – doi:10.1051/0004-6361:20021329. 15-03-24 3o4 OPD 8.22 45.53 0.41 , 0.03 Vasundhara, R., Selvakumar, G., Anbazhagan, P., 2017. Analysis of mu- 15-03-24 3o4 FOZ 8.22 46.49 0.41 , 0.03 tual events of Galilean satellites observed from VBO during 2014-2015. 15-03-25 2o1 FOZ 8.47 46.28 0.99 , 0.04 Monthly Notices of the Royal Astronomical Society 468, 501–508. 15-03-25 2o1 GOA 8.47 42.75 0.99 , 0.04 doi:10.1093/mnras/stx437, arXiv:1704.03518. Veiga, C.H., Vieira Martins, R., 2005. CCD astrometric observations of 15-03-26 2e1 OPD 8.48 52.32 – Amalthea and Thebe. Astronomy and Astrophysics 437, 1147–1150. 15-04-02 2o1 OPD 9.25 62.96 1.01 , 0.07 doi:10.1051/0004-6361:20042387. 15-04-02 2o1 FOZ 9.25 59.06 1.01 , 0.07 Zhang, X.L., Han, X.L., Arlot, J.E., 2019. Mutual events between Galilean 15-04-02 2e1 OPD 9.26 82.45 – satellites observed with SARA 0.9 m and 0.6 m telescopes during 2014- 15-04-03 1o3 FOZ 9.43 46.38 0.75 , 0.13 2015. Monthly Notices of the Royal Astronomical Society 483, 4518– 15-04-06 1e2 FOZ 9.70 46.72 – 4524. doi:10.1093/mnras/sty3030. 15-04-12 2e3 FOZ 10.09 65.99 – 15-04-14 1e2 OPD 10.22 69.43 – 15-04-17 4o1 OPD 10.43 54.20 3.90 , 0.08 15-04-17 4o1 GOA 10.43 55.82 3.90 , 0.08 15-04-18 4o3 OPD 10.48 44.23 2.32 , 0.05 15-04-18 1o3 OPD 10.44 72.62 0.64 , 0.11 15-04-25 1o3 OPD 10.73 58.34 0.69 , 0.07 15-04-25 1o3 FOZ 10.73 55.16 0.63 , 0.07 15-04-26 2o1 OPD 10.75 43.56 1.04 , 0.15 15-04-29 3o1 OPD 10.80 67.89 1.60 , 0.06 15-04-29 3o1 GOA 10.80 70.74 1.60 , 0.06 15-05-03 2o1 OPD 10.85 61.83 0.96 , 0.07 15-05-05 3o2 FOZ 10.85 47.24 1.61 , 0.06 15-05-13 3o2 OPD 10.77 86.09 1.48 , 0.07 15-06-04 2o1 FOZ 9.72 56.26 1.01 , 0.05 15-06-18 3o1 GOA 8.59 59.64 1.58 , 0.04 15-03-02 3e5 1.60 3.17 50.90 – Note: The solar phase angle, zenith distance and ratio of albedo in the sense S _S for each event. 2 1 B. Morgado, R. Vieira-Martins, M. Assafin et al.: Preprint submitted to Elsevier Page 11 of 11 Mutual Phenomena between the Jovian moons Table 3 Updated results of the 2009 mutual phenomena campaign observed in Brazil. Date Event Obs. t , t (UTC) t s , s s v , v v rms N 0 0 0 0 0 0 0 0 0 yy-mm-dd S xS hh:mm:ss.s , s.s mas mas mas mas/s mas/s 1 2 09-04-27 3o1 OPD 06:42:53.5 , 0.6 -00.8 121.0 , 20.0 -09.2 6.26 , 0.02 -0.01 0.010 0801 1.005 09-05-09 2o1 OPD 07:21:54.9 , 0.7 +01.5 538.0 , 14.3 -02.3 7.50 , 0.03 -0.03 0.010 0580 1.007 09-05-21 1o2 OPD 05:29:44.3 , 0.7 -00.2 033.6 , 21.1 +26.5 6.36 , 0.06 -0.00 0.014 0869 1.005 09-05-28 1o2 OPD 07:44:17.6 , 0.4 +00.0 219.1 , 17.5 +47.0 6.26 , 0.03 -0.01 0.007 0801 1.005 09-06-10 3e4 OPD 07:29:37.4 , 1.2 +14.2 315.0 , 17.9 +22.8 3.96 , 0.03 +0.01 0.009 1285 1.003 09-06-16 3e1 OPD 08:45:10.5 , 1.2 -25.6 928.4 , 13.3 +19.8 1.52 , 0.01 -0.00 0.008 2697 1.001 09-06-19 4e2 OPD 05:11:36.3 , 0.2 +25.1 489.1 , 13.1 +02.7 5.06 , 0.01 +0.02 0.007 0901 1.004 09-06-19 4e1 OPD 08:32:49.4 , 1.0 +15.6 930.1 , 13.5 +05.9 5.29 , 0.05 +0.02 0.006 1201 1.003 09-06-20 4e1 OPD 05:09:45.8 , 1.3 -22.6 530.5 , 12.8 +01.5 1.26 , 0.01 -0.01 0.006 2326 1.002 09-06-20 4e1 OPD 09:37:28.5 , 1.3 +31.6 412.3 , 18.5 -08.4 1.39 , 0.01 +0.00 0.013 1704 1.002 09-06-22 1o2 OPD 03:27:54.8 , 1.6 +12.0 576.1 , 14.0 +23.7 5.78 , 0.03 -0.03 0.010 0905 1.004 09-06-29 1o2 OPD 05:38:28.7 , 0.8 +00.5 606.4 , 13.3 +10.0 5.51 , 0.02 -0.01 0.004 0801 1.005 09-07-04 1e3 OPD 06:25:13.5 , 0.5 -12.7 395.5 , 16.3 -19.8 7.27 , 0.04 +0.02 0.008 1641 1.002 09-07-06 1e2 OPD 06:17:16.1 , 1.8 +06.2 718.5 , 14.2 -12.2 4.79 , 0.12 +0.02 0.006 2001 1.002 09-07-06 1o2 OPD 07:48:34.5 , 0.5 +17.4 603.5 , 12.9 -05.5 5.34 , 0.01 -0.03 0.005 1004 1.004 09-07-08 3e1 OPD 08:31:14.5 , 0.2 +30.1 223.9 , 14.7 -07.2 5.96 , 0.01 +0.02 0.007 1758 1.002 09-07-13 1e2 OPD 08:38:46.3 , 1.7 +10.3 623.4 , 15.3 -08.0 4.46 , 0.11 +0.02 0.012 2001 1.002 09-08-07 1e2 OPD 05:14:54.9 , 1.3 +12.4 444.6 , 15.1 +34.8 3.09 , 0.03 +0.01 0.021 1775 1.002 09-08-07 1o2 OPD 05:37:48.4 , 0.8 -08.8 283.1 , 20.3 -18.7 3.77 , 0.01 -0.01 0.008 1664 1.002 09-08-12 3o2 OPD 02:10:59.1 , 3.9 +01.0 1059.1 , 13.1 -24.7 2.71 , 0.02 -0.01 0.004 1296 1.003 09-08-22 1o2 OPD 04:07:54.9 , 2.2 +16.1 674.5 , 12.9 +11.2 1.80 , 0.01 +0.00 0.004 2454 1.002 09-09-16 1o2 OPD 00:46:04.4 , 0.7 -12.6 580.4 , 15.0 +06.5 3.71 , 0.02 +0.00 0.011 0976 1.004 09-09-16 1e2 OPD 02:15:11.0 , 0.4 +04.0 172.4 , 12.9 +00.8 3.53 , 0.02 +0.02 0.007 1095 1.004 09-10-24 3o2 OPD 00:35:33.9 , 1.5 +00.8 629.1 , 13.6 -48.8 4.17 , 0.02 -0.01 0.008 1032 1.004 09-10-25 1o2 OPD 01:21:30.8 , 3.9 +01.0 580.6 , 17.5 +08.5 5.35 , 0.07 -0.01 0.014 0252 1.016 Note: The results for the mutual phenomena campaign of 2009. t stand for the UTC central instant, s is the impact 0 0 parameter and v is the apparent relative velocity in the sky plane. Also contains the uncertainty in each parameter (t , s 0 0 0 and v ) and the difference between the fitted ones and the ones expected from the ephemeris jup310 and DE435 (t , s 0 0 0 and v ). In the last columns, we have the rms between the observed light fluxes and the fitted ones, the number of images utilised (N) and the normalised  of our fit. B. Morgado, R. Vieira-Martins, M. Assafin et al.: Preprint submitted to Elsevier Page 12 of 11 Mutual Phenomena between the Jovian moons Table 4 Comparison of the updated results for the 2009 mutual events with Arlot et al. (2014a), Dias-Oliveira et al. (2013) e Morgado et al. (2016). Date Event Obs. Central instant (t ) Impact parameter (s ) 0 0 yy-mm-dd S xS [2] - [1] [3] - [1] [4] - [1] [2] - [1] [3] - [1] [4] - [1] 1 2 09-04-27 3o1 OPD -0.65 -7.15 – +0.15 -0.74 – 09-05-09 2o1 OPD +0.78 -2.34 -0.12 +0.23 +0.97 – 09-05-21 1o2 OPD +0.28 -4.94 – +3.30 +1.05 – 09-05-28 1o2 OPD +0.98 -9.24 +0.57 -0.00 -2.36 – 09-06-10 3e4 OPD +0.13 +1.54 – +0.13 +1.56 – 09-06-16 3e1 OPD +0.99 -1.99 – -0.12 +8.81 – 09-06-19 4e2 OPD -1.87 +0.99 – +0.03 +4.45 – 09-06-19 4e1 OPD -0.83 +0.62 – +0.04 +8.26 – 09-06-20 4e1 OPD -0.70 +1.14 – +0.08 -4.27 – 09-06-20 4e1 OPD -0.03 -2.85 – -0.13 +9.77 – 09-06-22 1o2 OPD +0.27 -2.18 -0.07 +0.16 -0.79 – 09-06-29 1o2 OPD +0.31 -5.30 – +0.07 -0.86 – 09-07-04 1e3 OPD -0.03 +1.17 – -0.26 +4.57 – 09-07-06 1o2 OPD -0.13 +0.25 -0.30 +0.14 +5.74 – 09-07-06 1e2 OPD +0.50 -4.89 – +0.09 -0.01 – 09-07-08 3e1 OPD -0.61 +1.33 – -0.04 +1.84 – 09-07-13 1e2 OPD -0.31 +0.22 – +0.18 +8.36 – 09-08-07 1o2 OPD -0.22 +0.29 +0.44 +0.32 +7.38 – 09-08-07 1e2 OPD -0.34 -1.71 – -0.09 -0.17 – 09-08-12 3o2 OPD +0.83 +0.60 – -0.12 +0.50 – 09-08-22 1o2 OPD +0.31 +0.29 – +0.07 -0.99 – 09-09-16 1e2 OPD -0.64 +3.52 – -0.23 -1.43 – 09-09-16 1o2 OPD +0.51 -0.43 – -0.22 +1.02 – 09-10-24 3o2 OPD -0.25 +14.15 – -0.26 +2.78 – 09-10-25 1o2 OPD +0.62 +1.56 – +0.23 -0.98 – Mean -0.00 -0.61 +0.10 +0.15 +2.18 – Standard deviation 0.67 4.34 0.34 0.67 3.91 – Note: [1] This project, [2] Arlot et al. (2014a), [3] Dias-Oliveira et al. (2013), [4] Morgado et al. (2016). Comparison between the different reduction process divided by the uncertainty of each parameter. B. Morgado, R. Vieira-Martins, M. Assafin et al.: Preprint submitted to Elsevier Page 13 of 11 Mutual Phenomena between the Jovian moons Table 5 Results of the 2014-2015 mutual phenomena campaign observed in Brazil. Date Event Obs. t , t (UTC) t s , s s v , v v rms N 0 0 0 0 0 0 0 0 0 yy-mm-dd S xS hh:mm:ss.s , s.s mas mas mas mas/s mas/s 1 2 14-11-02 4o1 OPD 06:02:14.8 , 1.9 -07.2 288.0 , 18.6 -05.0 1.83 , 0.03 +0.00 0.041 230 1.018 14-11-19 4o2 OPD 07:37:43.9 , 0.2 +00.7 301.3 , 6.0 +00.2 5.04 , 0.03 +0.01 0.010 194 1.021 14-12-20 2o1 FOZ 05:41:38.7 , 0.8 -03.0 162.4 , 8.5 -17.5 1.62 , 0.02 -0.00 0.008 376 1.011 14-12-21 4e1 FOZ 03:22:39.2 , 1.5 +00.1 253.3 , 9.0 +26.2 2.58 , 0.02 +0.00 0.018 281 1.014 14-12-21 3o1 FOZ 04:17:20.7 , 3.2 -00.8 455.1 , 16.1 +14.0 3.66 , 0.03 -0.01 0.151 244 1.017 14-12-24 2e3 FOZ 06:35:06.5 , 1.1 -00.8 255.2 , 8.8 +09.5 1.93 , 0.02 +0.00 0.014 359 1.011 15-01-21 2e1 FOZ 03:52:29.3 , 1.0 +00.8 457.5 , 3.8 +12.8 3.08 , 0.04 +0.00 0.008 141 1.029 15-02-02 3o2 FOZ 07:56:42.3 , 1.9 -02.3 85.0 , 27.4 +08.0 7.78 , 0.05 -0.01 0.022 110 1.038 15-02-22 2o1 OPD 02:07:51.7 , 0.2 +06.5 125.0 , 12.9 -00.8 5.55 , 0.03 -0.02 0.007 201 1.020 15-02-22 2e1 OPD 02:45:11.0 , 0.3 +04.1 17.0 , 9.6 +07.5 4.62 , 0.03 +0.02 0.013 186 1.022 15-03-01 2o1 FOZ 04:12:38.3 , 0.7 -00.7 17.6 , 10.5 +16.8 5.73 , 0.04 -0.01 0.012 131 1.031 15-03-01 2e1 FOZ 05:05:06.7 , 2.4 +06.6 94.6 , 33.3 -06.8 4.94 , 0.04 +0.00 0.059 132 1.031 15-03-03 3o1 OPD 04:08:16.3 , 0.5 +00.3 95.0 , 14.4 +32.3 8.53 , 0.05 -0.03 0.010 100 1.042 15-03-03 3o1 FOZ 04:08:15.5 , 0.6 -07.0 86.8 , 10.8 +23.7 8.50 , 0.05 -0.02 0.010 107 1.039 15-03-06 1e2 OPD 01:16:16.6 , 0.4 +09.4 570.9 , 4.0 +06.2 7.50 , 0.04 +0.04 0.005 124 1.033 15-03-09 3e2 OPD 23:39:32.6 , 0.3 +08.1 67.7 , 7.2 +02.6 5.87 , 0.03 +0.03 0.012 213 1.019 15-03-13 1e2 FOZ 03:29:09.9 , 1.1 +10.8 445.2 , 20.0 +03.7 7.42 , 0.03 +0.03 0.020 198 1.021 15-03-13 1e3 FOZ 23:29:44.5 , 1.3 +04.0 236.3 , 4.4 +08.5 1.56 , 0.02 -0.00 0.015 408 1.010 15-03-16 4o2 FOZ 01:38:58.1 , 0.3 -00.5 373.0 , 6.3 +02.0 3.59 , 0.03 -0.00 0.013 246 1.017 15-03-17 3e2 FOZ 02:53:15.2 , 0.2 +08.2 225.4 , 14.4 -03.1 5.77 , 0.03 +0.02 0.012 234 1.017 15-03-18 2e1 GOA 22:50:43.6 , 1.9 +06.1 398.2 , 25.2 -02.7 5.55 , 0.02 +0.02 0.058 288 1.014 15-03-24 3o4 OPD 00:14:41.4 , 0.9 -38.1 499.0 , 13.9 -16.0 5.39 , 0.03 -0.01 0.006 182 1.022 15-03-24 3o4 FOZ 00:14:41.6 , 0.8 -37.5 519.1 , 9.9 +03.5 5.38 , 0.03 -0.01 0.008 207 1.020 15-03-25 2o1 FOZ 23:35:01.6 , 1.1 +00.8 400.2 , 13.2 +01.7 6.34 , 0.06 -0.02 0.011 083 1.051 15-03-25 2o1 GOA 23:35:01.3 , 0.7 -01.5 398.9 , 6.6 +00.7 6.40 , 0.05 -0.04 0.006 092 1.045 15-03-26 2e1 OPD 01:07:48.0 , 5.3 +10.3 516.7 , 44.7 -12.2 5.78 , 0.04 +0.02 0.048 151 1.027 15-04-02 2o1 OPD 01:43:55.9 , 0.7 +00.0 479.5 , 5.7 -00.5 6.46 , 0.06 -0.02 0.005 082 1.051 15-04-02 2o1 FOZ 01:43:55.6 , 1.3 -02.3 482.3 , 15.0 +02.5 6.50 , 0.05 -0.03 0.011 104 1.040 15-04-02 2e1 OPD 03:24:16.8 , 3.1 +07.8 658.7 , 19.1 -02.1 6.02 , 0.04 +0.01 0.018 129 1.032 15-04-03 1o3 FOZ 22:58:19.0 , 5.9 +05.7 737.3 , 6.0 -05.5 1.23 , 0.02 +0.01 0.021 383 1.011 15-04-06 1e2 FOZ 23:16:40.4 , 0.2 +11.3 54.2 , 10.8 +04.2 6.96 , 0.04 +0.02 0.007 152 1.027 15-04-12 2e3 FOZ 01:46:02.2 , 1.0 +13.6 142.6 , 8.3 +09.6 4.96 , 0.04 +0.00 0.013 121 1.034 15-04-14 1e2 OPD 01:30:58.3 , 0.4 +10.4 49.4 , 8.5 -01.5 6.83 , 0.05 +0.01 0.024 093 1.045 15-04-17 4o1 OPD 23:47:06.9 , 0.9 -01.3 711.5 , 4.5 +01.8 5.05 , 0.05 -0.01 0.007 101 1.041 15-04-17 4o1 GOA 23:47:06.9 , 0.9 -01.0 712.1 , 5.1 +02.7 5.06 , 0.04 -0.02 0.010 128 1.032 15-04-18 4o3 OPD 01:32:30.4 , 0.9 -02.3 69.5 , 15.1 +06.8 5.02 , 0.03 -0.00 0.021 161 1.025 15-04-18 1o3 OPD 20:54:45.6 , 3.7 -02.0 699.4 , 33.0 +03.3 5.50 , 0.04 +0.01 0.047 119 1.035 15-04-25 1o3 OPD 23:45:28.1 , 1.3 +00.5 679.2 , 7.5 -02.3 6.05 , 0.05 -0.02 0.007 100 1.042 15-04-25 1o3 FOZ 23:45:26.7 , 3.3 -08.0 685.0 , 28.8 +03.8 5.95 , 0.04 +0.01 0.014 133 1.031 15-04-26 2o1 OPD 21:25:00.0 , 3.4 +01.8 584.6 , 27.6 +00.5 7.08 , 0.08 -0.07 0.029 064 1.067 15-04-29 3o1 OPD 00:29:06.9 , 1.2 -15.6 661.3 , 19.2 -04.0 6.92 , 0.06 -0.01 0.009 087 1.048 15-04-29 3o1 GOA 00:29:07.6 , 2.4 -10.5 663.8 , 22.2 -01.7 7.00 , 0.05 -0.04 0.026 095 1.044 15-05-03 2o1 OPD 23:39:19.6 , 2.2 -08.1 571.1 , 25.6 +11.5 6.83 , 0.08 -0.04 0.033 063 1.068 15-05-05 3o2 FOZ 21:54:22.0 , 2.4 +00.3 780.5 , 21.9 -01.7 5.30 , 0.04 -0.01 0.010 117 1.035 15-05-13 3o2 OPD 01:13:50.0 , 3.0 +06.5 593.2 , 27.1 +03.5 5.04 , 0.04 -0.01 0.038 137 1.030 15-06-04 2o1 FOZ 21:55:27.9 , 0.2 -01.8 160.8 , 6.3 -06.8 6.89 , 0.05 -0.02 0.009 094 1.044 15-06-18 3o1 GOA 21:01:51.0 , 1.1 +01.0 237.6 , 18.6 +10.7 3.82 , 0.03 -0.02 0.017 215 1.019 15-03-02 3e5 1.60 23:17:06.0 , 2.3 -22.5 391.3 , 76.3 +20.9 8.44 , 0.53 -0.09 0.101 060 1.080 Note: Similar as the note in the Table 3 for the mutual phenomena campaign of the 2014-2015. B. Morgado, R. Vieira-Martins, M. Assafin et al.: Preprint submitted to Elsevier Page 14 of 11 http://www.deepdyve.com/assets/images/DeepDyve-Logo-lg.png Astrophysics arXiv (Cornell University)

Loading next page...
 
/lp/arxiv-cornell-university/the-2014-2015-brazilian-mutual-phenomena-campaign-for-the-jovian-Pm2XyHz4LW

References (46)

ISSN
0032-0633
eISSN
ARCH-3330
DOI
10.1016/j.pss.2019.104736
Publisher site
See Article on Publisher Site

Abstract

The 2014-2015 Brazilian Mutual Phenomena campaign for the Jovian satellites and improved results for the 2009 events. a,b,c,< a,b,c b,c a,b,d e,f B. Morgado , R. Vieira-Martins , M. Assafin , A. Dias-Oliveira , D. I. Machado , J. I. a,b g h h i,a,b a,b B. Camargo , M. Malacarne , R. Sfair , O. C. Winter , F. Braga-Ribas , G. Benedetti-Rossi , h h h a,b,c,h g L. A. Boldrin , B. C. B. Camargo , H. S. Gaspar , A. R. Gomes-Júnior , J. O. Miranda , T. h f de Santana and L. L. Trabuco Observatório Nacional/MCTIC, R. General José Cristino 77, Rio de Janeiro, RJ 20.921-400, Brazil Laboratório Interinstitucional de e-Astronomia - LIneA and INCT do e-Universo, Rua Gal. José Cristino 77, Rio de Janeiro, RJ 20921-400, Brazil Observatório do Valongo/UFRJ, Ladeira Pedro Antonio 43, Rio de Janeiro, RJ 20080-090, Brazil Escola SESC de Ensino Médio, Avenida Ayrton Senna, 5677, Rio de Janeiro - RJ, 22775-004, Brazil Universidade Estadual do Oeste do Paraná (Unioeste), Avenida Tarquínio Joslin dos Santos 1300, Foz do Iguaçu, PR 85870-650, Brazil Polo Astronômico Casimiro Montenegro Filho/FPTI-BR, Avenida Tancredo Neves 6731, Foz do Iguaçu, PR 85867-900, Brazil Universidade Federal do Espírito Santo, Av. Fernando Ferrari 514, Vitória, ES 29075-910, Brasil UNESP - São Paulo State University, Grupo de Dinâmica Orbital e Planetologia, CEP 12516-410, Guaratinguetá, SP 12516-410, Brazil Federal University of Technology - Paraná (UTFPR/DAFIS), Av. Sete de Setembro, 3165, Curitiba, PR 80230-901, Brazil A R T I C L E I N F O A B S T R A C T Keywords: Progress in astrometry and orbital modelling of planetary moons in the last decade enabled better Methods: Data analysis determinations of their orbits. These studies need accurate positions spread over extended periods. Astrometry We present the results of the 2014-2015 Brazilian campaign for 40 mutual events from 47 observed Planets and satellites: individual: Io, light curves by the Galilean satellites plus one eclipse of Amalthea by Ganymede. We also reanalysed Europa, Ganymede, Callisto, Amalthea. and updated results for 25 mutual events observed in the 2009 campaign. All telescopes were equipped with narrow-band filters centred at 889 nm with a width of 15 nm to eliminate the scattered light from Jupiter. The albedos’ ratio was determined using images before and after each event. We simulated images of moons, umbra, and penumbra in the sky plane, and integrated their fluxes to compute albedos, simulate light curves and fit them to the observed ones using a chi-square fitting procedure. For that, we used the complete version of the Oren-Nayer reflectance model. The relative satellite positions mean uncertainty was 11.2 mas (í35 km) and 10.1 mas (í31 km) for the 2014-2015 and 2009 campaigns respectively. The simulated and observed ASCII light curve files are freely available in electronic form at the Natural Satellites DataBase (NSDB). The 40/25 mutual events from our 2014-2015/2009 campaigns represent a significant contribu- tion of 17%/15% in comparison with the PHEMU campaigns lead by the IMCCE. Besides that, our th result for the eclipse of Amalthea is only the 4 such measurement ever published after the three ones observed by the 2014-2015 international PHEMU campaign. Our results are suitable for new orbital/ephemeris determinations for the Galilean moons and Amalthea. 2009, 2017). 1. Introduction The uncertainty of the positions obtained from mutual Mutual phenomena between natural satellites – occul- phenomena is usually smaller than the ones obtained by other tations and eclipses – have been successfully used to im- methods. For instance, classical CCD astrometry achieves prove the orbital studies of these moons. For the Galilean uncertainties around 100 mas (í 300 km) (Kiseleva et al., satellites, they have been systematically observed since 1976 2008). For satellite-pair distances, the uncertainties are at (Aksnes and Franklin, 1976). These phenomena occur as the the 30 mas level (í 90 km) (Peng et al., 2012). Mutual ap- Earth and the Sun cross the orbital plane of the satellites. For proximations, based in the same geometrical configuration Jupiter, they happen every six years. of mutual occultations, achieve uncertainties at the 10 mas The photometry of these events offers a reliable source of level (í 30 km) (Morgado et al., 2016, 2019). very precise relative positions between two satellites. They In this paper, we present results for 47 light curves, 31 often achieve uncertainties bellow 5 mas (í 15 km) (Emelyanov, occultations and 16 eclipses, representing 40 mutual events 2009; Dias-Oliveira et al., 2013; Arlot et al., 2014a; Saquet between the Galilean moons observed by three stations in et al., 2018). These relative positions can constrain the or- Brazil, during the 2014-2015 mutual phenomena campaign. bital studies of these moons and give us hints about their We also present one event, an eclipse involving the inner structure and formation processes (Lainey et al., 2004b,a, satellite Amalthea (J5). We also used our improved methods to re-analyse 25 light curves, 13 occultations and 12 eclipses, Based in part on observations made at the Laboratório Nacional de Astrofísica (LNA), Itajubá-MG, Brazil. of 25 mutual phenomena observed by our group during the Corresponding author 2009 mutual phenomena campaign. We compared the new Morgado.fis@gmail.com (B. Morgado) results with those by Dias-Oliveira et al. (2013), Arlot et al. ORCID(s): 0000-0003-0088-1808 (B. Morgado) B. Morgado, R. Vieira-Martins, M. Assafin et al.: Preprint submitted to Elsevier Page 1 of 11 arXiv:1909.05238v1 [astro-ph.EP] 11 Sep 2019 Mutual Phenomena between the Jovian moons (2014a) and Morgado et al. (2016). In Section 2 we detail the observational campaigns. In Section 3 we present the photometry used to produce the ob- served light curves and describe the new, improved light- curve fitting procedures developed and used in this work. Section 4.1 contains new results from the re-analysis of 25 mutual events observed in 2009 and the comparison with the older results. In Section 4.2, we present the results for the 47 light curves involving 40 mutual events observed from Brazil during the 2014-2015 campaign. In Section 4.3, we present the result for the eclipse involving Amalthea. Our conclusions are set on Section 5. Figure 1: Image of Jupiter, Io (J1), Europa (J2) and Callisto (J4) obtained with the 0.6 m diameter Zeiss telescope of the 2. Mutual phenomena campaign details Observatório Pico dos Dias, equipped with a methane filter Every six years, during Jupiter equinox, we can observe on 2014 November 19. The planet and the satellites present mutual occultations and eclipses between Jupiter’s regular about the same brightness due to the use of the narrow-band satellites. The results presented here come from the col- filter, centred at  = 889 nm with 15 nm width. laboration between five Brazilian institutes. The prediction of these events was provided by the Institut de Mécanique Céleste et de Calcul des Ephémérides (IMCCE) (Arlot et al., this filter because in this spectral region, Jupiter’s albedo 2014a,b). drops to 0.05 due to the absorption in the upper atmosphere The 2009 mutual phenomena campaign was the first large (Karkoschka, 1994, 1998). Figure 1 shows an example of an attempt of the kind carried out in Brazil for the Galilean image obtained with the 0.6 m Zeiss telescope from OPD. moons. Observations and instruments are described in de- This filter has been successfully used in the 2009 mutual tail in Dias-Oliveira et al. (2013). We re-analysed 25 mutual phenomena campaign (Dias-Oliveira et al., 2013) and the events encompassing 25 light curves, 13 from occultations mutual approximation campaigns started in 2014 (Morgado and 12 from eclipses, and obtained new results for this cam- et al., 2016, 2019). paign. Discrepancies between the results obtained by Dias- Oliveira et al. (2013) and by Arlot et al. (2014a) motivated 3. Light curve analysis this re-analysis, see more details in Section 4.1. The last Brazilian mutual phenomena campaign of 2014- In mutual phenomena, one can determine relative posi- 2015 obtained data from three telescopes spread on the South tions between the satellites through the analysis of the events’ and South-East of Brazil, with apertures ranging between light curves. In our procedure, we simulate theoretical light 28 and 60 cm. We obtained 47 light curves, 31 for occulta- curves and use them to fit the observed ones. tions and 16 for eclipses, from 40 events observed between The parameters of interest are: (i) the impact parameter November 2014 and June 2015. In Table 1 we present the (s ), the smallest apparent angular distance in the sky plane stations, observers and instrumental details of each station. between both satellite’s centres in the case of occultations It also contains the number of light curves obtained by each or between the eclipsed satellite centre and the centre of the observatory. Moreover, we added the Minor Planet Center eclipsing shadow in the sky plane for eclipses, both cases in (MPC) observatory code of the station (XXX for the station a topocentric frame; (ii) the central instant (t ), the instant of without a code). time that this smallest distance occurs; and (iii) the apparent In Table 2, we list the observational details of each ob- relative velocity (v ) between both satellites in the sky plane. served event. It contains the date of the event and the satel- In the supplementary material, we also provide the inter- lites’ pairs in the form "S oS " for occultation and "S eS " 1 2 1 2 satellite tangential coordinates (X, Y ). For occultations, these for eclipses, where 1 stands for Io, 2 for Europa, 3 for Ganymede coordinates between both satellitesâĂŹ centres are in a topocen- and 4 for Callisto. We furnish the sites involved in each ob- tric frame. For mutual eclipses, these coordinates are in servation (using the alias defined in Table 1). For each event, a topocentric frame and the mean difference between the we give the solar phase angle (i) and the zenith distance (z), eclipsed satellite centre and the centre of the eclipsing satel- both in degrees. In the last column, we list the instrumental liteâĂŹs shadow in the sky plane. albedos’ ratio of the involved satellites (and its uncertainty), determined by using images before and after the event. This 3.1. Obtaining the observed light-curves albedos’ ratio is only needed for occultations. Firstly all images were corrected by Bias and Flat-Field In all observations we used a narrow band filter cen- using standard procedures with the Image Reduction and tred at 889 nm with a width of 15 nm. This bandpass is in Analysis Facility (IRAF) (Butcher and Stevens, 1981). Then, the methane absorption region of the spectrum. We chose we determined the light flux of the targets in the images by differential aperture photometry using the PRAIA package Website: http://nsdb.imcce.fr/multisat/nsszph515he.htm B. Morgado, R. Vieira-Martins, M. Assafin et al.: Preprint submitted to Elsevier Page 2 of 11 Mutual Phenomena between the Jovian moons Table 1 2014-2015 mutual phenomena campaign observer list. City/ Country Longitude Observers Telescope N of positive Station alias Latitude CCD detections MPC code Height Itajubá/ MG, Brazil 45 34’ 57" W B. Morgado 60 cm 20 OPD 22 32’ 07" S H. S. Gaspar Andor/IKon-L 874 1.864 km R. Sfair B. C. B. Camargo T. Santana L. A. Boldrin M. E. Huaman G. Benedetti-Rossi A. R. Gomes-Júnior Foz do Iguaçu/ PR, Brazil 54 35’ 37" W D. I. Machado 28 cm 22 FOZ 25 26’ 05" S L. L. Trabuco SBIG/ST-7X-ME X57 0.184 km Vitória/ ES, Brazil 40 19’ 00" W M. Malacarne 35 cm 5 GOA 20 17’ 52" S J. O. Miranda SBIG/ST-8X-ME XXX 0.026 km (Assafin et al., 2011). During an occultation, both satellites the time and the dynamical parameters using equation (1) are measured together in the same aperture, and a third satel- (Assafin et al., 2009). lite is used as calibrator. In the case of eclipses, the eclipsed satellite is measured alone in the aperture and the eclipsing 2 2 2 s.t/ = s + v .t * t / (1) satellite (or any other) is used as calibrator. The light curve 0 0 is then normalised by a polynomial fit so that the flux ratio We also need a reflectance model to take into account the outside the flux drop gets equal to 1.0, and the flux drop can phase effect and how the surface of the satellite will reflect be adequately evaluated. the sunlight. For eclipses, we further need some information During the photometry of mutual events, one must take about the Sun, such as its radius and a model to consider the care with the possibility of a parasitic flux, as pointed out Sun’s limb darkening. by Emel’yanov (2017) and Arlot et al. (2017). The origin We could do simulations with triaxial bodies with vary- of this flux is likely to be the background (mostly Jupiter’s ing albedo, but not for practical purposes, due to photometry scattered light) or the CCD detector. In our case, we attenu- limitations. Thus, satellites are considered as spheres with a ate this parasitic flux with the Methane filter and a rigorous known appropriate radius, and the albedo is uniform along calibration process. Tests showed that the parasitic flux in the surface. The relative velocity between both satellites is our images is one order of magnitude below the noise of our constant during the mutual event (a few minutes only). For observations. notation, the occulting/eclipsing is denoted Sat and Sat is 1 2 the occulted/eclipsed one. 3.2. Simulating light curves In the simulations of occultations and eclipses, we used The procedures utilised here follow the same principles the same geometric relations described in detail in Dias-Oliveira outlined in Assafin et al. (2009) and Dias-Oliveira et al. (2013). et al. (2013). However, improvements were made in almost every step, as explained in sections 3.2.1, 3.2.2 and 3.3. 3.2.1. Occultation The parameters needed in the modelling of a mutual oc- The first step in simulating a light curve of a mutual oc- cultation or eclipse can be separated into two complemen- cultation is the production of a 2D satellite apparent pro- tary types. The first refers to physical characteristics: sizes file, simulating how the body reflects the sunlight as seen by and shapes of each satellite, and the satellites’ albedo for oc- an observer on Earth. As pointed out by Vasundhara et al. cultations. The second type relates to dynamics, and are the (2017), it is essential to use a realistic intensity distribution parameters of interest: s , t and v . 0 0 0 for the satellite. Physical parameters such as radius and shape are known However, this approach demands previous knowledge about from space probes’ data. Albedos are determined from aux- the satellite surface (albedo maps) that can change with time iliary observations made before and after events with the or even for different effective wavelengths, it is important to same instrument setup. The apparent separation in the sky highlight that these maps are not know for the wavelength plane between both satellites can be written as a function of of our observations ( = 889 nm). The same applies to the B. Morgado, R. Vieira-Martins, M. Assafin et al.: Preprint submitted to Elsevier Page 3 of 11 Mutual Phenomena between the Jovian moons Hapke scattering law (Hapke, 1981; Hapke and Wells, 1981; curve was then normalised by the sums of individual fluxes Hapke, 1984, 1986, 2002, 2008, 2012) used by Emelyanov (F + F ). Figure 2 shows a simulated light curve for the 1 2 (2009), Arlot et al. (2014a) and Saquet et al. (2018), which event when Europa occulted Io on February 22, 2015. The requires unknown parameters in the wavelength band of our black dots indicate seven instants for which the respective observations. 2D profiles are displayed on each corresponding box. For We successfully solved these issues by adopting a gener- this event the albedo ratio was 0.960, as determined before alisation of Lambert scattering law given by Oren and Nayar the event. (1994). The Oren-Nayer model takes into account the direc- 3.2.2. Eclipse tion of radiance and the roughness of the surface in a natural Following Dias-Oliveira et al. (2013), we developed a way, so that the reflectance depends only of the albedo and numerical 2D mask that incorporates the two regions of the in one more parameter that very smoothly tunes a wide range eclipsing satellite’s shadow, the umbra and the penumbra. of surface roughness, and most importantly, regardless of the This mask was then applied to the 2D profile of the eclipsed wavelength. This model realistically reproduces the illumi- satellite, considering the separation between both satellites nation of an object in modern computer graphic scenes for as seen from the heliocentre. Once again, the spatial resolu- movies and for the full Moon. (Oren and Nayar, 1994). In tion was set as 1 mas (í 3 km). Dias-Oliveira et al. (2013) a simplified version of the model For the penumbra region, the fraction of sunlight that was used. Here, we implemented the complete version in reaches the eclipsed satellite was computed by using a nu- Oren and Nayar (1994), taking into account the direct illu- merical method. The solar limb darkening was taken into mination and all inter-reflection components of the radiance. account by using Hestroffer and Magnan (1998) empirical The albedo ratio between the satellites is determined in- law, with parameters set for the 889 nm spectral region. dependently, by using observations right before and after The light flux was numerically integrated for a given in- the mutual event with the same instrument setup. In Dias- stant by using the profile after the mask was applied (F.t/). Oliveira et al. (2013), analytic expressions involving the ter- Then the light curve was normalised using the Light flux of minator were used to take solar phase angle effects into ac- the eclipsed satellite (F ). Figure 3 is a simulated light curve count in the evaluation of the flux measurements of albedo for the event when Europa eclipsed Ganymede on April 12, observations. Here, following a more rigorous approach, 2015. The black dots represent seven instants for which the we also simulated the 2D profiles of the satellites for these respective 2D profiles are displayed. observations, to better determine the effective area and re- flectance of the satellites to compute more accurate albedos. 3.3. Fitting procedure In fact, first we measured the light fluxes between both We took a somewhat different approach from that in Dias- satellites (F , F ) separately using images right before or af- 1 2 Oliveira et al. (2013). Here, the parameters of interest were ter the occultation. In the other hand, we compute the 2D determined by the minimisation of the Chi-square test, Eq. profile of each satellites for the given instants and obtained (3), where the simulated light curve is compared with the the simulated light flux for each satellite (F , F ), this S1 S2 observed one. simulations already take into consideration the size of the satellite and the phase angle. The ratio of albedo (A _A ) 1 2 can be determined using Equation (2). We considered each .LC * LC / obs model satellite’s albedo as uniform, due to the lack of information .N * P/ = (3) about albedo variations for the wavelength of our observa- tions. N is the number of observations used in the process, and P is the number of parameters fitted (P = 3).  is the light curve’s standard deviation outside the event (the noise). The F A F 1 1 s1 = : : (2) parameters (s ,t and v ) for which the chi-square is mini- 0 0 0 F A F 2 2 s2 2 mum ( =  ) were set as the solution for the obser- min vations. The normalised  is expected to be around 1 for Similar to Dias-Oliveira et al. (2013) we discretise the good fittings. satellite into a 2D profile with elements of a given spatial Initially, a large range of parameter values was tested, resolution. However, unlike Dias-Oliveira et al. (2013), we then the ranges were narrowed as the iterative process pro- created 2D satellite profiles with much better spatial reso- ceeded. For computing speed, the spatial resolution was de- lutions, 1 mas (í 3 km), avoiding eventual round off errors graded in the first steps and then is set to the nominal value in the simulated flux counts. The profiles were positioned of 1 mas as we approached the  minimum. for each instant of time (t) using the separation between the The uncertainty of each parameter (1 error bar) was de- satellites, overlapping the occulted satellite with the occult- termined by changing that parameter from its nominal solu- ing one when necessary, obtaining a new combined 2D pro- 2 2 2 tion value, so that  changes from  to  + 1. This file of both satellites. The light flux (F.t/) was numerically min min procedure is repeated for each parameter: t , s and v . integrated over the combined 2D profile for a given instant 0 0 0 The procedures described in Section 3 were developed (t). Then, this is repeated for every instant (t) of the event as a PYTHON software that analyses and fits observed light to produce the simulated light curve. This simulated light B. Morgado, R. Vieira-Martins, M. Assafin et al.: Preprint submitted to Elsevier Page 4 of 11 Mutual Phenomena between the Jovian moons Figure 2: Occultation simulation. Io was occulted by Europa on February 22, 2015, at 02:07 UTC as seen at OPD. The bottom panel shows the simulated light curve and seven instants (black dots) were highlighted. The 2D profiles for each instant are displayed in the top panel. The albedo ratio was 0.960, as determined before the event. The profile resolution is 1 mas (í 3.2 km). Figure 3: Eclipse simulation. Ganymede was eclipsed by Europa on April 12, 2015 at 01:46 UTC as seen at FOZ. The bottom panel shows the simulated light curve and seven instants (black dots) were highlighted. The 2D profiles for each instant are displayed in the top panel. The profile resolution was 1 mas (í 3.6 km). curves. This software uses functions from NUMPY, ASTROPY 4. Results of the mutual phenomena (Astropy Collaboration et al., 2013), SCIPY and MATPLOTLIB campaigns libraries. 4.1. Re-analysis of the mutual phenomena between the Galilean moons - 2009 From April to October of 2009, 25 light curves for 25 mutual events between the Galilean moons, 13 occultations and 12 eclipses, were observed by our group with the 60 cm Zeiss telescope of the Observatório Pico dos Dias (OPD) us- ing the methane filter. These events were analysed by Dias- B. Morgado, R. Vieira-Martins, M. Assafin et al.: Preprint submitted to Elsevier Page 5 of 11 Mutual Phenomena between the Jovian moons Oliveira et al. (2013) and by Arlot et al. (2014a). However, the results obtained by both presented a significant difference (higher than 3). More recently, Morgado et al. (2016) analysed 5 mutual approximations that were auxiliary observations designed for albedo determinations made immediately before and af- ter mutual occultations from the 2009 campaign. The mu- tual approximation technique is based on the same geometri- cal parameters that describe a mutual occultation (Morgado et al., 2016, 2019). From the measurements of the separation between both satellites, this technique allows for determin- ing the time of maximum apparent approach between these satellites, which can be directly compared with the central instant from the occultation. The results obtained with this independent technique agree within the errors with the re- sults obtained by Arlot et al. (2014a). The discrepancy of the results between Dias-Oliveira et al. Figure 4: Light curve of the event where Io was occulted by Europa on February 22 2015, observed with the 0.60 m tele- (2013) and Arlot et al. (2014a), the agreement between the scope at OPD. The measured normalised flux is denoted by later with Morgado et al. (2016) for 5 events, and the sig- black dots and the fitted model represented by the red line. In nificant improvements in our light curve fitting procedure the bottom panel, the red dots are the residuals of the fitting motivated us to re-analyse our 2009 campaign data. in the sense observation minus fit. Here, we present our updated results for the 2009 events in Table 3. Table 3 is organised as follows: the date of the event and the satellites’ pairs in the form "S oS " for occul- 1 2 (2014a) and Morgado et al. (2016) within 1. We have a rms tation and "S eS " for eclipses, where 1 stands for Io, 2 for 1 2 of 9.9 mas (í 30 km) and 14.8 mas (í 45 km) in comparison Europa, 3 for Ganymede and 4 for Callisto. We furnish the 3 4 with the JPL’s and the IMCCE’s ephemeris. sites involved in each observation (using the alias defined in Table 1). For each event, we give the obtained central in- 4.2. Mutual phenomena between the Galilean stant (t ) and its uncertainty (t ) in seconds of time, and the 0 0 moons - 2014-2015 ephemeris offset (t ) in mas; the impact parameter (s ), its 0 0 Here we present the results concerning the latest cam- uncertainty (s ) and its offset (s ) in mas; and the rela- 0 0 paign. We obtained new 47 light curves, 31 occultations tive velocity (v ), its uncertainty (v ) and its offset (v ) in 0 0 0 and 16 eclipses, from 40 events observed by 3 stations in mas per second. All times are UTC. In the last columns, we the South and South-East of Brazil. have the rms between the observed light fluxes and the fitted An example is the event where Europa occulted Io on ones, the number of images utilised (N) and the normalised February 22 2015. The observed light curve is illustrated of our fit. in Figure 4. In the upper panel, the black dots are the light The corresponding inter-satellite tangential coordinates flux observed and the red line the model fitted. The bottom (X and Y ) in the sense occulting/eclipsing satellite minus oc- panel contains the residuals in the sense observation minus culted/eclipsed satellite for the central instant can be found in model. For this event, the offset for the central instant was the supplementary material, such form is the same presented +6.5 mas (í 21 km) and for the impact parameter -0.8 mas by Emelyanov and Gilbert (2006); Emelyanov (2009); Ar- (í 3 km). The offsets regard to the JPL’s jup310 and the lot et al. (2014a); Saquet et al. (2018). The plots of the DE435 ephemeris. re-fitted light curves are available as online material in the A second example is the case when Europa eclipsed Ganymede supplementary material. The simulated and observed ASCII on April 12 2015. The observed light curve is illustrated in light curve files are freely available in electronic form at the Figure 5. For this event, the offset for the central instant was NSDB . +13.6 mas (í 49 km), and for the impact parameter +9.6 The re-analysis resulted in a mean uncertainty of 15.3 mas (í 35 km). mas (í 46 km) for the impact parameter and 4.9 mas (í 15 The multiple coverage observational strategy reduced the km) for the central instant. In Table 4, we compare the up- number of events lost by overcast weather or instrumental is- dated results with the ones from Arlot et al. (2014a), Dias- sues. An example was the Io occultation by Ganymede on Oliveira et al. (2013) and Morgado et al. (2016). The error March 03 2015, OPD and FOZ observed this event. Also, of each parameter normalises the differences. If the value is two other stations in the USA observed this event, one in less than one, both results agree within 1. At the bottom, Arnold (AAC) and another in Scottsdale (SCO). These ob- we have the mean difference and the standard deviation for servations were made in the context of the international mu- each parameter. The JPL ephemeris utilised was jup310 and DE435. The updated results now agree with those by Arlot et al. The IMCCE ephemeris utilised was NOE-5-2010-GAL and DE435. Website: http://nsdb.imcce.fr/nsdb/home.html B. Morgado, R. Vieira-Martins, M. Assafin et al.: Preprint submitted to Elsevier Page 6 of 11 Mutual Phenomena between the Jovian moons the satellites’ pairs in the form "S oS " for occultation and 1 2 "S eS " for eclipses, where 1 stands for Io, 2 for Europa, 3 1 2 for Ganymede and 4 for Callisto. We furnish the sites in- volved in each observation (using the alias defined in Ta- ble 1). For each event, we give the obtained central instant (t ) and its uncertainty (t ) in seconds of time, and the 0 0 ephemeris offset (t ) in mas; the impact parameter (s ), 0 0 its uncertainty (s ) and the offset (s ) in mas; and the rel- 0 0 ative velocity (v ), its uncertainty (v ) and the offset (v ) 0 0 0 in mas per second. All times are UTC. In the last columns, we have the rms between the observational curves and the fitted ones, the number of images utilised (N) and the nor- malised  of our fit. The corresponding inter-satellite tangential coordinates (X and Y ) in the sense occulting/eclipsing satellite minus oc- culted/eclipsed satellite for the central instant can be found in Figure 5: Light curve of the event when Ganymede was the supplementary material, such formalism is the same pre- eclipsed by Europa on April 12 2015, observed with the 0.28 sented by Emelyanov and Gilbert (2006); Emelyanov (2009); m telescope at FOZ. Arlot et al. (2014a); Saquet et al. (2018). The plots of the fit- ted light curves are available as online material in the supple- mentary material. The simulated and observed light curve ASCII files are freely available in electronic form at the NSDB. The mean uncertainty of our results is 14.8 mas (í 45 km) for the impact parameter and 7.5 mas (í 23 km) for the central instant. The rms relative to JPL ephemeris was 9.2 mas (í 28 km) and 13.5 mas (í 40 km) relative to IMCCE ephemeris. From the 2014-2015 events, 10 were also analysed using a different procedure and published by Saquet et al. (2018). In average, the comparison between this procedure and ours agrees within 1. 4.3. Amalthea eclipsed by Ganymede - 02 March One particular event in our 2014-2015 observational cam- paign was the eclipse of Amalthea by Ganymede. The as- trometry of this inner satellite is not easy to be done due to its proximity to Jupiter (major semi-axis equal to 2.54 Jupiter’s radius). Often, coronagraphy techniques are needed to sep- arate this object from Jupiter’s scattered light (Kulyk et al., 2002; Veiga and Vieira Martins, 2005; Robert et al., 2017). The positional uncertainty of classical astrometry for this Figure 6: Light curves of the event where Io was occulted by satellite is in the 120 mas level (í 360 km). Ganymede on March 03 2015, observed at AAC, SCO, OPD The observation of mutual eclipses involving Galilean and FOZ, respectively. The AAC and SCO light curves were reduced and analysed by Saquet et al. (2018). moons and inner satellite was strongly advocated by Vachier et al. (2002). The first registration of this kind of event was given by Christou et al. (2010), regarding three eclipses of Amalthea observed during the 2009 mutual phenomena cam- tual phenomena campaign PHEMU15, (Saquet et al., 2018; paign. More recently, Saquet et al. (2016) also analysed three Emel’yanov, 2017). Both light curves are available at the more eclipses of Amalthea and the first observation of an NSDB. In the Figure 6 we compare our light curves (OPD eclipse of Thebe during the 2014-2015 campaign. and FOZ) with the ones analysed by Saquet et al. (2018) Here we present the results of one eclipse of Amalthea (AAC and SCO). The central instant obtained by the obser- by Ganymede observed on March, 2 of 2015 at the 1.6 m vations agrees within 2. Notice that all curves present sim- Perkin-Elmer telescope of the Observatório Pico dos Dias ilar features and we highlight the small residual in our light (OPD, MPC code: 874). This observation was made using curves. the IKon-L CCD camera with the narrow Methane filter . The results for these events are presented in Table 5. Table 5 is organised as follows: the date of the event and 5 Centred at 889 nm with a width of 15 nm. B. Morgado, R. Vieira-Martins, M. Assafin et al.: Preprint submitted to Elsevier Page 7 of 11 Mutual Phenomena between the Jovian moons Figure 7: Digital coronagraphy on an image containing part of Jupiter and its inner satellite Amalthea, as obtained on 02 March, 2015, with the 1.6 m Perkin-Elmer telescope. The left panel (i) displays the original image. The central panel (ii) shows the clean bright object profile (Jupiter) obtained. The right panel (iii) displays the final coronagraphed image (see text). After correcting by Bias and Flat-Field using the same procedure described above, we applied a digital coronagra- phy technique to reduce the influence of Jupiter brightness in the images, this coronagraphy was done using the PRAIA package (Assafin et al., 2008, 2009). Briefly, the procedure is as follows. The centroid of the bright object is iteratively determined. Concentric rings with radius R are formed for each image pixel at a distance R to the centroid. Quartile statistics of weighted fluxes inside each ring are performed, and for each image pixel an average count is assigned. The result is an improved profile with cleaner pixel counts that better represent the bright object. The profile is then sub- tracted from the original image, resulting in the final coron- agraphed (science) image, see Figure 7. Aperture photometry was done using the PRAIA package, where the size of the aperture was manually determined to Figure 8: Light curve of the event when Amalthea was eclipsed maximise the signal to noise ratio. The light curve simula- by Ganymede on March 02 2015, observed with the 1.6 m tion and fitting procedure were the same described in Sec- telescope at OPD. tions 3.2.2 and 3.3. Notice that Amalthea’s triaxial shape is 125  73  64 km (uncertainty of 2 km in each axis; Thomas et al. (1998)) and its rotation phase during the event was un- 5. Conclusions known. However, without any loss of precision, in our sim- ulations, we considered Amalthea as an equivalent sphere We presented in this paper the results for 40 mutual events with a radius equal to 83.5 km. Because of the time resolu- from the observation and analysis of 47 light curves, 31 oc- tion of the observations (8 seconds), the spherical shape was cultations and 16 eclipses, obtained during the 2014-2015 indistinguishable from the elliptical one. mutual phenomena campaign between the Galilean satellites. The light curve of this event is illustrated in Fig. 8. We The observations were made at three stations in the South obtained a central instant with an uncertainty of 19.4 mas and South-East of Brazil, using telescopes with diameters (í 58.2 km) and an impact parameter with an uncertainty of ranging between 28 and 60 cm. We also obtained updated 76.3 mas (í 228 km). This corresponds to a mean uncer- results from the re-analysis of 25 mutual events, 13 occul- tainty of 47.8 mas (í 143 km). The result of this event is tations and 12 eclipses, observed in Brazil by our group in displayed in the last line of Table 5. The positions obtained 2009 with a 60 cm aperture telescope. In all observations, by Christou et al. (2010) had mean uncertainty of 82 mas (í we used a narrow band methane filter centred at 889 nm with 246 km) and the ones obtained by Saquet et al. (2016) had a a width of 15 nm, that eliminates Jupiter’s scattered light. mean uncertainty of 45 mas (í 135 km). We used the Oren-Nayer model (Oren and Nayar, 1994) to characterise the reflectance of the surface of the satellites. It is a generalisation of Lambert’s scattering law. The main advantage of this model is that it does not require previous B. Morgado, R. Vieira-Martins, M. Assafin et al.: Preprint submitted to Elsevier Page 8 of 11 Mutual Phenomena between the Jovian moons knowledge about the satellite surface, and delivers excellent sion Europa Clipper , scheduled to be launched in the next results. The albedo ratio was instrumentally obtained by us- decade (2020s). ing satellite observations before and after the mutual occulta- tions, and the same light curve simulation routines that take Acknowledgements into account solar phase angle and surface reflectance. In We thank our anonymous referee and N. Emelyanov for our procedure, the simulated light curves fitted to the ob- helpful comments. This study was financed by the Coorde- served ones had a normalised chi-square very close to 1.0, nação de Aperfeiçoamento de Pessoal de Nível Superior - indicating good agreement of our model to the data. Brasil (CAPES) - Finance Code 001. Part of this research The updated results for 2009 now agree within 1 with is suported by INCT do e-Universo, Brazil (CNPQ grants those from Arlot et al. (2014a) and from Morgado et al. (2016). 465376/2014-2). Based in part on observations made at the The internal mean uncertainty was 10.1 mas (í 31 km). For Laboratório Nacional de Astrofísica (LNA), Itajubá-MG, Brazil. the 2014-2015 campaign, the mean uncertainty was 11.2 mas BM thanks the CAPES/Cofecub-394/2016-05 grant. RVM (í 35 km). There is no significant difference between the acknowledges the grants: CNPq-304544/2017-5, 401903/2016- quality of the data for both campaigns. Our result is compa- 8, Capes/Cofecub-2506/2015, Faperj: PAPDRJ-45/2013 and rable with other independent observations during the same E-26/203.026/2015. MA thanks CNPq (Grants 427700/2018- campaign, such as the 23 light curves observed by Vasund- 3, 310683/2017-3 and 473002/2013-2) and FAPERJ (Grant hara et al. (2017) and the 21 light curves observed by Zhang E-26/111.488/2013). JIBC acknowledges CNPq grants 308489/2013- et al. (2019), which report uncertainties in the 10 and 20 mas 6 and 308150/2016-3. RS and OCW acknowledges Fapesp level (30 and 60 km), respectively. Saquet et al. (2018) pub- proc. 2016/24561-0 and 2011/08171-3, CNPq proc. 312813/2013- lished the light curves of the international campaign organ- 9 and 305737/2015-5. FBR acknowledges CNPq support, ised by the IMCCE. It reunites 609 light curves with a mean proc. 309578/2017-5. GBR thanks to the support of the internal uncertainty of 24 mas (í 75 km) . From these 609 CAPES and FAPERJ/PAPDRJ (E26/203.173/2016) grants. light curves, 10 were observed by our group, the parameters ARGJ thanks FAPESP proc. 2018/11239-8. This collabora- obtained by Saquet et al. (2018) agree, on average, within 1 tion, as part of the Encelade working group, has been sup- with the results presented here. ported by the International Space Sciences Institute (ISSI) Compared to the 236 mutual events covered by the 2014- in Bern, Switzerland. 2015 international PHEMU campaign lead by the IMCCE, Observatoire de Paris (Saquet et al., 2018), the 40 ones cov- ered by our 2014-2015 campaign (with only 5 events in com- References mon) represent a significant contribution of about 17%. No- Aksnes, K., Franklin, F.A., 1976. Mutual phenomena of the Galilean satel- tice that this campaign was favourable for the north hemi- lites in 1973. III. Final results from 91 light curves. Astronomical Journal sphere, enhancing the weight of our southern results due to 81, 464–481. doi:10.1086/111908. parallax effects. Equally, our new results for the 25 events Arlot, J.E., Cooper, N., Emelyanov, N., Lainey, V., Meunier, L.E., Mur- ray, C., Oberst, J., Pascu, D., Pasewaldt, A., Robert, V., Tajeddine, R., represent about 15% of the 172 mutual events covered by Willner, K., 2017. Natural satellites astrometric data from either space the international PHEMU campaign of 2009 (Arlot et al., probes and ground-based observatories produced by the European con- 2014a). In a similar way, our result for the eclipse by Amalthea sortium “ESPaCE”. Notes Scientifiques et Techniques de l’Institut de th is only the 4 such measurement ever published for the 2014- Mecanique Celeste 105. 2015 campaign, after the 3 ones observed by Saquet et al. Arlot, J.E., Emelyanov, N., 2019. Natural satellites mutual phenomena ob- (2016), representing a significant contribution to the orbit of servations: Achievements and future. Planetary and Space Science 169, 70–77. doi:10.1016/j.pss.2019.02.004. this inner satellite of Jupiter. All the data are freely available Arlot, J.E., Emelyanov, N., Varfolomeev, M.I., Amossé, A., Arena, C., As- to anyone at NSDB for further research and orbital fitting. safin, M., Barbieri, L., Bolzoni, S., Bragas-Ribas, F., Camargo, J.I.B., All these results can be used to improve the orbit and Casarramona, F., Casas, R., Christou, A., Colas, F., Collard, A., Combe, ephemeris of the Galilean satellites (plus Amalthea) taking S., Constantinescu, M., Dangl, G., De Cat, P., Degenhardt, S., Delcroix, M., Dias-Oliveira, A., Dourneau, G., Douvris, A., Druon, C., Ellington, into account the tidal forces, as pointed out by Lainey et al. C.K., Estraviz, G., Farissier, P., Farmakopoulos, A., Garlitz, J., Gault, (2009). D., George, T., Gorda, S.Y., Grismore, J., Guo, D.F., Herald, D., Ida, The next mutual phenomena events for the Galilean satel- M., Ishida, M., Ivanov, A.V., Klemt, B., Koshkin, N., Le Campion, J.F., lites will occur in 2021 and will favour the southern hemi- Liakos, A., Liao, S.L., Li, S.N., Loader, B., Lopresti, C., Lo Savio, E., sphere, due to Jupiter’s declination. An observational cam- Marchini, A., Marino, G., Masi, G., Massallé, A., Maulella, R., McFar- land, J., Miyashita, K., Napoli, C., Noyelles, B., Pauwels, T., Pavlov, paign such as this one will be organised in due time, the H., Peng, Q.Y., Perelló, C., Priban, V., Prost, J., Razemon, S., Rous- prediction of these events are already in the IMCCE web- selle, J.P., Rovira, J., Ruisi, R., Ruocco, N., Salvaggio, F., Sbarufatti, site (Arlot and Emelyanov, 2019). These campaigns can G., Shakun, L., Scheck, A., Sciuto, C., da Silva Neto, D.N., Sinyaeva, increase the accuracy and precision of ephemeris and can be N.V., Sofia, A., Sonka, A., Talbot, J., Tang, Z.H., Tejfel, V.G., Thuillot, helpful to space missions aimed at the Jovian system. For ex- W., Tigani, K., Timerson, B., Tontodonati, E., Tsamis, V., Unwin, M., Venable, R., Vieira-Martins, R., Vilar, J., Vingerhoets, P., Watanabe, H., ample, we have the ESA mission JUICE and NASA’s mis- Yin, H.X., Yu, Y., Zambelli, R., 2014a. The PHEMU09 catalogue and The standard deviation after fitting the light curves. astrometric results of the observations of the mutual occultations and Website: http://nsdb.imcce.fr/multisat/nssephme.htm Website: https://www.nasa.gov/europa/. Website: http://sci.esa.int/juice/. B. Morgado, R. Vieira-Martins, M. Assafin et al.: Preprint submitted to Elsevier Page 9 of 11 Mutual Phenomena between the Jovian moons eclipses of the Galilean satellites of Jupiter made in 2009. Astronomy particle phase function hockey stick relation. Icarus 221, 1079–1083. and Astrophysics 572, A120. doi:10.1051/0004-6361/201423854. doi:10.1016/j.icarus.2012.10.022. Arlot, J.E., Saquet, E., Robert, V., Lainey, V., 2014b. The Phemu 2015 Hapke, B., Wells, E., 1981. Bidirectional reflectance spectroscopy. 2. Ex- campaign of observations of the mutual events of the Galilean satellites periments and observations. Journal of Geophysical Research 86, 3055– of Jupiter, in: European Planetary Science Congress, pp. EPSC2014–57. 3060. doi:10.1029/JB086iB04p03055. Assafin, M., Campos, R.P., Vieira Martins, R., da Silva Neto, D.N., Ca- Hestroffer, D., Magnan, C., 1998. Wavelength dependency of the Solar margo, J.I.B., Andrei, A.H., 2008. Instrumental and digital coronagra- limb darkening. Astronomy and Astrophysics 333, 338–342. phy for the observation of the Uranus satellites’ upcoming mutual events. Karkoschka, E., 1994. Spectrophotometry of the Jovian Planets and Titan Planetary and Space Science 56, 1882–1887. doi:10.1016/j.pss.2007. at 300- to 1000-nm Wavelength: The Methane Spectrum. Icarus 111, 05.030. 174–192. doi:10.1006/icar.1994.1139. Assafin, M., Vieira-Martins, R., Braga-Ribas, F., Camargo, J.I.B., da Silva Karkoschka, E., 1998. Methane, Ammonia, and Temperature Measure- Neto, D.N., Andrei, A.H., 2009. Observations and Analysis of Mutual ments of the Jovian Planets and Titan from CCD-Spectrophotometry. Events between the Uranus Main Satellites. Astronomical Journal 137, Icarus 133, 134–146. doi:10.1006/icar.1998.5913. 4046–4053. doi:10.1088/0004-6256/137/4/4046. Kiseleva, T.P., Kiselev, A.A., Kalinichenko, O.A., Vasilyeva, N.A., Assafin, M., Vieira Martins, R., Camargo, J.I.B., Andrei, A.H., Da Silva Khovricheva, M.L., 2008. Results of astrometric observations of Neto, D.N., Braga-Ribas, F., 2011. PRAIA - Platform for Reduction Jupiter’s Galilean satellites at the Pulkovo Observatory from 1986 of Astronomical Images Automatically, in: Gaia follow-up network for to 2005. Solar System Research 42, 414–433. doi:10.1134/ the solar system objects : Gaia FUN-SSO workshop proceedings, pp. S0038094608050055. 85–88. Kulyk, I., Jockers, K., Karpov, N., Sergeev, A., 2002. Astrometric CCD Astropy Collaboration, Robitaille, T.P., Tollerud, E.J., Greenfield, P., observations of the inner Jovian satellites in 1999-2000. Astronomy and Droettboom, M., Bray, E., Aldcroft, T., Davis, M., Ginsburg, A., Price- Astrophysics 383, 724–728. doi:10.1051/0004-6361:20011770. Whelan, A.M., Kerzendorf, W.E., Conley, A., Crighton, N., Barbary, Lainey, V., Arlot, J.E., Karatekin, Ö., van Hoolst, T., 2009. Strong tidal K., Muna, D., Ferguson, H., Grollier, F., Parikh, M.M., Nair, P.H., Un- dissipation in Io and Jupiter from astrometric observations. Nature 459, ther, H.M., Deil, C., Woillez, J., Conseil, S., Kramer, R., Turner, J.E.H., 957–959. doi:10.1038/nature08108. Singer, L., Fox, R., Weaver, B.A., Zabalza, V., Edwards, Z.I., Azalee Lainey, V., Arlot, J.E., Vienne, A., 2004a. New accurate ephemerides for Bostroem, K., Burke, D.J., Casey, A.R., Crawford, S.M., Dencheva, N., the Galilean satellites of Jupiter. II. Fitting the observations. Astronomy Ely, J., Jenness, T., Labrie, K., Lim, P.L., Pierfederici, F., Pontzen, A., and Astrophysics 427, 371–376. doi:10.1051/0004-6361:20041271. Ptak, A., Refsdal, B., Servillat, M., Streicher, O., 2013. Astropy: A Lainey, V., Duriez, L., Vienne, A., 2004b. New accurate ephemerides for community Python package for astronomy. Astronomy and Astrophysics the Galilean satellites of Jupiter. I. Numerical integration of elaborated 558, A33. doi:10.1051/0004-6361/201322068, arXiv:1307.6212. equations of motion. Astronomy and Astrophysics 420, 1171–1183. Butcher, H., Stevens, R., 1981. Image Reduction and Analysis Facility doi:10.1051/0004-6361:20034565. Development. Kitt Peak National Observatory Newsletter 16, 6. Lainey, V., Jacobson, R.A., Tajeddine, R., Cooper, N.J., Murray, C., Robert, Christou, A.A., Lewis, F., Roche, P., Hidas, M.G., Brown, T.M., 2010. Ob- V., Tobie, G., Guillot, T., Mathis, S., Remus, F., Desmars, J., Arlot, J.E., servational detection of eclipses of J5 Amalthea by the Galilean satel- De Cuyper, J.P., Dehant, V., Pascu, D., Thuillot, W., Le Poncin-Lafitte, lites. Astronomy and Astrophysics 522, A6. doi:10.1051/0004-6361/ C., Zahn, J.P., 2017. New constraints on Saturn’s interior from Cassini 201014822, arXiv:1104.0042. astrometric data. Icarus 281, 286–296. doi:10.1016/j.icarus.2016.07. Dias-Oliveira, A., Vieira-Martins, R., Assafin, M., Camargo, J.I.B., Braga- 014, arXiv:1510.05870. Ribas, F., da Silva Neto, D.N., Gaspar, H.S., Pires dos Santos, P.M., Morgado, B., Assafin, M., Vieira-Martins, R., Camargo, J.I.B., Dias- Domingos, R.C., Boldrin, L.A.G., Izidoro, A., Carvalho, J.P.S., Sfair, Oliveira, A., Gomes-Júnior, A.R., 2016. Astrometry of mutual approx- R., Sampaio, J.C., Winter, O.C., 2013. Analysis of 25 mutual eclipses imations between natural satellites. Application to the Galilean moons. and occultations between the Galilean satellites observed from Brazil in Monthly Notices of the Royal Astronomical Society 460, 4086–4097. 2009. Monthly Notices of the Royal Astronomical Society 432, 225– doi:10.1093/mnras/stw1244, arXiv:1605.06573. 242. doi:10.1093/mnras/stt447. Morgado, B., Vieira-Martins, R., Assafin, M., Machado, D.I., Camargo, Emelyanov, N.V., 2009. Mutual occultations and eclipses of the Galilean J.I.B., Sfair, R., Malacarne, M., Braga-Ribas, F., Robert, V., Bassallo, satellites of Jupiter in 2002-2003: final astrometric results. Monthly T., Benedetti-Rossi, G., Boldrin, L.A., Borderes-Motta, G., Camargo, Notices of the Royal Astronomical Society 394, 1037–1044. doi:10. B.C.B., Crispim, A., Dias-Oliveira, A., Gomes-Júnior, A.R., Lainey, V., 1111/j.1365-2966.2009.14398.x. Miranda, J.O., Moura, T.S., Ribeiro, F.K., Santana, T., Santos-Filho, Emel’yanov, N.V., 2017. Current problems of dynamics of moons of planets S., Trabuco, L.L., Winter, O.C., Yamashita, T.A.R., 2019. APPROX and binary asteroids based on observations. Solar System Research 51, - mutual approximations between the Galilean moons: the 2016-2018 20–37. doi:10.1134/S0038094617010014. observational campaign. Monthly Notices of the Royal Astronomical Emelyanov, N.V., Gilbert, R., 2006. Astrometric results of observations Society 482, 5190–5200. doi:10.1093/mnras/sty3040, arXiv:1811.02913. of mutual occultations and eclipses of the Galilean satellites of Jupiter Oren, M., Nayar, S.K., 1994. Generalization of lambert’s reflectance model, in 2003. Astronomy and Astrophysics 453, 1141–1149. doi:10.1051/ in: Proceedings of the 21st Annual Conference on Computer Graph- 0004-6361:20064810. ics and Interactive Techniques, ACM, New York, NY, USA. pp. 239– Hapke, B., 1981. Bidirectional reflectance spectroscopy. 1. Theory. Journal 246. URL: http://doi.acm.org/10.1145/192161.192213, doi:10.1145/ of Geophysical Research 86, 4571–4586. 192161.192213. Hapke, B., 1984. Bidirectional reflectance spectroscopy 3. Correction for Peng, Q.Y., He, H.F., Lainey, V., Vienne, A., 2012. Precise CCD positions macroscopic roughness. Icarus 59, 41–59. doi:10.1016/0019-1035(84) of Galilean satellite-pairs. Monthly Notices of the Royal Astronomical 90054-X. Society 419, 1977–1982. doi:10.1111/j.1365-2966.2011.19852.x. Hapke, B., 1986. Bidirectional reflectance spectroscopy 4. The extinction Robert, V., Saquet, E., Colas, F., Arlot, J.E., 2017. CCD astrometric obser- coefficient and the opposition effect. Icarus 67, 264–280. doi:10.1016/ vations of Amalthea and Thebe in the Gaia era. Monthly Notices of the 0019-1035(86)90108-9. Royal Astronomical Society 467, 694–698. doi:10.1093/mnras/stx123. Hapke, B., 2002. Bidirectional Reflectance Spectroscopy. 5. The Coherent Saquet, E., Emelyanov, N., Colas, F., Arlot, J.E., Robert, V., Christophe, B., Backscatter Opposition Effect and Anisotropic Scattering. Icarus 157, Dechambre, O., 2016. Eclipses of the inner satellites of Jupiter observed 523–534. doi:10.1006/icar.2002.6853. in 2015. Astronomy and Astrophysics 591, A42. doi:10.1051/0004-6361/ Hapke, B., 2008. Bidirectional reflectance spectroscopy. 6. Effects of poros- 201628246, arXiv:1605.06935. ity. Icarus 195, 918–926. doi:10.1016/j.icarus.2008.01.003. Saquet, E., Emelyanov, N., Robert, V., Arlot, J.E., Anbazhagan, P., Bail- Hapke, B., 2012. Bidirectional reflectance spectroscopy 7. The single lié, K., Bardecker, J., Berezhnoy, A.A., Bretton, M., Campos, F., Ca- B. Morgado, R. Vieira-Martins, M. Assafin et al.: Preprint submitted to Elsevier Page 10 of 11 Mutual Phenomena between the Jovian moons pannoli, L., Carry, B., Castet, M., Charbonnier, Y., Chernikov, M.M., Christou, A., Colas, F., Coliac, J.F., Dangl, G., Dechambre, O., Del- croix, M., Dias-Oliveira, A., Drillaud, C., Duchemin, Y., Dunford, R., Dupouy, P., Ellington, C., Fabre, P., Filippov, V.A., Finnegan, J., Foglia, Table 2 S., Font, D., Gaillard, B., Galli, G., Garlitz, J., Gasmi, A., Gaspar, H.S., Mutual events and observation conditions. Gault, D., Gazeas, K., George, T., Gorda, S.Y., Gorshanov, D.L., Gual- doni, C., Guhl, K., Halir, K., Hanna, W., Henry, X., Herald, D., Houdin, Date Event Obs. i z Ratio of o o G., Ito, Y., Izmailov, I.S., Jacobsen, J., Jones, A., Kamoun, S., Karda- yy-mm-dd S xS ( ) ( ) Albedo 1 2 sis, E., Karimov, A.M., Khovritchev, M.Y., Kulikova, A.M., Laborde, 14-11-02 4o1 OPD 10.60 73.49 3.56 , 0.03 J., Lainey, V., Lavayssiere, M., Le Guen, P., Leroy, A., Loader, B., 14-11-19 4o2 OPD 10.69 43.44 3.59 , 0.02 Lopez, O.C., Lyashenko, A.Y., Lyssenko, P.G., Machado, D.I., Maig- 14-12-20 2o1 FOZ 8.63 50.72 0.97 , 0.03 urova, N., Manek, J., Marchini, A., Midavaine, T., Montier, J., Mor- 14-12-21 4e1 FOZ 8.52 77.19 – gado, B.E., Naumov, K.N., Nedelcu, A., Newman, J., Ohlert, J.M., Ok- 14-12-21 3o1 FOZ 8.52 65.78 1.65 , 0.19 sanen, A., Pavlov, H., Petrescu, E., Pomazan, A., Popescu, M., Pratt, A., Raskhozhev, V.N., Resch, J.M., Robilliard, D., Roschina, E., Rothen- 14-12-24 2e3 FOZ 8.13 40.85 – berg, E., Rottenborn, M., Rusov, S.A., Saby, F., Saya, L.F., Selvaku- 15-01-21 2e1 FOZ 3.45 46.09 – mar, G., Signoret, F., Slesarenko, V.Y., Sokov, E.N., Soldateschi, J., 15-02-02 3o2 FOZ 0.95 67.15 1.53 , 0.09 Sonka, A., Soulie, G., Talbot, J., Tejfel, V.G., Thuillot, W., Timerson, 15-02-22 2o1 OPD 3.20 41.50 0.99 , 0.04 B., Toma, R., Torsellini, S., Trabuco, L.L., Traverse, P., Tsamis, V., Un- 15-02-22 2e1 OPD 3.20 43.44 – win, M., Abbeel, F.V.D., Vand enbruaene, H., Vasundhara, R., Velikod- 15-03-01 2o1 FOZ 4.57 57.20 1.01 , 0.06 sky, Y.I., Vienne, A., Vilar, J., Vugnon, J.M., Wuensche, N., Zeleny, P., 15-03-01 2e1 FOZ 4.58 65.89 – 2018. The PHEMU15 catalogue and astrometric results of the Jupiter’s 15-03-03 3o1 OPD 4.94 62.37 1.56 , 0.03 Galilean satellite mutual occultation and eclipse observations made in 15-03-03 3o1 FOZ 4.94 58.30 1.56 , 0.03 2014-2015. Monthly Notices of the Royal Astronomical Society 474, 15-03-06 1e2 OPD 5.46 43.08 – 4730–4739. doi:10.1093/mnras/stx2957. Thomas, P.C., Burns, J.A., Rossier, L., Simonelli, D., Veverka, J., Chap- 15-03-09 3e2 OPD 6.14 44.04 – man, C.R., Klaasen, K., Johnson, T.V., Belton, M.J.S., Galileo Solid 15-03-13 1e2 FOZ 6.79 48.89 – State Imaging Team, 1998. The Small Inner Satellites of Jupiter. Icarus 15-03-13 1e3 FOZ 6.66 60.76 – 135, 360–371. doi:10.1006/icar.1998.5976. 15-03-16 4o2 FOZ 7.11 48.84 3.58 , 0.03 Vachier, F., Arlot, J.E., Thuillot, W., 2002. Mutual phenomena involving J5 15-03-17 3e2 FOZ 7.27 58.36 – Amalthea in 2002-2003. Astronomy and Astrophysics 394, L19–L21. 15-03-18 2e1 GOA 7.54 41.46 – doi:10.1051/0004-6361:20021329. 15-03-24 3o4 OPD 8.22 45.53 0.41 , 0.03 Vasundhara, R., Selvakumar, G., Anbazhagan, P., 2017. Analysis of mu- 15-03-24 3o4 FOZ 8.22 46.49 0.41 , 0.03 tual events of Galilean satellites observed from VBO during 2014-2015. 15-03-25 2o1 FOZ 8.47 46.28 0.99 , 0.04 Monthly Notices of the Royal Astronomical Society 468, 501–508. 15-03-25 2o1 GOA 8.47 42.75 0.99 , 0.04 doi:10.1093/mnras/stx437, arXiv:1704.03518. Veiga, C.H., Vieira Martins, R., 2005. CCD astrometric observations of 15-03-26 2e1 OPD 8.48 52.32 – Amalthea and Thebe. Astronomy and Astrophysics 437, 1147–1150. 15-04-02 2o1 OPD 9.25 62.96 1.01 , 0.07 doi:10.1051/0004-6361:20042387. 15-04-02 2o1 FOZ 9.25 59.06 1.01 , 0.07 Zhang, X.L., Han, X.L., Arlot, J.E., 2019. Mutual events between Galilean 15-04-02 2e1 OPD 9.26 82.45 – satellites observed with SARA 0.9 m and 0.6 m telescopes during 2014- 15-04-03 1o3 FOZ 9.43 46.38 0.75 , 0.13 2015. Monthly Notices of the Royal Astronomical Society 483, 4518– 15-04-06 1e2 FOZ 9.70 46.72 – 4524. doi:10.1093/mnras/sty3030. 15-04-12 2e3 FOZ 10.09 65.99 – 15-04-14 1e2 OPD 10.22 69.43 – 15-04-17 4o1 OPD 10.43 54.20 3.90 , 0.08 15-04-17 4o1 GOA 10.43 55.82 3.90 , 0.08 15-04-18 4o3 OPD 10.48 44.23 2.32 , 0.05 15-04-18 1o3 OPD 10.44 72.62 0.64 , 0.11 15-04-25 1o3 OPD 10.73 58.34 0.69 , 0.07 15-04-25 1o3 FOZ 10.73 55.16 0.63 , 0.07 15-04-26 2o1 OPD 10.75 43.56 1.04 , 0.15 15-04-29 3o1 OPD 10.80 67.89 1.60 , 0.06 15-04-29 3o1 GOA 10.80 70.74 1.60 , 0.06 15-05-03 2o1 OPD 10.85 61.83 0.96 , 0.07 15-05-05 3o2 FOZ 10.85 47.24 1.61 , 0.06 15-05-13 3o2 OPD 10.77 86.09 1.48 , 0.07 15-06-04 2o1 FOZ 9.72 56.26 1.01 , 0.05 15-06-18 3o1 GOA 8.59 59.64 1.58 , 0.04 15-03-02 3e5 1.60 3.17 50.90 – Note: The solar phase angle, zenith distance and ratio of albedo in the sense S _S for each event. 2 1 B. Morgado, R. Vieira-Martins, M. Assafin et al.: Preprint submitted to Elsevier Page 11 of 11 Mutual Phenomena between the Jovian moons Table 3 Updated results of the 2009 mutual phenomena campaign observed in Brazil. Date Event Obs. t , t (UTC) t s , s s v , v v rms N 0 0 0 0 0 0 0 0 0 yy-mm-dd S xS hh:mm:ss.s , s.s mas mas mas mas/s mas/s 1 2 09-04-27 3o1 OPD 06:42:53.5 , 0.6 -00.8 121.0 , 20.0 -09.2 6.26 , 0.02 -0.01 0.010 0801 1.005 09-05-09 2o1 OPD 07:21:54.9 , 0.7 +01.5 538.0 , 14.3 -02.3 7.50 , 0.03 -0.03 0.010 0580 1.007 09-05-21 1o2 OPD 05:29:44.3 , 0.7 -00.2 033.6 , 21.1 +26.5 6.36 , 0.06 -0.00 0.014 0869 1.005 09-05-28 1o2 OPD 07:44:17.6 , 0.4 +00.0 219.1 , 17.5 +47.0 6.26 , 0.03 -0.01 0.007 0801 1.005 09-06-10 3e4 OPD 07:29:37.4 , 1.2 +14.2 315.0 , 17.9 +22.8 3.96 , 0.03 +0.01 0.009 1285 1.003 09-06-16 3e1 OPD 08:45:10.5 , 1.2 -25.6 928.4 , 13.3 +19.8 1.52 , 0.01 -0.00 0.008 2697 1.001 09-06-19 4e2 OPD 05:11:36.3 , 0.2 +25.1 489.1 , 13.1 +02.7 5.06 , 0.01 +0.02 0.007 0901 1.004 09-06-19 4e1 OPD 08:32:49.4 , 1.0 +15.6 930.1 , 13.5 +05.9 5.29 , 0.05 +0.02 0.006 1201 1.003 09-06-20 4e1 OPD 05:09:45.8 , 1.3 -22.6 530.5 , 12.8 +01.5 1.26 , 0.01 -0.01 0.006 2326 1.002 09-06-20 4e1 OPD 09:37:28.5 , 1.3 +31.6 412.3 , 18.5 -08.4 1.39 , 0.01 +0.00 0.013 1704 1.002 09-06-22 1o2 OPD 03:27:54.8 , 1.6 +12.0 576.1 , 14.0 +23.7 5.78 , 0.03 -0.03 0.010 0905 1.004 09-06-29 1o2 OPD 05:38:28.7 , 0.8 +00.5 606.4 , 13.3 +10.0 5.51 , 0.02 -0.01 0.004 0801 1.005 09-07-04 1e3 OPD 06:25:13.5 , 0.5 -12.7 395.5 , 16.3 -19.8 7.27 , 0.04 +0.02 0.008 1641 1.002 09-07-06 1e2 OPD 06:17:16.1 , 1.8 +06.2 718.5 , 14.2 -12.2 4.79 , 0.12 +0.02 0.006 2001 1.002 09-07-06 1o2 OPD 07:48:34.5 , 0.5 +17.4 603.5 , 12.9 -05.5 5.34 , 0.01 -0.03 0.005 1004 1.004 09-07-08 3e1 OPD 08:31:14.5 , 0.2 +30.1 223.9 , 14.7 -07.2 5.96 , 0.01 +0.02 0.007 1758 1.002 09-07-13 1e2 OPD 08:38:46.3 , 1.7 +10.3 623.4 , 15.3 -08.0 4.46 , 0.11 +0.02 0.012 2001 1.002 09-08-07 1e2 OPD 05:14:54.9 , 1.3 +12.4 444.6 , 15.1 +34.8 3.09 , 0.03 +0.01 0.021 1775 1.002 09-08-07 1o2 OPD 05:37:48.4 , 0.8 -08.8 283.1 , 20.3 -18.7 3.77 , 0.01 -0.01 0.008 1664 1.002 09-08-12 3o2 OPD 02:10:59.1 , 3.9 +01.0 1059.1 , 13.1 -24.7 2.71 , 0.02 -0.01 0.004 1296 1.003 09-08-22 1o2 OPD 04:07:54.9 , 2.2 +16.1 674.5 , 12.9 +11.2 1.80 , 0.01 +0.00 0.004 2454 1.002 09-09-16 1o2 OPD 00:46:04.4 , 0.7 -12.6 580.4 , 15.0 +06.5 3.71 , 0.02 +0.00 0.011 0976 1.004 09-09-16 1e2 OPD 02:15:11.0 , 0.4 +04.0 172.4 , 12.9 +00.8 3.53 , 0.02 +0.02 0.007 1095 1.004 09-10-24 3o2 OPD 00:35:33.9 , 1.5 +00.8 629.1 , 13.6 -48.8 4.17 , 0.02 -0.01 0.008 1032 1.004 09-10-25 1o2 OPD 01:21:30.8 , 3.9 +01.0 580.6 , 17.5 +08.5 5.35 , 0.07 -0.01 0.014 0252 1.016 Note: The results for the mutual phenomena campaign of 2009. t stand for the UTC central instant, s is the impact 0 0 parameter and v is the apparent relative velocity in the sky plane. Also contains the uncertainty in each parameter (t , s 0 0 0 and v ) and the difference between the fitted ones and the ones expected from the ephemeris jup310 and DE435 (t , s 0 0 0 and v ). In the last columns, we have the rms between the observed light fluxes and the fitted ones, the number of images utilised (N) and the normalised  of our fit. B. Morgado, R. Vieira-Martins, M. Assafin et al.: Preprint submitted to Elsevier Page 12 of 11 Mutual Phenomena between the Jovian moons Table 4 Comparison of the updated results for the 2009 mutual events with Arlot et al. (2014a), Dias-Oliveira et al. (2013) e Morgado et al. (2016). Date Event Obs. Central instant (t ) Impact parameter (s ) 0 0 yy-mm-dd S xS [2] - [1] [3] - [1] [4] - [1] [2] - [1] [3] - [1] [4] - [1] 1 2 09-04-27 3o1 OPD -0.65 -7.15 – +0.15 -0.74 – 09-05-09 2o1 OPD +0.78 -2.34 -0.12 +0.23 +0.97 – 09-05-21 1o2 OPD +0.28 -4.94 – +3.30 +1.05 – 09-05-28 1o2 OPD +0.98 -9.24 +0.57 -0.00 -2.36 – 09-06-10 3e4 OPD +0.13 +1.54 – +0.13 +1.56 – 09-06-16 3e1 OPD +0.99 -1.99 – -0.12 +8.81 – 09-06-19 4e2 OPD -1.87 +0.99 – +0.03 +4.45 – 09-06-19 4e1 OPD -0.83 +0.62 – +0.04 +8.26 – 09-06-20 4e1 OPD -0.70 +1.14 – +0.08 -4.27 – 09-06-20 4e1 OPD -0.03 -2.85 – -0.13 +9.77 – 09-06-22 1o2 OPD +0.27 -2.18 -0.07 +0.16 -0.79 – 09-06-29 1o2 OPD +0.31 -5.30 – +0.07 -0.86 – 09-07-04 1e3 OPD -0.03 +1.17 – -0.26 +4.57 – 09-07-06 1o2 OPD -0.13 +0.25 -0.30 +0.14 +5.74 – 09-07-06 1e2 OPD +0.50 -4.89 – +0.09 -0.01 – 09-07-08 3e1 OPD -0.61 +1.33 – -0.04 +1.84 – 09-07-13 1e2 OPD -0.31 +0.22 – +0.18 +8.36 – 09-08-07 1o2 OPD -0.22 +0.29 +0.44 +0.32 +7.38 – 09-08-07 1e2 OPD -0.34 -1.71 – -0.09 -0.17 – 09-08-12 3o2 OPD +0.83 +0.60 – -0.12 +0.50 – 09-08-22 1o2 OPD +0.31 +0.29 – +0.07 -0.99 – 09-09-16 1e2 OPD -0.64 +3.52 – -0.23 -1.43 – 09-09-16 1o2 OPD +0.51 -0.43 – -0.22 +1.02 – 09-10-24 3o2 OPD -0.25 +14.15 – -0.26 +2.78 – 09-10-25 1o2 OPD +0.62 +1.56 – +0.23 -0.98 – Mean -0.00 -0.61 +0.10 +0.15 +2.18 – Standard deviation 0.67 4.34 0.34 0.67 3.91 – Note: [1] This project, [2] Arlot et al. (2014a), [3] Dias-Oliveira et al. (2013), [4] Morgado et al. (2016). Comparison between the different reduction process divided by the uncertainty of each parameter. B. Morgado, R. Vieira-Martins, M. Assafin et al.: Preprint submitted to Elsevier Page 13 of 11 Mutual Phenomena between the Jovian moons Table 5 Results of the 2014-2015 mutual phenomena campaign observed in Brazil. Date Event Obs. t , t (UTC) t s , s s v , v v rms N 0 0 0 0 0 0 0 0 0 yy-mm-dd S xS hh:mm:ss.s , s.s mas mas mas mas/s mas/s 1 2 14-11-02 4o1 OPD 06:02:14.8 , 1.9 -07.2 288.0 , 18.6 -05.0 1.83 , 0.03 +0.00 0.041 230 1.018 14-11-19 4o2 OPD 07:37:43.9 , 0.2 +00.7 301.3 , 6.0 +00.2 5.04 , 0.03 +0.01 0.010 194 1.021 14-12-20 2o1 FOZ 05:41:38.7 , 0.8 -03.0 162.4 , 8.5 -17.5 1.62 , 0.02 -0.00 0.008 376 1.011 14-12-21 4e1 FOZ 03:22:39.2 , 1.5 +00.1 253.3 , 9.0 +26.2 2.58 , 0.02 +0.00 0.018 281 1.014 14-12-21 3o1 FOZ 04:17:20.7 , 3.2 -00.8 455.1 , 16.1 +14.0 3.66 , 0.03 -0.01 0.151 244 1.017 14-12-24 2e3 FOZ 06:35:06.5 , 1.1 -00.8 255.2 , 8.8 +09.5 1.93 , 0.02 +0.00 0.014 359 1.011 15-01-21 2e1 FOZ 03:52:29.3 , 1.0 +00.8 457.5 , 3.8 +12.8 3.08 , 0.04 +0.00 0.008 141 1.029 15-02-02 3o2 FOZ 07:56:42.3 , 1.9 -02.3 85.0 , 27.4 +08.0 7.78 , 0.05 -0.01 0.022 110 1.038 15-02-22 2o1 OPD 02:07:51.7 , 0.2 +06.5 125.0 , 12.9 -00.8 5.55 , 0.03 -0.02 0.007 201 1.020 15-02-22 2e1 OPD 02:45:11.0 , 0.3 +04.1 17.0 , 9.6 +07.5 4.62 , 0.03 +0.02 0.013 186 1.022 15-03-01 2o1 FOZ 04:12:38.3 , 0.7 -00.7 17.6 , 10.5 +16.8 5.73 , 0.04 -0.01 0.012 131 1.031 15-03-01 2e1 FOZ 05:05:06.7 , 2.4 +06.6 94.6 , 33.3 -06.8 4.94 , 0.04 +0.00 0.059 132 1.031 15-03-03 3o1 OPD 04:08:16.3 , 0.5 +00.3 95.0 , 14.4 +32.3 8.53 , 0.05 -0.03 0.010 100 1.042 15-03-03 3o1 FOZ 04:08:15.5 , 0.6 -07.0 86.8 , 10.8 +23.7 8.50 , 0.05 -0.02 0.010 107 1.039 15-03-06 1e2 OPD 01:16:16.6 , 0.4 +09.4 570.9 , 4.0 +06.2 7.50 , 0.04 +0.04 0.005 124 1.033 15-03-09 3e2 OPD 23:39:32.6 , 0.3 +08.1 67.7 , 7.2 +02.6 5.87 , 0.03 +0.03 0.012 213 1.019 15-03-13 1e2 FOZ 03:29:09.9 , 1.1 +10.8 445.2 , 20.0 +03.7 7.42 , 0.03 +0.03 0.020 198 1.021 15-03-13 1e3 FOZ 23:29:44.5 , 1.3 +04.0 236.3 , 4.4 +08.5 1.56 , 0.02 -0.00 0.015 408 1.010 15-03-16 4o2 FOZ 01:38:58.1 , 0.3 -00.5 373.0 , 6.3 +02.0 3.59 , 0.03 -0.00 0.013 246 1.017 15-03-17 3e2 FOZ 02:53:15.2 , 0.2 +08.2 225.4 , 14.4 -03.1 5.77 , 0.03 +0.02 0.012 234 1.017 15-03-18 2e1 GOA 22:50:43.6 , 1.9 +06.1 398.2 , 25.2 -02.7 5.55 , 0.02 +0.02 0.058 288 1.014 15-03-24 3o4 OPD 00:14:41.4 , 0.9 -38.1 499.0 , 13.9 -16.0 5.39 , 0.03 -0.01 0.006 182 1.022 15-03-24 3o4 FOZ 00:14:41.6 , 0.8 -37.5 519.1 , 9.9 +03.5 5.38 , 0.03 -0.01 0.008 207 1.020 15-03-25 2o1 FOZ 23:35:01.6 , 1.1 +00.8 400.2 , 13.2 +01.7 6.34 , 0.06 -0.02 0.011 083 1.051 15-03-25 2o1 GOA 23:35:01.3 , 0.7 -01.5 398.9 , 6.6 +00.7 6.40 , 0.05 -0.04 0.006 092 1.045 15-03-26 2e1 OPD 01:07:48.0 , 5.3 +10.3 516.7 , 44.7 -12.2 5.78 , 0.04 +0.02 0.048 151 1.027 15-04-02 2o1 OPD 01:43:55.9 , 0.7 +00.0 479.5 , 5.7 -00.5 6.46 , 0.06 -0.02 0.005 082 1.051 15-04-02 2o1 FOZ 01:43:55.6 , 1.3 -02.3 482.3 , 15.0 +02.5 6.50 , 0.05 -0.03 0.011 104 1.040 15-04-02 2e1 OPD 03:24:16.8 , 3.1 +07.8 658.7 , 19.1 -02.1 6.02 , 0.04 +0.01 0.018 129 1.032 15-04-03 1o3 FOZ 22:58:19.0 , 5.9 +05.7 737.3 , 6.0 -05.5 1.23 , 0.02 +0.01 0.021 383 1.011 15-04-06 1e2 FOZ 23:16:40.4 , 0.2 +11.3 54.2 , 10.8 +04.2 6.96 , 0.04 +0.02 0.007 152 1.027 15-04-12 2e3 FOZ 01:46:02.2 , 1.0 +13.6 142.6 , 8.3 +09.6 4.96 , 0.04 +0.00 0.013 121 1.034 15-04-14 1e2 OPD 01:30:58.3 , 0.4 +10.4 49.4 , 8.5 -01.5 6.83 , 0.05 +0.01 0.024 093 1.045 15-04-17 4o1 OPD 23:47:06.9 , 0.9 -01.3 711.5 , 4.5 +01.8 5.05 , 0.05 -0.01 0.007 101 1.041 15-04-17 4o1 GOA 23:47:06.9 , 0.9 -01.0 712.1 , 5.1 +02.7 5.06 , 0.04 -0.02 0.010 128 1.032 15-04-18 4o3 OPD 01:32:30.4 , 0.9 -02.3 69.5 , 15.1 +06.8 5.02 , 0.03 -0.00 0.021 161 1.025 15-04-18 1o3 OPD 20:54:45.6 , 3.7 -02.0 699.4 , 33.0 +03.3 5.50 , 0.04 +0.01 0.047 119 1.035 15-04-25 1o3 OPD 23:45:28.1 , 1.3 +00.5 679.2 , 7.5 -02.3 6.05 , 0.05 -0.02 0.007 100 1.042 15-04-25 1o3 FOZ 23:45:26.7 , 3.3 -08.0 685.0 , 28.8 +03.8 5.95 , 0.04 +0.01 0.014 133 1.031 15-04-26 2o1 OPD 21:25:00.0 , 3.4 +01.8 584.6 , 27.6 +00.5 7.08 , 0.08 -0.07 0.029 064 1.067 15-04-29 3o1 OPD 00:29:06.9 , 1.2 -15.6 661.3 , 19.2 -04.0 6.92 , 0.06 -0.01 0.009 087 1.048 15-04-29 3o1 GOA 00:29:07.6 , 2.4 -10.5 663.8 , 22.2 -01.7 7.00 , 0.05 -0.04 0.026 095 1.044 15-05-03 2o1 OPD 23:39:19.6 , 2.2 -08.1 571.1 , 25.6 +11.5 6.83 , 0.08 -0.04 0.033 063 1.068 15-05-05 3o2 FOZ 21:54:22.0 , 2.4 +00.3 780.5 , 21.9 -01.7 5.30 , 0.04 -0.01 0.010 117 1.035 15-05-13 3o2 OPD 01:13:50.0 , 3.0 +06.5 593.2 , 27.1 +03.5 5.04 , 0.04 -0.01 0.038 137 1.030 15-06-04 2o1 FOZ 21:55:27.9 , 0.2 -01.8 160.8 , 6.3 -06.8 6.89 , 0.05 -0.02 0.009 094 1.044 15-06-18 3o1 GOA 21:01:51.0 , 1.1 +01.0 237.6 , 18.6 +10.7 3.82 , 0.03 -0.02 0.017 215 1.019 15-03-02 3e5 1.60 23:17:06.0 , 2.3 -22.5 391.3 , 76.3 +20.9 8.44 , 0.53 -0.09 0.101 060 1.080 Note: Similar as the note in the Table 3 for the mutual phenomena campaign of the 2014-2015. B. Morgado, R. Vieira-Martins, M. Assafin et al.: Preprint submitted to Elsevier Page 14 of 11

Journal

AstrophysicsarXiv (Cornell University)

Published: Sep 11, 2019

There are no references for this article.