A novel method based on ultrasonic solvent extraction and stir bar sorptive extraction (SBSE) for the analysis of 51 persistent organic pollutants including organochlorine pesticides (OCPs), polychlorinated biphenyls (PCBs), polycyclic aromatic hydrocarbons (PAHs), and polybrominated diphenylethers (PBDEs) in soil samples was developed. The different parameters that affect both the extraction of analytes from the soil samples, such as solvent selection, solvent volume, mass of soil, and extraction time, and the partitioning from the solvent/water mixture to the PDMS were studied. The final selected conditions consisted of the extraction of 1 g of soil with 15 mL methanol by sonication for 30 min. The methanol extract was mixed with 85 mL of Milli-Q water and extracted by means of SBSE for 14 h at 900 rpm. The stir bars were analyzed by thermal desorption-GC-mass spectometry (TD-GC-MS). The effects of the matrix on the recovery of the various pollutants under the developed method were studied using two soils with very different physicochemical properties. Method sensitivity, linearity, repeatability, and reproducibility were also studied. Validation and accuracy of the method were conducted by analyzing two commercial certified reference materials (CRMs). The main advantage of this method resides in the fact that a small amount of a nontoxic solvent (methanol) is needed for the extraction of only 1 g of solid sample allowing LODs ranging from 0.01 to 2.0 microg/kg. Repeatability and reproducibility variations were lower than 20% for all investigated compounds. Results of the CRMs verify the high accuracy of this method.