Abstract Survival of larval fish is often linked to production of preferred prey such as copepods, both inter- and intra-annually. In turn, copepod production depends not only the quantity of food, but also on the nutritional quality, edibility and/or toxicity of their micro-algal food. Hence, larval fish survival can become de-coupled from levels of nutrient input depending on the resulting composition of the plankton. Here we use a plankton dynamics model to study nutrient input, phytoplankton composition and copepod, Paracalanus, production in relation to interannual variation in recruitment of snapper, Chrysophrys auratus, in Port Phillip Bay, Australia. The model was able to simulate the ratio of diatoms to flagellates in the plume of the main river entering Port Phillip Bay. Interannual variability in the copepod, Paracalanus, abundance during the C. auratus spawning period over 5 years was accurately predicted. The seasonal peak in Paracalanus production depended on the timing and magnitude (match-mismatch) of nutrient inputs and how these were reflected in temporal change in the diatom:flagellate ratio. In turn, the model-predicted Paracalanus abundance was strongly related to inter-annaul variability in abundance of snapper, C. auratus, larvae over 7 years. Years of highest larval C. auratus abundance coincided with a matching of the spawning period with the peak in Paracalanus abundance. High freshwater flows and nutrient inputs led to an early seasonal dominance of diatoms, and consequently reduced abundances of copepods over the C. auratus spawning period with correspondingly low abundances of larvae. Conversely years of very low rainfall and nutrient input also led to low phytoplankton and copepod concentrations and larval C. auratus abundances. Highest abundances of larval C. auratus occurred in years of low to intermediate rainfall and nutrient inputs, particularly when pulses of nutrients occurred in the spring period, the latter supporting the match-mismatch hypothesis.