On the synergy between Ariel and ground-based high-resolution spectroscopy

Since the first discovery of an extra-solar planet around a main-sequence star, in 1995, the number of detected exoplanets has increased enormously. Over the past two decades, observational instruments (both onboard and on ground-based facilities) have revealed an astonishing diversity in planetary physical features (i. e. mass and radius), and orbital parameters (e.g. period, semi-major axis, inclination). Exoplanetary atmospheres provide direct clues to understand the origin of these differences through their observable spectral imprints. In the near future, upcoming ground and space-based telescopes will shift the focus of exoplanetary science from an era of “species discovery” to one of “atmospheric characterization”. In this context, the Atmospheric Remote-sensing Infrared Exoplanet Large (Ariel) survey, will play a key role. As it is designed to observe and characterize a large and diverse sample of exoplanets, Ariel will provide constraints on a wide gamut of atmospheric properties allowing us to extract much more information than has been possible so far (e.g. insights into the planetary formation and evolution processes). The low resolution spectra obtained with Ariel will probe layers different from those observed by ground-based high resolution spectroscopy, therefore the synergy between these two techniques offers a unique opportunity to understanding the physics of planetary atmospheres. In this paper, we set the basis for building up a framework to effectively utilise, at near-infrared wavelengths, high-resolution datasets (analyzed via the cross-correlation technique) with spectral retrieval analyses based on Ariel low-resolution spectroscopy. We show preliminary results, using a benchmark object, namely HD 209458 b, addressing the possibility of providing improved constraints on the temperature structure and molecular/atomic abundances.

from the telluric and stellar photosphere contamination, which are nearly stationary signals during one observing night. However, sometimes strong telluric residuals can remain in the data biasing the interpretation of the atmospheric composition and properties (e.g., [12]). On the contrary, LRS observations, as they are mainly gathered with space telescopes, are not affected by telluric lines.
HRS from the ground and LRS from space are complementary, but they are difficult to combine. For instance, as the HRS data are self-calibrated (i.e. by fitting the flux in each spectra channel as a function of time and dividing by the fit), any variation due to the planetary atmosphere is being measured with respect to a local 'pseudo-continuum', and the absolute level of the planetary absorption is lost. The lack of a robust continuum can thus introduce degeneracies, because model templates with different abundances and T/P profiles may look similar [7]. Consequently, the removal of the continuum makes the typical retrieval framework for the analysis of LRS data unusable, and, therefore, putting strong constraints on both the molecular abundances and the temperature/pressure structure is harder. Furthermore, another trouble in using retrieval algorithms with HRS data is converting the cross-correlation values into a goodness-of-fit estimator. The first attempt at combining LRS and HRS was performed by [11], who employed a new retrieval algorithm able to compute the joint probability distribution of LRS and HSR data. Their results showed that the combination of these two techniques improves constraints on both the vertical thermal structure and the retrieved molecular abundances of a planetary atmosphere. However, this algorithm requires significant computational power, thus its application is limited to the evaluation of a few thousand model HRS spectra sampled from an LRS posterior. More recently, [9] introduced a new approach to solving these restrictions. They built a new robust and unbiased framework to perform Bayesian retrieval analyses on HRS data. This framework is based on the cross-correlation between models to extract the planetary spectral signal. In this way, it permits to combine and explore the synergies between HRS and LRS. However, [9] applied this framework to a narrow spectral coverage -i.e. the VLT CRIRES K-band -while, as current spectrographs like GIANO-B, CARMENES, SPIROU, and upcoming such as NIRPS and CRIRES+, have broader wavelength ranges, they would require major computational efforts.
In this work, we present a case study of a representative hot Jupiter, namely HD 209458 b, with the objective of analyzing both HRS data acquired with the spectrograph GIANO-B mounted on the Telescopio Nazionale Galileo, TNG, (Section 2), and LRS data gathered with the HST/WFC3 instrument (Section 3). More precisely, we carried out the analysis at low resolution by taking as reference the results obtained from our simulation at high resolution. After comparing the results obtained with the two different methods, in Section 4 we perform a simulation to test the Atmospheric Remote-sensing Infrared Exoplanet Large-survey (hereafter, Ariel, [52]) capability to probe exoplanetary atmospheres. Finally, we conclude (Section 5) by highlighting the importance of a future synergy between space-borne low-resolution telescope, like Ariel, and ground-based high-resolution facilities.

