Polarized reflectance (Rp) and degree of linear polarization (DOLP) provide essential information about polarized characteristics of land surfaces. For a given target, DOLP determines the magnitude of Rp. It has been proved that DOLP can be used for some remote monitoring cases that cannot be well detected with either non-polarized or polarized reflectance. Several bidirectional polarization distribution function (BPDF) models have been proposed in the last several decades to reproduce the angular distribution of Rp, but much less attention has been devoted to modeling and analyzing of DOLP. In this study, the Nadal-Bré on BPDF model was transferred for calculating the DOLP of earth targets, and characteristics of DOLP were analyzed based on the modeling results. To evaluate the model's feasibility, two experiments were executed: a fitting and a a priori modeling. The results showed good correlations (r>0.9) between estimated and measured DOLP when the model was fitted with POLDER/PARASOL (a space-borne multiangle multi-spectral polarimetric sensor) measurements. An increase of accuracy from 490 nm to 865 nm for fitting modeling was achieved and the highest accuracy was found at 865 nm for both experiments, with overall relative root mean square errors of 1.1 and 1.3 for fitting and a priori modeling, respectively. Classbased free parameters can be used for the a priori model of DOLP. The dispersion of the target-based free parameters controls the correlation of the a priori modeling results. Moreover, the maximum DOLP was found to be strongly determined by the corresponding bidirectional reflectance factor for every surface type (R 2 =0.86). This study provides an additional approach for obtaining DOLP from remote sensing platform and is helpful for studies of typical land surfaces.