Mushroom picking has become a widespread autumn recreational activity in the Central Pyrenees and other regions of Spain. Predictive models that relate mushroom production or fungal species richness with forest stand and site characteristics are not available. This study used mushroom production data from 24 Scots pine plots over 3 years to develop a predictive model that could facilitate forest management decisions when comparing silvicultural options in terms of mushroom production. Mixed modelling was used to model the dependence of mushroom production on stand and site factors. The results showed that productions were greatest when stand basal area was approximately 20 m ha. Increasing elevation and northern aspect increased total mushroom production as well as the production of edible and marketed mushrooms. Increasing slope decreased productions. Marketed Lactarius spp., the most important group collected in the region, showed similar relationships. The annual variation in mushroom production correlated with autumn rainfall. Mushroom species richness was highest when the total production was highest.