Observations and analysis of GIANO-B high-resolution spectra
The HRS dataset has been analysed in [24]. It encompassed 4 transit events gathered with the TNG telescope in the GIARPS observing mode [16], that allows for a simultaneous acquisition of high-resolution spectra in the optical (0.39-0.69 μm) and in the nIR (0.95-2.45 μm) with the HARPS-N (resolving power R∼115,000), and the GIANO-B (resolving power R∼50,000) spectrographs, respectively. We observed the system HD 209458 as part of the Large Program "GAPS2: the origin of planetary systems diversity" [PI: G. Micela]. For the work presented here, we focused on the nIR observations collected with GIANO-B. This nIR spectrograph acquires images with the nodding acquisition mode ABAB to enable an optimal subtraction of the background and detector noise. The detailed description of the analysis we carried out on these spectroscopic observations of HD 209458 b is described in [24]. In this manuscript, for completeness, we summarize the most important steps we performed to extract the planetary signal from the GIANO-B spectra. After the extraction and the wavelength calibration of the 50 GIANO-B orders with the GOFIO instrument pipeline [47], we proceeded to eliminate the stellar and telluric contamination. To do it, we took advantage of the fact that the planetary signal is Doppler shifted during one night of observation -the planet is moving around its host star-, whereas the stellar and telluric spectra are stationary (or quasi-stationary) signals in wavelength. In this work, we employed the Principal Component Analysis (PCA) and linear regression techniques to model and remove the telluric and stellar signals. (For more details on this new analysis framework see [24]). Successively, we looked for HD 209458 b's atmospheric signatures by cross-correlating our GIANO-B spectra (now free from the telluric and stellar contamination) with model templates. The theoretical spectra have been generated with the GENESIS model [22]. Our models spanned a wide range of pressures (10 2 -10 −8 bar) and wavelengths (0.9-2.6 μm). Collision-induced absorption from H 2 -H 2 and H 2 -He was included. An isothermal atmosphere and constant Volume Mixing Ratios (VMR) were assumed within the range 1,000<T< 1,500 K and 10 −5 <VMR< 10 −2 . We utilised the ExoMol database for H 2 O, NH 3 , HCN, and C 2 H 2 [4,14,17,46], the HITEMP database for CH 4 and CO [25,37], and the Ames database for CO 2 [28]. For each theoretical model, to maximize the planetary signal, we chose to perform the cross-correlation technique on a subset of GIANO-B orders, that is those that did not exhibit strong telluric residuals and contained the strongest spectral lines of the planet spectrum. The cross-correlation functions of these selected orders were then co-added in time as a function of the planetary radial velocity in the Earth's rest frame and maximum radial velocity semi-amplitude (K p ). For all species, we calculated the detection significance by performing a Welch t-test [61] , see e.g.
[10], on the cross-correlation values by assuming as null hypothesis that out-of-trail (far from the planet radial velocity) and in-trail (around the planet radial velocity) values have the same mean. We assumed a confidence level limit of 3σ . Table 1 lists the molecules we detected (H 2 O, NH 3 , HCN, C 2 H 2 , CH 4 , and CO) and the corresponding abundances used for the cross-correlation between atmospheric templates and our GIANO-B spectra. We need to underline that these abundances do not match From left to right, we report the investigated molecules, the theoretical abundances employed for the isothermal models, and the significance. These abundances do not match any specific chemical scenarios, but they were used to maximise the detection significance any specific chemical scenarios, but they were used to maximise the detection significance. We did not find any evidence of CO 2 . Thus, we discarded this molecule from the rest of the analysis carried out in this paper.
To better characterize the atmosphere of HD 209458 b, in the work presented in [24], we computed two other sets of non-isothermal atmospheric models. For the first set of models, we used input temperature-pressure abundance (T-p-VMR) profiles calculated under the assumptions of a cloud-free atmosphere in chemical and radiative equilibrium. The second set of models accounts for the presence of clouds/aerosol by adding a grey cloud deck with a top-deck pressure of 10 −5.5 bar and a cloud fraction of 0.4 [5]. Following the receipt proposed by [9], we thus converted the cross-correlation values into likelihood values, and we used the likelihood-ratio test to compare the different models. Our results statistically favour the presence of aerosols in the atmosphere of HD 209458 b which dampen the amplitude of the molecular lines but do not evidently hamper their detection [23,26]. Furthermore, the atmospheric models we tested in thermochemical equilibrium seem to prefer a carbon-to-oxygen ratio close to or greater than 1, higher than the solar value (0.55).

