Abstract: We revisit the cosmological viability of the Hu-Sawicki modified gravity scenario. The impact of such a modification on the different cosmological observables, including gravitational waves, is carefully described. The most recent cosmological data, as well as constraints on the relationship between the clustering parameter sigma(8) and the current matter mass-energy density Omega(m) from cluster number counts and weak lensing tomography, are considered in our numerical calculations. The strongest bound we find is vertical bar f(R0)vertical bar < 3.7 x 10(-6) at 95% C.L. Forthcoming cluster surveys covering 10 000 deg(2) in the sky, with galaxy surface densities of O(10) arcmin(-2) could improve the precision in the sigma(8)-Omega(m) relationship, tightening the above constraint.