Is It Possible to Predict a Driver's Distortion from Thiele Small Parameters?

Title says it all. I'm shopping for subwoofer drivers, and I would like to have some idea how much harmonic distortion they will produce.
Is there any way to do this mathematically? Given the advanced state of speaker simulation software, I would hope there is.
Any thoughts and ideas would be appreciated!
 
There are other parameters of the loudspeaker electromechanical model in addition to the Thiele-Small parameters. Like Xmax, Bl variation over the gap lenth and beyond, Rdc variation as a function of temperature which is the function of instantenous power, spider and surround nonlinearity as the function of excursion (simplified in a single Xmax parameter), cone breakup vs. frequency vs. amplitude, and there should be more, I believe. Also the same T-S parameter triplet can be achived by many different construction. For example ferrite and AlNiCo magnet could yield same parameters, yet they sound different.
Also just like for amplifiers, harmonic distortion does not characterize the quality of a loudspeaker exclusively.