Observations and analysis of HST/WFClow-resolution data
The LRS observations employed in this work have been acquired with the HST/WFC3/G141 grism as part of the 12181 program, 'The Atmospheric Structure of Giant Hot Exoplanets' [PI: D. Deming]. We downloaded publicly available observations of HD 209458 b from the Mikulski Archive for Space Telescopes (MAST). We used the publicy available Python package Iraclis [56] to analyze the raw HST spatially scanned spectroscopic images and to extract the planetary spectrum. We then carried out the modeling of the extracted spectrum using the publicly available spectral retrieval algorithm TauREx3 [3,59,60].
Step (i) consists in bias-level and zero-read corrections; non-linearity correction; dark current subtraction; gain variations correction; sky background subtraction; calibration; flat-field correction; bad pixels and cosmic ray correction [54][55][56]. After these initial operations, the stellar flux was extracted from the raw images to create the wavelength-dependent light curves. Two types of light curves were extracted, a broad wavelength band white light curve covering the entire G141 grism wavelength range (1.088 -1.68 μm, see Fig. 1) and spectral light curves with a resolving power of 70 at 1.4 μm (see Fig. 2). When extracting the spectral light curves, Iraclis accounts for two "ramps", i.e. time-dependent systematics introduced by the WFC3/IR detector. The first ramp has a linear behavior and affects each HST visit, while the second one alters each orbit and has an exponential behavior. To correct for these effects, the white light curve is fitted by Iraclis with a model for the systematics simultaneously with the transit model. Since, in this data set, the long-term ramp can be approximated by a linear function only after the third orbit, following [55], we decided to discard the first two orbits. We fitted for T 0 (the mid-transit time) and R P /R (the ratio of planet to stellar radius) as free parameters, we fixed the orbital period (P), the inclination (i) and a/R to the values of Table 2, and assumed a circular orbit given that the eccentricity is consistent with zero (e.g., [8]). We modeled the stellar limb darkening effect by employing the non-linear formula with four terms  [55], we note additional residuals in the white light curve fit. The model fails to fit the egress. In this previous work, they attributed this behavior to either non-optimal values used for i and a/R or remaining systematics by [15]. The coefficients are calculated by fitting the stellar profile from an ATLAS model [27,35] and by using the stellar parameters presented in Table 2. We then fitted the spectral light curves by using the dividing method proposed by [32]. This technique considers the white light curve as a comparison source, indeed each spectral light curve is fitted with a model that includes the white light curve and its best-fit model. This has as consequence that the residuals from fitting one of the spectral light curves do not show trends similar to those in the white light curve. In this fitting procedure, the only free parameter is R P /R , while the other parameters are the same employed for the white light curve fitting. The transmission spectrum is then constructed from the spectral light curves by determining the planet-star radii ratio as a function of wavelength (see Table 3).
Successively, we used the Bayesian atmospheric retrieval framework TauREx3 [3,59,60] to fit the WFC3 spectrum. Thanks to the implementation of the nested sampling code Multinest [20], TauREx allows us to explore the parameter space and find the best fit for the transmission spectrum extracted with Iraclis. In our retrieval analysis, we used 1000 live points and an evidence tolerance of 0.5. We used the plane-parallel approximation to model the atmospheres, with pressures ranging from 10 −2 to 10 6 Pa sampled uniformly in log-space by 100 atmospheric layers. The atmosphere is simulated by assuming an isothermal T/P profile (T=T eq =1484 K, see Table 2) and constant molecular abundances as a function of altitude. These assumptions are reasonable because, due to the limited wavelength range of the HST/WFC3/G141 grism, we are probing a narrow range of the planetary atmosphere [56]. Considering the results obtained at HRS with GIANO-B (see Section 2), we decided to consider in our fit the following active-gases: H 2 O [46], CH 4 [62], CO [37], NH 3 [17], C 2 H 2 [14], HCN [4]. Each molecular abundance was allowed to vary between 10 −12 and 10 −2 in volume mixing ratios as a log uniform prior. We assumed the bulk composition of the atmosphere, similar to that of Jupiter, i.e. made up of a mixture of 85% hydrogen and 15% helium with a ratio H 2 /He=0.17647, 85% H 2 and 15% He -the Jupiter Chemical composition. We used absorption crosssections at 15000. Furthermore, we assumed uniform priors for the 10 bar radius (R = 1.3−1.4 R Jup ), and the temperature (T = 1000−1800 K). Rayleigh scattering and collision-induced absorption of H 2 -H 2 and H 2 -He [1,2,21] were also included. We tested two different scenarios by assuming a cloudy atmosphere (model [1] in Table 4) and a cloud-free atmosphere (model [2] in Table 4). Clouds were fitted assuming a grey opacity model and cloud top pressures -i.e., the pressure at which the cloud starts to be opaque-bounds were set between 10 −2 and 10 6 Pa. All the priors we assumed are listed in Table 4. In our analysis, the fitted parameters were the molecular abundances, the temperature, the mean molecular weight, the radius at 10 bar, and, in the cloudy scenario, the cloud top pressure. In agreement with [56], we quantified the significance of our detections with the Atmospheric Detection Index (ADI) a positively defined Bayes Factor between the nominal atmospheric model and a flat-line model (i.e. a model representing a fully cloudy atmosphere, which contains no active trace gases, Rayleigh scattering or collision-induced absorption). This value was then translated into a statistical significance [30] by using Table 2 of [6]. Our TauREx3 retrieval results are listed in Table 4. Figure 3 shows the best-fit models and the contribution plots for the two scenarios tested. The posterior distributions are shown in Fig. 4. In both cases, we retrieved a significant atmosphere around HD 209458 b with an ADI of 22.1-17.4 for the cloudy and clear model, respectively. Accordingly to the Bayesian evidences the cloudy scenario seems to be strongly, but not decisively, favorite. For both scenarios the retrieved temperature is lower than the predicted equilibrium temperature. This could be explained by the fact that we are probing the atmosphere in the terminator area, and we modeled the atmosphere in 1D using an isothermal profile (e.g. [13,40]). This is in agreement with what was highlighted by [49], who found evidence of a global trend between the equilibrium and the retrieved temperatures, with the latter almost always showing lower values. The retrieved radii are compatible with the theoretical value (1.359±0.019 R J ) within ∼ 1σ . In the cloudy-model, we note a correlation between the H 2 O abundance, the radius, and the cloud pressure. For less H 2 O, the model requires deeper clouds and a Table 4 List of the retrieved parameters, their uniform prior bounds, the scaling used and the retrieved value for both the two scenarios tested in this paper (model [1] and model [2]) Model  In the last two columns the ADI index and the log of the Bayesian Evidence are also reported  The cloudy and the cloud-free scenarios are represented in green and in pink, respectively. For consistency, we do not plot the clouds pressure posterior distribution by TauREx for the best fit solution. The signal is theoretically blocked by this layer and nothing can be observed at higher pressures. Molecules found below this line are unconstrained. On the contrary, HRS is able to detect molecular absorbers even if there is a high altitude clouds deck, given that it is most sensitive to the spectral lines cores that are formed above the clouds at lower pressure (see [23,26]).

