The petrochemical industry is a major contributor to air pollution emissions, so determining the variations of these pollutants in the petrochemical industrial area is needed to reduce health risks due to PM_(10) exposures. In this study, a specific industrial monitoring database was used to collect PM_(10) data from May 2015 to September 2018 in Kaohsiung Petrochemical Parks area (including Linyuan Industrial Park and Linhai Industrial Park), Taiwan. Conventional land use regression and hybrid kriging approach which consider both the spatial and temporal heterogeneity of air pollutants was then conducted to develop the prediction model for estimating the spatial-temporal variability of PM_(10) concentrations. The resultant models showed a high-estimate performance level with the value of R^2, Adjusted R^2, RMSE and cross-validation for conventional approach was 0.89, 0.89, 7.29 μg/m^3 and 0.83, 0.83, 9.0 μg/m^3, respectively; and for hybrid approach was 0.95, 0.95, 5.2 μg/m^3 and 0.94, 0.94, 5.41 μg/m^3, respectively. The hybrid approach showed better explanatory power than conventional approach. Moreover, the results of hybrid model simulation showed that the highest concentrations of PM_(10) formed in the pure residential area, the manufacturing areas and agricultural fields; a similar decreasing trend was observed as well.