The Positive Matrix Factorization (PMF) receptor model and the Observation Based Model (OBM) were combined to analyze volatile organic compound (VOC) data collected at a suburban site (WQS) in the PRD region. The purposes are to estimate the VOC source apportionment and investigate the contributions of these sources and species of these sources to the O(3) formation in PRD. Ten VOC sources were identified. We further applied the PMF-extracted concentrations of these 10 sources into the OBM and found "solvent usage 1", "diesel vehicular emissions" and "biomass/biofuel burning" contributed most to the O(3) formation at WQS. Among these three sources, higher Relative Incremental Reactivity (RIR)-weighted values of ethene, toluene and m/p-xylene indicated that they were mainly responsible for local O(3) formation in the region. Sensitivity analysis revealed that the sources of "diesel vehicular emissions", "biomass/biofuel burning" and "solvent usage 1" had low uncertainties whereas "gasoline evaporation" showed the highest uncertainty.