This work describes the development, optimization, and validation of a new method for the simultaneous determination of a wide range of pharmaceuticals (beta-blockers, lipid regulators…) and personal care products (fragrances, UV filters, phthalates…) in both aqueous and solid environmental matrices. Target compounds were extracted from sediments using pressurized hot water extraction followed by stir bar sorptive extraction. The first stage was performed at 1,500 psi during three static extraction cycles of 5 min each after optimizing the extraction temperature (50-150 °C) and addition of organic modifiers (% methanol) to water, the extraction solvent. Next, aqueous extracts and water samples were processed using polydimethylsiloxane bars. Several parameters were optimized for this technique, including extraction and desorption time, ionic strength, presence of organic modifiers, and pH. Finally, analytes were extracted from the bars by ultrasonic irradiation using a reduced amount of solvent (0.2 mL) prior to derivatization and gas chromatography-mass spectrometry analysis. The optimized protocol uses minimal amounts of organic solvents (<10 mL/sample) and time (≈8 h/sample) compared to previous existing methodologies. Low standard deviation (usually below 10%) and limits of detection (sub-ppb) vouch for the applicability of the methodology for the analysis of target compounds at trace levels. Once developed, the method was applied to determine concentrations of these compounds in several types of sample (wastewater, seawater, pore water, and sediment) from Cadiz Bay (SW Spain). To our knowledge, these findings represent the first information available on the presence of some of the target compounds in the marine environment.