A total number of 471 soil samples collected during the period of 1996 2006 from the agricultural and forest areas of the Czech Republic were analyzed for their content of persistent organic pollutants (POPs). Spatial variability of the POP concentrations was assessed using an IDW spatial GIS model analysis. For every grid of the network, resulting modeled levels of contamination allowed for estimation of the total burden of POPs in soils. Potential risks associated with contaminated soils were assessed as well. Database of the old ecological burdens counting 3061 sampling sites was used to adjust the model and incorporate the risks of nesvily contaminated sites. The high levels of health risks were only found at less than 1% of the area of interest. The IDW modeling proved to be a useful tool for screening of the health risks in the large areas with scarce monitoring data. Presented approach can be applied in the risk management, to support an efficient targeting of the risk reduction measures, or to improve a design of the national monitoring.