Abstract Objective analysis of large-scale simulated anomalies is used to estimate statistically the effectiveness of different sampling arrays. The methodology, which includes a new use of Bayesian inference, is of fairly general applicability, but is illustrated by a specific requirement to measure five-year changes in the North Atlantic to an accuracy equivalent to ±10 W m−2. The numerical procedure require 1) specification of an overall measure of accuracy on an appropriate resolution, 2) a reference field approximating the long-term mean, 3) an assumed ensemble of large-scale anomalies, 4) the distribution of mesoscale eddy noise as derived from previous analyses of historical data and 5) a set of ship tracks and a sampling, in space and time along them. The numerical output is the expected accuracy and other diagnostic information. The ship tracks and sampling are varied by trial and error until the expected accuracy is within requirements in an economical manner. The study indicates that the optimi...