Atmospheric composition and circulation of the ultra-hot Jupiter WASP-121b with joint NIRPS, HARPS and CRIRES+ transit spectroscopy

Ultra-hot gas giants such as WASP-121b provide unique laboratories for exploring atmospheric chemistry and dynamics under extreme irradiation conditions. Uncovering their chemical composition and atmospheric circulation is critical for tracing planet formation pathways. Here, we present a comprehensive atmospheric characterisation of WASP-121b using high-resolution transit spectroscopy across the optical to infrared with HARPS, NIRPS, and CRIRES+ spanning nine transit events. These observations are complemented with five TESS photometric sectors, two EulerCam light curves simultaneous to the HARPS and NIRPS transits, and an extensive radial velocity dataset in order to refine WASP-121b's orbital parameters. A cross-correlation analysis detected iron (Fe), carbon monoxide (CO) and vanadium (V) absorption signals with SNR of 5.8, 5.0, and 4.7, respectively. Our retrieval analysis constrains the water (H2O) abundance to −6.52−0.68+0.49 dex, although its absorption signal is effectively muted by the hydride (H−) continuum. We constrained the relative abundances of the volatile and refractory elements - which represents a crucial diagnostic of atmospheric chemistry, evolution, and planet formation pathways. The retrieved abundance ratios are broadly consistent with expected values of a solar composition atmosphere in chemical equilibrium, likely indicating minimal disequilibrium chemistry alterations at the probed pressures (∼10−4−10−3 bar). We update the orbital parameters of WASP-121b with its largest radial velocity dataset to date. By comparing orbital velocities derived from both the radial velocity analysis and the atmospheric retrieval, we determined a non-zero velocity offset caused by atmospheric circulation, ΔKp = −15 ± 3 km s−1 (assuming M⋆ = 1.38 ± 0.02 M⊙), consistent with predictions from either drag-free or weak-drag 3D global circulation models, while we caution the non-negligible dependence on the assumed stellar mass. These results place new constraints on the thermal structure, dynamics, and chemical inventory of WASP-121b, highlighting the power of multi-wavelength high-resolution spectroscopy to probe exoplanetary atmospheres.