# Photodynamical analysis of the nearly resonant planetary system WASP-148. Accurate transit-timing variations and mutual orbital inclination

J. M. Almenara, G. Hébrard, R. F. Díaz, J. Laskar, A. C. M. Correia, D. R.

Anderson, I. Boisse, X. Bonfils, D. J. A. Brown, V. Casanova, et al.

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

► To cite this version:

J. M. Almenara, G. Hébrard, R. F. Díaz, J. Laskar, A. C. M. Correia, et al. Photodynamical analysis of the nearly resonant planetary system WASP-148. Accurate transit-timing variations and mutual orbital inclination. Astronomy and Astrophysics - A&A, 2022, 663, 10.1051/0004-6361/202142964. insu-03748255

HAL Id: insu-03748255

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

Submitted on 9 Aug 2022

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

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

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

Distributed under a Creative Commons Attribution 4.0 International License
# Photodynamical analysis of the nearly resonant planetary system WASP-148

## Accurate transit-timing variations and mutual orbital inclination

J. M. Almenara<fnref n="1" />,², G. Hébrard³,⁴, R. F. Díaz⁵, J. Laskar⁶, A. C. M. Correia⁷,⁶, D. R. Anderson⁸,⁹,¹⁰, I. Boisse¹¹, X. Bonfils², D. J. A. Brown⁸,⁹, V. Casanova¹², A. Collier Cameron¹³, M. Fernández¹², J. M. Jenkins¹⁴, F. Kiefer³,¹⁵, A. Lecavelier des Étangs³, J. J. Lissauer¹⁴, G. Maciejewski¹⁶, J. McCormac⁸,⁹, H. Osborn⁸,⁹,¹¹,¹⁷, D. Pollacco⁸,⁹, G. Ricker¹⁸, J. Sánchez¹², S. Seager¹⁸,¹⁹,²⁰, S. Udry<fnref n="1" />, D. Verilhac²¹, and J. Winn²²

*(Affiliations can be found after the references)*

Received 21 December 2021 / Accepted 6 April 2022

## ABSTRACT

WASP-148 is a recently announced extra-solar system harbouring at least two giant planets. The inner planet transits its host star. The planets travel on eccentric orbits and are near the 4:1 mean-motion resonance, which implies significant mutual gravitational interactions. In particular, this causes transit-timing variations of a few minutes, which were detected based on ground-based photometry. This made WASP-148 one of the few cases where such a phenomenon was detected without space-based photometry. Here, we present a self-consistent model of WASP-148 that takes into account the gravitational interactions between all known bodies in the system. Our analysis simultaneously fits the available radial velocities and transit light curves. In particular, we used the photometry secured by the Transiting Exoplanet Survey Satellite (TESS) and made public after the WASP-148 discovery announcement. The TESS data confirm the transit-timing variations, but only in combination with previously measured transit times. The system parameters we derived agree with those previously reported and have a significantly improved precision, including the mass of the non-transiting planet. We found a significant mutual inclination between the orbital planes of the two planets: $I = 41.0^{+6.2}_{-7.6}^{\circ}$ based on the modelling of the observations, although we found $I = 20.8 \pm 4.6^{\circ}$ when we imposed a constraint on the model enforcing long-term dynamical stability. When a third planet was added to the model – based on a candidate signal in the radial velocity – the mutual inclination between planets b and c changed significantly allowing solutions closer to coplanar. We conclude that more data are needed to establish the true architecture of the system. If the significant mutual inclination is confirmed, WASP-148 would become one of the only few candidate non-coplanar planetary systems. We discuss possible origins for this misalignment.

**Key words.** stars: individual: WASP-148 – planetary systems – techniques: photometric – techniques: radial velocities

## 1. Introduction

While the orbit of a single planet around its host star is well reproduced by a Keplerian model with a constant orbital period, multi-planetary systems have more complex orbits. Indeed, the mutual gravitational interactions between planets imply small deviations from the Keplerian orbits and in particular slight variations of the orbital periods. Such effects are negligible for most systems, but they are amplified when orbital periods are exactly or nearly commensurable. When such a resonant or nearly resonant system includes transiting planets, their orbital periods can be accurately measured, allowing their variations to be detected. This was predicted in particular by Holman & Murray (2005) and Agol et al. (2005), who showed how interactions in multi-planetary systems might cause transit-timing variations (TTVs).

When they are measured, TTVs are a powerful tool to characterise planetary systems. In particular, they allow constraints to be put on masses as well as on the presence of additional planets. This was done on the first system in which TTVs were detected (Holman et al. 2010) as well as on dozens of other detections reported thereafter (e.g. Lissauer et al. 2011; Nesvorný et al. 2013; Wang et al. 2014; Hadden & Lithwick 2014; Gillon et al.

2017; Freudenthal et al. 2019; Jontof-Hutter et al. 2021). Among different methods employed to analyse TTV observations, the photodynamical modelling<fnref n="1" /> (Carter et al. 2011) of light curves can be used together with radial velocities to constrain planetary system parameters without using external inputs on the masses or radii of the host stars often derived from stellar evolution models (Agol et al. 2005). One of the strengths of this technique is that it fits the full transit light curves (including transit durations) instead of only using the transit timings. Photodynamical analyses have been used to characterise several TTV systems (e.g. Almenara et al. 2018a,b). Here, we apply it to the system WASP-148 recently reported by Hébrard et al. (2020).

Most of the confirmed TTV detections to date have been discovered using light curves obtained from space telescopes (namely *Kepler* or *Spitzer*). On the other hand, the TTVs in the WASP-148 system were detected using ground-based telescopes only. The WASP-148 system includes (at least) two giant planets. The inner one, WASP-148 b, is a hot Saturn of $0.72 \pm 0.06 R_J$.

<footnote n="1">¹ A photodynamical model is a light curve model for more than a two-body system that includes the gravitational interactions between the assumed bodies in the system, as the planet-planet interactions (and not just the planet-star interaction) in a multi-planetary system.</footnote>
and $0.29 \pm 0.03 M_J$ that transits its host with an orbital period of 8.80 days. The outer planet, WASP-148 c, has an orbital period of 34.5 days and a minimum mass of $0.40 \pm 0.05 M_J$ and a true mass <0.60 $M_J$. It was discovered and characterised using radial velocities obtained with the SOPHIE spectrograph. No transits of this planet were detected. The orbits of both planets have significant eccentricities ($e_b = 0.22 \pm 0.06$ and $e_c = 0.36 \pm 0.09$) and their orbital periods fall near the 1:4 mean-motion resonance. This particular configuration induces amplified dynamical effects, and TTVs were detected with an amplitude of about $\pm 15$ min and a period of roughly 460 days (Hébrard et al. 2020).

Several analyses of the available data sets have been presented by Hébrard et al. (2020). The first one simultaneously fitted transit light curves and radial velocities, but it did not include any mutual interactions. The model allowed for small, artificial, ad hoc shifts in the transit times in order to reproduce the TTVs. The second model was fitted to the radial velocities taking mutual interactions into account, constrained by the average period and phase derived from the transit light curves. Both analyses gave similar results and the derived system properties imply large TTVs for both planets, whereas only those of WASP-148 b could actually be observed as no transits of WASP-148 c were detected. However, Hébrard et al. (2020) did not present a complete model fitted to both radial velocity and transit light curve data sets that takes the gravitational interactions between the planets into account. Finally, Hébrard et al. (2020) indicate that there is a hint of a possible third planet in the system with a period near 150 days and a minimum mass around $0.25 M_J$, without confirming that however.

Following the detection and characterisation of the WASP-148 system by Hébrard et al. (2020), observations of that star secured with the Transiting Exoplanet Survey Satellite (TESS) were released. TESS provides continuous, high-quality photometry over 28 days (Ricker et al. 2015). In the case of WASP-148, it observed seven new transits of WASP-148 b, and provided a unique opportunity to search for possible transits of WASP-148 c or potential additional planets in the system. Maciejewski et al. (2020) published an analysis of the TTVs and radial velocities of WASP-148 presented by Hébrard et al. (2020) in addition to TESS data, as well as two new ground-based transit observations. They took into account the gravitational interactions between the planets, but they used only the transit timings instead of the whole transit light curves. Recently, Wang et al. (2022) have shown the orbit of WASP-148 b is aligned and prograde.

Here, we present new analyses of the WASP-148 system, applying the full photodynamical approach on the available transit light curves and radial velocities. The article is organised as follows: in Sect. 2 we describe the data used; in Sect. 3 we determine the stellar parameters; in Sects. 4 and 5 we analyse the radial velocity and transits without accounting for the gravitational interactions between the planets, respectively; in Sect. 6 we detail the photodynamical modelling; and in Sect. 7 we present the results. Finally, we discuss the results of our work in Sect. 8.

## 2. Observations

We used the ground-based photometry and SOPHIE radial velocities presented in Hébrard et al. (2020). We converted the time of the photometry from $BJD_{UTC}$ to $BJD_{TDB}^2$ (Eastman et al. 2010), except for the Telescopio Carlos Sánchez

transit, and NITES observation on June 13, 2016, which were already in $BJD_{TDB}$. In addition to the data in Hébrard et al. (2020), we added the photometry from TESS and four transits observed with the 1.5 m Ritchey-Chrétien Telescope at the Sierra Nevada Observatory (OSN150), the first two of which have been presented in Maciejewski et al. (2020).

### 2.1. TESS

TESS observed WASP-148 in sectors 24, 25, and 26, with a total time span of 79.3 days (TIC 115524421, TOI-2064). Two transits were lost during the interruption of the observations that occur every TESS orbit (around 14 days) at the middle of sectors 25 and 26. Each TESS sector lasts two orbits of the satellite. Around the perigee of the TESS orbit, data collection is paused. In total, seven new transits of WASP-148 b were observed. The photometry is available in the full-frame images (FFIs) at 30-min cadence. The FFIs were calibrated by the Science Processing Operations Center (SPOC) at NASA Ames Research Center (Jenkins et al. 2016). We used *eleanor* (Feinstein et al. 2019) to extract the light curve of WASP-148 from the TESS FFIs. We chose the point spread function photometry that has the lower dispersion for this object. *eleanor* corrects the times for the object coordinates, which otherwise are set to the centre of the CCD in FFIs. There is a star 4.5 mag fainter in the Gaia G-band at 26" of WASP-148 that partially contaminates the TESS photometry<fnref n="3" />. This is taken into account in our modelling (Sects. 5 and 6). The TESS data show no evidence for any transit of WASP-148 c, which is in agreement with Maciejewski et al. (2020).

### 2.2. OSN150

Two new precise photometric time series for transits of WASP-148 b were acquired in March and June 2021 using the 1.5 m Ritchey-Chrétien telescope (OSN150) at the Sierra Nevada Observatory (OSN, Spain). The instrument was equipped with a Roper Scientific VersArray 2048B CCD camera with a $2048 \times 2048 \times 13.5 \ \mu m$ back-illuminated matrix. The field of view was $7.92 \times 7.92$ with the pixel scale of 0.232 per pixel. The instrument was mildly defocussed to allow for longer exposure times and a lower fraction of time lost for CCD readout. The observations were gathered without any filter to increase the signal-to-noise ratio for transit-timing purposes. The observing runs were scheduled following an ephemeris from Maciejewski et al. (2020). The telescope was auto-guided to keep the star at the same position in the CCD matrix. The details on the individual runs are given in Table 1. The observations started about 90 min before a transit ingress and lasted about 90 min after an egress. This out-of-transit monitoring was secured for detrending purposes. On March 22, 2021 the observations were stopped about 12 min after the transit due to dawn.

Photometric data reduction was performed with the AstroImageJ software (Collins et al. 2017) following a standard calibration procedure. The science frames were de-biased and flat-field calibrated using sky flat frames. The light curves were generated with the differential aperture photometry method. The aperture size and a collection of comparison stars were set after

<footnote n="2">² $BJD_{TDB} = BJD_{UTC} + \delta$, with $\delta$ being between 65.184 and 69.184 s depending on the date of the transit.</footnote>

<footnote n="3">³ Assuming the flux of the contaminant star is completely inside the TESS aperture and that the difference in magnitude in the *Gaia* band is the same as in the TESS band, the dilution is about 1.6%, which can change the planet to star radius ratio by 0.8%, and the final planet radius by ~0.4–$\sigma$.</footnote>
<table><caption>**Table 1.** Details on the observing runs.</caption><thead><tr><th>Date UT</th><th>UT start-end</th><th>X</th><th>N<sub>obs</sub></th><th>t<sub>exp</sub> (s)</th><th>Γ</th><th>pnr (ppth)</th></tr></thead><tbody><tr><td>2021 March 22</td><td>00:20–05:34</td><td>1.70 → 1.01</td><td>515</td><td>30</td><td>1.69</td><td>1.10</td></tr><tr><td>2021 June 26</td><td>21:19–03:00</td><td>1.09 → 1.01 → 1.01 → 1.47</td><td>575</td><td>30</td><td>1.69</td><td>0.76</td></tr></tbody></table>

**Notes.** Date UT is given for the beginning of an observing run. X shows the changes of the target’s airmass during a run. N<sub>obs</sub> is the number of useful scientific exposures. t<sub>exp</sub> is the exposure time. Γ is the median number of exposures per minute. pnr is the photometric noise rate (Fulton et al. 2011) in parts per thousand (ppth) of the normalised flux per minute of the observation.

a series of test runs to minimise the data point scatter. The fluxes were de-trended against airmass and time using out-of-transit data only and then they were normalised to unity outside the transits. Timestamps were converted into barycentric Julian dates in barycentric dynamical time BJD<sub>TDB</sub>. For the homogeneity of our analysis, we also reprocessed the two transit light curves from Maciejewski et al. (2020) following the same de-trending procedure.

### 3. Stellar parameters

To determine the stellar parameters of WASP-148, we used the precise parallax determination by Gaia (Gaia Collaboration 2016; Gaia Collaboration 2021), which was not exploited by Hébrard et al. (2020). Stellar atmosphere models and stellar evolution models are also required to model the observed spectral energy distribution (SED). We constructed the SED of WASP-148 using the magnitudes from Gaia Early Data Release 3 (Gaia EDR3, Riello et al. 2021), the 2-Micron All-Sky Survey (2MASS, Skrutskie et al. 2006; Cutri et al. 2003), and the Wide-field Infrared Survey Explorer (WISE, Wright et al. 2010; Cutri & et al. 2013). The measurements are listed in Table A.1. We modelled these magnitude measurements using the procedure described by Díaz et al. (2014), with informative priors for the effective temperature (T<sub>eff</sub>), surface gravity (log g), and metallicity ([Fe/H]) from Hébrard et al. (2020), and for the distance from Gaia EDR3 (Lindegren et al. 2021). We used non-informative priors for the rest of parameters. The priors are listed in Table 2.

We decided to use an additive jitter for each set of photometric bands (Gaia, 2MASS, and WISE), which had the effect of slightly broadening the posteriors of the stellar parameters. We used the two stellar atmosphere models, PHOENIX/BT-Settl (Allard et al. 2012) and ATLAS/Castelli & Kurucz (Castelli & Kurucz 2003), and two stellar evolution models, Dartmouth (Dotter et al. 2008) and PARSEC (Chen et al. 2014). We obtained posterior samples for the four combinations of stellar atmosphere models and stellar evolution models using the Markov chain Monte Carlo (MCMC) algorithm from Díaz et al. (2014). The posteriors of the stellar parameters of the individual combinations agree within 1σ (see Fig. 1). We merged the results assuming an equal probability for each model combination (labelled as ‘merged’ in Fig. 1). The posteriors’ median and 68.3% credible intervals (CI) for model parameters and for derived physical quantities of interest are listed in Table 2. Those results agree with the ones reported by Hébrard et al. (2020), but they are more accurate, mainly due to the use of the parallax from Gaia.

The data with the maximum a posteriori (MAP) stellar atmosphere model is shown in Fig. 2. Before Gaia EDR3, we carried out the same analysis with Gaia DR2 (Gaia Collaboration 2018), obtaining similar results. In particular, we obtained the same stellar radius error, despite a factor ~3 increase in the precision

<table><caption>**Table 2.** Modelling of the spectral energy distribution.</caption><thead><tr><th>Parameter</th><th>Prior</th><th>Posterior median and 68.3% CI</th></tr></thead><tbody><tr><td>Effective temperature, T<sub>eff</sub> [K]</td><td>N(5460, 130)</td><td>5555 ± 90</td></tr><tr><td>Surface gravity, log g [cgs]</td><td>N(4.40, 0.15)</td><td>4.490<sup>+0.027</sup><sub>-0.034</sub></td></tr><tr><td>Metallicity, [Fe/H] [dex]</td><td>N(0.11, 0.08)</td><td>0.099 ± 0.078</td></tr><tr><td>Distance [pc]</td><td>N(247.73, 0.60)</td><td>247.73 ± 0.60</td></tr><tr><td>E<sub>(B-V)</sub> [mag]</td><td>U(0, 3)</td><td>0.026<sup>+0.041</sup><sub>-0.019</sub></td></tr><tr><td>Jitter Gaia [mag]</td><td>U(0, 1)</td><td>0.140<sup>+0.22</sup><sub>-0.080</sub></td></tr><tr><td>Jitter 2MASS [mag]</td><td>U(0, 1)</td><td>0.056<sup>+0.072</sup><sub>-0.030</sub></td></tr><tr><td>Jitter WISE [mag]</td><td>U(0, 1)</td><td>0.033<sup>+0.076</sup><sub>-0.024</sub></td></tr><tr><td>Radius, R<sub>⋆</sub> [R<sub>☉</sub>]</td><td></td><td>0.921 ± 0.016</td></tr><tr><td>Mass, M<sub>⋆</sub> [M<sub>☉</sub>]</td><td></td><td>0.958 ± 0.048</td></tr><tr><td>Density, ρ<sub>⋆</sub> [g cm<sup>-3</sup>]</td><td></td><td>1.73 ± 0.15</td></tr><tr><td>Isochronal age [Gyr]</td><td></td><td>3.6<sup>+4.0</sup><sub>-2.7</sub></td></tr><tr><td>Luminosity [L<sub>☉</sub>]</td><td></td><td>0.726 ± 0.036</td></tr></tbody></table>

**Notes.** N(μ,σ): Normal distribution prior with mean μ, and standard deviation σ. U(l,u): Uniform distribution prior in the range [l, u].

<figure><img src="image_3.png" /><figcaption>**Fig. 1.** Stellar parameters derived in this work compared with Hébrard et al. (2020), the TESS Input Catalog (TIC) version 8 (Stassun et al. 2019), and Gaia DR2 (Andrae et al. 2018). We finally adopt the values labelled as ‘merged’.</figcaption></figure>

of the parallax. Thus, at least for WASP-148, the precision in the parallax is not the limiting factor to improve the stellar radius determination with the SED technique.

We used the stellar rotation period, P<sub>rot</sub> = 26.2 ± 1.3 days, derived in Hébrard et al. (2020) and the stellar mass derived in this section (Table 2) to estimate a gyrochronological age, neglecting the influence of the planets, of 4.0<sup>+0.9</sup><sub>-0.7</sub> Gyr (Barnes 2010; Barnes & Kim 2010, using a P<sub>0</sub> between 0.12 and 3.4 days), where we added a systematic 10% error to the statistical one (Meibom et al. 2015). The isochronal age (Table 2) agrees with the gyrochronological age within the uncertainties, but it is less precise.
<figure><img src="image_4.png" /><figcaption>**Fig. 2.** Spectral energy distribution of WASP-148. Top panel: solid line is the MAP PHOENIX/BT-Settl interpolated synthetic spectrum, red circles are the absolute photometric observations, and grey open circles are the result of integrating the synthetic spectrum in the observed band-passes. Bottom panel: residuals of the MAP model (the jitter has been added quadratically to the data error bars).</figcaption></figure>

## 4. SOPHIE-only analysis

When radial velocity and transit observations are modelled simultaneously (as was done by Hébrard et al. 2020), both data sets constrain the eccentricity. In this section, we aim for a planetary eccentricity determination using only radial velocity data. We analysed the SOPHIE radial velocities with a two-planet Keplerian model, without taking mutual gravitational interactions into account and using juliet (Espinoza et al. 2019), radvel (Fulton et al. 2018), and dynasty (Speagle 2020). We used normal priors from Hébrard et al. (2020) for the period and the time of inferior conjunction of planet b, and non-informative priors for the rest of the parameters. We included a jitter parameter that was added quadratically to the velocity uncertainties. We performed a second analysis using a Gaussian process (GP) regression model with a quasi-periodic (QP) kernel to model the error terms. We employed the QP kernel included in celerite (Foreman-Mackey et al. 2017), with non-informative priors for the hyperparameters. We found an odds ratio of $ln(Z_{jitter}/Z_{jitter+QP}) = 0.5 \pm 0.4$. Neither of the models is strongly favoured over the other (Kass & Raftery 1995), and we therefore adopted the simpler model of uncorrelated error terms with an additive jitter. The posterior median and 68.3% credible interval are shown in Table 3. They agree with the results presented by Hébrard et al. (2020).

