Quantitative study of urban-rural gradients for persistent organic pollutants (POPs) concentrations is extremely important for understanding POPs behavior and for ecological risk assessment and management. However, little has been studied about quantitative calculation for soil concentrations of POPs from urban to rural areas. In this paper, 153 soil samples were collected in the coastal region of Bohai Sea and Yellow Sea, and the distributions of polychlorinated biphenyls, Polyaromatic hydrocarbons, Perfluoroalkyl and polyfluoroalkyl substances were examined. Urbanrural gradient model (URGM) was derived by using atmospheric point source diffusion model combined with fugacity model to prove potential mathematical relationships among urban and rural soils, and the URGM was proved suitable for simulating urbanrural soil concentrations through validation of measured data. Significantly linear correlation was found between POPs amount in surface soil and city population, between POP concentrations and the artificial surface area. Urban-rural POPs concentrations were simulated by URGM model, and calibrated by urban population and land-cover data. The results showed that it was suitable for simulating urban-rural POPs concentrations at a regional scale. This study could provide a new method for quantifying urban-rural gradients for POPs concentrations, and promote quantitative research on coupling between land cover, socioeconomic data and POPs concentrations.