Starting from the joint statistics of the orthogonal components of the broadband thermal light viberation, which are normal to the propagation direction of light wave, we derive the concept of polarization degree and the joint probability density functions of those two orthogonal components in various polarization states, so that an overall description of polarization of the broadband thermal light is obtained.