A new set of cross-over marine data has been used to generate a regional model for the secular variation of the total geomagnetic field, showing the potential of the suggested approach for gaining a better knowledge of the field over oceanic regions. The model, which is valid for the Northern Atlantic region during the temporal interval 1960-2000, was obtained using spherical cap harmonic analysis (SCHA) in space and penalized splines in time. The maximum spatial expansion is equivalent to degree 9 in ordinary spherical harmonic analysis. Annual mean intensity data from different geomagnetic observatories have been used to improve the spatial and temporal resolution of the original dataset. Results indicate that the regional model improves, in terms of the root mean square error, the prediction given by the 11th generation of IGRF and CM4 global models, especially for the geomagnetic observatories considered. We also provide the uncertainty of the model coefficients and the secular variation prediction given by a bootstrap algorithm. The model is available in the EarthRef. org Digital Archive at http://earthref.org/ERDA/1728/.