We investigated the significance of the main peak in the residuals of the two-planet model (Fig. 3 in Hébrard et al. 2020), thus we repeated the analysis adding a third Keplerian with a period of ~150 days. We found a minimum mass of $0.24 \pm 0.05 M_J$, similar to the $0.25 M_J$ value reported by Hébrard et al. (2020). The model comparison gives an odds ratio of $ln(Z_{3 planets}/Z_{2 planets}) = 1.6 \pm 0.4$, which is positive evidence in favour of the three-planet model but below the strong evidence cutoff to be preferred over the two-planet model, so we continued with the adopted two-planet model.

## 5. Transit-only analysis

We derived the transit times of planet b without considering planet c or the radial velocities. This modelling is independent

<table><caption>**Table 3.** Two-Keplerian fit to the SOPHIE data.</caption><thead><tr><th>Keplerian</th><th></th><th>WASP-148 b</th><th>WASP-148 c</th></tr></thead><tbody><tr><td>P</td><td>[d]</td><td>8.803809 ± 0.000043</td><td>34.524 ± 0.029</td></tr><tr><td>T<sub>c</sub></td><td>[BJD]</td><td>2 457 957.4876 ± 0.0060</td><td>2 457 935.6 ± 1.2</td></tr><tr><td>e</td><td></td><td>0.183 ± 0.070</td><td>0.352 ± 0.085</td></tr><tr><td>ω</td><td>[°]</td><td>59 <sup>+15</sup><sub>-20</sub></td><td>9 <sup>+17</sup><sub>-14</sub></td></tr><tr><td>K</td><td>[m s<sup>-1</sup>]</td><td>28.8 ± 2.0</td><td>27.0 ± 2.9</td></tr><tr><td>γ<sub>SOPHIE</sub></td><td>[km s<sup>-1</sup>]</td><td colspan="2">-5.6174 ± 0.0014</td></tr><tr><td>jitter</td><td>[m s<sup>-1</sup>]</td><td colspan="2">10.9 ± 1.3</td></tr></tbody></table>

**Notes.** The parameters are: orbital period, time of conjunction, eccentricity, argument of pericentre, radial velocity semi-amplitude, systemic velocity, and jitter.

of the two-planet system hypothesis adopted in the photodynamical modelling (Sect. 6), and it can be used to verify if it is an appropriate hypothesis given the data.

In Hébrard et al. (2020), the error term of the transit photometry data were modelled using a simple additive jitter term. However, if systematics are not correctly accounted for, the posteriors of the transit modelling can be biased (Barros et al. 2013). This should be particularly severe for incomplete transits. Here, we reanalyse the transits presented in Hébrard et al. (2020) and the TESS observations using a more sophisticated error model.

To model the transits of planet b simultaneously with the systematics, we used juliet (Espinoza et al. 2019), using batman (Kreidberg 2015) for the transit model, and we chose the approximate Matern kernel GP included in celerite (Foreman-Mackey et al. 2017). We used different GPs for each ground-based transit, and for each TESS sector. The timing of each individual transit is a free parameter. We used a prior for the stellar density from Sect. 3, and for $\sqrt{e} \cos \omega$ and $\sqrt{e} \sin \omega$ from Hébrard et al. (2020). Otherwise, we adopted non-informative or large priors for the rest of the parameters. We oversampled (Kipping 2010a) the model of the MARS transit by a factor 3, which has an 88 s cadence, and the TESS data by a factor 30, which have about a 30 min cadence. The remaining transit observations have cadences between 0.9 and 37 s, and we deemed them unnecessary to oversample their model. We set a dilution factor for the wide-field observations: WASP (one dilution for the three WASP transits), and TESS (one dilution per sector). To reduce the number of free parameters, we made the choice to use one set of limb darkening parameters for all transit observations without a filter ('clear'), but they certainly do not correspond to exactly the same instrument bandpass. In total, there are 76 free parameters. To sample from the posterior, we used the nested sampling code dynasty (Speagle 2020). The complete list of parameters, priors, and posteriors are shown in Table A.2. Figure 3 shows the data and the model posterior. We compare the timing of the ground-based observations with the results in Hébrard et al. (2020), and found differences within $1\sigma$, except for the partial transit observed with MARS where it is $12.0 \pm 2.5$ min later (see Fig. 4). The transit normalisation can affect the transit-timing determination (Barros et al. 2013). The MARS transit can be more affected because it is a partial transit and the baseline is not well-defined.

We repeat this analysis for the four OSN150 transits, which were included at a later stage in this work. The results are presented in Table A.3 and Fig. 5.

This analysis could be used to de-trend the transits for the photodynamical modelling. However, in this process, one set of
<figure><img src="image_5.png" /><figcaption>**Fig. 3.** Modelling of the planet b transits presented in Hébrard et al. (2020) and observed by TESS (blue error bars). Each panel is labelled with the observatory or instrument name and sequential night of observation if there were more than one. Black lines and intervals in grey show the model median and 68.3% credible interval computed with 1000 random samples of the posterior. The mean of the predictive distribution of the kernel model is shown in orange. On top of each panel, the transit-timing posterior is shown for Hébrard et al. (2020) (red line histogram), the modelling with juliet (Sect. 5, grey histogram), and the photodynamical modelling (blue points with black error bars, although the errors are barely visible).</figcaption></figure>

parameters of the transit model needs to be fixed, including a transit time. Therefore, if the de-trended transits are then modelled, the results can be biased. Instead, we decided to include the GP in the photodynamical modelling.

**6. Photodynamical modelling**

While Sects. 4 and 5 do not take mutual interactions into account, here we report our fits of the observed photometry and radial velocity measurements accounting for the gravitational interactions between the three bodies known in the system using a photodynamical model. Its positions and velocities in time were obtained through an n-body integration. The sky-projected positions were used to compute the light curve (Mandel & Agol 2002) using a quadratic limb-darkening law (Manduca et al. 1977), which we parametrised following Kipping (2013). To account for the integration time, the model was oversampled by a factor of 30 and 3 for the TESS and MARS data, respectively, and then binned back to match the cadence of the data points (Kipping 2010a). The line-of-sight projected velocity of the star issued from the n-body integration was used to model the radial velocity measurements. We used the n-body code REBOUND (Rein & Liu 2012) with the WHFast integrator (Rein & Tamayo 2015) and an integration step of 0.01 days, which results

<figure><img src="image_6.png" /><figcaption>**Fig. 4.** Comparison between transit-time posteriors of Hébrard et al. (2020) (red line histograms), the modelling with juliet (Sect. 5, grey histograms), and the photodynamical modelling (Sect. 6, blue points with black error bars and blue histograms).</figcaption></figure>
<figure><img src="image_7.png" /><figcaption>**Fig. 5.** Idem as Fig. 3, but for the OSN150 transits.</figcaption></figure>

in a maximum error of ~20 ppm for the photometric model. The light-time effect (Irwin 1952) is included, although with an amplitude of ~0.05 s in the transit timing; this is a negligible effect for this system. The model was parametrised using the stellar mass and radius, planet-to-star mass ratios, planet b-to-star radius ratio, and Jacobi orbital elements (Table 4) at the reference time, $t_{\text{ref}} = 2\,457\,957.48167$ BJD$_{\text{TDB}}$, close to the RISE transit centre. Due to the symmetry of the problem, we fixed the longitude of the ascending node of the interior planet $\Omega_b$ at $t_{\text{ref}}$ to 180°, and we limited the inclination of the outer one $i_c > 90°$. We used a GP regression model, with an approximate Matern kernel (celerite, Foreman-Mackey et al. 2017) for the model of the error terms of the transit light curves. We used different kernel hyperparameters for each transit, except for the TESS data for which we used different kernel hyperparameters for each sector. We added one dilution factor<fnref n="4" /> for each TESS sector, and another one for the WASP transits. For each photometric data set, we added a transit normalisation factor and an additive jitter parameter. For the radial velocity, we added a systemic radial velocity and an additive jitter parameter. In total, the model has 90 free parameters. We used normal priors for the stellar mass and radius from Sect. 3, a non-informative sinusoidal prior for the orbital inclinations, and non-informative uniform prior distributions for the rest of the parameters. The joint posterior distribution was sampled using the emcee algorithm (Goodman & Weare 2010; Foreman-Mackey et al. 2013) with 1000 walkers with starting points based on the results of Hébrard et al. (2020) and Sect. 5.

## 7. Results

In Table 4 we list the prior, the median, and the 68% credible interval of the inferred system parameters' marginal distributions. The one- and two-dimensional projections of the posterior sample are shown in Fig. A.1. The MAP model is plotted in Figs. A.2 and A.3. Figure 6 shows the posterior of the planet orbits. With an inferred impact parameter of $21.2 \pm 4.5$, the transits of planet c are highly disfavoured, which is in agreement with the null result of the transit search in Maciejewski et al. (2020). The derived parameters agree with the ones reported by Hébrard et al. (2020), but they have a significantly improved precision. A difference, however, is the slightly larger radius ratio $R_p/R_\star$.

<figure><img src="image_8.png" /><figcaption>**Fig. 6.** Orbital projections. Top: orbital projections around the time of the RISE transit for planet b (blue) and planet c (orange). The origin is the system barycentre, and the orbits are projected in the sky plane seen by the observer (X–Y, left) and X–Z plane (right, system top view, the movement is clockwise, the positive Z-axis points towards the observer). A thousand random orbits were drawn from the posterior samples, and the MAP is shown as a black orbit. The black points mark the position of the star (size to scale) and the planets at the central time of the RISE transit for the MAP (the size of planet b is enlarged by a factor of 10, and the size of planet c is not known). Bottom: Idem as for 1000 stable solutions (Sect. 7.5).</figcaption></figure>

Instead of $M_p \sin i$, the photodynamical modelling allows the true mass of planet c to be measured at $M_c = 0.424 \pm 0.046 M_J$, which is in agreement with the dynamical upper limit of $0.60 M_J$ reported by Hébrard et al. (2020). Most of the posterior samples have the apsces of the orbits that librate around alignment. From tests we determined that the improved precision in model posteriors over the earlier study of Hébrard et al. (2020) is due to both the additional data and the photodynamical modelling that account for the interactions between the planets.

### 7.1. Transit-timing variations

Figures 7 and A.4 show the posterior TTVs of planet b, obtained from the time of conjunction. The posterior transit times, which

<footnote n="4">⁴ The definition of the dilution factor and transit normalisation factor is different from the one in *juliet*. For the photodynamical model, we used $M = f_0(f_i T + f_c)/(f_i + f_c)$, with M being the model to be compared with the observations, $f_0$ being the transit normalisation factor (the light curve flux level out of transit), $f_i$ being the flux of the target star in the aperture, T being the transit model, and $f_c$ being the flux of the contaminant star in the aperture. The dilution factor is defined as $d = f_c/(f_i + f_c)$.</footnote>
<table><caption>**Table 4.** Inferred system parameters.</caption><thead><tr><th>Parameter</th><th>Units</th><th>Prior</th><th>Median and 68.3% CI</th><th>Stable MAP</th><th>Stable median and 68.3% CI</th></tr></thead><tbody><tr><td colspan="6"><strong>Star</strong></td></tr><tr><td>Stellar mass, M<sub>★</sub></td><td>[M<sub>☉</sub>]</td><td>N(0.958, 0.048)</td><td>0.967 ± 0.049</td><td>0.954</td><td>0.954<sup>+0.030</sup><sub>-0.052</sub></td></tr><tr><td>Stellar radius, R<sub>★</sub></td><td>[R<sup>N</sup><sub>☉</sub>]</td><td>N(0.921, 0.016)</td><td>0.920 ± 0.016</td><td>0.913</td><td>0.912 ± 0.014</td></tr><tr><td>Stellar mean density, ρ<sub>★</sub></td><td>[g cm<sup>-3</sup>]</td><td></td><td>1.76 ± 0.12</td><td>1.77</td><td>1.76 ± 0.12</td></tr><tr><td>Surface gravity, log g</td><td>[cgs]</td><td></td><td>4.497 ± 0.025</td><td>4.497</td><td>4.496<sup>+0.022</sup><sub>-0.028</sub></td></tr><tr><td>Kipping (2013) q<sub>1</sub>, q<sub>2</sub> CLEAR</td><td></td><td>U(0, 1)</td><td>0.42<sup>+0.13</sup><sub>-0.11</sub>, 0.240<sup>+0.11</sup><sub>-0.094</sub></td><td>0.304, 0.421</td><td>0.388<sup>+0.15</sup><sub>-0.083</sub>, 0.279<sup>+0.093</sup><sub>-0.16</sub></td></tr><tr><td>Kipping (2013) q<sub>1</sub>, q<sub>2</sub> Johnson-R</td><td></td><td>U(0, 1)</td><td>0.230<sup>+0.13</sup><sub>-0.096</sub>, 0.48<sup>+0.32</sup><sub>-0.26</sub></td><td>0.106, 0.98</td><td>0.162<sup>+0.16</sup><sub>-0.056</sub>, 0.59<sup>+0.28</sup><sub>-0.32</sub></td></tr><tr><td>Kipping (2013) q<sub>1</sub>, q<sub>2</sub> RISE</td><td></td><td>U(0, 1)</td><td>0.80 ± 0.11, 0.269 ± 0.062</td><td>0.802, 0.272</td><td>0.802<sup>+0.091</sup><sub>-0.13</sub>, 0.278<sup>+0.056</sup><sub>-0.063</sub></td></tr><tr><td>Kipping (2013) q<sub>1</sub>, q<sub>2</sub> TESS</td><td></td><td>U(0, 1)</td><td>0.38<sup>+0.21</sup><sub>-0.15</sub>, 0.40<sup>+0.33</sup><sub>-0.24</sub></td><td>0.26, 0.72</td><td>0.38<sup>+0.17</sup><sub>-0.12</sub>, 0.40<sup>+0.31</sup><sub>-0.20</sub></td></tr><tr><td colspan="6"><strong>Planet b</strong></td></tr><tr><td>Semi-major axis, a</td><td>[au]</td><td></td><td>0.0825 ± 0.0014</td><td>0.08214</td><td>0.08215<sup>+0.00086</sup><sub>-0.0015</sub></td></tr><tr><td>Eccentricity, e</td><td></td><td></td><td>0.214<sup>+0.021</sup><sub>-0.018</sub></td><td>0.208</td><td>0.208<sup>+0.020</sup><sub>-0.025</sub></td></tr><tr><td>Argument of pericentre, ω</td><td>[°]</td><td></td><td>63.7 ± 8.3</td><td>62.6</td><td>60.6<sup>+5.0</sup><sub>-6.0</sub></td></tr><tr><td>Inclination, i</td><td>[°]</td><td>S(0, 180)</td><td>89.61<sup>+0.21</sup><sub>-0.27</sub></td><td>89.32</td><td>89.30 ± 0.24</td></tr><tr><td>Longitude of the ascending node, Ω</td><td>[°]</td><td></td><td>180 (fixed at t<sub>ref</sub>)</td><td>180</td><td>180</td></tr><tr><td>Mean anomaly, M<sub>0</sub></td><td>[°]</td><td></td><td>16.7<sup>+5.8</sup><sub>-5.1</sub></td><td>17.8</td><td>19.2<sup>+4.0</sup><sub>-3.5</sub></td></tr><tr><td>√e cos ω</td><td></td><td>U(-1, 1)</td><td>0.204 ± 0.056</td><td>0.210</td><td>0.219 ± 0.041</td></tr><tr><td>√e sin ω</td><td></td><td>U(-1, 1)</td><td>0.413 ± 0.039</td><td>0.405</td><td>0.394<sup>+0.028</sup><sub>-0.031</sub></td></tr><tr><td>Mass ratio, M<sub>p</sub>/M<sub>★</sub></td><td></td><td>U(0, 1)</td><td>0.000284 ± 0.000019</td><td>0.000309</td><td>0.000292<sup>+0.000012</sup><sub>-0.000014</sub></td></tr><tr><td>Radius ratio, R<sub>p</sub>/R<sub>★</sub></td><td></td><td>U(0, 1)</td><td>0.08436 ± 0.00058</td><td>0.08532</td><td>0.08498<sup>+0.00064</sup><sub>-0.00080</sub></td></tr><tr><td>Scaled semi-major axis, a/R<sub>★</sub></td><td></td><td></td><td>19.31<sup>+0.41</sup><sub>-0.46</sub></td><td>19.34</td><td>19.33<sup>+0.39</sup><sub>-0.44</sub></td></tr><tr><td>Impact parameter, b</td><td></td><td></td><td>0.105<sup>+0.073</sup><sub>-0.056</sub></td><td>0.185</td><td>0.199<sup>+0.050</sup><sub>-0.070</sub></td></tr><tr><td>T'<sub>0</sub> - 2 450 000</td><td>[BJD<sub>TDB</sub>]</td><td>U(6957, 8957)</td><td>7957.48172 ± 0.00022</td><td>7957.48190</td><td>7957.48166<sup>+0.00025</sup><sub>-0.00015</sub></td></tr><tr><td>P'</td><td>[d]</td><td>U(0, 1000)</td><td>8.80369<sup>+0.00019</sup><sub>-0.00017</sub></td><td>8.80366</td><td>8.80354<sup>+0.00011</sup><sub>-0.00016</sub></td></tr><tr><td>K'</td><td>[m s<sup>-1</sup>]</td><td></td><td>29.7 ± 2.0</td><td>32.0</td><td>30.1<sup>+1.6</sup><sub>-1.5</sub></td></tr><tr><td>Planet mass, M<sub>p</sub></td><td>[M<sub>J</sub>]</td><td></td><td>0.288 ± 0.021</td><td>0.308</td><td>0.287<sup>+0.022</sup><sub>-0.016</sub></td></tr><tr><td>Planet radius, R<sub>p</sub></td><td>[R<sup>N</sup><sub>e</sub>]</td><td></td><td>0.756 ± 0.014</td><td>0.758</td><td>0.756<sup>+0.013</sup><sub>-0.017</sub></td></tr><tr><td>Planet mean density, ρ<sub>p</sub></td><td>[g cm<sup>-3</sup>]</td><td></td><td>0.829 ± 0.077</td><td>0.878</td><td>0.831<sup>+0.087</sup><sub>-0.069</sub></td></tr><tr><td>Planet surface gravity, log g<sub>p</sub></td><td>[cgs]</td><td></td><td>3.096<sup>+0.032</sup><sub>-0.036</sub></td><td>3.124</td><td>3.098<sup>+0.037</sup><sub>-0.033</sub></td></tr><tr><td colspan="6"><strong>Planet c</strong></td></tr><tr><td>Semi-major axis, a</td><td>[au]</td><td></td><td>0.2053 ± 0.0034</td><td>0.2044</td><td>0.2044<sup>+0.0021</sup><sub>-0.0038</sub></td></tr><tr><td>Eccentricity, e</td><td></td><td></td><td>0.228<sup>+0.014</sup><sub>-0.019</sub></td><td>0.1791</td><td>0.1809<sup>+0.018</sup><sub>-0.0072</sub></td></tr><tr><td>Argument of pericentre, ω</td><td>[°]</td><td></td><td>22.9<sup>+6.8</sup><sub>-5.1</sub></td><td>25.7</td><td>26.1 ± 9.4</td></tr><tr><td>Inclination, i</td><td>[°]</td><td>S(90, 180)</td><td>120.6 ± 7.3</td><td>106.7</td><td>104.9<sup>+4.6</sup><sub>-7.3</sub></td></tr><tr><td>Longitude of the ascending node, Ω</td><td>[°]</td><td>U(90, 270)</td><td>207.4<sup>+4.4</sup><sub>-6.4</sub></td><td>191.6</td><td>192.2<sup>+6.0</sup><sub>-2.7</sub></td></tr><tr><td>Mean anomaly, M<sub>0</sub></td><td>[°]</td><td></td><td>260.5 ± 3.5</td><td>260.4</td><td>258.0<sup>+5.0</sup><sub>-4.0</sub></td></tr><tr><td>√e cos ω</td><td></td><td>U(-1, 1)</td><td>0.437<sup>+0.025</sup><sub>-0.032</sub></td><td>0.381</td><td>0.382 ± 0.030</td></tr><tr><td>√e sin ω</td><td></td><td>U(-1, 1)</td><td>0.185<sup>+0.051</sup><sub>-0.037</sub></td><td>0.183</td><td>0.187 ± 0.062</td></tr><tr><td>Mass ratio, M<sub>p</sub>/M<sub>★</sub></td><td></td><td>U(0, 1)</td><td>0.000417 ± 0.000040</td><td>0.000406</td><td>0.000394 ± 0.000027</td></tr><tr><td>Scaled semi-major axis, a/R<sub>★</sub></td><td></td><td></td><td>48.0 ± 1.1</td><td>48.12</td><td>48.10<sup>+0.98</sup><sub>-1.1</sub></td></tr><tr><td>Impact parameter, b</td><td></td><td></td><td>21.2 ± 4.5</td><td>12.4</td><td>11.0<sup>+3.4</sup><sub>-5.3</sub></td></tr><tr><td>T'<sub>0</sub> - 2 450 000</td><td>[BJD<sub>TDB</sub>]</td><td>U(6957, 8957)</td><td>7971.27<sup>+0.33</sup><sub>-0.28</sub></td><td>7971.54</td><td>7971.62<sup>+0.31</sup><sub>-0.26</sub></td></tr><tr><td>P'</td><td>[d]</td><td>U(0, 1000)</td><td>34.5412 ± 0.0028</td><td>34.54626</td><td>34.54619<sup>+0.00094</sup><sub>-0.0016</sub></td></tr><tr><td>K'</td><td>[m s<sup>-1</sup>]</td><td></td><td>23.7 ± 2.2</td><td>25.47</td><td>24.85<sup>+1.5</sup><sub>-0.97</sub></td></tr><tr><td>Planet mass, M<sub>p</sub></td><td>[M<sub>J</sub>]</td><td></td><td>0.424 ± 0.046</td><td>0.406</td><td>0.392<sup>+0.023</sup><sub>-0.027</sub></td></tr><tr><td>M<sub>p</sub> sin i</td><td>[M<sub>J</sub>]</td><td></td><td>0.361 ± 0.036</td><td>0.389</td><td>0.378<sup>+0.021</sup><sub>-0.024</sub></td></tr><tr><td>Mutual inclination, I</td><td>[°]</td><td></td><td>41.0<sup>+6.2</sup><sub>-7.6</sub></td><td>20.7</td><td>20.8 ± 4.6</td></tr></tbody></table>

