Abstract: We present 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, beta-decay, and neutrinoless double-beta-decay (0 nu beta beta) experiments. We directly constrain the full 3 + 1 active-sterile mixing matrix elements vertical bar U-alpha 4 vertical bar(2) , with alpha is an element of (e,mu,tau), and the mass-squared splitting Delta m(41)(2) (math) m(4)(2) – m(1)(2). We find that results for a 3 + 1 case differ from previously studied 1 + 1 scenarios where the sterile is coupled to only 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 datasets. The exact results are slightly prior dependent, but we reliably find all matrix elements to be constrained below vertical bar U-alpha 4 vertical bar(2) less than or similar to 10(-3) . Short-baseline neutrino oscillation hints in favor of eV-scale sterile neutrinos arc 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(beta) or m(beta)(beta), the effective mass parameters probed by beta-decay and 0 nu beta beta searches, respectively. When allowing for mixing with a light sterile neutrino, cosmology leads to upper bounds of m(beta) < 0.09 eV and m(beta)(beta )< 0.07 eV at 95% CL, more stringent than the limits from current laboratory experiments.