A review of compressibility and rarefaction effects on spherical particle drag was conducted based on existing experimental data, theoretical limits, and direct simulation Monte Carlo method results. The data indicated a nexus point with respect to effects of Mach number and Knudsen number. In particular, it was found that a single drag coefficient (of about 1.63) is obtained for all particle conditions when the particle Reynolds number is about 45, and is independent of compressibility or rarefaction effects. At lower Reynolds numbers, the drag is dominated by rarefaction, and at higher Reynolds numbers, it is dominated by compressibility. The nexus, therefore, allows construction of two separate models for these two regimes. The compression-dominated regime is obtained using a modification of the Clift-Gauvin model to specifically incorporate Mach number effects. The resulting model was based on a wide range of experimental data and showed superior prediction robustness compared with previous models. For the rarefaction-dominated regime, the present model was constructed to directly integrate the theoretical creeping flow limits, including the incompressible continuum flow limit (Stokes drag), the incompressible weak rarefaction limit (Basset-Knudsen correction), and the incompressible free-molecular flow limit (Epstein theory). Empirical correlations are used to extend this model to finite particle Reynolds numbers within the rarefaction-dominated regime.