**Notes.** The table lists: Prior, posterior median, and 68.3% credible interval (CI) for the photodynamical analysis (Sect. 6), maximum a posteriori (MAP) and median, and 68.3% CI for the stable samples (Sect. 7.5). The Jacobi orbital elements are given for the reference time $t_{\text{ref}} = 2 457 957.48167$ BJD$_{TDB}$. Additional parameters are in Table A.4. The ‘stable MAP’ solution is the nominal solution used in Sects. 7.5 and 7.7. We note that $P'$ and $T'_0$ should not be confused with the period or the time of conjunction, respectively, and they were only used to reduce the correlations between jump parameters, replacing the semi-major axis and the mean anomaly at $t_{\text{ref}}$. $T'_0 \equiv t_{\text{ref}} - \frac{P'}{2\pi}(M_0 - E + e \sin E)$ with $E = 2 \arctan\left\{\frac{\sqrt{1-e^2+e^2\tan[\frac{1}{2}(\frac{\pi}{2}-\omega)]}}{2}\right\}$, $P' \equiv \sqrt{\frac{4\pi^2 a^2}{GM_*}}$, $K' \equiv \frac{M_p \sin i}{M_*^2 \sqrt{1-e^2}})^{1/3}$. CODATA 2018: $G = 6.674 \times 10^{-11} \text{m}^3 \text{kg}^{-1} \text{s}^{-2}$. IAU 2012: au = 149 597 870 700 m. IAU 2015: $R_{\odot}^N = 6.957 \times 10^8 \text{m}$, $(GM)^N_{\odot} = 1.327 124 4 \times 10^{20} \text{m}^3 \text{s}^{-2}$, $R_{eJ}^N = 7.149 2 \times 10^7 \text{m}$, $(GM)^N_{eJ} = 1.266 865 3 \times 10^{17} \text{m}^3 \text{s}^{-2}$. $M_{\odot} = (GM)^N_{\odot}/G$, $M_{eJ} = (GM)^N_{eJ}/G$, $K^2 = (GM)^N_{\odot} (86 400 \text{s})^2/\text{au}^3$. $N(\mu, \sigma)$: normal distribution with mean $\mu$ and standard deviation $\sigma$. $U(a, b)$: a uniform distribution defined between a lower $a$ and upper $b$ limit. $S(a, b)$: A sinusoidal distribution defined between a lower $a$ and upper $b$ limit.
rely on the three-body system hypothesis, agree with the individually derived transit times (Sect. 5). In addition, the periodicity of the TTVs agrees with the ≈450 days super-period (Lithwick et al. 2012) for two planets near the 4:1 mean-motion resonance<fnref n="5" />. The periodicity of the TTVs also agrees with the one reported by Hébrard et al. (2020).

The posterior timing from the photodynamical modelling is quite wide at some epochs (lower panel of Figs. 7 and A.4). The uncertainty in the posterior transit times is related with the knowledge of the system parameters. Future transit observations should favour the epochs where the posterior transit uncertainty is large to further improve the characterisation of the system. There is room for improvement with ~1 min transit-timing precision observations. Predictions of transit times up to 2026 are listed in Table A.5.

If only the TESS observations are considered, TTVs would not have been detected (Fig. 7). TESS nearly continuous photometry observations of WASP-148, with a time span of approximately nine planet b periods, are insufficient to detect the TTVs of planet b. This is due to the particular configuration of the TESS observations, which accidentally only cover a part of increasing TTVs (see upper panel Fig. 7). Thus, it is possible that a similar situation is occurring in other systems on which TTVs are also missed. The period derived using only the TESS data lasts 8.80604 ± 0.00014 days, and from all the observed transit times it is 8.8038083 ± 0.0000026 days (Sect. 7.2). This is a difference of 193 ± 12 s. Using the TESS ephemeris to predict future transit timing induces an offset of ~2.2 h yr⁻¹.

The posterior TTVs' evolution in the near future is shown in Fig. A.4. The posterior transit duration variations are shown in Fig. A.5.

## 7.2. Model-independent linear ephemeris

The predictions issued from the photodynamical model and presented in Table A.5 have the disadvantage of having been obtained under the assumption that the two-planet model is correct. The results from the transit-only analysis presented in Sect. 5 can be used to provide a model independent ephemeris which should provide valid predictions, albeit less precise ones, even if the system is discovered to contain additional planets in the future.

A value of a mean period and time of transit can be straightforwardly produced by fitting a slope to the transit times in Table A.2. However, in the presence of transit-timing variations, a slope is not a flexible enough model to describe the transit times. Therefore, the parameters inferred from such a model are likely biased (e.g. Bishop 2007). Hébrard et al. (2020) dealt with this by inflating the error bars of the individual transit times to reach χ² ~ 1. Here we decided to use a non-parametric model to describe the variation of the transit times over the linear ephemeris model.

More precisely, we chose a GP regression model whose mean function was specified to be a linear function with parameters – slope and intercept – to be inferred from the data (see Rasmussen & Williams 2005, Sect. 2.7). If normal priors are chosen for these parameters, the computation of the marginal likelihood can be performed analytically. We chose diffuse priors with widths of 100 days, centred at 8.804 days for the slope (period) and at the observed time of the second TESS transit observed in sector 25 for the intercept.

<figure><img src="image_9.png" /><figcaption>Fig. 7. Posterior TTV predictions of planet b computed relative to a linear ephemeris⁷ are shown in blue. A thousand random draws from the posterior distribution were used to estimate the TTV median value and its uncertainty (68.3% credible interval). In the upper panel, the median TTV values are shown and compared with individual transit-time determinations (Hébrard et al. 2020 in red, and Sect. 5 in black). The thick grey line represents a linear ephemeris computed using only the transits observed by TESS (whose residuals are shown in the small panel in the upper right). In the lower panel, the posterior median transit-timing value was subtracted to visualise the uncertainty of the distribution. The posterior median transit time was also subtracted from each observed epoch for the individual transit-time determinations to allow for better comparison with the photodynamical modelling. The orange curve in the upper panel represents the variation in the times of inferior conjunction for planet c.</figcaption></figure>

We tried several kernel functions to define the covariance function of the GP. We used a modified version of the implementation in scikit-learn (Pedregosa et al. 2011) and optimised the hyperparameters using the L-BFGS-B (Byrd et al. 1995; Zhu et al. 1997) and the Sequential Least Squares Programming (SLSQP) algorithms (Kraft et al. 1988), implemented in the scipy package (Virtanen et al. 2020). All tested kernels produced almost identical results for the parameters of the linear ephemeris. We discuss below the results issued from the kernel choice that provided the largest value of optimised marginal likelihood.

The kernel function, $k(x, x')$, producing the largest marginal likelihood value was an exp-sine-squared (ESS) kernel, without exponential decay:

$$k(x, x') = A_{ESS} \exp \left\{ -\frac{2 \sin^2 \left[ \pi(x-x')^2 / P \right]}{\epsilon^2} \right\}, \quad (1)$$

where P and ε are the kernel hyperparameters that correspond to the period and length-scale, respectively, and AESS is the amplitude of the covariance function. The posterior mean and 68.3% intervals for the TTVs are shown in Fig. 8. A kernel function

<footnote n="6">⁶ Computed using the posterior median of the mid-transit times between the first and last transit observation, 2 457 957.493798 + 8.803808 × Epoch [BJD$_{TDB}$]</footnote>

<footnote n="5">⁵ $P_{sup} = \frac{P_c}{4|\Delta|}$, $\Delta = \frac{P_c}{4P_b} - 1$, with $P_{b,c}$ being the planet orbital period of the corresponding subscript.</footnote>
<figure><img src="image_10.png" /><figcaption>Fig. 8. Model-independent analysis of the timings of the observed transits (black circles). The black line is the mean prediction of the posterior GP model, and the blue band is the 68.3% credible interval.</figcaption></figure>

with an additional decay term produced similar results, with a decay timescale far exceeding the time span of the observations.

Under the assumption of normal priors, the linear model parameter posterior is also normal. We find the posterior ephemeris:

$$ P = 8.8038083 \pm 0.0000026 \text{ d} \\ T_c = 2\,459\,005.148 \pm 0.016 \text{ [BJD}_{\text{TDB}}\text{],} $$

with covariance term, cov(P, T_c) = 6.2 × 10⁻¹⁰. The large uncertainties in T_c are misleading. The predicted transit times issued from this model have a covariance which includes additional terms coming from the interaction between the mean function and the non-parametric part of the model. In fact, the predictions of the transit times remain precise to better than 10 min for over 2640 transits of planet b, that is to say over 60 years. Comparison with the transit times reported in Table A.5 shows that this method predicts transit times that are in agreement with the fully photodynamical one to better than 5 minutes up to the end of 2026.

In comparison, an ordinary least squares (OLS) fit to the transit times, with uncertainties scaled to have a reduced χ² of one as presented in Hébrard et al. (2020), produces a period of P = 8.803824 ± 0.000029 days, which is in agreement with the result from above and with the values presented in Hébrard et al. (2020). The OLS intercept estimator has a standard deviation six times smaller than the one in our model, which is probably unrealistic. However, the predictions from the OLS model remain precise to better than 10 min for less than four years.

The optimised covariance amplitude is A_ESS = 25.3 min, and the period P is 49.7 orbits, corresponding to 437 days. The length-scale ϵ was fixed to 2.0. The fact that the model with the largest marginal likelihood does not include a long-term evolution of the transit times means that with the current data, such a trend is not detected.

### 7.3. Mutual inclination

We inferred a mutual inclination between the planets<fnref n="7" />, I, at t ref of 41.0+6.2°/−7.6° (median and 68.3% credible interval), and a 95% highest density interval (HDI) of [22.8°, 54.8°] (Fig. 9). With these values, a coplanar system, as assumed in Maciejewski et al. (2020), was discarded. Maciejewski et al. (2020) tried a non-coplanar model and found a best fit solution (with Ω_c − Ω_b ≈

−17° and i_c = 47°, or Ω_c − Ω_b ≈ +17° and i_c = 133°) which corresponds to a mutual inclination of ~46°, although they also found that the Bayesian information criterion (BIC) disfavours the non-coplanar solution. However, for high-dimensional models such as these ones, the BIC is known to provide unreliable results (Díaz et al. 2016; Nelson et al. 2020). The coplanar solution in Maciejewski et al. (2020) has a lower eccentricity for planet b (Fig. 10). From a stability analysis of the orbital solution they derived, Hébrard et al. (2020) found an upper limit of 35° for the mutual inclination. We note that the reduced-χ² level curves of the Newtonian fit plotted in Fig. 9 of Hébrard et al. (2020) favoured mutual inclinations around 30°.

We tried to investigate which observable favours the significant mutual inclination. For this, we repeated the photodynamical analysis (Sect. 6) assuming coplanar orbits; we fixed the longitude of the ascending node of both planets to the same value, and we matched the orbital inclination of planet c to the one of planet b, with the latter still being a free parameter. The results are presented in Table A.6, and Figs. 10, A.4, and A.5. The TTVs' posteriors of coplanar and inclined orbits are roughly the same for the observed transits (Fig. A.4). On the other hand, the transit duration posteriors are different, the transit duration of the observed transits remains almost constant for the coplanar model, whereas it decreased for the model with inclined orbits (Fig. A.5). However, we found that the precision of individually determined transit durations<fnref n="8" /> is not enough to confirm the results of the modelling with inclined orbits independently. The photodynamical modelling is in principle more sensible than the analysis on individual transits (Almenara et al. 2015), but the heterogeneity of the transit observations analysed in this work call for caution. Small variations in transit shape or timing intrinsic to the different observations could be wrongly interpreted by photodynamical modelling as an evolution in the orbital parameters<fnref n="9" />. If not correctly taken into account, the limb darkening dependence with the observation bandpass could be misinterpreted as changes in the impact parameter, and therefore high mutual inclinations between the planets. However, no correlation is seen between the mutual inclination and the limb-darkening parameters. In addition, the posterior distributions of the limb-darkening parameters agree with those expected from theoretical computation (Claret & Bloemen 2011; Claret 2017), although they are much wider.

In addition, models starting with a co-planar configuration were run. We found that the co-planar region of parameter space is left quickly by the MCMC walkers. The detection of a non-negative mutual inclination does not seem to be produced by an inadequate exploration of a space parameter, or by bad mixing or lack of convergence of the MCMCs. Future observations focussed on distinguishing between the increasingly and differently predicted transit duration could conclude about the mutual inclination.

<footnote n="8">⁸ We analysed each transit observation individually to determine the transit duration (with the techniques described in Sect. 5), excluding partial transits, as well as transits observed by SuperWASP (poor precision) and TESS (30-min. cadence). The transit duration was computed numerically from 1000 models from the posterior predictive distribution sampled to one-second cadence. This values differs from the transit duration computed in Hébrard et al. (2020) using the approximations in Tingley & Sackett (2005), whose limitations were discussed by Kipping (2010b).</footnote>

<footnote n="9">⁹ For example, transit depth variations, which have not been discussed in this work, can be due to changes in the impact parameter, but also due to a different type of contamination in the photometric aperture.</footnote>

<footnote n="7">⁷ cos I = cos i_b cos i_c + sin i_b sin i_c cos (Ω_b − Ω_c).</footnote>
<figure><img src="image_11.png" /></figure>

**Fig. 9.** Posterior of the mutual inclination at $t_{\text{ref}}$ from the photodynamical modelling (Sect. 6), the stable samples (Sect. 7.5, the stable MAP value is marked by a vertical grey line), and the three-planet model (Sect. 7.8).

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

**Fig. 10.** Planet b and c eccentricity, as well as the mutual inclination for the different analysis presented in this work, and in previous literature. Maciejewski et al. (2020) used a $t_{\text{ref}} = 2\,459\,048.747$ BJD$_{TDB}$, which is different from the one used in this work, but the difference in the inferred orbital parameters cannot be explained by their evolution between the two epochs (Fig. A.6). We adopted the values labelled as ‘stable’.

## 7.4. Photodynamical modelling without stellar priors

The photodynamical modelling of photometry and radial velocity data allows one to measure the absolute radius and mass in multi-planetary systems (Agol et al. 2005; Almenara et al. 2015). To test the precision of this determination in this system, we ran the same analysis as in Sect. 6, but with non-informative priors for the stellar mass and radius. Masses, radii, densities, and orbital parameters are listed in Table A.7. The absolute parameters’ precision is poor. The radii were determined with a precision of 32% relative uncertainty for the star and planet b. The bulk densities of the star and planet b were determined with a precision of 20 and 39%, respectively. The masses of the star, planet b, and planet c were determined with a precision of 89, 58, and 82%, respectively. The precision on planet masses outperform the one on the star (Almenara et al. 2018b). The posterior eccentricities are also more precise for the photodynamical modelling with the stellar priors (Fig. 10). This means that the photoeccentric effect (Dawson & Johnson 2012) puts constraints on the eccentricities in addition to the ones coming from the TTVs. The same is true for the mutual inclination.

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

**Fig. 11.** Stability analysis of 50 000 samples of the posterior distribution obtained from the photodynamical modelling (Sect. 6), projected in a ($\Omega_c, i_c$) diagram. The colour scale corresponds to values between -8.5 (black) and 0.9 (red) for the decimal logarithm of the stability index $D$ used in Correia et al. (2010). The red zones correspond to highly unstable orbits, while the dark blue region can be assumed to be stable on a billion-year timescale.

## 7.5. Stability analysis

The dynamical analysis of the WASP-148 orbital solution reported in Hébrard et al. (2020) has shown that the system is stable, despite significant mutual gravitational interactions between the planets. For this study, we repeated a similar stability study, but instead we performed a global frequency analysis (Laskar 1990, 1993) on 50 000 samples of the posterior distribution obtained from our new photodynamical model. We used the symplectic integrator SABA1064 of Farrés et al. (2013), with a step size of $5 \times 10^{-3}$ yr and general relativity corrections. Each initial condition was integrated over 50 kyr, and a stability indicator $D$ was derived with the frequency analysis of the mean longitude, that is the variation in the measured mean motion over the two consecutive 25 kyr intervals of time (for more details, see Couetdic et al. 2010). For regular motion there is no significant variation in the mean motion along the trajectory, while it can significantly vary for chaotic trajectories. In Fig. 11 we show the distribution of the 50 000 samples projected in a ($\Omega_c, i_c$) diagram, which corresponds to the two less constrained parameters. The colour index gives the value of $D$ for each solution. The values of $\log_{10} D < -6$ for both WASP-148 b and c correspond to stable systems on scales of billions of years (Correia et al. 2010). Only a small region from this diagram is stable for $\Omega_c < 205^\circ$ and $i_c < 120^\circ$, corresponding to 1239 solutions (2.5% of the total). In Fig. 9 we show the probability density function (PDF) of this stable subset of solutions. The solutions cluster around a mutual inclination of $20.8 \pm 4.6^\circ$. We hence conclude from that three-body analysis that WASP-148 can only be stable for mutual inclinations below about 30°. In the last two columns in Table 4, we provide the system parameters corresponding to the MAP, the median, and the 68% credible interval of the stable solutions.

As in Hébrard et al. (2020), we also explored the stability around the stable MAP solution (Table 4), which we refer to as the nominal solution henceforward. As expected, in the ($a_c, e_c$) domain, we confirm that it lies in a stable area, with the orbits close to the 4:1 mean-motion resonance. In the ($i_c, \Omega_c$) domain (Fig. 12), we confirm that stable orbits must have a mutual inclination of $I \le 30^\circ$. In comparison to the stability map shown in
Fig. 9 in Hébrard et al. (2020) with $i_c < 90°$, here we only zoom into the most stable regions, and imposed $i_c > 90°$: this is equivalent from a dynamical point of view, but now the stable solutions are centred around $Ω_c = 180°$ rather than $Ω_c = 0°$. The main difference in the new analysis is the fact that the nodes and the inclination are now constrained by the observations.

The two stability analyses presented above use different approaches: Fig. 11 presents a stability analysis carried out on different solutions allowing all parameters to vary, whereas Fig. 12 studies the stability with only two free parameters ($i_c$, $Ω_c$) and the other parameters are fixed at their MAP values. We show here that they provide similar results in terms of the derived parameters for stability, which supports the reliability of both approaches and their results.

## 7.6. Orbital evolution

