Fisheries management has historically focused on the population elasticity of target fish based primarily on demographic modeling, with the key assumptions of stability in environmental conditions and static trophic relationships. The predictive capacity of this fisheries framework is poor, especially in closed systems where the benthic diversity and boundary effects are important and the stock levels are low. Here, we present a probabilistic model that couples key fish populations with a complex suite of trophic, environmental, and geomorphological factors. Using 41 years of observations we model the changes in eastern Baltic cod (Gadus morhua), herring (Clupea harengus), and Baltic sprat (Sprattus sprattus balticus) for the Baltic Sea within a Bayesian network. The model predictions are spatially explicit and show the changes of the central Baltic Sea from cod- to sprat-dominated ecology over the 41 years. This also highlights how the years 2004 to 2014 deviate in terms of the typical cod–environment relationship, with environmental factors such as salinity being less influential on cod population abundance than in previous periods. The role of macrozoobenthos abundance, biotopic rugosity, and flatfish biomass showed an increased influence in predicting cod biomass in the last decade of the study. Fisheries management that is able to accommodate shifting ecological and environmental conditions relevant to biotopic information will be more effective and realistic. Non-stationary modelling for all of the homogeneous biotope regions, while acknowledging that each has a specific ecology relevant to understanding the fish population dynamics, is essential for fisheries science and sustainable management of fish stocks.