A 3D hydrodynamic model (FVCOM) coupled to a carbonate system (ERSEM) has been used to model a number of seabed CO2 release scenarios ranging from 3 to 3000 t d -1 for the Goldeneye complex in the northern North Sea. The results of the scenario runs were used to characterise the fate of CO2 in the water column in space and time. A new approach to designing monitoring networks has been implemented and compared with a simple approach. A weighted greedy set algorithm is used to identify the positions within the model domain which yield the greatest combined coverage for the smallest number of sampling stations, further limited by selecting only a feasible number of sample sites. The weighted greedy set algorithm incorporates the effect of the unstructured grid in FVCOM as well as the proximity of the candidate sample locations to the Goldeneye complex. For the range of release rates simulated, the design of the optimal sampling strategy changes depending on the magnitude of the release. The role of the tides discriminates the four release scenarios into two categories: for the lower release rates (3 and 30 t/d), the effect of the tide is relatively unimportant in the distribution; for the larger release rates (300 and 3000 t/d), the direction of the principal tidal axis controls the distribution of the sampling stations more strongly. Comparison of the weighted greedy set approach shows it is able to identify releases sooner and with a stronger signal than a simple regular sampling approach.