Since 1950, increase in nitrogen (N) and phosphorus (P) river loadings in the North-East Atlantic (NEA) continental seas has induced a deep change in the marine coastal ecosystems, leading to eutrophication symptoms in some areas. In order to recover a Good Ecological Status (GES) in the NEA, as required by European Water Framework Directive (WFD) and Marine Strategy Framework Directive (MSFD), reductions in N- and P-river loadings are necessary but they need to be minimal due to their economic impact on the farming industry. In the frame of the “EMoSEM” European project, we used two marine 3D ecological models (ECO-MARS3D, MIRO&CO) covering the Bay of Biscay, the English Channel and the southern North Sea to estimate the contributions of various sources (riverine, oceanic and atmospheric) to the winter nitrate and phosphate marine concentrations. The various distributed descriptors provided by the simulations allowed also to find a log-linear relationship between the 90th percentile of satellite-derived chlorophyll concentrations and the “fully bioavailable” nutrients, i.e. simulated nutrient concentrations weighted by light and stoichiometric limitation factors. Any GES threshold on the 90th percentile of marine chlorophyll concentration can then be translated in maximum admissible ‘fully bioavailable’ DIN and DIP concentrations, from which an iterative linear optimization method can compute river-specific minimal abatements of N and P loadings. The method has been applied to four major river groups, assuming either a conservative (8 μg Chl L−1) or a more socially acceptable (15 μg Chl L−1) GES chlorophyll concentration threshold. In the conservative case, maximum admissible winter concentrations for nutrients correspond to marine background values, whereas in the lenient case, they are close to values recommended by the WFD/MSFD. Both models suggest that to reach chlorophyll GES, strong reductions of DIN and DIP are required in the Eastern French and Belgian-Dutch river groups.