Artificial Neural Networks are used to find the influence of habitat types on the quality of the environment expressed by the concentrations of toxic and harmful elements in avian tissue. The main habitat types were described according to the Corine Land Cover CLC2012 model. Eggs of free-living species of a colonial waterbird, the grey heron Ardea cinerea, were used as a biological data storing media for biomonitoring. For modeling purposes, pollution indices expressing the sum of the concentration of harmful and toxic elements (multi-contamination rank index) and indices for single elements were created. In the case of all the examined indices apart from Cd, the generated topologies were a multi-layer perceptron (MLP) with 1 hidden layer. Interestingly, in the case of Cd, the generated optimal topology was a network with a radial basis function (RBF). The data analysis showed that the increase in environmental pollution was mainly influenced by human industrial activity. The increase in Hg, Cd, and Pb content correlated mainly with the increase in the areas characterized by human activity (industrial, commercial, and transport units) in the vicinity of a grey heron breeding colony. The decrease in the above elements was conditioned by relative areas of farmland and inland waters. Pollution with Fe, Mn, Zn, and As was associated mainly with areas affected by industrial activities. As the location variable did not affect the quality of the obtained networks, it was removed from the models making them more universal