Sorption linearity and reversibility are implicit in models for the fate and transport of per- and polyfluoroalkyl substances (PFAS). In this study, however, we found that the sorption of cationic and zwitterionic PFAS in natural soils was highly nonlinear. The nonlinearity was so severe that it led to a variation in the coefficient of sorption by several orders of magnitude over the experimental concentration range. This implies a considerable increase in sorption as concentration falls in the natural environment. Sorption of cationic PFAS correlated strongly with the soil organic matter (SOM) content and was reversible in all soils. Sorption of zwitterionic PFAS, on the other hand, displayed concentration-dependent hysteresis in soils with a low SOM content. The irreversibility, which was associated with neither SOM, pore deformation, nor surface complexation, was likely caused by the entrapment of molecules in porous structures within inorganic components of soil aggregates. Furthermore, electrostatic interactions with negatively charged soil constituents and the hydrophobic effect were found to be major sorption driving forces for cationic/zwitterionic PFAS at low and high concentrations, respectively. The maximum electrostatic potential of PFAS ions, computed using density functional theory, was found to be a useful predictor of the sorption of ionic PFAS species.