Ariel spectra simulation
Upcoming observatories in space and on the ground, thanks to their broader spectral coverage and higher signal-to-noise ratio, will enable the analysis of a large number of planets. In this contest, during its 4-yr mission, the ESA/Ariel mission, scheduled to launch in 2028, will allow for the atmospheric characterisation of a statistically significant population of exoplanets (∼1000) -ranging from Jupiters and Neptunes down to super-Earths. It will thus permit to constrain the abundances of atmospheric constituents which are strictly linked to the planetary formation/evolution environment.
In this section, we describe a simulation we performed to explore the potential of using Ariel with the approach described above. As in the rest of the paper our benchmark object is HD 209458 b. Firstly, we simulated an high-resolution transmission spectrum of HD 209458 b by using the TauREx3 algorithm in forward mode. We simulated the atmosphere of HD 209458 b composed of H 2 and He with a ratio H 2 /He=0.17647. We modeled the atmosphere with pressures ranging from 10 −2 to 10 6 Pa, uniformly sampled in log-space with 100 atmospheric layers. We assumed an isothermal T/P profile, with T=1080 K, i.e. the maximum of the temperature posterior distribution obtained from the retrieval of HST/WFC3 observations 2 (see Fig. 3), and constant chemistry profile. The planetary mass and stellar parameters (i.e. radius and effective temperature) were fixed to the values reported in Table 2, whilst the 10 bar radius was assumed to be equal to that found in the retrieval analysis of real HST/WFC3 observation, i.e. R P =1.331 R jup (see Table 4). We included Rayleigh scattering and collision-induced absorption (CIA) of H 2 -H 2 and H 2 -He. We assumed a grey opaque cloud deck at P clouds ∼ 118 Pa (value obtained from the retrieval analysis of real HST observations -see Table 4). The trace gases we considered were the same included in the HST retrieval analysis (i.e. H 2 O, CO, CH 4 , NH 3 , C 2 H 2 , and HCN). We performed two different simulations: ). (b) We then tried a second experiment -a more unrealistic example-where we imposed the tracers abundances to be equal to those values that maximise the cross-correlation analysis at HRS (see Table 1). We must emphasize that this simulation is only an exercise, it does not claim to reproduce a real observed HST spectrum. We are in fact using the abundances reported in Table 1 which, as previously pointed out, do not correspond to any specific chemico-physical scenario of the atmosphere. The aim of this experiment is to assess Ariel's improvement in abundance retrievals with respect to HST/WFC3 in case of higher abundances.
Next, we binned the high-resolution transmission spectrum to the resolutions of the Ariel spectrometers (i.e. NIRSpec, 1.1-1.95 μm at a resolving power of R=20, AIRS-Ch0, 1.95-3.9 μm at R=100, and AIRS-Ch1, 3.9-7.8 μm at R=30), and we used the instrument noise simulator Ariel Radiometric Model (ArielRad) [45] to provide a realistic noise model. Finally, we performed atmospheric retrievals using TauREx3 in fitting mode. Table 5 lists the TauREx retrieval results we obtained, whilst the posterior distribution and the best-fit spectrum are plotted in turquoise in Fig. 5 (for the (a) scenario) and in Fig. 6 (for the (b) simulation). As a comparison, Figs. 5, and 6 show also the posterior distribution and the best-fit transmission spectrum we obtained for simulated HST data. More precisely, we binned our simulated atmospheric spectrum at the WFC3/HST resolution and we perturbed it with the noise obtained in the previous section for the real HST data (right middle box).
-Simulation (a) Figure 5 shows the retrieval results we obtained for both HST and Ariel observations for the (a) scenario. Only the water vapour posterior distribution is well constrained, whereas we are not able to constrain the CO abundance, and we can only put an upper limit of ∼ 10 −3 on the amount of the other molecules. The 10 bar radius correlates with the temperature distribution: the higher the temperature is, the smaller radius. The Ariel retrieval seems to prefer a slightly higher temperature -and thus a lower 10 bar radius-than the HST simulation.
The results obtained for the scenario (b) are instead shown in Fig. 6. Now that the simulated molecular abundances are higher than before, we can appreciate some differences between Ariel and HST/WFC3. As in the previous simulation, we are not able to constrain the CO abundance distribution. CH 4 , C 2 H 2 , NH 3 and HCN are better constrained in the Ariel simulation.
These simulations highlight several aspects, which are important to discuss. We can appreciate the improvement in putting constraints on the chemical abundances that we can obtain with Ariel. Indeed, while there is still considerable degeneracy for retrievals from HST/WFC3 spectra, Ariel simulations in Fig. 6 show constrained posterior distributions for most of the active gases. This aspect has been also pointed out by [52]. Our simulation highlights how the three Ariel spectrometers represent a step forward compared to the current HST/WFC3. On one side, thanks to its wide spectral coverage, Ariel will permit the detection of several molecular species that do not have strong absorption bands in the WFC3 wavelength range. On the other side, Ariel observations will probe a much wider pressure range than WFC3, from approximately 1 to 10 3 mbar [52].
However, when dealing with low molecular abundances, such as those tested in simulation (a), our experiments revealed that also Ariel has some limitations to well constrain the investigated trace-gas abundances and we have to use a combination of HRS results and Ariel data. Indeed, if we consider the Ariel results alone, we could not claim the presence of CO and we can only put an upper limit on the abundance of

