Ongoing eutrophication frequently causes toxic phytoplankton blooms. This induces huge worldwide challenges for drinking water quality, food security and public health. Of crucial importance in avoiding and reducing blooms is to determine the maximum nutrient load ecosystems can absorb, while remaining in a good ecological state. These so called critical nutrient loads for lakes depend on the shape of the load-response curve. Due to spatial variation within lakes, load-response curves and therefore critical nutrient loads could vary throughout the lake. In this study we determine spatial patterns in critical nutrient loads for Lake Taihu (China) with a novel modelling approach called Spatial Ecosystem Bifurcation Analysis (SEBA). SEBA evaluates the impact of the lake's total external nutrient load on the local lake dynamics, resulting in a map of critical nutrient loads for different locations throughout the lake. Our analysis shows that the largest part of Lake Taihu follows a nonlinear load-response curve without hysteresis. The corresponding critical nutrient loads vary within the lake and depend on management goals, i.e. the maximum allowable chlorophyll concentration. According to our model, total nutrient loads need to be more than halved to reach chlorophyll-a concentrations of 30–40 μg L⁻¹ in most sections of the lake. To prevent phytoplankton blooms with 20 μg L⁻¹ chlorophyll-a throughout Lake Taihu, both phosphorus and nitrogen loads need a nearly 90% reduction. We conclude that our approach is of great value to determine critical nutrient loads of lake ecosystems such as Taihu and likely of spatially heterogeneous ecosystems in general.