Economic development in an area attracts more people to live in it. This increment drives the necessity to improve available infrastructure, like roads for instance, to satisfy a higher demand. Bigger roads and higher number of vehicles have raised the concern about possible pollution coming from these sources In this thesis, copper coming from road runoff in a peri-urban catchment in Portugal was analyzed. The catchment is located in the Coimbra region, center of Portugal. In order to model copper pollution in road runoff, it is necessary to couple a hydrological model and a water quality Based on a previous study (Kalantari, Ferreira, Walsh, Ferreira, & Destouni, 2017) a physical based hydrological model  MIKE SHE  coupled with the hydraulic model MIKE 11was updated and further coupled with  ECO Lab to simulate water quality and ecological processes. The results show an improvement of the hydrological model compared with the original one, nash-sutcliffe efficiency was raised from 0.59 to 0.77 and the coefficient of determination varied from 0.64 to 0.79. For copper the model behavior for punctual and distributed sources was analyzed. For punctual sources, highest concentrations were present in the grid points where the incoming sources were located in the tributaries, and these concentrations are rapidly reduced downstream. On the other hand, distributed sources approach gives higher concentrations near the end of the river than in the tributaries upstream. Comparing time-averaged model results along the river, with fresh water quality criteria according to U.S. EPA (2004), for punctual sources an extension of 978 meters (7,6% of the river) presents a concentration above CCC or CMC, on the contrary, for distributed sources the extension is lower with 494 meters (3,8%). The organic carbon partitioning coefficient have bigger influence on the results than other factors, nevertheless this influence is not marked. Modeled copper values do not agree well with the mesured values specially for periods with higher discharge as the model simulates lower concentration with higher discharge and viceversa. Copper boundary values for the model represent a big challenge considering limited data available. This thesis gives good overview about the coupling process between MIKE 11 and ECO Lab, as well as analyzes the importance of some factors as well as model limitations and uncertainties.