To explore the dynamics of the system, we analyse 1000 stable samples from the joint parameter posterior distribution of the photodynamical model and performed numerical integrations<fnref n="10" /> for 1 kyr after $t_{\text{ref}}$. The results for the selected parameters are plotted in Figs. A.6 and A.7. The posterior median and 68.3% credible interval of the mutual inclination over 1 kyr is 20.3 ± 4.5°, and for the eccentricities it is $e_b = 0.126 \pm 0.076$ and $e_c = 0.217^{+0.018}_{-0.027}$. The mutual inclination remains above ~10° at the 95.4% credible interval. Over the 1 kyr integration, the orbital inclination of planet c is too low for transits to occur for most of the samples. Interestingly, planet b only transits for a small fraction of that time (Fig. A.6); in particular, it will not transit anymore after about 200 yr, then it will transit again in about 600 yr.

## 7.7. Tidal evolution

The semi-major axis of the innermost planet is only 0.082 au, which means that the planet is close enough to the star to undergo some tidal evolution. As a result, the eccentricity can be damped and the inner orbit circularised (e.g. Hut 1981). Adopting a value identical to Jupiter's value (Lainey et al. 2009) for the tidal quality factor $Q = 10^5$, we get a characteristic timescale ~10 Gyr (e.g. Correia & Laskar 2010) for the circularisation, which is comparable to the lifetime of the system. Therefore, it is not surprising that at present the innermost planet still shows a significant eccentricity. Moreover, the orbits of the two planets strongly interact, and secular or resonant effects can also excite the eccentricity of the innermost planet (e.g. Correia et al. 2012, 2013). To check this scenario, using a direct three-body model with linear tides (Correia 2018), and general relativity corrections, we ran a simulation over 800 Myr, starting with the initial conditions of the nominal solution from Table 4. We adopted a Love number of $k_2 = 0.5$ and a time lag of $\Delta t = 1$ s (equivalent to $Q \sim 10^5$). We observed that the eccentricity of the innermost planet undergoes large oscillations owing to secular interactions. Moreover, after 130 Myr, the system crosses a resonance which pumps the inner planet eccentricity and damps the mutual inclination (Fig. 13). As a result, at the end of the simulation, the average inner planet eccentricity is higher than the initial one. We conclude that the presently observed non-zero value is compatible with the tidal evolution of the WASP-148 system. As expected, the semi-major axis of the innermost planet also slightly decreases, and we have $a_b \approx 0.072$ au after 800 Myr. Since the age of the star at present is already 4 Gyr, we can

<figure><img src="image_14.png" /><figcaption>**Fig. 12.** Global stability analysis of the WASP-148 planetary system. We fixed all orbital parameters of the stable MAP solution (Table 4), and we varied the inclination, $i_c$, and the longitude of the ascending node, $Ω_c$, of planet-c. The step size is 0.25° in both axes. For each initial condition, the system was integrated over 50 kyr and a stability criterion was derived with the frequency analysis of the mean longitude. White dashed curves give the isolines of constant mutual inclination $I = 10°, 20°, 30°, \text{and } 40°$. The white dot marks the position of the stable MAP solution from Table 4. The colour bar corresponds to the one in Fig. 11.</figcaption></figure>

assume that the initial semi-major axis was slightly larger than the present value, and so the two planets could even be trapped in the 4:1 mean motion resonance. This is an interesting formation scenario that deserves more attention in future work on the system.

## 7.8. Three-planet model

All these results rely on the three-body system hypothesis. If the dynamics of planet b is affected by additional planets other than planet c, increasing the mutual inclination could provide the additional variability in the model required to fit the data. The photodynamical model is flexible and could overfit the data.

<footnote n="10">¹⁰ With the same n-body integrator and time step used in Sect. 6.</footnote>
<figure><img src="image_15.png" /><figcaption>**Fig. 13.** Tidal evolution of the WASP-148 system. Starting with the nominal solution from Table 4, we show the evolution of the mutual inclination (top), the eccentricities (middle), and semi-major axes (bottom) over 800 Myr. We used a direct three-body model with linear tides (Correia 2018) with general relativity corrections, and adopted $k_2 = 0.5$ and $\Delta t = 1$ s for both planets.</figcaption></figure>

The argument of an inaccurately specified model is discussed by Petit et al. (2020) in the context of the determination of the eccentricity of the K2-19 planets.

To test the influence in the system parameters of an unaccounted for additional planet in the system, we repeated the photodynamical analysis (Sect. 6) with a third planet, which we started at the period of 150 days, which is the peak in the residuals of a 2 Keplerian fit of the radial velocities (Sect. 4.1 and Fig. 3

<figure><img src="image_16.png" /><figcaption>**Fig. 14.** Planet mass versus orbital semi-major axis in logarithmic scale. Error bars mark the position of planet b (blue) and planet c (orange). Grey dots are planets listed in the NASA Exoplanet Archive accessed through DACE API (dace.unige.ch). A Gaussian kernel density estimate is shown in different intensities of green. Solar System planets (black points, from left to right: Venus, Earth, Jupiter, Saturn, Uranus, and Neptune) from NASA.</figcaption></figure>

in Hébrard et al. 2020). The results are presented in Table A.8 and Figs. 9, 10, A.4, and A.5. The third planet converges to a period of $151.2^{+2.1}_{-1.7}$ days and a mass of $0.262^{+0.10}_{-0.076} M_J$, which is in agreement with the estimations reported in Hébrard et al. (2020). Its mutual inclination relative to the planets b and c is not well-constrained. The mutual inclination between planet b and c is reduced from $41.0^{+6.2}_{-7.6}^\circ$ in the two-planet model (without imposing long-term stability) to $21 \pm 12^\circ$ in the three-planet model. The masses of planets b and c as well as the radius of planet b are compatible within $1\sigma$ to the two-planet model. The transit timing and duration variations do not allow one to distinguish between the two- and three-planet models with the current data (Figs. A.4 and A.5). More data are needed to assert the presence of additional planets in the system.

## 8. Discussion

The WASP-148 system is composed of a G5V star orbited by a hot Saturn ($0.287^{+0.022}_{-0.016} M_J$, $0.756^{+0.013}_{-0.017} R_J$, and $0.08215^{+0.00086}_{-0.0015}$ au) and a warm Jupiter ($0.392^{+0.023}_{-0.027} M_J$ and $0.2044^{+0.0021}_{-0.0038}$ au) near a 4:1 mean-motion resonance$^{10}$. The planets have eccentricities $e_b = 0.208^{+0.020}_{-0.025}$ and $e_c = 0.1809^{+0.018}_{-0.0072}$, and a mass ratio of $M_c/M_b = 1.351^{+0.11}_{-0.076}$. Also, assuming only two planets in the system, their orbits have a mutual inclination of $20.8 \pm 4.6^\circ$.

Both planets are located in a scarcely populated part of the mass-separation diagram (Fig. 14), which cannot be fully explained by a detection bias$^{11}$, thus indicating a low probability outcome for the planetary formation process. The scarcity

<footnote>$^{10}$ The period ratio at $t_{\text{ref}}$ is $3.92411^{+0.00014}_{-0.00017}$ and computed as $(a_c/a_b)^{3/2}$, with $a_c/a_b = 2.487868^{+0.000061}_{-0.000073}$. For the stable MAP solution, the period ratio (computed over a 1000 orbits of the planet c) is 3.922.</footnote>

<footnote>$^{11}$ While the detection of systems similar to WASP-148 is accessible for high-precision stable radial velocities' instruments, they are difficult to detect with transit searches. Due to the mutually inclined orbits and the secular oscillation of the orbital inclination, the probability that both planets transit for a given observer is low (Sect. 7.6).</footnote>
of planets such as WASP-148 b can be partially explained by the low efficiency of the Lidov-Kozai mechanism (Lidov 1962; Kozai 1962) to form hot Saturns because most migrating planets are tidally disrupted (Anderson et al. 2016). WASP-148 c was placed in the 'period valley' (Udry et al. 2003) between the hot and warm Jupiters clumps.

The WASP-148 giant planets are expected to form in a protoplanetary disk beyond the snow line, which is located at a few au (Lecar et al. 2006), and then migrate inwards (Goldreich & Tremaine 1980). Convergent migration leads to planets being captured in mean-motion resonances, but usually of first order. High order resonances such as the 4:1 one require high initial orbital eccentricities (Rein & Papaloizou 2010). However, disk migration does not favour significant orbital eccentricity excitation because planet-disk interactions tend to damp eccentricities (Bitsch et al. 2013; Dunhill et al. 2013). After the disk disappears, planet-planet scattering can cause high eccentricity migration for the inner planet (Rasio & Ford 1996), or the eccentricities can be excited by chaotic secular interactions (Wu & Lithwick 2011). Both mechanisms are compatible with the non-coplanarity of the observed orbits.

Another possibility that does not require planet-planet scattering or secular chaos is explored in Lee & Peale (2002): the planets are captured in the mean-motion resonance and are massive enough to open gaps in the disk. With the mechanism for damping the eccentricity being reduced, the orbits can become elliptical during the inward migration within the resonance (Artymowicz 1992). Capture in the 4:1 mean-motion resonance is then able to produce some inclination excitation given that the inner planet is not too massive (Thommes & Lissauer 2003; Libert & Tsiganis 2009).

Finally, at a later stage, tidal interactions with the star should shrink the inner planet orbit afterwards and pull the system out of resonance. Tidal interactions usually also damp the eccentricity, but the presence of the outer planet delays this process due to secular and resonant interactions.

There are only a few other planetary systems for which a large mutual inclination has been reported: 102.3°$^{+7.4°}_{-8.0°}$ in HD 3167 (Dalal et al. 2019; Bourrier et al. 2021), 30±1° in *υ* Andromedae (McArthur et al. 2010), 24$^{+11°}_{-8°}$ in Kepler-108 (Mills & Fabrycky 2017), [34.5, 140.6]° (95% credible interval) in $\pi$ Mensae (De Rosa et al. 2020), and some systems with ultra short period planets with mutual inclinations larger than 10° reported in Dai et al. (2018). HD 3167 is a system with three small-sized planets, including two transiting that allowed Rossiter-McLaughlin measurements. $\pi$ Mensae is a hierarchical system with a period ratio of ~330, and a large difference in planet masses. *υ* Andromedae is composed of three planets, and the mutual inclination is measured for the outer pair of super Jovian planets with a period ratio of 5.3. Kepler-108, composed of two Saturn-mass planets with a period ratio of 3.88, is the one that resembles WASP-148 the most, but with longer period planets in lower eccentricity orbits. Interestingly, the period ratio is also close to the 4:1 mean-motion resonance.

As shown in Sect. 7.8, additional planets not accounted for in the modelling can affect the determination of the system parameters, in particular the mutual inclination. Thus more data are needed to conclude on the true architecture of the system.

**Acknowledgements.** This paper includes data collected by the TESS mission. Funding for the TESS mission is provided by the NASA Explorer Program. Resources supporting this work were provided by the NASA High-End Computing (HEC) Program through the NASA Advanced Supercomputing (NAS) Division at Ames Research Center for the production of the SPOC data products. This work has made use of data from the European Space Agency (ESA)

