We study finite-temperature properties of strongly correlated fermions in two-dimensional optical lattices by means of numerical linked cluster expansions, a computational technique that allows one to obtain exact results in the thermodynamic limit. We focus our analysis on the strongly interacting regime, where the on-site repulsion is of the order of or greater than the band width. We compute the equation of state, double occupancy, entropy, uniform susceptibility, and spin correlations for temperatures that are similar to or below the ones achieved in current optical lattice experiments. We provide a quantitative analysis of adiabatic cooling of trapped fermions in two dimensions, by means of both flattening the trapping potential and increasing the interaction strength.