Analysing the estuarine use patterns of juveniles of marine migrant fish species is vital for identifying important sites for juveniles as well as the basic environmental features that characterize these sites for different species. This is a key aspect towards understanding nursery function. Various estuarine systems along the Portuguese coast (Minho, Douro, Ria de Aveiro, Mondego, Tejo, Sado, Mira, Ria Formosa and Guadiana) were sampled during Spring and Summer 2005 and 2006. Juveniles of commercially important marine fish species Solea solea, Solea senegalensis, Platichthys flesus, Diplodus vulgaris and Dicentrarchus labrax, predominantly 0-group individuals, were amongst the most abundant species and had distinct patterns of estuarine use as well as conspicuous associations with several environmental features. Juvenile occurrence and density varied amongst estuaries and sites within them, and differed with species. Sites with consistently high juvenile densities were identified as important juvenile sites (i.e. putative nursery grounds). Through generalized linear models (GLM), intra-estuarine variation in occurrence and density of each of the individual species was largely explained by environmental variables (temperature; salinity; depth; percentage of mud in the sediment; presence of seagrass; importance of intertidal areas; relative distance to estuary mouth; macrozoobenthos densities; and latitude). Decisive environmental factors defining important sites for juveniles varied depending on the system as a result of different environmental gradients, though there were common dominant features for each species regardless of the estuary considered. Analysed environmental variables in the GLM also accounted for inter-estuarine variation in species' occurrence and density. In several estuaries, the identified important juvenile sites were used by many of these species simultaneously and may be of increased value to both management and conservation. Overall, the variability in site features amongst estuaries highlighted the tolerance of these species to different available environmental conditions and provided fundamental information for future spatially explicit modelling of their distribution. This should ultimately enable the prediction of species response to habitat alterations.