Known since the ancient times, cotton continues to be one of the essential materials for the human civilization. Cotton fibers are almost pure cellulose and contain both crystalline and amorphous nanodomains with different physicochemical properties. While understanding of interactions between the individual cellulose chains within the crystalline phase is important from a perspective of mechanical properties, studies of the amorphous phase lead to characterization of the essential transport parameters, such as solvent diffusion, dyeing, drug release, and toxin absorption, as well as more complex processes of enzymatic degradation. Here, we describe the use of spin probe electron paramagnetic resonance methods to study local polarity and heterogeneous viscosity of two types of unprocessed cotton fibers, G. hirsutum and G. barbadense, harvested in the State of North Carolina, USA. These fibers were loaded with two small molecule nitroxide probes that differ in polarity-Tempo and its more hydrophilic derivative Tempol-using a series of polar and non-polar solvents. The electron paramagnetic resonance spectra of the nitroxide-loaded cotton fibers were analyzed both semi-empirically and by least-squares simulations using a rigorous stochastic theory of electron paramagnetic resonance spectra developed by Freed and coworkers. A software package and least-squares fitting protocols were developed to carry out automatic simulations of multi-component electron paramagnetic resonance spectra in both first-derivative and the absorption forms at multiple resonance frequencies such as X-band (9.5 GHz) and W-band (94.3 GHz). The results are compared with the preceding electron paramagnetic resonance spin probe studies of a commercial bleached cotton sheeting carried out by Batchelor and coworkers. One of the results of this study is a demonstration of a co-existence of cellulose nanodomains with different physicochemical properties such as polarity and microviscosity that are affected by solvents and temperature. Spin labeling studies also revealed a macroscopic heterogeneity in the domain distribution along the cotton fibers and a critical role the cuticular layer is playing as a barrier for spin probe penetration. Finally but not lastly, the simultaneous multi-component least-squares simulation method of electron paramagnetic resonance spectra acquired at different resonant frequencies and the display forms (e.g., absorption and first-derivative displays) and the strategy of spectral parameter sharing could be potentially applicable to other heterogeneous biological systems in addition to the cotton fibers studies here.