We provide a consistent framework to set limits on properties of light sterile neutrinos coupled to all three active neutrinos using a combination of the latest cosmological data and terrestrial measurements from oscillations, β-decay and neutrinoless double-β decay (0νββ) experiments. We directly constrain the full 3 + 1 active-sterile mixing matrix elements |U α4 | 2 , with α ∈ (e, µ, τ ), and the mass-squared splitting ∆m 2 41 ≡ m 2 4 − m 2 1 . We find that results for a 3 + 1 case differ from previously studied 1 + 1 scenarios where the sterile is only coupled to one of the neutrinos, which is largely explained by parameter space volume effects. Limits on the mass splitting and the mixing matrix elements are currently dominated by the cosmological data sets. The exact results are slightly prior dependent, but we reliably find all matrix elements to be constrained below |U α4 | 2 10 −3 .Short-baseline neutrino oscillation hints in favor of eV-scale sterile neutrinos are in serious tension with these bounds, irrespective of prior assumptions. We also translate the bounds from the cosmological analysis into constraints on the parameters probed by laboratory searches, such as m β or m ββ , the effective mass parameters probed by β-decay and 0νββ searches, respectively. When allowing for mixing with a light sterile neutrino, cosmology leads to upper bounds of m β < 0.09 eV and m ββ < 0.07 eV at 95% C.L, more stringent than the limits from current laboratory experiments.