This paper presents a novel numerical environmental multimedia modeling system (RNEMM) for assessing the environmental fate of emerging organic contaminants and their relative health risk at a regional scale. The RNEMM is developed based on an integrated numerical algorithm that comprises four sub-models: a river network simulation module, a gaseous phase simulation module, a mass balance based simulation module for soil compartment, and a food web analysis module. This RNEMM has been applied to simulate the spatial distribution of PFOS and assess the consequent health risks for a central water basin region of the Pearl River in China. The study region includes the urban areas of Guangzhou, Foshan, and Dongguan Cities with emission sources of PFOS, which was detected in local water, sediments, and air environment. The spatial concentration distributions of PFOS in water, sediment, air, soil, and various fish species are examined based on RNEMM and compared with the measured data. With a focus on water environment, it shows that the simulated results essentially agree well with measured concentrations. Comparing the simulated results and the measured data collected in 2013, the relative errors are mostly less than 40 % in the surface water and sediment zones for this regional scale field study. Whereas the relative error in the atmosphere zone is less than 5%. In addition, the health risk assessment for children and adults is conducted based on the RNEMM approach. The hazard quotient (HQ) values for the 95th percentile in most subareas of the study region are higher than 0.1, showing a low-risk level for the study period. The results indicate that the RNEMM is a useful modeling tool to manage the environmental and health risks associated with emerging contaminants on regional air, water, soil, and ecosystem at an adequate spatial-temporal resolution.