AbstractThe evolution of environmentally‐induced changes in phenotype or reaction norm implies both the existence at some time of genetic variation within a population for that plasticity measured by the presence of genotype x environment interaction (G x E), and that phenotypic variation affects fitness. Otherwise, the genetic structure of polygenic traits may restrict the evolution of the reaction norm by the lack of independent evolution of a given trait in different environments or by genetic trade‐offs with other traits that affect fitness. In this paper, we analyze the existence of G x E in metamorphic traits to two environmental factors, larval density and pond duration in a factorial experiment with Bufo calamita tadpoles in semi‐natural conditions and in the laboratory.Results showed no plastic temporal response in metamorphosis to pond durability at low larval density. The rank of genotypes did not change across different hydroperiods, implying a high genetic correlation that may constrain the evolution of the reaction norm. At high larval density a significant G x E interaction was found, suggesting the potential for the evolution of the reaction norm. A sibship (#1) attained the presumed “optimal” reaction norm by accelerating developmental rate in short duration ponds and delaying it in longer ponds. This could be translated in fitness by an increment in metamorphic survival and size at metamorphosis in short and long ponds respectively with respect to non‐plastic sibships. However, genetic variability for plasticity suggests that optimal reaction norm for developmental rates may be variable and hard to achieve in the heterogeneous pond environment.Mass at metamorphosis was not plastic across different pond durations but decreased at high larval density. Significant adaptive plasticity for growth rates appeared in environments that differed drastically in level of crowding conditions, both in the field and in the laboratory. The fact that survival of juveniles metamorphosed at high density ponds was a monotonic function of metamorphic size, implies that response to selection may occur in this population of natterjacks and that genetic variability in plasticity may be a reliable mechanism maintaining adaptive genetic variation in growth rates in the highly variable pond environment.