Exoplanet mass estimation for a sample of targets for the Ariel mission

Ariel’s ambitious goal to survey a quarter of known exoplanets will transform our knowledge of planetary atmospheres. Masses measured directly with the radial velocity technique are essential for well determined planetary bulk properties. Radial velocity masses will provide important checks of masses derived from atmospheric fits or alternatively can be treated as a fixed input parameter to reduce possible degeneracies in atmospheric retrievals. We quantify the impact of stellar activity on planet mass recovery for the Ariel mission sample using Sun-like spot models scaled for active stars combined with other noise sources. Planets with necessarily well-determined ephemerides will be selected for characterisation with Ariel. With this prior requirement, we simulate the derived planet mass precision as a function of the number of observations for a prospective sample of Ariel targets. We find that quadrature sampling can significantly reduce the time commitment required for follow-up RVs, and is most effective when the planetary RV signature is larger than the RV noise. For a typical radial velocity instrument operating on a 4 m class telescope and achieving 1 m s−1 precision, between ~17% and ~ 37% of the time commitment is spent on the 7% of planets with mass Mp < 10 M⊕. In many low activity cases, the time required is limited by asteroseismic and photon noise. For low mass or faint systems, we can recover masses with the same precision up to ~3 times more quickly with an instrumental precision of ~10 cm s−1.

Why simple quadrature is just as good as Monte Carlo

We motive and calculate Newton–Cotes quadrature integration variance and compare it directly with Monte Carlo (MC) integration variance. We find an equivalence between deterministic quadrature sampling and random MC sampling by noting that MC random sampling is statistically indistinguishable from a method that uses deterministic sampling on a randomly shuffled (permuted) function. We use this statistical equivalence to regularize the form of permissible Bayesian quadrature integration priors such that they are guaranteed to be objectively comparable with MC. This leads to the proof that simple quadrature methods have expected variances that are less than or equal to their corresponding theoretical MC integration variances. Separately, using Bayesian probability theory, we find that the theoretical standard deviations of the unbiased errors of simple Newton–Cotes composite quadrature integrations improve over their worst case errors by an extra dimension independent factor {\propto N^{-\frac{1}{2}}}. This dimension independent factor is validated in our simulations.

New Sensitivity Indices of a 2D Flood Inundation Model Using Gauss Quadrature Sampling

A new method for sensitivity analysis of water depths is presented based on a two-dimensional hydraulic model as a convenient and cost-effective alternative to Monte Carlo simulations. The method involves perturbation of the probability distribution of input variables. A relative sensitivity index is calculated for each variable, using the Gauss quadrature sampling, thus limiting the number of runs of the hydraulic model. The variable-related highest variation of the expected water depths is considered to be the most influential. The proposed method proved particularly efficient, requiring less information to describe model inputs and fewer model executions to calculate the sensitivity index. It was tested over a 45 km long reach of the Richelieu River, Canada. A 2D hydraulic model was used to solve the shallow water equations (SWE). Three input variables were considered: Flow rate, Manning’s coefficient, and topography of a shoal within the considered reach. Four flow scenarios were simulated with discharge rates of 759, 824, 936, and 1113 m 3 / s . The results show that the predicted water depths were most sensitive to the topography of the shoal, whereas the sensitivity indices of Manning’s coefficient and the flow rate were comparatively lower. These results are important for making better hydraulic models, taking into account the sensitivity analysis.