2023MNRAS.525..607F

With its near-to-mid-infrared high-contrast imaging capabilities, JWST is ushering us into a golden age of directly imaging Jupiter-like planets. As the two closest cold Jupiters, ε Ind A b and ε Eridani b have sufficiently wide orbits and adequate infrared emissions to be detected by JWST. To detect more Jupiter-like planets for direct imaging, we develop a GOST-based method to analyse radial velocity data and multiple Gaia data releases simultaneously. Without approximating instantaneous astrometry by catalogue astrometry, this approach enables the use of multiple Gaia data releases for detection of both short-period and long-period planets. We determine a mass of $2.96_{-0.38}^{+0.41}$ M_Jup and a period of $42.92_{-4.09}^{+6.38}$ yr for ε Ind A b. We also find a mass of $0.76_{-0.11}^{+0.14}$ M_Jup , a period of $7.36_{-0.05}^{+0.04}$ yr, and an eccentricity of 0.26$_{-0.04}^{+0.04}$ M_Jup, for ε Eridani b. The eccentricity differs from that given by some previous solutions, probably due to the sensitivity of orbital eccentricity to noise modelling. Our work refines the constraints on orbits and masses of the two nearest Jupiters and demonstrate the feasibility of using multiple Gaia data releases to constrain Jupiter-like planets.