mission Gaia (https://www.cosmos.esa.int/gaia), processed by the Gaia Data Processing and Analysis Consortium (DPAC, https://www.cosmos.esa.int/web/gaia/dpac/consortium). Funding for the DPAC has been provided by national institutions, in particular the institutions participating in the Gaia Multilateral Agreement. Simulations in this paper made use of the REBOUND code which can be downloaded freely at http://github.com/hannorein/rebound. Part of these simulations have been run on the Lesta cluster kindly provided by the Observatoire de Genève. This research has made use of the NASA Exoplanet Archive, which is operated by the California Institute of Technology, under contract with the National Aeronautics and Space Administration under the Exoplanet Exploration Program. A.C. acknowledges support by CFIS/UC projects (UIDB/04564/2020 and UIDP/04564/2020), GRAVITY (PTDC/FIS-AST/7002/2020), ENGAGE SKA (POCI-01-0145-FEDER-022217), and PHO-BOS (POCI-01-0145-FEDER-029932), funded by COMPETE 2020 and FCT, Portugal. G.M. acknowledges the financial support from the National Science Centre, Poland through grant no. 2016/23/B/ST9/00579. MF acknowledges financial support from grant PID2019-109522GB-C5X/AEI/10.13039/501100011033 of the Spanish Ministry of Science and Innovation (MICINN). M.F., V.C. and J.S. acknowledge financial support from the State Agency for Research of the Spanish MCIU through the Center of Excellence Severo Ochoa award to the Instituto de Astrofísica de Andalucía (SEV-2017-0709). J.M.A. and X.B. acknowledge funding from the European Research Council under the ERC Grant Agreement n. 337591-ExTrA. We thank L. Kreidberg for her Mandel & Agol code.

## References

Agol, E., Steffen, J., Sari, R., & Clarkson, W. 2005, MNRAS, 359, 567

Allard, F., Homeier, D., & Freytag, B. 2012, Phil. Trans. R. Soc. London Ser. A, 370, 2765

Almenara, J. M., Díaz, R. F., Mardling, R., et al. 2015, MNRAS, 453, 2644

Almenara, J. M., Díaz, R. F., Dorn, C., Bonfils, X., & Udry, S. 2018a, MNRAS, 478, 460

Almenara, J. M., Díaz, R. F., Hébrard, G., et al. 2018b, A&A, 615, A90

Anderson, K. R., Storch, N. I., & Lai, D. 2016, MNRAS, 456, 3671

Andrae, R., Fouesneau, M., Creevey, O., et al. 2018, A&A, 616, A8

Artymowicz, P. 1992, PASP, 104, 769

Barnes, S. A. 2010, ApJ, 722, 222

Barnes, S. A., & Kim, Y.-C. 2010, ApJ, 721, 675

Barros, S. C. C., Boué, G., Gibson, N. P., et al. 2013, MNRAS, 430, 3032

Bishop, C. M. 2007, Pattern Recognition and Machine Learning, Information Science and Statistics, 1st edn. (Berlin: Springer)

Bitsch, B., Crida, A., Libert, A. S., & Lega, E. 2013, A&A, 555, A124

Bourrier, V., Lovis, C., Cretignier, M., et al. 2021, A&A, 654, A152

Byrd, R. H., Lu, P., Nocedal, J., & Zhu, C. 1995, SIAM J. Sci. Comput., 16, 1190

Carter, J. A., Fabrycky, D. C., Ragozzine, D., et al. 2011, Science, 331, 562

Castelli, F., & Kurucz, R. L. 2003, IAU Symp., 210, A20

Chen, Y., Girardi, L., Bressan, A., et al. 2014, MNRAS, 444, 2525

Claret, A. 2017, A&A, 600, A30

Claret, A., & Bloemen, S. 2011, A&A, 529, A75

Collins, K. A., Kielkopf, J. F., Stassun, K. G., & Hessman, F. V. 2017, AJ, 153, 77

Correia, A. C. M. 2018, Icarus, 305, 250

Correia, A. C. M., & Laskar, J. 2010, Tidal Evolution of Exoplanets, ed. S. Seager (Tucson: University of Arizona Press), 239

Correia, A. C. M., Couetdic, J., Laskar, J., et al. 2010, A&A, 511, A21

Correia, A. C. M., Boué, G., & Laskar, J. 2012, ApJ, 744, L23

Correia, A. C. M., Boué, G., Laskar, J., & Morais, M. H. M. 2013, A&A, 553, A39

Couetdic, J., Laskar, J., Correia, A. C. M., Mayor, M., & Udry, S. 2010, A&A, 519, A10

Cutri, R. M., Wright, E. L., Conrow, T., et al. 2013, VizieR Online Data Catalog: II/328

Cutri, R. M., Skrutskie, M. F., van Dyk, S., et al. 2003, VizieR Online Data Catalog: II/246

Dai, F., Masuda, K., & Winn, J. N. 2018, ApJ, 864, L38

Dalal, S., Hébrard, G., Lecavelier des Étangs, A., et al. 2019, A&A, 631, A28

Dawson, R. I., & Johnson, J. A. 2012, ApJ, 756, 122

De Rosa, R. J., Dawson, R., & Nielsen, E. L. 2020, A&A, 640, A73

Díaz, R. F., Almenara, J. M., Santerne, A., et al. 2014, MNRAS, 441, 983

Díaz, R. F., Ségransan, D., Udry, S., et al. 2016, A&A, 585, A134

Dotter, A., Chaboyer, B., Jevremović, D., et al. 2008, ApJS, 178, 89

Dunhill, A. C., Alexander, R. D., & Armitage, P. J. 2013, MNRAS, 428, 3072

Eastman, J., Siverd, R., & Gaudí, B. S. 2010, PASP, 122, 935

Espinoza, N. 2018, Res. Notes Am. Astron Soc, 2, 209

Espinoza, N., Kossakowski, D., & Brahm, R. 2019, MNRAS, 490, 2262

Farrés, A., Laskar, J., Blanes, S., et al. 2013, Celest. Mech. Dyn. Astron., 116, 141
Feinstein, A. D., Montet, B. T., Foreman-Mackey, D., et al. 2019, PASP, 131, 094502

Foreman-Mackey, D., Hogg, D. W., Lang, D., & Goodman, J. 2013, PASP, 125, 306

Foreman-Mackey, D., Agol, E., Ambikasaran, S., & Angus, R. 2017, AJ, 154, 220

Freudenthal, J., von Essen, C., Ofir, A., et al. 2019, A&A, 628, A108

Fulton, B. J., Shporer, A., Winn, J. N., et al. 2011, AJ, 142, 84

Fulton, B. J., Petigura, E. A., Blunt, S., & Sinukoff, E. 2018, PASP, 130, 044504

Gaia Collaboration (Prusti, T., et al.) 2016, A&A, 595, A1

Gaia Collaboration (Brown, A. G. A., et al.) 2018, A&A, 616, A1

Gaia Collaboration (Brown, A. G. A., et al.) 2021, A&A, 649, A1

Gillon, M., Triaud, A. H. M. J., Demory, B.-O., et al. 2017, Nature, 542, 456

Goldreich, P., & Tremaine, S. 1980, ApJ, 241, 425

Goodman, J., & Weare, J. 2010, Comm. Appl. Math. Comput. Sci., 5, 65

Hadden, S., & Lithwick, Y. 2014, ApJ, 787, 80

Hébrard, G., Díaz, R. F., Correia, A. C. M., et al. 2020, A&A, 640, A32

Holman, M. J., & Murray, N. W. 2005, Science, 307, 1288

Holman, M. J., Fabrycky, D. C., Ragozzine, D., et al. 2010, Science, 330, 51

Hut, P. 1981, A&A, 99, 126

Irwin, J. B. 1952, ApJ, 116, 211

Jenkins, J. M., Twicken, J. D., McCauliff, S., et al. 2016, SPIE Conf. Ser., 9913, 99133E

Jontof-Hutter, D., Wolfgang, A., Ford, E. B., et al. 2021, AJ, 161, 246

Kass, R. E., & Raftery, A. E. 1995, J. Am. Stat. Assoc., 90, 773

Kipping, D. M. 2010a, MNRAS, 408, 1758

Kipping, D. M. 2010b, MNRAS, 407, 301

Kipping, D. M. 2013, MNRAS, 435, 2152

Kozai, Y. 1962, AJ, 67, 591

Kraft, D. 1988, A software package for sequential quadratic programming, DFVLR-FB 88-28

Kreidberg, L. 2015, PASP, 127, 1161

Lainey, V., Arlot, J.-E., Karatekin, Ö., & van Hoolst, T. 2009, Nature, 459, 957

Laskar, J. 1990, Icarus, 88, 266

Laskar, J. 1993, Phys. D: Nonlinear Phenomena, 67, 257

Lecar, M., Podolak, M., Sasselov, D., & Chiang, E. 2006, ApJ, 640, 1115

Lee, M. H., & Peale, S. J. 2002, ApJ, 567, 596

Libert, A. S., & Tsiganis, K. 2009, MNRAS, 400, 1373

Lidov, M. L. 1962, Planet. Space Sci., 9, 719

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

Lissauer, J. J., Fabrycky, D. C., Ford, E. B., et al. 2011, Nature, 470, 53

Lithwick, Y., Xie, J., & Wu, Y. 2012, ApJ, 761, 122

Maciejewski, G., Fernández, M., Sota, A., & García Segura, A. J. 2020, Acta Astron., 70, 203

Mandel, K., & Agol, E. 2002, ApJ, 580, L171

Manduca, A., Bell, R. A., & Gustafsson, B. 1977, A&A, 61, 809

McArthur, B. E., Benedict, G. F., Barnes, R., et al. 2010, ApJ, 715, 1203

Meibom, S., Barnes, S. A., Platais, I., et al. 2015, Nature, 517, 589

Mills, S. M. & Fabrycky, D. C. 2017, AJ, 153, 45

Nelson, B. E., Ford, E. B., Buchner, J., et al. 2020, AJ, 159, 73

Nesvorný, D., Kipping, D., Terrell, D., et al. 2013, ApJ, 777, 3

Pedregosa, F., Varoquaux, G., Gramfort, A., et al. 2011, J. Mach. Learn. Res., 12, 2825

Petit, A. C., Petigura, E. A., Davies, M. B., & Johansen, A. 2020, MNRAS, 496, 3101

Rasio, F. A., & Ford, E. B. 1996, Science, 274, 954

Rasmussen, C. E., & Williams, C. K. I. 2005, Gaussian Processes for Machine Learning, Adaptive Computation and Machine Learning (Cambridge: The MIT Press)

Rein, H., & Liu, S.-F. 2012, A&A, 537, A128

Rein, H., & Papaloizou, J. C. B. 2010, EAS Pub. Ser., 42, 299

Rein, H., & Tamayo, D. 2015, MNRAS, 452, 376

Ricker, G. R., Winn, J. N., Vanderspek, R., et al. 2015, J. Astrono. Teles. Instrume. Syst., 1, 014003

Riello, M., De Angeli, F., Evans, D. W., et al. 2021, A&A, 649, A3

Skrutskie, M. F., Cutri, R. M., Stiening, R., et al. 2006, AJ, 131, 1163

Speagle, J. S. 2020, MNRAS, 493, 3132

Stassun, K. G., Oelkers, R. J., Paegert, M., et al. 2019, AJ, 158, 138

Thommes, E. W., & Lissauer, J. J. 2003, ApJ, 597, 566

Tingley, B., & Sackett, P. D. 2005, ApJ, 627, 1011

Udry, S., Mayor, M., & Santos, N. C. 2003, A&A, 407, 369

Virtanen, P., Gommers, R., Oliphant, T. E., et al. 2020, Nat. Methods, 17, 261

Wang, J., Xie, J.-W., Barclay, T., & Fischer, D. A. 2014, ApJ, 783, 4

Wang, X.-Y., Rice, M., Wang, S., et al. 2022, ApJ, 926, L8

Wright, E. L., Eisenhardt, P. R. M., Mainzer, A. K., et al. 2010, AJ, 140, 1868

Wu, Y., & Lithwick, Y. 2011, ApJ, 735, 109

Zhu, C., Byrd, R. H., Lu, P., & Nocedal, J. 1997, ACM Trans. Math. Softw., 23, 550

<footnote n="1">1 Observatoire de Genève, Département d'Astronomie, Université de Genève, Chemin des Maillettes 51, 1290 Versoix, Switzerland</footnote>

<footnote n="2">2 Univ. Grenoble Alpes, CNRS, IPAG, 38000 Grenoble, France e-mail: Jose-Manuel.Almenara-Villa@univ-grenoble-alpes.fr</footnote>

<footnote n="3">3 Institut d'astrophysique de Paris, UMR7095 CNRS, Université Pierre & Marie Curie, 98bis boulevard Arago, 75014 Paris, France</footnote>

<footnote n="4">4 Observatoire de Haute-Provence, 04670 Saint Michel l'Observatoire, France</footnote>

<footnote n="5">5 International Center for Advanced Studies (ICAS) and ICIFI (CON-ICET), ECyT-UNSAM, Campus Miguelete, 25 de Mayo y Francia, (1650) Buenos Aires, Argentina</footnote>

<footnote n="6">6 IMCC, UMR8028 CNRS, Observatoire de Paris, PSL University, Sorbonne Univ., 77 av. Denfert-Rochereau, 75014 Paris, France</footnote>

<footnote n="7">7 CFisUC, Departamento de Física, Universidade de Coimbra, 3004-516 Coimbra, Portugal</footnote>

<footnote n="8">8 Centre for Exoplanets and Habitability, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, UK</footnote>

<footnote n="9">9 Department of Physics, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, UK</footnote>

<footnote n="10">10 Astrophysics Group, Keele University, Staffordshire, ST5 5BG, UK</footnote>

<footnote n="11">11 Laboratoire d'Astrophysique de Marseille, Univ. de Provence, UMR6110 CNRS, 38 r. F. Joliot Curie, 13388 Marseille cedex 13, France</footnote>

<footnote n="12">12 Instituto de Astrofísica de Andalucía (IAA-CSIC), Glorieta de la Astronomía 3, 18008 Granada, Spain</footnote>

<footnote n="13">13 School of Physics and Astronomy, Physical Science Building, North Haugh, St Andrews, UK</footnote>

<footnote n="14">14 NASA Ames Research Center, Moffett Field, CA 94035, USA</footnote>

<footnote n="15">15 LESIA, Observatoire de Paris, Université PSL, CNRS, Sorbonne Université, Université de Paris, 92195 Meudon, France</footnote>

<footnote n="16">16 Institute of Astronomy, Faculty of Physics, Astronomy and Informat- ics, Nicolaus Copernicus University, Grudziadzka 5, 87-100 Toruń, Poland</footnote>

<footnote n="17">17 Center for Space and Habitability, University of Bern, Gesellschaftsstrasse 6, 3012 Bern, Switzerland</footnote>

<footnote n="18">18 Department of Physics and Kavli Institute for Astrophysics and Space Research, Massachusetts Institute of Technology, Cambridge, MA 02139, USA</footnote>

<footnote n="19">19 Department of Earth, Atmospheric and Planetary Sciences, Massachusetts Institute of Technology, Cambridge, MA 02139, USA</footnote>

<footnote n="20">20 Department of Aeronautics and Astronautics, MIT, 77 Massachusetts Avenue, Cambridge, MA 02139, USA</footnote>

<footnote n="21">21 Observatoire Hubert-Reeves, 07320 Mars, France</footnote>

<footnote n="22">22 Department of Astrophysical Sciences, Princeton University, NJ 08544, USA</footnote>
**Appendix A: Additional figures and tables**

<table><caption>**Table A.1.** Photometric measurements used for the SED analysis of WASP-148.</caption><thead><tr><th>Filter</th><th>Magnitude</th><th>±1σ</th></tr></thead><tbody><tr><td>Gaia-G</td><td>12.0790</td><td>0.0028</td></tr><tr><td>Gaia-BP</td><td>12.4636</td><td>0.0030</td></tr><tr><td>Gaia-RP</td><td>11.5315</td><td>0.0029</td></tr><tr><td>2MASS-J</td><td>10.938</td><td>0.024</td></tr><tr><td>2MASS-H</td><td>10.585</td><td>0.018</td></tr><tr><td>2MASS-Ks</td><td>10.506</td><td>0.017</td></tr><tr><td>WISE-W1</td><td>10.466</td><td>0.022</td></tr><tr><td>WISE-W2</td><td>10.519</td><td>0.020</td></tr><tr><td>WISE-W3<sup>(†)</sup></td><td>10.476</td><td>0.064</td></tr></tbody></table>

**Notes.** (†) Not covered by the Castelli & Kurucz models.
<table><caption>**Table A.2.** Inferred system parameters for the transit-only analysis.</caption><thead><tr><th>Parameter</th><th>Units</th><th>Prior</th><th>Posterior median and 68.3% CI</th></tr></thead><tbody><tr><td colspan="4"><strong>Star</strong></td></tr><tr><td>Stellar mean density</td><td>[g cm<sup>-3</sup>]</td><td>N(1.75, 0.15)</td><td>1.77 ± 0.13</td></tr><tr><td>Kipping (2013) q<sub>1</sub> for Clear, Johnson-R, RISE, TESS</td><td></td><td>U(0, 1)</td><td>0.38<sup>+0.25</sup><sub>-0.18</sub>, 0.204<sup>+0.13</sup><sub>-0.090</sub>, 0.78 ± 0.16, 0.35<sup>+0.22</sup><sub>-0.15</sub></td></tr><tr><td>Kipping (2013) q<sub>2</sub> for Clear, Johnson-R, RISE, TESS</td><td></td><td>U(0, 1)</td><td>0.29<sup>+0.26</sup><sub>-0.19</sub>, 0.57 ± 0.28, 0.40<sup>+0.14</sup><sub>-0.12</sub>, 0.39<sup>+0.29</sup><sub>-0.23</sub></td></tr><tr><td colspan="4"><strong>Planet b</strong></td></tr><tr><td>Espinoza (2018) r<sub>1</sub></td><td></td><td>U(0, 1)</td><td>0.432<sup>+0.080</sup><sub>-0.065</sub></td></tr><tr><td>Espinoza (2018) r<sub>2</sub></td><td></td><td>U(0, 1)</td><td>0.0854 ± 0.0013</td></tr><tr><td>√e cos ω</td><td></td><td>N(0.24, 0.09)</td><td>0.243 ± 0.080</td></tr><tr><td>√e sin ω</td><td></td><td>N(0.40, 0.12)</td><td>0.384<sup>+0.043</sup><sub>-0.049</sub></td></tr><tr><td colspan="4"><strong>Transit timings (instrument, band, epoch)</strong></td></tr><tr><td>WASP1 Clear</td><td>-377</td><td>[BJD<sub>TDB</sub>] U(2454638.265, 2454638.773)</td><td>2454638.461<sup>+0.21</sup><sub>-0.097</sub></td></tr><tr><td>WASP2 Clear</td><td>-296</td><td>[BJD<sub>TDB</sub>] U(2455351.288, 2455351.835)</td><td>2455351.517<sup>+0.033</sup><sub>-0.041</sub></td></tr><tr><td>WASP3 Clear</td><td>-255</td><td>[BJD<sub>TDB</sub>] U(2455712.42, 2455712.65)</td><td>2455712.5208<sup>+0.0080</sup><sub>-0.0067</sub></td></tr><tr><td>NITES1 Clear</td><td>-122</td><td>[BJD<sub>TDB</sub>] U(2456883.440, 2456883.449)</td><td>2456883.4442 ± 0.0010</td></tr><tr><td>SANCHEZ Johnson-R</td><td>-87</td><td>[BJD<sub>TDB</sub>] U(2457191.55, 2457191.56)</td><td>2457191.5555 ± 0.0010</td></tr><tr><td>NITES2 Johnson-R</td><td>-46</td><td>[BJD<sub>TDB</sub>] U(2457552.5005, 2457552.5082)</td><td>2457552.50421 ± 0.00088</td></tr><tr><td>MARS Clear</td><td>-11</td><td>[BJD<sub>TDB</sub>] U(2457860.648, 2457860.660)</td><td>2457860.6550 ± 0.0012</td></tr><tr><td>RISE RISE</td><td>0</td><td>[BJD<sub>TDB</sub>] U(2457957.480, 2457957.484)</td><td>2457957.48202 ± 0.00042</td></tr><tr><td>TESS sector 24 TESS</td><td>114</td><td>[BJD<sub>TDB</sub>] U(2458961.1150, 2458961.1331)</td><td>2458961.1226 ± 0.0019</td></tr><tr><td>TESS sector 24 TESS</td><td>115</td><td>[BJD<sub>TDB</sub>] U(2458969.9220, 2458969.9358)</td><td>2458969.9285 ± 0.0015</td></tr><tr><td>TESS sector 24 TESS</td><td>116</td><td>[BJD<sub>TDB</sub>] U(2458978.7245, 2458978.7438)</td><td>2458978.7346 ± 0.0022</td></tr><tr><td>TESS sector 25 TESS</td><td>117</td><td>[BJD<sub>TDB</sub>] U(2458987.531, 2458987.5476)</td><td>2458987.5403 ± 0.0019</td></tr><tr><td>TESS sector 25 TESS</td><td>119</td><td>[BJD<sub>TDB</sub>] U(2459005.1466, 2459005.1581)</td><td>2459005.1521 ± 0.0016</td></tr><tr><td>TESS sector 26 TESS</td><td>120</td><td>[BJD<sub>TDB</sub>] U(2459013.9516, 2459013.970)</td><td>2459013.9600 ± 0.0021</td></tr><tr><td>TESS sector 26 TESS</td><td>122</td><td>[BJD<sub>TDB</sub>] U(2459031.5623, 2459031.5792)</td><td>2459031.5701 ± 0.0018</td></tr><tr><td colspan="4"><strong>Instruments</strong></td></tr><tr><td>Dilution factor for WASP</td><td></td><td>TN(1.0, 0.1, 0.0, 1.0)</td><td>0.923<sup>+0.047</sup><sub>-0.069</sub></td></tr><tr><td>Dilution factor for TESS sector 24, 25, 26</td><td></td><td>TN(1.0, 0.1, 0.0, 1.0)</td><td>0.949 ± 0.033, 0.936 ± 0.035, 0.945<sup>+0.032</sup><sub>-0.041</sub></td></tr><tr><td>Offset relative flux for WASP1, WASP2, WASP3</td><td>[Relative flux]</td><td>N(0.0, 0.1)</td><td>0.0030<sup>+0.0038</sup><sub>-0.0033</sub>, -0.0006<sup>+0.0037</sup><sub>-0.0030</sub>, 0.0008 ± 0.0022</td></tr><tr><td>Offset relative flux for NITES1, NITES2</td><td>[Relative flux]</td><td>N(0.0, 0.1)</td><td>0.00061<sup>+0.00058</sup><sub>-0.00086</sub>, -0.00015<sup>+0.0016</sup><sub>-0.00068</sub></td></tr><tr><td>Offset relative flux for MARS</td><td>[Relative flux]</td><td>N(0.0, 0.1)</td><td>-0.0019<sup>+0.0062</sup><sub>-0.0060</sub></td></tr><tr><td>Offset relative flux for SANCHEZ</td><td>[Relative flux]</td><td>N(0.0, 0.1)</td><td>-0.00201<sup>+0.00081</sup><sub>-0.00090</sub></td></tr><tr><td>Offset relative flux for RISE</td><td>[Relative flux]</td><td>N(0.0, 0.1)</td><td>0.0034<sup>+0.035</sup><sub>-0.0059</sub></td></tr><tr><td>Offset relative flux for TESS sector 24, 25, 26</td><td>[Relative flux]</td><td>N(0.0, 0.1)</td><td>0.0004 ± 0.0013, 0.0005 ± 0.0019, -0.0003 ± 0.0014</td></tr><tr><td>Additive jitter for WASP1, WASP2, WASP3</td><td>[ppm]</td><td>J(1, 10000)</td><td>82<sup>+980</sup><sub>-76</sub>, 61<sup>+630</sup><sub>-56</sub>, 210<sup>+2800</sup><sub>-200</sub></td></tr><tr><td>Additive jitter for NITES1, NITES2</td><td>[ppm]</td><td>J(1, 4000)</td><td>18<sup>+100</sup><sub>-15</sub>, 3290 ± 130</td></tr><tr><td>Additive jitter for MARS</td><td>[ppm]</td><td>J(1, 1600)</td><td>860 ± 120</td></tr><tr><td>Additive jitter for SANCHEZ</td><td>[ppm]</td><td>J(1, 1600)</td><td>27<sup>+150</sup><sub>-24</sub></td></tr><tr><td>Additive jitter for RISE</td><td>[ppm]</td><td>J(1, 300)</td><td>10.5<sup>+36</sup><sub>-81</sub></td></tr><tr><td>Additive jitter for TESS sector 24, 25, 26</td><td>[ppm]</td><td>J(1, 900)</td><td>27<sup>+110</sup><sub>-23</sub>, 94<sup>+120</sup><sub>-85</sub>, 744 ± 27</td></tr><tr><td>Timescale of the GP for WASP1, WASP2, WASP3</td><td>[days]</td><td>J(0.001, 1000)</td><td>1.0<sup>+91</sup><sub>-1.0</sub>, 1.8<sup>+75</sup><sub>-1.8</sub>, 0.021<sup>+0.034</sup><sub>-0.012</sub></td></tr><tr><td>Timescale of the GP for NITES1, NITES2</td><td>[days]</td><td>J(0.001, 1000)</td><td>6.0<sup>+160</sup><sub>-5.9</sub>, 12<sup>+220</sup><sub>-12</sub></td></tr><tr><td>Timescale of the GP for MARS</td><td>[days]</td><td>J(0.001, 1000)</td><td>1.3<sup>+15</sup><sub>-1.1</sub></td></tr><tr><td>Timescale of the GP for SANCHEZ</td><td>[days]</td><td>J(0.001, 1000)</td><td>0.0233<sup>+0.016</sup><sub>-0.0098</sub></td></tr><tr><td>Timescale of the GP for RISE</td><td>[days]</td><td>J(0.001, 1000)</td><td>0.34<sup>+0.75</sup><sub>-0.22</sub></td></tr><tr><td>Timescale of the GP for TESS sector 24, 25, 26</td><td>[days]</td><td>J(0.001, 1000)</td><td>1.81<sup>+0.38</sup><sub>-0.30</sub>, 5.1<sup>+2.4</sup><sub>-1.4</sub>, 1.04<sup>+0.16</sup><sub>-0.13</sub></td></tr><tr><td>Amplitude of the GP for WASP1, WASP2, WASP3</td><td>[Relative flux]</td><td>J(10<sup>-6</sup>, 1)</td><td>0.00022<sup>+0.0046</sup><sub>-0.00021</sub>, 0.00072<sup>+0.013</sup><sub>-0.00070</sub>, 0.0040<sup>+0.0020</sup><sub>-0.0013</sub></td></tr><tr><td>Amplitude of the GP for NITES1, NITES2</td><td>[Relative flux]</td><td>J(10<sup>-6</sup>, 1)</td><td>0.00020<sup>+0.0044</sup><sub>-0.00019</sub>, 0.00040<sup>+0.0088</sup><sub>-0.00038</sub></td></tr><tr><td>Amplitude of the GP for MARS</td><td>[Relative flux]</td><td>J(10<sup>-6</sup>, 1)</td><td>0.014<sup>+0.10</sup><sub>-0.012</sub></td></tr><tr><td>Amplitude of the GP for SANCHEZ</td><td>[Relative flux]</td><td>J(10<sup>-6</sup>, 1)</td><td>0.00177<sup>+0.00085</sup><sub>-0.00045</sub></td></tr><tr><td>Amplitude of the GP for RISE</td><td>[Relative flux]</td><td>J(10<sup>-6</sup>, 1)</td><td>0.0104<sup>+0.041</sup><sub>-0.0074</sub></td></tr><tr><td>Amplitude of the GP for TESS sector 24, 25, 26</td><td>[Relative flux]</td><td>J(10<sup>-6</sup>, 1)</td><td>0.00312<sup>+0.00077</sup><sub>-0.00050</sub>, 0.00285<sup>+0.0019</sup><sub>-0.00090</sub>, 0.00430<sup>+0.00077</sup><sub>-0.00058</sub></td></tr><tr><td colspan="4"><strong>Derived</strong></td></tr><tr><td>Radius ratio, R<sub>p</sub>/R<sub>*</sub></td><td></td><td></td><td>0.0854 ± 0.0013</td></tr><tr><td>Impact parameter</td><td></td><td></td><td>0.148<sup>+0.12</sup><sub>-0.098</sub></td></tr><tr><td>Eccentricity, e</td><td></td><td></td><td>0.213 ± 0.035</td></tr><tr><td>Argument of pericentre, ω</td><td>[°]</td><td></td><td>58 ± 11</td></tr></tbody></table>

**Notes.** The table lists: Priors and posterior median as well as 68.3% CI for the transit-only analysis with *juliet* (Section 5). *N*(*μ*, *σ*): Normal distribution with mean *μ* and standard deviation *σ*. *TN*(*μ*, *σ*, *a*, *b*): Normal distribution with mean *μ* and standard deviation *σ*, truncated between a lower *a* and upper *b* limit. *U*(*a*, *b*): A uniform distribution defined between a lower *a* and upper *b* limit. *J*(*a*, *b*): Freyje (or log-uniform) distribution defined between a lower *a* and upper *b* limit.
<table><caption>**Table A.3.** Idem. as Figure 3, but for the analysis of OSN150 transits.</caption><thead><tr><th>Parameter</th><th>Units</th><th>Prior</th><th>Posterior median and 68.3% CI</th></tr></thead><tbody><tr><td><strong>Star</strong></td><td></td><td></td><td></td></tr><tr><td>Stellar mean density</td><td>[g cm<sup>-3</sup>]</td><td>N(1.75, 0.15)</td><td>1.73 ± 0.15</td></tr><tr><td>Kipping (2013) q<sub>1</sub> for Clear</td><td></td><td>U(0, 1)</td><td>0.45 <sup>+0.25</sup><sub>-0.16</sub></td></tr><tr><td>Kipping (2013) q<sub>2</sub> for Clear</td><td></td><td>U(0, 1)</td><td>0.17 <sup>+0.17</sup><sub>-0.11</sub></td></tr><tr><td><strong>Planet b</strong></td><td></td><td></td><td></td></tr><tr><td>Espinoza (2018) r<sub>1</sub></td><td></td><td>U(0, 1)</td><td>0.550 <sup>+0.064</sup><sub>-0.11</sub></td></tr><tr><td>Espinoza (2018) r<sub>2</sub></td><td></td><td>U(0, 1)</td><td>0.0851 <sup>+0.0020</sup><sub>-0.0016</sub></td></tr><tr><td>√e cos ω</td><td></td><td>N(0.24, 0.09)</td><td>0.238 ± 0.091</td></tr><tr><td>√e sin ω</td><td></td><td>N(0.40, 0.12)</td><td>0.366 <sup>+0.053</sup><sub>-0.062</sub></td></tr><tr><td><strong>Transit timings (instrument, band, epoch)</strong></td><td></td><td></td><td></td></tr><tr><td>OSN1 Clear 123</td><td>[BJDTDB]</td><td>U(2459040.365914, 2459040.49543)</td><td>2459040.37337 +0.00052<br>-0.00047</td></tr><tr><td>OSN2 Clear 128</td><td>[BJDTDB]</td><td>U(2459084.327292, 2459084.50194)</td><td>2459084.39721 +0.00058<br>-0.00067</td></tr><tr><td>OSN3 Clear 152</td><td>[BJDTDB]</td><td>U(2459295.515766, 2459295.733722)</td><td>2459295.65861 +0.00053<br>-0.00053</td></tr><tr><td>OSN4 Clear 163</td><td>[BJDTDB]</td><td>U(2459392.390943, 2459392.628101)</td><td>2459392.51001 +0.00036<br>-0.00036</td></tr><tr><td><strong>Instruments</strong></td><td></td><td></td><td></td></tr><tr><td>Offset relative flux for OSN1</td><td>[Relative flux]</td><td>U(-0.001, 0.001)</td><td>-0.00035 <sup>+0.00067</sup><sub>-0.00041</sub></td></tr><tr><td>Offset relative flux for OSN2</td><td>[Relative flux]</td><td>U(-0.001, 0.001)</td><td>0.00037 <sup>+0.00043</sup><sub>-0.00068</sub></td></tr><tr><td>Offset relative flux for OSN3</td><td>[Relative flux]</td><td>U(-0.001, 0.001)</td><td>0.00040 <sup>+0.00042</sup><sub>-0.00066</sub></td></tr><tr><td>Offset relative flux for OSN4</td><td>[Relative flux]</td><td>U(-0.001, 0.001)</td><td>0.00006 <sup>+0.00047</sup><sub>-0.00055</sub></td></tr><tr><td>Additive jitter for OSN1</td><td>[ppm]</td><td>J(1, 2000)</td><td>979 ± 63</td></tr><tr><td>Additive jitter for OSN2</td><td>[ppm]</td><td>J(1, 2000)</td><td>1374 ± 77</td></tr><tr><td>Additive jitter for OSN3</td><td>[ppm]</td><td>J(1, 2000)</td><td>1051 ± 56</td></tr><tr><td>Additive jitter for OSN4</td><td>[ppm]</td><td>J(1, 2000)</td><td>19 <sup>+130</sup><sub>-17</sub></td></tr><tr><td>Timescale of the GP for OSN1</td><td>[days]</td><td>J(0.001, 1000)</td><td>0.118 <sup>+0.33</sup><sub>-0.075</sub></td></tr><tr><td>Timescale of the GP for OSN2</td><td>[days]</td><td>J(0.001, 1000)</td><td>0.063 <sup>+0.068</sup><sub>-0.032</sub></td></tr><tr><td>Timescale of the GP for OSN3</td><td>[days]</td><td>J(0.001, 1000)</td><td>0.048 <sup>+0.060</sup><sub>-0.024</sub></td></tr><tr><td>Timescale of the GP for OSN4</td><td>[days]</td><td>J(0.001, 1000)</td><td>0.066 <sup>+0.076</sup><sub>-0.035</sub></td></tr><tr><td>Amplitude of the GP for OSN1</td><td>[Relative flux]</td><td>J(10<sup>-6</sup>, 0.05)</td><td>0.00106 <sup>+0.0023</sup><sub>-0.00056</sub></td></tr><tr><td>Amplitude of the GP for OSN2</td><td>[Relative flux]</td><td>J(10<sup>-6</sup>, 0.05)</td><td>0.00120 <sup>+0.0010</sup><sub>-0.00056</sub></td></tr><tr><td>Amplitude of the GP for OSN3</td><td>[Relative flux]</td><td>J(10<sup>-6</sup>, 0.05)</td><td>0.00117 <sup>+0.0010</sup><sub>-0.00044</sub></td></tr><tr><td>Amplitude of the GP for OSN4</td><td>[Relative flux]</td><td>J(10<sup>-6</sup>, 0.05)</td><td>0.00090 <sup>+0.00095</sup><sub>-0.00034</sub></td></tr><tr><td><strong>Derived</strong></td><td></td><td></td><td></td></tr><tr><td>Radius ratio, R<sub>p</sub>/R<sub>★</sub></td><td></td><td></td><td>0.0851 <sup>+0.0020</sup><sub>-0.0016</sub></td></tr><tr><td>Impact parameter</td><td></td><td></td><td>0.324 <sup>+0.096</sup><sub>-0.17</sub></td></tr><tr><td>Eccentricity, e</td><td></td><td></td><td>0.197 ± 0.046</td></tr><tr><td>Argument of pericentre, ω</td><td>[°]</td><td></td><td>57 ± 13</td></tr></tbody></table>

**Notes.** N(μ, σ): Normal distribution with mean μ and standard deviation σ. TN(μ, σ, a, b): Normal distribution with mean μ and standard deviation σ, truncated between a lower a and upper b limit. U(a, b): A uniform distribution defined between a lower a and upper b limit. J(a, b): Jeffreys (or log-uniform) distribution defined between a lower a and upper b limit.
<table><caption>**Table A.4.** Continuation of Table 4 for the photodynamical modelling.</caption><thead><tr><th>Parameter</th><th>Units</th><th>Prior</th><th>Median and 68.3% CI</th></tr></thead><tbody><tr><td>SOPHIE offset</td><td>[km s<sup>-1</sup>]</td><td>U(-100, 100)</td><td>5.6179 ± 0.0014</td></tr><tr><td>SOPHIE multiplicative jitter</td><td></td><td>U(0, 10)</td><td>1.53<sup>+0.12</sup><sub>-0.10</sub></td></tr><tr><td>Dilution factor for WASP</td><td></td><td>U(0, 1)</td><td>0.53<sup>+0.25</sup><sub>-0.31</sub></td></tr><tr><td>Dilution factor for TESS sector 24</td><td></td><td>U(0, 1)</td><td>0.104<sup>+0.035</sup><sub>-0.041</sub></td></tr><tr><td>Dilution factor for TESS sector 25</td><td></td><td>U(0, 1)</td><td>0.122<sup>+0.037</sup><sub>-0.044</sub></td></tr><tr><td>Dilution factor for TESS sector 26</td><td></td><td>U(0, 1)</td><td>0.116 ± 0.053</td></tr><tr><td>Relative flux for WASP1</td><td>[Relative flux]</td><td>U(0.9, 1.1)</td><td>1.0002<sup>+0.0023</sup><sub>-0.0027</sub></td></tr><tr><td>Relative flux for WASP2</td><td>[Relative flux]</td><td>U(0.9, 1.1)</td><td>1.0022 ± 0.0015</td></tr><tr><td>Relative flux for WASP3</td><td>[Relative flux]</td><td>U(0.9, 1.1)</td><td>0.9995 ± 0.0011</td></tr><tr><td>Relative flux for NITES1</td><td>[Relative flux]</td><td>U(0.9, 1.1)</td><td>1.00014 ± 0.00018</td></tr><tr><td>Relative flux for NITES2</td><td>[Relative flux]</td><td>U(0.9, 1.1)</td><td>0.99994 ± 0.00016</td></tr><tr><td>Relative flux for MARS</td><td>[Relative flux]</td><td>U(0.9, 1.1)</td><td>0.99994 ± 0.00013</td></tr><tr><td>Relative flux for SANCHEZ</td><td>[Relative flux]</td><td>U(0.9, 1.1)</td><td>1.00006<sup>+0.00015</sup><sub>-0.00013</sub></td></tr><tr><td>Relative flux for RISE</td><td>[Relative flux]</td><td>U(0.9, 1.1)</td><td>1.000018<sup>+0.000065</sup><sub>-0.000057</sub></td></tr><tr><td>Relative flux for TESS sector 24</td><td>[Relative flux]</td><td>U(0.9, 1.1)</td><td>0.999921 ± 0.00011</td></tr><tr><td>Relative flux for TESS sector 25</td><td>[Relative flux]</td><td>U(0.9, 1.1)</td><td>1.00001 ± 0.00012</td></tr><tr><td>Relative flux for TESS sector 26</td><td>[Relative flux]</td><td>U(0.9, 1.1)</td><td>0.99991 ± 0.00014</td></tr><tr><td>Relative flux for OSN1</td><td>[Relative flux]</td><td>U(0.9, 1.1)</td><td>1.000054 ± 0.000086</td></tr><tr><td>Relative flux for OSN2</td><td>[Relative flux]</td><td>U(0.9, 1.1)</td><td>1.000009 ± 0.000097</td></tr><tr><td>Relative flux for OSN3</td><td>[Relative flux]</td><td>U(0.9, 1.1)</td><td>0.999996 ± 0.000072</td></tr><tr><td>Relative flux for OSN4</td><td>[Relative flux]</td><td>U(0.9, 1.1)</td><td>0.999998 ± 0.000052</td></tr><tr><td>Additive jitter for WASP1</td><td>[ppm]</td><td>J(1, 10000)</td><td>69<sup>+1000</sup><sub>-65</sub></td></tr><tr><td>Additive jitter for WASP2</td><td>[ppm]</td><td>J(1, 10000)</td><td>53<sup>+540</sup><sub>-49</sub></td></tr><tr><td>Additive jitter for WASP3</td><td>[ppm]</td><td>J(1, 10000)</td><td>320<sup>+2500</sup><sub>-310</sub></td></tr><tr><td>Additive jitter for NITES1</td><td>[ppm]</td><td>J(1, 10000)</td><td>25<sup>+130</sup><sub>-22</sub></td></tr><tr><td>Additive jitter for NITES2</td><td>[ppm]</td><td>J(1, 10000)</td><td>3300 ± 140</td></tr><tr><td>Additive jitter for MARS</td><td>[ppm]</td><td>J(1, 10000)</td><td>870 ± 140</td></tr><tr><td>Additive jitter for SANCHEZ</td><td>[ppm]</td><td>J(1, 10000)</td><td>31<sup>+180</sup><sub>-29</sub></td></tr><tr><td>Additive jitter for RISE</td><td>[ppm]</td><td>J(1, 10000)</td><td>15<sup>+43</sup><sub>-12</sub></td></tr><tr><td>Additive jitter for TESS sector 24</td><td>[ppm]</td><td>J(1, 10000)</td><td>26<sup>+100</sup><sub>-23</sub></td></tr><tr><td>Additive jitter for TESS sector 25</td><td>[ppm]</td><td>J(1, 10000)</td><td>110<sup>+120</sup><sub>-100</sub></td></tr><tr><td>Additive jitter for TESS sector 26</td><td>[ppm]</td><td>J(1, 10000)</td><td>745 ± 27</td></tr><tr><td>Additive jitter for OSN1</td><td>[ppm]</td><td>J(1, 10000)</td><td>972 ± 64</td></tr><tr><td>Additive jitter for OSN2</td><td>[ppm]</td><td>J(1, 10000)</td><td>1372 ± 77</td></tr><tr><td>Additive jitter for OSN3</td><td>[ppm]</td><td>J(1, 10000)</td><td>1052 ± 61</td></tr><tr><td>Additive jitter for OSN4</td><td>[ppm]</td><td>J(1, 10000)</td><td>20<sup>+110</sup><sub>-18</sub></td></tr><tr><td>Timescale of the GP for WASP1</td><td>[days]</td><td>J(0.001, 1000)</td><td>3.0<sup>+100</sup><sub>-3.0</sub></td></tr><tr><td>Timescale of the GP for WASP2</td><td>[days]</td><td>J(0.001, 1000)</td><td>2.0<sup>+63</sup><sub>-1.9</sub></td></tr><tr><td>Timescale of the GP for WASP3</td><td>[days]</td><td>J(0.001, 1000)</td><td>0.040<sup>+0.061</sup><sub>-0.024</sub></td></tr><tr><td>Timescale of the GP for NITES1</td><td>[days]</td><td>J(0.001, 1000)</td><td>2.7<sup>+72</sup><sub>-2.7</sub></td></tr><tr><td>Timescale of the GP for NITES2</td><td>[days]</td><td>J(0.001, 1000)</td><td>13<sup>+160</sup><sub>-13</sub></td></tr><tr><td>Timescale of the GP for MARS</td><td>[days]</td><td>J(0.001, 1000)</td><td>1.2<sup>+12</sup><sub>-1.1</sub></td></tr><tr><td>Timescale of the GP for SANCHEZ</td><td>[days]</td><td>J(0.001, 1000)</td><td>0.0230<sup>+0.017</sup><sub>-0.0076</sub></td></tr><tr><td>Timescale of the GP for RISE</td><td>[days]</td><td>J(0.001, 1000)</td><td>1.03<sup>+1.5</sup><sub>-0.91</sub></td></tr><tr><td>Timescale of the GP for TESS sector 24</td><td>[days]</td><td>J(0.001, 1000)</td><td>1.79<sup>+0.40</sup><sub>-0.30</sub></td></tr><tr><td>Timescale of the GP for TESS sector 25</td><td>[days]</td><td>J(0.001, 1000)</td><td>4.7<sup>+2.4</sup><sub>-1.3</sub></td></tr><tr><td>Timescale of the GP for TESS sector 26</td><td>[days]</td><td>J(0.001, 1000)</td><td>1.05<sup>+0.17</sup><sub>-0.14</sub></td></tr><tr><td>Timescale of the GP for OSN1</td><td>[days]</td><td>J(0.001, 1000)</td><td>0.130<sup>+0.16</sup><sub>-0.083</sub></td></tr><tr><td>Timescale of the GP for OSN2</td><td>[days]</td><td>J(0.001, 1000)</td><td>0.081<sup>+3.6</sup><sub>-0.041</sub></td></tr><tr><td>Timescale of the GP for OSN3</td><td>[days]</td><td>J(0.001, 1000)</td><td>0.062<sup>+0.080</sup><sub>-0.028</sub></td></tr><tr><td>Timescale of the GP for OSN4</td><td>[days]</td><td>J(0.001, 1000)</td><td>0.046<sup>+0.040</sup><sub>-0.018</sub></td></tr><tr><td>Amplitude of the GP for WASP1</td><td>[Relative flux]</td><td>J(10<sup>-6</sup>, 1)</td><td>0.00036<sup>+0.0050</sup><sub>-0.00035</sub></td></tr><tr><td>Amplitude of the GP for WASP2</td><td>[Relative flux]</td><td>J(10<sup>-6</sup>, 1)</td><td>0.0019<sup>+0.0053</sup><sub>-0.0019</sub></td></tr><tr><td>Amplitude of the GP for WASP3</td><td>[Relative flux]</td><td>J(10<sup>-6</sup>, 1)</td><td>0.0058<sup>+0.0038</sup><sub>-0.0024</sub></td></tr><tr><td>Amplitude of the GP for NITES1</td><td>[Relative flux]</td><td>J(10<sup>-6</sup>, 1)</td><td>0.00051<sup>+0.0011</sup><sub>-0.00049</sub></td></tr><tr><td>Amplitude of the GP for NITES2</td><td>[Relative flux]</td><td>J(10<sup>-6</sup>, 1)</td><td>0.00023<sup>+0.0049</sup><sub>-0.00022</sub></td></tr><tr><td>Amplitude of the GP for MARS</td><td>[Relative flux]</td><td>J(10<sup>-6</sup>, 1)</td><td>0.0052<sup>+0.0053</sup><sub>-0.0040</sub></td></tr><tr><td>Amplitude of the GP for SANCHEZ</td><td>[Relative flux]</td><td>J(10<sup>-6</sup>, 1)</td><td>0.00181<sup>+0.00082</sup><sub>-0.00040</sub></td></tr><tr><td>Amplitude of the GP for RISE</td><td>[Relative flux]</td><td>J(10<sup>-6</sup>, 1)</td><td>0.051<sup>+0.11</sup><sub>-0.049</sub></td></tr><tr><td>Amplitude of the GP for TESS sector 24</td><td>[Relative flux]</td><td>J(10<sup>-6</sup>, 1)</td><td>0.00312<sup>+0.00084</sup><sub>-0.00056</sub></td></tr><tr><td>Amplitude of the GP for TESS sector 25</td><td>[Relative flux]</td><td>J(10<sup>-6</sup>, 1)</td><td>0.00266<sup>+0.0018</sup><sub>-0.00087</sub></td></tr><tr><td>Amplitude of the GP for TESS sector 26</td><td>[Relative flux]</td><td>J(10<sup>-6</sup>, 1)</td><td>0.00436<sup>+0.00076</sup><sub>-0.00062</sub></td></tr><tr><td>Amplitude of the GP for OSN1</td><td>[Relative flux]</td><td>J(10<sup>-6</sup>, 1)</td><td>0.0019<sup>+0.0039</sup><sub>-0.0013</sub></td></tr><tr><td>Amplitude of the GP for OSN2</td><td>[Relative flux]</td><td>J(10<sup>-6</sup>, 1)</td><td>0.00108<sup>+0.0013</sup><sub>-0.00037</sub></td></tr><tr><td>Amplitude of the GP for OSN3</td><td>[Relative flux]</td><td>J(10<sup>-6</sup>, 1)</td><td>0.00102<sup>+0.00076</sup><sub>-0.00022</sub></td></tr><tr><td>Amplitude of the GP for OSN4</td><td>[Relative flux]</td><td>J(10<sup>-6</sup>, 1)</td><td>0.00076<sup>+0.00058</sup><sub>-0.00022</sub></td></tr></tbody></table>

**Notes.** *U(a, b)*: A uniform distribution defined between a lower *a* and upper *b* limit. *J(a, b)*: Jeffreys (or log-uniform) distribution defined between a lower *a* and upper *b* limit.
<figure><img src="image_17.png" /><figcaption>**Fig. A.1.** Two-parameter joint posterior distributions for the most relevant model parameters from the photodynamical modelling (Section 6). The 39.3, 86.5, and 98.9% two-variable joint confidence regions are denoted by three different grey levels; in the case of a Gaussian posterior, these regions project on to the one-dimensional 1, 2, and 3 σ intervals. The histogram of the marginal distribution for each parameter is shown at the top of each column, except for the parameter on the last line, which is shown at the end of the line. Units are the same as in Table 4.</figcaption></figure>
<figure><img src="image_18.png" /></figure>

**Fig. A.2.** Transits of WASP-148 b (blue points) and the MAP model (black line) from the photodynamical modelling (Section 6). Each panel is centred at the linear ephemeris (indicated by the vertical grey lines, and reported in the caption of Fig. 7). For the SANCHEZ transit, 90-second binned data are shown in grey in addition to the observed data points. Each panel is labelled with the epoch; zero is the transit at $t_{\text{ref}}$. In the lower part of each panel, the residuals after subtracting the MAP model are shown.

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

**Fig. A.3.** SOPHIE radial velocities of WASP-148 (green error bars) and the MAP model (black line) from the photodynamical modelling (Section 6). In the lower panel, the residuals after subtracting the MAP model are shown.
<table><caption>Table A.5. Planet b transit-time predictions from 2021 to 2026.</caption><thead><tr><th colspan="3">Date</th><th colspan="3">68.3% CI</th><th colspan="3">68.3% CI</th><th colspan="3">68.3% CI</th></tr><tr><th></th><th>Mid-transit time<br>[BJD<sub>TDB</sub>]</th><th>[days]</th><th></th><th>Mid-transit time<br>[BJD<sub>TDB</sub>]</th><th>[days]</th><th></th><th>Mid-transit time<br>[BJD<sub>TDB</sub>]</th><th>[days]</th><th></th><th>Mid-transit time<br>[BJD<sub>TDB</sub>]</th><th>[days]</th></tr></thead><tbody><tr><td>2021-01-01</td><td>2459216.4351</td><td>-0.0014</td><td>+0.0012</td><td>2023-01-02</td><td>2459947.1690</td><td>-0.0005</td><td>+0.0005</td><td>2025-01-02</td><td>2460677.8592</td><td>-0.0008</td><td>+0.0007</td></tr><tr><td>2021-01-10</td><td>2459225.2363</td><td>-0.0012</td><td>+0.0011</td><td>2023-01-11</td><td>2459955.9730</td><td>-0.0006</td><td>+0.0006</td><td>2025-01-11</td><td>2460686.6651</td><td>-0.0007</td><td>+0.0008</td></tr><tr><td>2021-01-19</td><td>2459234.0382</td><td>-0.0011</td><td>+0.0011</td><td>2023-01-20</td><td>2459964.7781</td><td>-0.0006</td><td>+0.0006</td><td>2025-01-19</td><td>2460695.4686</td><td>-0.0008</td><td>+0.0008</td></tr><tr><td>2021-01-28</td><td>2459242.8399</td><td>-0.0012</td><td>+0.0011</td><td>2023-01-29</td><td>2459973.5816</td><td>-0.0007</td><td>+0.0006</td><td>2025-01-28</td><td>2460704.2750</td><td>-0.0008</td><td>+0.0008</td></tr><tr><td>2021-02-06</td><td>2459251.6440</td><td>-0.0010</td><td>+0.0009</td><td>2023-02-06</td><td>2459982.3843</td><td>-0.0008</td><td>+0.0007</td><td>2025-02-06</td><td>2460713.0804</td><td>-0.0008</td><td>+0.0008</td></tr><tr><td>2021-02-14</td><td>2459260.4460</td><td>-0.0008</td><td>+0.0007</td><td>2023-02-15</td><td>2459991.1872</td><td>-0.0009</td><td>+0.0009</td><td>2025-02-15</td><td>2460721.8870</td><td>-0.0008</td><td>+0.0008</td></tr><tr><td>2021-02-23</td><td>2459269.2488</td><td>-0.0006</td><td>+0.0007</td><td>2023-02-24</td><td>2459999.9909</td><td>-0.0011</td><td>+0.0011</td><td>2025-02-24</td><td>2460730.6909</td><td>-0.0008</td><td>+0.0009</td></tr><tr><td>2021-03-04</td><td>2459278.0509</td><td>-0.0006</td><td>+0.0006</td><td>2023-03-05</td><td>2460008.7940</td><td>-0.0012</td><td>+0.0011</td><td>2025-03-04</td><td>2460739.4972</td><td>-0.0008</td><td>+0.0008</td></tr><tr><td>2021-03-13</td><td>2459286.8557</td><td>-0.0005</td><td>+0.0005</td><td>2023-03-14</td><td>2460017.5955</td><td>-0.0014</td><td>+0.0012</td><td>2025-03-13</td><td>2460748.3030</td><td>-0.0008</td><td>+0.0009</td></tr><tr><td>2021-03-22</td><td>2459295.6588</td><td>-0.0004</td><td>+0.0004</td><td>2023-03-22</td><td>2460026.3974</td><td>-0.0016</td><td>+0.0014</td><td>2025-03-22</td><td>2460757.1099</td><td>-0.0008</td><td>+0.0009</td></tr><tr><td>2021-03-30</td><td>2459304.4627</td><td>-0.0004</td><td>+0.0004</td><td>2023-03-31</td><td>2460035.1999</td><td>-0.0018</td><td>+0.0017</td><td>2025-03-31</td><td>2460765.9140</td><td>-0.0009</td><td>+0.0010</td></tr><tr><td>2021-04-08</td><td>2459313.2654</td><td>-0.0004</td><td>+0.0005</td><td>2023-04-09</td><td>2460044.0029</td><td>-0.0018</td><td>+0.0017</td><td>2025-04-09</td><td>2460774.7198</td><td>-0.0008</td><td>+0.0009</td></tr><tr><td>2021-04-17</td><td>2459322.0709</td><td>-0.0004</td><td>+0.0004</td><td>2023-04-18</td><td>2460052.8039</td><td>-0.0019</td><td>+0.0017</td><td>2025-04-18</td><td>2460783.5252</td><td>-0.0008</td><td>+0.0009</td></tr><tr><td>2021-04-26</td><td>2459330.8749</td><td>-0.0004</td><td>+0.0005</td><td>2023-04-27</td><td>2460061.6053</td><td>-0.0020</td><td>+0.0018</td><td>2025-04-26</td><td>2460792.3318</td><td>-0.0008</td><td>+0.0008</td></tr><tr><td>2021-05-05</td><td>2459339.6799</td><td>-0.0004</td><td>+0.0005</td><td>2023-05-05</td><td>2460070.4070</td><td>-0.0022</td><td>+0.0020</td><td>2025-05-05</td><td>2460801.1358</td><td>-0.0008</td><td>+0.0009</td></tr><tr><td>2021-05-13</td><td>2459348.4830</td><td>-0.0004</td><td>+0.0005</td><td>2023-05-14</td><td>2460079.2103</td><td>-0.0021</td><td>+0.0019</td><td>2025-05-14</td><td>2460809.9403</td><td>-0.0008</td><td>+0.0008</td></tr><tr><td>2021-05-22</td><td>2459357.2891</td><td>-0.0005</td><td>+0.0004</td><td>2023-05-23</td><td>2460088.0113</td><td>-0.0020</td><td>+0.0019</td><td>2025-05-23</td><td>2460818.7448</td><td>-0.0008</td><td>+0.0008</td></tr><tr><td>2021-05-31</td><td>2459366.0938</td><td>-0.0004</td><td>+0.0004</td><td>2023-06-01</td><td>2460096.8130</td><td>-0.0019</td><td>+0.0019</td><td>2025-06-01</td><td>2460827.5504</td><td>-0.0007</td><td>+0.0008</td></tr><tr><td>2021-06-09</td><td>2459374.8997</td><td>-0.0004</td><td>+0.0004</td><td>2023-06-10</td><td>2460105.6146</td><td>-0.0020</td><td>+0.0020</td><td>2025-06-09</td><td>2460836.3541</td><td>-0.0008</td><td>+0.0008</td></tr><tr><td>2021-06-18</td><td>2459383.7032</td><td>-0.0004</td><td>+0.0004</td><td>2023-06-18</td><td>2460114.4184</td><td>-0.0017</td><td>+0.0018</td><td>2025-06-18</td><td>2460845.1572</td><td>-0.0009</td><td>+0.0009</td></tr><tr><td>2021-06-27</td><td>2459392.5095</td><td>-0.0003</td><td>+0.0003</td><td>2023-06-27</td><td>2460123.2200</td><td>-0.0015</td><td>+0.0016</td><td>2025-06-27</td><td>2460853.9604</td><td>-0.0011</td><td>+0.0010</td></tr><tr><td>2021-07-05</td><td>2459401.3150</td><td>-0.0003</td><td>+0.0003</td><td>2023-07-06</td><td>2460132.0224</td><td>-0.0013</td><td>+0.0016</td><td>2025-07-06</td><td>2460862.7646</td><td>-0.0013</td><td>+0.0013</td></tr><tr><td>2021-07-14</td><td>2459410.1215</td><td>-0.0003</td><td>+0.0004</td><td>2023-07-15</td><td>2460140.8243</td><td>-0.0013</td><td>+0.0015</td><td>2025-07-15</td><td>2460871.5678</td><td>-0.0014</td><td>+0.0014</td></tr><tr><td>2021-07-23</td><td>2459418.9253</td><td>-0.0004</td><td>+0.0004</td><td>2023-07-24</td><td>2460149.6289</td><td>-0.0011</td><td>+0.0013</td><td>2025-07-23</td><td>2460880.3697</td><td>-0.0016</td><td>+0.0016</td></tr><tr><td>2021-08-01</td><td>2459427.7316</td><td>-0.0004</td><td>+0.0004</td><td>2023-08-01</td><td>2460158.4316</td><td>-0.0009</td><td>+0.0010</td><td>2025-08-01</td><td>2460889.1719</td><td>-0.0019</td><td>+0.0018</td></tr><tr><td>2021-08-10</td><td>2459436.5372</td><td>-0.0004</td><td>+0.0004</td><td>2023-08-10</td><td>2460167.2351</td><td>-0.0009</td><td>+0.0009</td><td>2025-08-10</td><td>2460897.9747</td><td>-0.0023</td><td>+0.0021</td></tr><tr><td>2021-08-18</td><td>2459445.3440</td><td>-0.0004</td><td>+0.0005</td><td>2023-08-19</td><td>2460176.0375</td><td>-0.0009</td><td>+0.0009</td><td>2025-08-19</td><td>2460906.7777</td><td>-0.0023</td><td>+0.0022</td></tr><tr><td>2021-08-27</td><td>2459454.1481</td><td>-0.0004</td><td>+0.0005</td><td>2023-08-28</td><td>2460184.8428</td><td>-0.0007</td><td>+0.0007</td><td>2025-08-28</td><td>2460915.5787</td><td>-0.0024</td><td>+0.0023</td></tr><tr><td>2021-09-05</td><td>2459462.9537</td><td>-0.0004</td><td>+0.0004</td><td>2023-09-06</td><td>2460193.6464</td><td>-0.0006</td><td>+0.0006</td><td>2025-09-05</td><td>2460924.3803</td><td>-0.0026</td><td>+0.0025</td></tr><tr><td>2021-09-14</td><td>2459471.7591</td><td>-0.0004</td><td>+0.0005</td><td>2023-09-14</td><td>2460202.4510</td><td>-0.0006</td><td>+0.0006</td><td>2025-09-14</td><td>2460933.1821</td><td>-0.0029</td><td>+0.0029</td></tr><tr><td>2021-09-23</td><td>2459480.5656</td><td>-0.0004</td><td>+0.0005</td><td>2023-09-23</td><td>2460211.2539</td><td>-0.0007</td><td>+0.0007</td><td>2025-09-23</td><td>2460941.9852</td><td>-0.0028</td><td>+0.0028</td></tr><tr><td>2021-10-01</td><td>2459489.3695</td><td>-0.0004</td><td>+0.0005</td><td>2023-10-02</td><td>2460220.0598</td><td>-0.0006</td><td>+0.0006</td><td>2025-10-02</td><td>2460950.7861</td><td>-0.0027</td><td>+0.0028</td></tr><tr><td>2021-10-10</td><td>2459498.1741</td><td>-0.0004</td><td>+0.0005</td><td>2023-10-11</td><td>2460228.8643</td><td>-0.0006</td><td>+0.0006</td><td>2025-10-11</td><td>2460959.5876</td><td>-0.0027</td><td>+0.0028</td></tr><tr><td>2021-10-19</td><td>2459506.9785</td><td>-0.0005</td><td>+0.0005</td><td>2023-10-20</td><td>2460237.6699</td><td>-0.0006</td><td>+0.0006</td><td>2025-10-19</td><td>2460968.3892</td><td>-0.0028</td><td>+0.0030</td></tr><tr><td>2021-10-28</td><td>2459515.7840</td><td>-0.0005</td><td>+0.0005</td><td>2023-10-28</td><td>2460246.4732</td><td>-0.0007</td><td>+0.0007</td><td>2025-10-28</td><td>2460977.1927</td><td>-0.0026</td><td>+0.0029</td></tr><tr><td>2021-11-06</td><td>2459524.5876</td><td>-0.0005</td><td>+0.0005</td><td>2023-11-06</td><td>2460255.2795</td><td>-0.0006</td><td>+0.0006</td><td>2025-11-06</td><td>2460985.9941</td><td>-0.0024</td><td>+0.0026</td></tr><tr><td>2021-11-14</td><td>2459533.3907</td><td>-0.0005</td><td>+0.0005</td><td>2023-11-15</td><td>2460264.0848</td><td>-0.0006</td><td>+0.0006</td><td>2025-11-15</td><td>2460994.7961</td><td>-0.0021</td><td>+0.0025</td></tr><tr><td>2021-11-23</td><td>2459542.1940</td><td>-0.0006</td><td>+0.0006</td><td>2023-11-24</td><td>2460272.8911</td><td>-0.0006</td><td>+0.0006</td><td>2025-11-24</td><td>2461003.5978</td><td>-0.0022</td><td>+0.0026</td></tr><tr><td>2021-12-02</td><td>2459550.9982</td><td>-0.0008</td><td>+0.0008</td><td>2023-12-03</td><td>2460281.6949</td><td>-0.0006</td><td>+0.0007</td><td>2025-12-02</td><td>2461012.4021</td><td>-0.0020</td><td>+0.0023</td></tr><tr><td>2021-12-11</td><td>2459559.8014</td><td>-0.0008</td><td>+0.0008</td><td>2023-12-12</td><td>2460290.5013</td><td>-0.0006</td><td>+0.0006</td><td>2025-12-11</td><td>2461021.2043</td><td>-0.0018</td><td>+0.0020</td></tr><tr><td>2021-12-20</td><td>2459568.6032</td><td>-0.0009</td><td>+0.0009</td><td>2023-12-20</td><td>2460299.3070</td><td>-0.0006</td><td>+0.0007</td><td>2025-12-20</td><td>2461030.0074</td><td>-0.0016</td><td>+0.0018</td></tr><tr><td>2021-12-28</td><td>2459577.4054</td><td>-0.0011</td><td>+0.0010</td><td>2023-12-29</td><td>2460308.1138</td><td>-0.0006</td><td>+0.0007</td><td>2025-12-29</td><td>2461038.8096</td><td>-0.0015</td><td>+0.0018</td></tr><tr><td>2022-01-06</td><td>2459586.2083</td><td>-0.0013</td><td>+0.0012</td><td>2024-01-07</td><td>2460316.9718</td><td>-0.0006</td><td>+0.0008</td><td>2026-01-07</td><td>2461047.6146</td><td>-0.0014</td><td>+0.0015</td></tr><tr><td>2022-01-15</td><td>2459595.0113</td><td>-0.0014</td><td>+0.0012</td><td>2024-01-16</td><td>2460325.7238</td><td>-0.0006</td><td>+0.0007</td><td>2026-01-15</td><td>2461056.4179</td><td>-0.0013</td><td>+0.0013</td></tr><tr><td>2022-01-24</td><td>2459603.8124</td><td>-0.0014</td><td>+0.0013</td><td>2024-01-25</td><td>2460334.5294</td><td>-0.0007</td><td>+0.0007</td><td>2026-01-24</td><td>2461065.2221</td><td>-0.0012</td><td>+0.0012</td></tr><tr><td>2022-02-02</td><td>2459612.6139</td><td>-0.0015</td><td>+0.0014</td><td>2024-02-02</td><td>2460343.3361</td><td>-0.0007</td><td>+0.0007</td><td>2026-02-02</td><td>2461074.0248</td><td>-0.0012</td><td>+0.0012</td></tr><tr><td>2022-02-10</td><td>2459621.4159</td><td>-0.0018</td><td>+0.0016</td><td>2024-02-11</td><td>2460352.1402</td><td>-0.0007</td><td>+0.0008</td><td>2026-02-11</td><td>2461082.8305</td><td>-0.0011</td><td>+0.0010</td></tr><tr><td>2022-02-19</td><td>2459630.2190</td><td>-0.0017</td><td>+0.0015</td><td>2024-02-20</td><td>2460360.9452</td><td>-0.0007</td><td>+0.0007</td><td>2026-02-20</td><td>2461091.6347</td><td>-0.0010</td><td>+0.0010</td></tr><tr><td>2022-02-28</td><td>2459639.0200</td><td>-0.0017</td><td>+0.0015</td><td>2024-02-29</td><td>2460369.7500</td><td>-0.0007</td><td>+0.0007</td><td>2026-02-28</td><td>2461100.4399</td><td>-0.0010</td><td>+0.0010</td></tr><tr><td>2022-03-09</td><td>2459647.8215</td><td>-0.0017</td><td>+0.0015</td><td>2024-03-09</td><td>2460378.5560</td><td>-0.0006</td><td>+0.0006</td><td>2026-03-09</td><td>2461109.2431</td><td>-0.0011</td><td>+0.0010</td></tr><tr><td>2022-03-18</td><td>2459656.6231</td><td>-0.0018</td><td>+0.0017</td><td>2024-03-17</td><td>2460387.3598</td><td>-0.0007</td><td>+0.0007</td><td>2026-03-18</td><td>2461118.0493</td><td>-0.0010</td><td>+0.0010</td></tr><tr><td>2022-03-26</td><td>2459665.4267</td><td>-0.0016</td><td>+0.0015</td><td>2024-03-26</td><td>2460396.1633</td><td>-0.0007</td><td>+0.0007</td><td>2026-03-27</td><td>2461126.8543</td><td>-0.0010</td><td>+0.0010</td></tr><tr><td>2022-04-04</td><td>2459674.2281</td><td>-0.0014</td><td>+0.0014</td><td>2024-04-04</td><td>2460404.9670</td><td>-0.0008</td><td>+0.0007</td><td>2026-04-05</td><td>2461135.6604</td><td>-0.0009</td><td>+0.0010</td></tr><tr><td>2022-04-13</td><td>2459683.0302</td><td>-0.0012</td><td>+0.0013</td><td>2024-04-13</td><td>2460413.7717</td><td>-0.0009</td><td>+0.0009</td><td>2026-04-13</td><td>2461144.4641</td><td>-0.0010</td><td>+0.0011</td></tr><tr><td>2022-04-22</td><td>2459691.8320</td><td>-0.0013</td><td>+0.0013</td><td>2024-04-22</td><td>2460422.5750</td><td>-0.0010</td><td>+0.0009</td><td>2026-04-22</td><td>2461153.2706</td><td>-0.0010</td><td>+0.0010</td></tr><tr><td>2022-05-01</td><td>2459700.6363</td><td>-0.0010</td><td>+0.0012</td><td>2024-04-30</td><td>2460431.3772</td><td>-0.0011</td><td>+0.0011</td><td>2026-05-01</td><td>2461162.0762</td><td>-0.0010</td><td>+0.0010</td></tr><tr><td>2022-05-09</td><td>2459709.4387</td><td>-0.0008</td><td>+0.0009</td><td>2024-05-09</td><td>2460440.1798</td><td>-0.0014</td><td>+0.0013</td><td>2026-05-10</td><td>2461170.8830</td><td>-0.0011</td><td>+0.0010</td></tr><tr><td>2022-05-18</td><td>2459718.2418</td><td>-0.0007</td><td>+0.0008</td><td>2024-05-18</td><td>2460448.9831</td><td>-0.0017</td><td>+0.0016</td><td>2026-05-19</td><td>2461179.6870</td><td>-0.0011</td><td>+0.0012</td></tr><tr><td>2022-05-27</td><td>2459727.0441</td><td>-0.0007</td><td>+0.0008</td><td>2024-05-27</td><td>2460457.7861</td><td>-0.0017</td><td>+0.0016</td><td>2026-05-27</td><td>2461188.4932</td><td>-0.0010</td><td>+0.0011</td></tr><tr><td>2022-06-05</td><td>2459735.8492</td><td>-0.0006</td><td>+0.0006</td><td>2024-06-05</td><td>2460466.5873</td><td>-0.0019</td><td>+0.0017</td><td>2026-06-05</td><td>2461197.2990</td><td>-0.0011</td><td>+0.0011</td></tr><tr><td>2022-06-14</td><td>2459744.6525</td><td>-0.0005</td><td>+0.0005</td><td>2024-06-13</td><td>2460475.3890</td><td>-0.0021</td><td>+0.0020</td><td>2026-06-14</td><td>2461206.1059</td><td>-0.0010</td><td>+0.0011</td></tr><tr><td>2022-06-22</td><td>2459753.4568</td><td>-0.0005</td><td>+0.0005</td><td>2024-06-22</td><td>2460484.1912</td><td>-0.0024</td><td>+0.0022</td><td>2026-06-23</td><td>2461214.9101</td><td>-0.0011</td><td>+0.0012</td></tr><tr><td>2022-07-01</td><td>2459762.2596</td><td>-0.0005</td><td>+0.0005</td><td>2024-07-01</td><td>2460492.9942</td><td>-0.0023</td><td>+0.0021</td><td>2026-07-02</td><td>2461223.7155</td><td>-0.0010</td><td>+0.0010</td></tr><tr><td>2022-07-10</td><td>2459771.0653</td><td>-0.0005</td><td>+0.0005</td><td>2024-07-10</td><td>2460501.7951</td><td>-0.0023</td><td>+0.0022</td><td>2026-07-11</td><td>2461232.5207</td><td>-0.0010</td><td>+0.0010</td></tr><tr><td>2022-07-19</td><td>2459779.8695</td><td>-0.0005</td><td>+0.0005</td><td>2024-07-19</td><td>2460510.5965</td><td>-0.0024</td><td>+0.0023</td><td>2026-07-19</td><td>2461241.3270</td><td>-0.0009</td><td>+0.0009</td></tr><tr><td>2022-07-28</td><td>2459788.6748</td><td>-0.0005</td><td>+0.0005</td><td>2024-07-27</td><td>2460519.3981</td><td>-0.0025</td><td>+0.0026</td><td>2026-07-28</td><td>2461250.1309</td><td>-0.0010</td><td>+0.0009</td></tr><tr><td>2022-08-05</td><td>2459797.4781</td><td>-0.0005</td><td>+0.0005</td><td>2024-08-05</td><td>2460528.014</td><td>-0.0023</td><td>+0.0024</td><td>2026-08-06</td><td>2461258.9350</td><td>-0.0009</td><td>+0.0009</td></tr><tr><td>2022-08-14</td><td>2459806.2842</td><td>-0.0005</td><td>+0.0004</td><td>2024-08-14</td><td>2460537.0027</td><td>-0.0022</td><td>+0.0023</td><td>2026-08-15</td><td>2461267.7391</td><td>-0.0009</td><td>+0.0009</td></tr><tr><td>2022-08-23</td><td>2459815.0892</td><td>-0.0005</td><td>+0.0005</td><td>2024-08-23</td><td>2460545.8045</td><td>-0.0020</td><td>+0.0022</td><td>2026-08-24</td><td>2461276.5442</td><td>-0.0011</td><td>+0.0010</td></tr><tr><td>2022-09-01</td><td>2459823.8953</td><td>-0.0004</td><td>+0.0004</td><td>2024-09-01</td><td>2460554.6061</td><td>-0.0021</td><td>+0.0023</td><td>2026-09-01</td><td>2461285.3478</td><td>-0.0012</td><td>+0.0010</td></tr><tr><td>2022-09-10</td><td>2459832.6990</td><td>-0.0004</td><td>+0.0005</td><td>2024-09-09</td><td>2460563.4101</td><td>-0.0018</td><td>+0.0021</td><td>2026-09-10</td><td>2461294.1504</td><td>-0.0013</td><td>+0.0013</td></tr><tr><td>2022-09-19</td><td>2459841.5053</td><td>-0.0004</td><td>+0.0004</td><td>2024-09-18</td><td>2460572.2121</td><td>-0.0017</td><td>+0.0018</td><td>2026-09-19</td><td>2461302.9533</td><td>-0.0016</td><td>+0.0015</td></tr><tr><td>2022-09-27</td><td>2459850.3109</td><td>-0.0004</td><td>+0.0004</td><td>2024-09-27</td><td>2460581.0148</td><td>-0.0014</td><td>+0.0017</td><td>2026-09-28</td><td>2461311.7569</td><td>-0.0019</td><td>+0.0019</td></tr><tr><td>2022-10-06</td><td>2459859.1176</td><td>-0.0004</td><td>+0.0005</td><td>2024-10-06</td><td>2460589.8169</td><td>-0.0014</td><td>+0.0016</td><td>2026-10-07</td><td>2461320.5600</td><td>-0.0020</td><td>+0.0020</td></tr><tr><td>2022-10-15</td><td>2459867.9216</td><td>-0.0005</td><td>+0.0006</td><td>2024-10-15</td><td>2460598.6217</td><td>-0.0013</td><td>+0.0014</td><td>2026-10-15</td><td>2461329.3616</td><td>-0.0022</td><td>+0.0022</td></tr><tr><td>2022-10-24</td><td>2459876.7277</td><td>-0.0005</td><td>+0.0005</td><td>2024-10-23</td><td>2460607.4246</td><td>-0.0011</td><td>+0.0011</td><td>2026-10-24</td><td>2461338.1634</td><td>-0.0025</td><td>+0.0025</td></tr><tr><td>2022-11-02</td><td>2459885.5334</td><td>-0.0005</td><td>+0.0005</td><td>2024-11-01</td><td>2460616.2285</td><td>-0.0010</td><td>+0.0010</td><td>2026-11-02</td><td>2461346.9659</td><td>-0.0029</td><td>+0.0029</td></tr><tr><td>2022-11-10</td><td>2459894.3402</td><td>-0.0005</td><td>+0.0005</td><td>2024-11-10</td><td>2460625.0311</td><td>-0.0010</td><td>+0.0010</td><td>2026-11-11</td><td>2461355.7689</td><td>-0.0029</td><td>+0.0028</td></tr><tr><td>2022-11-19</td><td>2459903.1443</td><td>-0.0006</td><td>+0.0006</td><td>2024-11-19</td><td>2460633.8365</td><td>-0.0009</td><td>+0.0009</td><td>2026-11-20</td><td>2461364.5698</td><td>-0.0030</td><td>+0.0029</td></tr><tr><td>2022-11-28</td><td>2459911.9496</td><td>-0.0005</td><td>+0.0006</td><td>2024-11-28</td><td>2460642.6405</td><td>-0.0008</td><td>+0.0008</td><td>2026-11-28</td><td>2461373.3712</td><td>-0.0030</td><td>+0.0032</td></tr><tr><td>2022-12-07</td><td>2459920.7548</td><td>-0.0005</td><td>+0.0006</td><td>2024-12-06</td><td>2460651.4453</td><td>-0.0008</td><td>+0.0008</td><td>2026-12-07</td><td>2461382.1729</td><td>-0.0033</td><td>+0.0035</td></tr><tr><td>2022-12-16</td><td>2459929.5610</td><td>-0.0005</td><td>+0.0005</td><td>2024-12-15</td><td>2460660.2484</td><td>-0.0009</td><td>+0.0009</td><td>2026-12-16</td><td>2461390.9761</td><td>-0.0032</td><td>+0.0034</td></tr><tr><td>2022-12-24</td><td>2459938.3649</td><td>-0.0005</td><td>+0.0006</td><td>2024-12-24</td><td>2460669.0544</td><td>-0.0008</td><td>+0.0007</td><td>2026-12-25</td><td>2461399.7771</td><td>-0.0030</td><td>+0.0033</td></tr></tbody></table>

Notes. The table lists: Calendar UT date (YYYY-MM-DD), mid-transit time, and 68.3% CI.
<table><caption>**Table A.6.** Inferred system parameters with the photodynamical modelling assuming coplanar orbits.</caption><thead><tr><th>Parameter</th><th></th><th>Median<br>and 68.3% CI</th></tr></thead><tbody><tr><td><strong>Star</strong></td><td></td><td></td></tr><tr><td>Mass, M<sub>★</sub></td><td>(M<sub>☉</sub>)</td><td>0.968 ± 0.050</td></tr><tr><td>Radius, R<sub>★</sub></td><td>(R<sub>☉</sub><sup>N</sup>)</td><td>0.918 ± 0.017</td></tr><tr><td>Stellar mean density, ρ<sub>★</sub></td><td>[g cm<sup>-3</sup>]</td><td>1.76 <sup>+0.14</sup><sub>-0.12</sub></td></tr><tr><td><strong>Planet b</strong></td><td></td><td></td></tr><tr><td>Semi-major axis, a</td><td>[au]</td><td>0.0825 ± 0.0015</td></tr><tr><td>Eccentricity, e</td><td></td><td>0.194 ± 0.035</td></tr><tr><td>Argument of pericentre, ω</td><td>[°]</td><td>61.4 ± 7.2</td></tr><tr><td>Inclination, i</td><td>[°]</td><td>89.06 <sup>+0.43</sup><sub>-0.31</sub></td></tr><tr><td>Longitude of the ascending node, Ω</td><td>[°]</td><td>180 (fixed at t<sub>ref</sub>)</td></tr><tr><td>Mean anomaly, M<sub>0</sub></td><td>[°]</td><td>19.0 <sup>+6.1</sup><sub>-5.0</sub></td></tr><tr><td>Impact parameter, b</td><td></td><td>0.259 <sup>+0.096</sup><sub>-0.12</sub></td></tr><tr><td>Mass, M<sub>p</sub></td><td>(M<sub>1</sub>)</td><td>0.289 ± 0.023</td></tr><tr><td>Radius, R<sub>p</sub></td><td>(R<sub>☉</sub><sup>N</sup>)</td><td>0.756 ± 0.017</td></tr><tr><td>Planet mean density, ρ<sub>p</sub></td><td>[g cm<sup>-3</sup>]</td><td>0.825 <sup>+0.094</sup><sub>-0.078</sub></td></tr><tr><td><strong>Planet c</strong></td><td></td><td></td></tr><tr><td>Semi-major axis, a</td><td>[au]</td><td>0.2054 ± 0.0036</td></tr><tr><td>Eccentricity, e</td><td></td><td>0.168 <sup>+0.016</sup><sub>-0.010</sub></td></tr><tr><td>Argument of pericentre, ω</td><td>[°]</td><td>26.3 ± 6.0</td></tr><tr><td>Inclination, i</td><td>[°]</td><td>≡ Planet b inclination</td></tr><tr><td>Longitude of the ascending node, Ω</td><td>[°]</td><td>180 (fixed at t<sub>ref</sub>)</td></tr><tr><td>Mean anomaly, M<sub>0</sub></td><td>[°]</td><td>258.1 ± 3.6</td></tr><tr><td>Impact parameter, b</td><td></td><td>0.71 <sup>+0.23</sup><sub>-0.33</sub></td></tr><tr><td>Mass, M<sub>p</sub></td><td>(M<sub>1</sub>)</td><td>0.375 ± 0.040</td></tr></tbody></table>

**Notes.** The table lists: Posterior median and 68.3% CI for the photodynamical modelling assuming coplanar orbits (Section 7.3). Only masses, radii, densities, and orbital parameters (at $t_{ref} = 2 457 957.48167 \text{ BJD}_{\text{TDB}}$) are listed.

<table><caption>**Table A.7.** Inferred system parameters with the photodynamical modelling without stellar priors.</caption><thead><tr><th>Parameter</th><th></th><th>Median<br>and 68.3% CI</th></tr></thead><tbody><tr><td><strong>Star</strong></td><td></td><td></td></tr><tr><td>Mass, M<sub>★</sub></td><td>(M<sub>☉</sub>)</td><td>9.7 <sup>+12</sup><sub>-5.3</sub></td></tr><tr><td>Radius, R<sub>★</sub></td><td>(R<sub>☉</sub><sup>N</sup>)</td><td>2.14 <sup>+0.69</sup><sub>-0.54</sub></td></tr><tr><td>Stellar mean density, ρ<sub>★</sub></td><td>[g cm<sup>-3</sup>]</td><td>1.41 <sup>+0.29</sup><sub>-0.22</sub></td></tr><tr><td><strong>Planet b</strong></td><td></td><td></td></tr><tr><td>Semi-major axis, a</td><td>[au]</td><td>0.178 <sup>+0.053</sup><sub>-0.041</sub></td></tr><tr><td>Eccentricity, e</td><td></td><td>0.254 ± 0.045</td></tr><tr><td>Argument of pericentre, ω</td><td>[°]</td><td>65.6 ± 9.6</td></tr><tr><td>Inclination, i</td><td>[°]</td><td>88.90 <sup>+0.56</sup><sub>-0.45</sub></td></tr><tr><td>Longitude of the ascending node, Ω</td><td>[°]</td><td>180 (fixed at t<sub>ref</sub>)</td></tr><tr><td>Mean anomaly, M<sub>0</sub></td><td>[°]</td><td>13.9 <sup>+7.0</sup><sub>-5.6</sub></td></tr><tr><td>Mass, M<sub>p</sub></td><td>(M<sub>1</sub>)</td><td>1.32 <sup>+0.88</sup><sub>-0.53</sub></td></tr><tr><td>Radius, R<sub>p</sub></td><td>(R<sub>☉</sub><sup>N</sup>)</td><td>1.77 <sup>+0.58</sup><sub>-0.45</sub></td></tr><tr><td>Planet mean density, ρ<sub>p</sub></td><td>[g cm<sup>-3</sup>]</td><td>0.293 <sup>+0.14</sup><sub>-0.089</sub></td></tr><tr><td><strong>Planet c</strong></td><td></td><td></td></tr><tr><td>Semi-major axis, a</td><td>[au]</td><td>0.44 <sup>+0.13</sup><sub>-0.10</sub></td></tr><tr><td>Eccentricity, e</td><td></td><td>0.280 ± 0.023</td></tr><tr><td>Argument of pericentre, ω</td><td>[°]</td><td>23.6 ± 7.0</td></tr><tr><td>Inclination, i</td><td>[°]</td><td>129 ± 13</td></tr><tr><td>Longitude of the ascending node, Ω</td><td>[°]</td><td>214.4 ± 8.2</td></tr><tr><td>Mean anomaly, M<sub>0</sub></td><td>[°]</td><td>262.1 ± 4.4</td></tr><tr><td>Mass, M<sub>p</sub></td><td>(M<sub>1</sub>)</td><td>2.4 <sup>+2.1</sup><sub>-1.2</sub></td></tr><tr><td>Mutual inclination, I</td><td>[°]</td><td>52 <sup>+10</sup><sub>-13</sub></td></tr></tbody></table>

**Notes.** The table lists: Posterior median and 68.3% CI for the photodynamical modelling without stellar priors (Section 7.4). Only masses, radii, densities, and orbital parameters (at $t_{ref} = 2 457 957.48167 \text{ BJD}_{\text{TDB}}$) are listed.
<figure><img src="image_20.png" /></figure>

**Fig. A.4.** Comparison of TTVs for different models. Left: Idem. as Figure 7, but up to the year 2100 with the two-planet photodynamical modelling (Section 6), three-planet photodynamical modelling (Section 7.8), and coplanar two-planet photodynamical modelling (Section 7.3). Right: Zoom of the plot on the left.

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

**Fig. A.5.** Comparison of transit duration variations for different models. Left: Posterior transit duration variations of planet b up to the end of the year 2100 from 1,000 random draws from the posterior distribution of the two-planet photodynamical modelling (Section 6), three-planet photodynamical modelling (Section 7.8), and coplanar two-planet photodynamical modelling (Section 7.3). The colours correspond with the ones used in Figure A.4. The points with error bars are the durations measured on individual transits (Section 7.3). Right: Zoom of the plot on the left.
<table><caption>**Table A.8.** Inferred system parameters with the photodynamical modelling with three planets.</caption><thead><tr><th>Parameter</th><th></th><th>Median<br>and 68.3% CI</th></tr></thead><tbody><tr><td><strong>Star</strong></td><td></td><td></td></tr><tr><td>Mass, M<sub>★</sub></td><td>(M<sub>☉</sub>)</td><td>0.976 ± 0.053</td></tr><tr><td>Radius, R<sub>★</sub></td><td>(R<sub>☉</sub><sup>N</sup>)</td><td>0.918 ± 0.017</td></tr><tr><td>Stellar mean density, ρ<sub>★</sub></td><td>[g cm<sup>-3</sup>]</td><td>1.79 ± 0.13</td></tr><tr><td><strong>Planet b</strong></td><td></td><td></td></tr><tr><td>Semi-major axis, a</td><td>[au]</td><td>0.0828 ± 0.0015</td></tr><tr><td>Eccentricity, e</td><td></td><td>0.183 <sup>+0.028</sup><sub>-0.031</sub></td></tr><tr><td>Argument of pericentre, ω</td><td>[°]</td><td>50.7 <sup>+9.8</sup><sub>-8.6</sub></td></tr><tr><td>Inclination, i</td><td>[°]</td><td>88.87 ± 0.21</td></tr><tr><td>Longitude of the ascending node, Ω</td><td>[°]</td><td>180 (fixed at t<sub>ref</sub>)</td></tr><tr><td>Mean anomaly, M<sub>0</sub></td><td>[°]</td><td>27.3 <sup>+6.4</sup><sub>-7.1</sub></td></tr><tr><td>Impact parameter, b</td><td></td><td>0.323 ± 0.064</td></tr><tr><td>T'<sub>0</sub> - 2 450 000</td><td>[BJD<sub>TDB</sub>]</td><td>7957.48176 ± 0.00025</td></tr><tr><td>P'</td><td>[d]</td><td>8.80358 ± 0.00018</td></tr><tr><td>K'</td><td>[m s<sup>-1</sup>]</td><td>29.1 ± 1.9</td></tr><tr><td>Mass, M<sub>p</sub></td><td>(M<sub>J</sub>)</td><td>0.284 <sup>+0.021</sup><sub>-0.018</sub></td></tr><tr><td>Radius, R<sub>p</sub></td><td>(R<sub>e</sub><sup>N</sup>)</td><td>0.765 ± 0.016</td></tr><tr><td>Planet mean density, ρ<sub>p</sub></td><td>[g cm<sup>-3</sup>]</td><td>0.788 ± 0.071</td></tr><tr><td><strong>Planet c</strong></td><td></td><td></td></tr><tr><td>Semi-major axis, a</td><td>[au]</td><td>0.2059 ± 0.0037</td></tr><tr><td>Eccentricity, e</td><td></td><td>0.196 <sup>+0.023</sup><sub>-0.019</sub></td></tr><tr><td>Argument of pericentre, ω</td><td>[°]</td><td>25.4 ± 5.6</td></tr><tr><td>Inclination, i</td><td>[°]</td><td>108 <sup>+11</sup><sub>-14</sub></td></tr><tr><td>Longitude of the ascending node, Ω</td><td>[°]</td><td>188.4 <sup>+4.5</sup><sub>-2.9</sub></td></tr><tr><td>Mean anomaly, M<sub>0</sub></td><td>[°]</td><td>259.0 ± 3.8</td></tr><tr><td>Impact parameter, b</td><td></td><td>13.3 <sup>+7.5</sup><sub>-9.6</sub></td></tr><tr><td>T'<sub>0</sub> - 2 450 000</td><td>[BJD<sub>TDB</sub>]</td><td>7971.56 ± 0.31</td></tr><tr><td>P'</td><td>[d]</td><td>34.5441 ± 0.0026</td></tr><tr><td>K'</td><td>[m s<sup>-1</sup>]</td><td>25.0 <sup>+2.5</sup><sub>-2.2</sub></td></tr><tr><td>Mass, M<sub>p</sub></td><td>(M<sub>J</sub>)</td><td>0.411 <sup>+0.054</sup><sub>-0.044</sub></td></tr><tr><td><strong>Planet (d)</strong></td><td></td><td></td></tr><tr><td>Semi-major axis, a</td><td>[au]</td><td>0.551 <sup>+0.010</sup><sub>-0.012</sub></td></tr><tr><td>Eccentricity, e</td><td></td><td>0.19 <sup>+0.19</sup><sub>-0.13</sub></td></tr><tr><td>Argument of pericentre, ω</td><td>[°]</td><td>148 ± 96</td></tr><tr><td>Inclination, i</td><td>[°]</td><td>82 <sup>+39</sup><sub>-33</sub></td></tr><tr><td>Longitude of the ascending node, Ω</td><td>[°]</td><td>199 <sup>+110</sup><sub>-150</sub></td></tr><tr><td>Mean anomaly, M<sub>0</sub></td><td>[°]</td><td>134 <sup>+110</sup><sub>-85</sub></td></tr><tr><td>Impact parameter, b</td><td></td><td>48 <sup>+40</sup><sub>-34</sub></td></tr><tr><td>T'<sub>0</sub> - 2 450 000</td><td>[BJD<sub>TDB</sub>]</td><td>7891 <sup>+20</sup><sub>-11</sub></td></tr><tr><td>P'</td><td>[d]</td><td>151.2 <sup>+2.1</sup><sub>-1.7</sub></td></tr><tr><td>K'</td><td>[m s<sup>-1</sup>]</td><td>8.9 ± 2.2</td></tr><tr><td>Mass, M<sub>p</sub></td><td>(M<sub>J</sub>)</td><td>0.262 <sup>+0.10</sup><sub>-0.076</sub></td></tr><tr><td>Mutual inclination b, c, I<sub>b,c</sub></td><td>[°]</td><td>21 ± 12</td></tr><tr><td>Mutual inclination b, d, I<sub>b,d</sub></td><td>[°]</td><td>93 ± 46</td></tr><tr><td>Mutual inclination c, d, I<sub>c,d</sub></td><td>[°]</td><td>95 ± 43</td></tr></tbody></table>

**Notes.** The table lists: Posterior median and 68.3% CI for the photodynamical modelling with three planets (Section 7.8). The orbital parameters are given for the reference time $t_{ref} = 2 457 957.48167$ BJD<sub>TDB</sub>.
<figure><img src="image_22.png" /><figcaption>**Fig. A.6.** Integrations over 1 kyr since the RISE transit from 1000 stable samples of the photodynamical modelling posterior. The 68.3% and 95.4% Bayesian CIs are plotted in different intensities. The solid colour curve marks the median of the posterior distribution. The solid grey curves correspond to the simulation based on the stable MAP values. The Lomb-Scargle periodogram of the median curve is shown, and the position of the most prominent peak is annotated.</figcaption></figure>
<figure><img src="image_23.png" /><figcaption>Fig. A.7. Idem. as Figure A.6, but for correlations of selected parameters ($\omega = \Omega + \omega$ is the longitude of the periapsis).</figcaption></figure>