Discussion
In the near future, Ariel will allow for detailed characterisation of thousands of exoplanetary atmospheres. In this perspective, exploiting the synergies between Ariel and current (e.g. GIANO-B, SPIROU, CARMENES) and upcoming ground-based high-resolution spectrographs (NIRPS, CRIRES+, and ELTs) represents the future of exoplanetary science. The combination of HRS and LRS can break the degeneracies that arise when we apply the two techniques individually. On one hand, the simultaneous spectra coverage of Ariel will permit us to determine important atmospheric properties as the molecular abundances, the clouds coverage, and the T/P profile. Moreover, it could provide a local pseudo-continuum to HRS observations. On the other hand, the contributions of multiple species are expected to overlap at the low resolution of Ariel, and high-resolution instruments like GIANO-B can help in understanding which atmospheric constituents are expected. Furthermore, in presence of a deep cloud deck, some spectral features that could be masked with LRS, may be accessible with HRS, since it probes higher atmospheric levels. It is thus clear that only a combination of low-and high-resolution data in a consistent retrieval framework can provide absolute molecular abundances, avoiding confusion between species. The power of having a joint retrieval of multi-spectral resolution data was shown by [9] who provided a simple framework for combining HRS data and LRS data within a unified likelihood function. However, these calculation were restricted to an extremely narrow wavelength range at high-resolution (2270-2350 nm), and limited to emission spectroscopy data. No such framework has ever been applied to transmission spectroscopy. An efficient LRS+HRS retrieval framework could allow to achieve unprecedented atmospheric characterization and thus to (i) better understand what are the main parameters that govern the chemistry and physics of exoplanet's atmospheres, and (ii) derive the C/O and O/H ratios which provide insights into formation and evolution mechanisms (e.g. [43]).
In this work, we took HD 209458 b as a reference, but the same study can be performed on several exoplanets. To understand the number of potential objects on which we can perform a similar analysis, it would be useful to know how many of the future Ariel targets could be observed with ground-based HRS. We downloaded Fig. 6 Posterior distributions for the simulated transmission spectrum of HD 209458 b with molecular abundances that maximise the cross-correlation analysis at HRS (simulation (b)) as observed by Ariel (in turquoise), with overplotted the posterior distribution obtained for simulated HST data (in orange). The 'true values' we assumed for each parameter are shown in red. Insets: simulated transmission spectra (black points) for Ariel (upper panel) and HST (bottom panel) with overplotted the best fit solution found by TauREx and the correspondent 1σ and 2σ error bars the list of transiting planets from the TEPCat Catalog [51], and for each of them we estimated the atmospheric signal to noise ratio in transmission (S/N) under the following assumptions: (i) giant planets with helium-hydrogen dominated atmospheres (R P ≥ 3.5 R ⊕ ), (ii) no clouds, (iii) and pure photon noise. We calculated the expected S/N as follows: where we account for the transit duration (T ), the stellar magnitude in K band (F = 10 − magK 2.5 ), the atmospheric scale height (H ), the planetary (R P ) and the stellar radius (R ). We divided each of the estimated S/N to the expected S/N of HD 209458 b.
In Section 2 we reported a detection of water vapour with a statistical confidence of 10 σ by considering 4 transits, therefore around 5 σ (σ H D209 = 10 σ/ √ 4.) for each transit. If we assume a statistical significance threshold of 4 σ for water detection and to gather 5 transit observations for each target, ∼25 targets in our sample of transiting planets could be studied with both Ariel and HRS instruments such as GIANO-B at 4-m class telescope. The number of suitable targets is limited by the 4-m aperture of the TNG, which restricts the observable sample to relatively bright stars. Indeed, most of the planet-hosting stars are remarkably fainter than HD209458b, and thus we are not able to observe them with a sufficiently high S/N. In the near future, as soon as high-resolution instrumentation at ELTs is available, this sample is expected to increase. Even though the analysis presented in this work is focused on the nIR band, we have to highlight also the benefits that can be obtained by combining HRS and Ariel data in the VIS. Studies of hot/ultra-hot Jupiters with the Hubble STIS/WFC3 (e.g. [19]) have recently shown indications of excess absorption possibly due to TiO, VO, and FeH or H-opacity. Ariel could measure these excesses but is not able to distinguish between the different opacities due to the spectral bins' width. Thus, HRS could be employed to place upper limits on the molecular abundances. The HRS prior information could then be incorporated into the retrieval of the Ariel data -in a similar way to that proposed in this study-to help break the degeneracy that will be seen. This approach has recently been taken for the VO molecule. Indeed, HST observations of WASP-121 b seemed to show evidence of VO [19], whilst high-resolution data only succeeded in putting an upper limit on this opacity [44]. Thus a future combination between optical high-resolution spectrographs (such as HARPS-N@TNG) and Ariel's VIS photometric channels could put constraints also on the TiO/VO/Na opacities.

