Mineral resource exploitation by human societies throughout history led to the deposit of mining and smelting wastes and the subsequent contamination of surrounding soils by trace metals. After several centuries, the impact of these legacy hazardous wastes may remain a cause of environmental concern, especially for indigenous soil invertebrate populations such as earthworms. Therefore, we conducted a passive biomonitoring campaign in a former metallurgical district (Vosges Mountains, eastern France). According to community descriptors, we evidenced a significant decrease of anecic and endogeic earthworm density in the former mining stations. To link these results to soil contamination and bioaccumulation levels in earthworm tissues, we propose an original modelling approach using nonlinear mixed-effects regression models. Beyond a dose-response relationship between metal internal concentrations and their levels in soils, we highlighted contrasted behaviors according to ecological groups (epianecics and endogeics most impacted). We interpreted these results in relation to some eco-physiological features without completely exclude the influence of textural characteristics of soil, especially for deep-burrowing species such as anecic strict. Nonetheless, the presence of earthworm populations currently living in highly contaminated sites and handling elevated internal concentrations raises the question of the acquisition of genetic adaptive traits and the trophic transfers of metals.