One of the best ways to understand hadronization in QCD is to study the production of quarkonium. The color evaporation model (CEM) and Nonrelativistic QCD (NRQCD) can describe production yields rather well but spin-related measurements like the polarization are stronger tests. In this talk, I will outline the quarkonium polarization puzzle and present recent attempts to use the color evaporation model to describe the polarization of quarkonium production.