Conclusions
In this paper, we compared results between HRS and LRS gathered in transmission for HD 209458 b. We applied different methods to analyse the two different datasets. On the one hand, we disentangle the planetary spectrum from the stellar and telluric contamination and the cross-correlation technique to extract the planet's signal from our GIANO-B data. On the other one, we used the Python package Iraclis to extract the transmission spectrum of HD 209458 b from the HST/WFC3 raw images and the algorithm TauRex3 to perform retrieval analysis. We noted a considerable degeneracy in the posteriors distributions of many molecular species obtained at low resolution.
Successively, we performed a simulation to test the ability of the Ariel space mission to give precise constraints on the atmospheric chemical abundances. By using TauRex3 in the forward mode, we simulated a transmission spectrum for HD 209458 b containing H 2 O, CO, CH 4 , NH 3 , C 2 H 2 ,and HCN. We added instrumental noise to the model by using the ArielRad simulator and we interpreted it with TauREx3 in retrieval mode. Our simulation showed how, unlike HST/WFC3, the retrieved atmospheric information we could obtain from Ariel is much better constrained, and how much the exoplanetary community could benefit from exploiting the synergy between Ariel, and current or upcoming HRS ground-based telescopes.
The work presented in this paper set only the basis to combine LRS and HRS. In this respect, we need to further stress that in this work we did not carry out a true combination of LRS and HRS results such as that advocated in recent literature (e.g., [9,11]) and based on a statistically supported method to combine datasets, but rather we selected the chemical absorbers to be investigated at low-resolutions based on HRS results, and for each molecule we used as starting value of the TauREx retrieval the abundance that maximises the detection significance at HRS. However, to perform a true combination of these two techniques, work still needs to be done. The natural starting point will be the effective combination of HST/WFC3 and GIANO-B spectra by maximizing a total likelihood which is the combination of a low-resolution likelihood [42] and an high-resolution likelihood [9], over much broader wavelength coverage than that used in the pioneering work of [9]. Such a framework will i) yield unprecedented constraints on both the atmospheric composition of the studied exoplanets and their formation/migration mechanisms and, even more importantly, ii) constitute a reference for combined analyses of future low-resolution data gathered with JWST and/or ARIEL and ground-based high-resolution spectra, including those that will be acquired with future ELTs.