Articles | Volume 19, issue 5
Research article
21 Sep 2023
Research article |  | 21 Sep 2023

Evaluating altimetry-derived surface currents on the south Greenland shelf with surface drifters

Arthur Coquereau and Nicholas P. Foukal

The pathways and fate of freshwater in the East Greenland Coastal Current (EGCC) are crucial to the climate system. The EGCC transports large amounts of freshwater in close proximity to sites of deep open-ocean convection in the Labrador and Irminger seas. Many studies have attempted to analyze this system from models and various observational platforms, but the modeling results largely disagree with one another, and observations are limited due to the harsh conditions typical of the region. Altimetry-derived surface currents, constructed from remote-sensing observations and applying geostrophic equations, provide a continuous observational data set beginning in 1993. However, these products have historically encountered difficulties in coastal regions, and thus their validity must be checked. In this work, we use a comprehensive methodology to compare these Eulerian data to a Lagrangian data set of 34 surface drifter trajectories and demonstrate that the altimetry-derived surface currents are surprisingly capable of recovering the spatial structure of the flow field on the south Greenland shelf and can mimic the Lagrangian nature of the flow as observed from surface drifters.

1 Introduction

Along the east Greenland shelf, meltwater from both the Arctic pack ice and the Greenland ice sheet flow southward toward the southern tip of Greenland (Foukal et al.2020). This southward flow is concentrated in two current cores that are each supported by bathymetric gradients: the East Greenland Current (EGC) at the shelf break and the East Greenland Coastal Current (EGCC) at the coastline (Håvik et al.2017). Around Greenland, the deep shelves (150–500 m), steep gradient at the coastline, and strong along-shelf barrier winds sustain this two-core system. If the buoyant meltwater water masses in the EGCC are mixed offshore, they may stratify regions of deepwater formation and water mass transformation, inhibit air–sea heat exchange, and slow the Atlantic meridional overturning circulation (AMOC; Rahmstorf et al.2015; Hansen et al.2016; Glikson2019).

As these currents approach Cape Farewell at the southern tip of Greenland, strong northerly barrier winds on the southeast Greenland shelf push the fresh, surface waters toward the coastline and into the EGCC (Duyck and De Jong2021). At Cape Farewell, the downwelling-favorable barrier winds relax as they meet the northerly, upwelling-favorable winds from the southwest Greenland shelf (Pacini and Pickart2023). The confluence of these wind regimes at Cape Farewell leads to the mean winds crossing the shelf and potentially allowing for large amounts of freshwater to leave the shelf here.

Many modeling studies have analyzed this potential for shelf–basin exchange around Cape Farewell (Fichefet et al.2003; Marsh et al.2010; Dukhovskoy et al.2016; Böning et al.2016; Gillard et al.2016; Luo et al.2016; Schulze Chretien and Frajka-Williams2018; Oliver et al.2018; Dukhovskoy et al.2019; Pennelly et al.2019; Castelao et al.2019; Garcia-Quintana et al.2019; Tagklis et al.2020; Dukhovskoy et al.2021; Gou et al.2021, 2022; Duyck et al.2022; Morrison et al.2023). There are many conflicting results between these various publications; notably some detect a very clear signal of Greenland meltwater influencing the large-scale circulation (Garcia-Quintana et al.2019), while others document a more muted oceanic response to the freshwater forcing (Böning et al.2016). Dukhovskoy et al. (2019) note the importance of the vertical mixing rate in the effect of this freshwater and the rapid dilution of the meltwater signature if mixed below the upper 100 m. In contrast, Tagklis et al. (2020) and Pennelly et al. (2019) find a clear relationship with freshwater input and shoaling of the mixed-layer depths in the Labrador Sea. Luo et al. (2016) and Gillard et al. (2016) both note the strong sensitivity to where the meltwater is introduced from Greenland to the shelf, with meltwater introduced on the east Greenland shelf being more likely to leave the shelf in the Labrador Sea than meltwater introduced on the west Greenland shelf. Schulze Chretien and Frajka-Williams (2018) and Duyck et al. (2022) report an Ekman, wind-driven mechanism to transport freshwater off the Greenland shelf, while Tagklis et al. (2020) note the importance of boundary current instabilities in driving the shelf–basin exchange.

Reconciling these various model results requires an observational data set of ocean currents that covers geographic extents similar to the models. Satellite altimetry is an example of such a data set. Though each individual pass of a nadir-altimeter satellite measures just a single line of sea surface height (SSH) directly below the satellite’s path, the combination of multiple satellites that cross a region can yield significant data coverage, even over smaller regions such as the southern tip of Greenland. These individual passes are then gridded onto a 14 global grid (Rio et al.2011) to yield a spatially contiguous and temporally continuous data set of SSH. The pressure gradients resulting from the structure in the SSH field are then used to calculate surface geostrophic velocities (Rio et al.2011). In addition, Ekman velocities derived from atmospheric reanalyses are combined with the geostrophic velocities to yield a surface velocity that more accurately accounts for the large-scale wind-driven flow (Rio et al.2014). It is often assumed that the gridded SSH data do not perform well on the continental shelves due to the smaller-scale features in the coastal circulation, as well as errors in the altimetry associated with sampling close to coastlines (Pujol et al.2023). The Greenland shelf is unique in this regard though because it is deep (150–500 m), the bathymetric gradient from the coast to the shelf is steep, and the shelf circulation is vigorous (velocities on the order of 0.4–1 m s−1). Thus determining whether altimetry-derived surface currents (ADSCs) in the region are valid is of interest.

In August 2021, we deployed 38 surface drifters on the southeast Greenland shelf to determine the pathways of the shelf circulation around Cape Farewell. This observational campaign provides an unprecedented view into the surface circulation of the region, though it is limited to about 2 weeks during which the drifters passed through the area. Here, we use these drifters to evaluate the performance of the ADSCs and assess whether the gridded altimetry data product is suitable for longer-term studies. The paper is structured as follows: in Sect. 2, we describe the data sets and a novel three-step methodology to robustly compare Lagrangian and gridded Eulerian data. In Sect. 3, we present the results of this methodology. And in Sect. 4 we discuss the results in light of other published literature and implications for future work.

Figure 1(a) The study region, with Cape Farewell at 44 W, Julianehaab Bight at 46 W, and Cape Desolation at 48 W. Contour lines represent isobaths. (b–c) Deployment location of drifters and floats, respectively, for the (b) first (15 August) and (c) second (20 August) batch. (d–e) Drifter trajectories, respectively, for the (d) first and (e) second batch. Colors represent lines of drifters. (f) Péclet number derived from the lines of drifters, representing the ratio of advection to dispersion. The trajectories of the centroids of each line are represented; overlapping data are averaged. The calculation of the Péclet number is detailed in Sect. 2.1. Please note that the date format used in the figures is month/day/year.

2 Data and methods

2.1 Description of surface drifters

In August 2021, 38 Surface Velocity Program (SVP) drifters were released on the southeastern shelf of Greenland (Fig. 1a) from R/V Neil Armstrong. These surface drifters consist of a surface buoy attached to a drogue that is centered at 15 m depth and a drogue-to-buoy drag ratio of 40:1 to ensure drifters follow the 15 m current and are not directly pushed by the wind. Two distinct programs provided these drifters: the (US) National Science Foundation-funded Greenland Fresh Water Experiment (GFWE) with 12 drifters and the (EU) European Research Council-funded Targeted Experiment to Reconcile Increased Freshwater with Increased Convection (TERIFIC) with 26 drifters. The drifters were separated evenly into two batches of 19 drifters that were deployed 5 d (d denotes days) apart (15 and 20 August). Each batch consisted of one line of four drifters deployed at the shelf break and three zonal lines of five drifters deployed across the inner shelf. The three lines were separated meridionally by 25–30 km between one another, and the drifters along each line were separated by about 5 km (Fig. 1b, c).

A number of parameters were then checked. The continuous presence of the drogue was verified using two methods: (1) the buoy parameters and (2) the coefficient of determination R2 between complex drifter velocity and complex wind velocity (from ERA5; Hersbach et al.2020) or similarly the least-squares complex linear regression of drifter velocity from wind velocity (Kundu1976; Poulain et al.2009). Both methods confirmed that all drifters maintained their drogues throughout the period of analysis described in this paper. We also tested whether any of the data points qualified as outliers according to the Elipot et al. (2016) definitions, but we did not identify any. Four of the drifters were trapped in fjords for prolonged periods; thus we did not consider them further. These steps left us with 34 trajectories to analyze. We then filtered out high-frequency variability using a 24 h cutoff low-pass filter. More details on drifter data processing are available in Sect. S1 of the Supplement.

To assess the relative importance of advection compared to dispersion, we calculated the Péclet number (Fig. 1f) from the drifter data set. This is calculated by considering each line of drifters as a group. For each line and time step, the dispersion of the drifters is computed as the root mean square distance of drifters from the center of gravity of the line. The advection is estimated from the distance traveled by the center of mass of the line since the previous time step. The advection is then divided by the dispersion to obtain the Péclet number. Values greater than 1 indicate the advective component is dominant.

2.2 Description of satellite altimetry data sets

Two data sets of ADSCs have been evaluated in this study. The first is the Copernicus Marine Environment Monitoring Service (CMEMS) near-real-time product denominated Global Ocean Gridded L4 Sea Surface Heights And Derived Variables Nrt. The processing used is the DUACS (Data Unification and Altimeter Combination System) multi-mission altimeter data processing system DT-2021 (Faugère et al.2022) provided by the CNES CLS (Centre National d'Études Spatiales Collecte Localisation Satellites) with a methodology detailed in Pujol et al. (2016). The altimetry data are merged from all available altimetry missions (Fig. S1 in the Supplement) and interpolated on a 14 grid with a daily resolution. The geostrophic currents provided in this data set benefit from the CNES-CLS18 mean dynamic topography (MDT) product (Mulet et al.2021). This MDT gathers altimetry, gravity, and drifter data. It shows better results in all regions around the globe compared to the previous product, CNES-CLS15, especially in coastal areas. The geostrophic currents are computed using a nine-point stencil width methodology (Arbic et al.2012). In the following work, we refer to this data set as “Geostrophy”.

The second ADSC data set is the Global Total Surface and 15 m Current (COPERNICUS-GLOBCURRENT) from Altimetric Geostrophic Current and Modeled Ekman Current Processing. The Ekman velocities at 15 m depth are computed from ERA5 (Hersbach et al.2020) 3 h wind data following the methodology developed in Rio and Hernandez (2003), Rio et al. (2011), Rio (2012), and Rio et al. (2014). An empirical Ekman spiral-like model is estimated based on two parameters determined from a least-squares regression from SVP drifter data, Argo float data, and wind stress measurements from ECMWF ERA5. We refer to this data set as “Geostrophy + Ekman” because it is the addition of the Ekman contribution to the previous geostrophic-current data set. The final product provides a 6 h frequency data set with zonal and meridional components of surface velocity. For fair comparison with the daily Geostrophy product, here we primarily evaluate daily velocities but have also looked at the 6 h frequency product with higher-frequency wind data to see the impact of higher temporal resolution. The 15 m depth velocities are used for consistency with a 15 m depth drogue drifter.

In addition to the two ADSC data sets, we also compare the drifter velocities to the Arctic Ocean Physics Analysis and Forecast product from the operational TOPAZ4 Arctic Ocean system (Sakov et al.2012) with the same 14 resolution. This is a completely different data set in that it is a model that assimilates all possible data streams (e.g., Argo, satellites) rather than only relying on satellite altimetry and wind products. The goal of including TOPAZ in the first step of this methodology is to provide context to the results from the ADSC data sets. This data assimilation product is based on HyCOM (Hybrid Coordinate Ocean Model) and a 100-member ensemble Kalman filter (EnKF) assimilation scheme.

2.3 Motivation for a framework to compare Lagrangian and gridded Eulerian velocity fields

The goal of this paper is to directly compare the velocity fields derived from the surface drifters with those derived from altimetry. But direct comparisons between Eulerian and Lagrangian data are difficult because the reference frames provide fundamentally different information about the flow field. Eulerian data inform how the flow field evolves through time, while Lagrangian data provide information on the origin, pathways, and fate of fluid particles. The spatial scales also typically vary between these two reference frames. While Eulerian data can come from point measurements or data mapped onto a regular grid that typically cover larger spatial scales, Lagrangian data only come from point measurements, typically on smaller scales. Comparing point measurements from Eulerian and Lagrangian reference frames is not useful in most applications – these two data streams are only directly comparable when they physically intersect, and, when they do, it is straightforward to compare them. Thus, we do not discuss this comparison any further. In contrast, the need to compare gridded Eulerian data with Lagrangian data arises quite often in oceanography. The proliferation of gridded satellite and reanalysis products, as well as numerical-model output, has produced a large number of gridded Eulerian data sets (e.g., Haine et al.2021). These data sets are often compared to surface drifters from the Global Drifter Program (GDP; Lumpkin and Johnson2013) and profiling floats from the Argo network (Johnson et al.2022) to yield a velocity field that utilizes the accuracy of the in situ Lagrangian data with the spatial perspective of the gridded Eulerian data.

Figure 2 Representation of the point-wise comparison. The hourly Lagrangian trajectory (gray) is averaged to daily resolution (black) to match the resolution of the altimetry-derived surface current (ADSC) field. The Lagrangian velocity at each daily location is calculated as the distance between the location 12 h before the daily location and the location 12 h after, divided by 24 h. The nearest grid cells of the Eulerian product (blue dots with red edges) are determined along the Lagrangian trajectory, and the velocity at those points is compared to the Lagrangian velocity. The gray areas represent the coastline, and the contours represent the bathymetry.

Though the need for a robust framework for comparing gridded Eulerian and Lagrangian velocities exists, there is not yet a well-accepted, published methodology for these comparisons. Previous work on this topic typically converts one of the data sets into the other reference frame and then applies standard statistical methods such as correlation or variance metrics (Liu and Weisberg2011; Liu et al.2014; Rio et al.2011; Pujol et al.2016; Rio and Santoleri2018; Mulet et al.2021). This can be done by either gridding Lagrangian velocities onto a Eulerian grid or simulating Lagrangian trajectories through a Eulerian velocity field. While both methods yield data sets that can be directly compared, the process of transforming the data between reference frames inevitably degrades the converted data by interpolating or extrapolating when/where data are not available. If this comparison is only done in one direction, this process likely biases the results toward one reference frame. Thus, to robustly compare Eulerian and Lagrangian data, this conversion must proceed in both directions.

2.4 Comparison framework

In this work, we propose a framework from which Eulerian and Lagrangian velocity fields can be directly compared. The framework directly addresses three questions that leverage the relative strengths of each perspective. How well do the gridded Eulerian velocities resolve the velocities directly observed by the Lagrangian platform (step no. 1)? How well do the Lagrangian data recover the spatial structure of the gridded Eulerian data (step no. 2)? How well do the Eulerian data capture the origins, pathways, fate, and connectivity of water masses (step no. 3)? To answer these questions, we use three analysis steps: directly compare the velocity fields as measured at collocated and contemporaneous points (step no. 1), map the Lagrangian velocities onto a Eulerian grid (step no. 2), and simulate Lagrangian trajectories in the Eulerian flow fields (step no. 3). The first step directly compares the two data sets without any transformation between reference frames, while the second two steps convert one data set into the other's reference frame. All three steps utilize the strengths of the data sets: the accuracy of the Lagrangian data (step no. 1); the spatial perspective of the Eulerian data (step no. 2); and the origins, pathways, fate, and connectivity of the Lagrangian data (step no. 3).

2.4.1 Point-wise comparison between gridded Eulerian and Lagrangian velocity fields

The point-wise method consists of extracting the velocities of the gridded Eulerian field along the Lagrangian trajectories (Fig. 2) and comparing them to the Lagrangian velocities. To select the collocated and contemporaneous node on the ADSC grid that matches a specific spatiotemporal drifter data point, we extracted the nearest grid point from the drifter location corresponding at the same time step. It is also possible to use other interpolation method such as linear interpolation or optimal interpolation, though we found that using the nearest grid cell retained more of the variance in the ADSCs and was thus preferable to interpolation schemes. Here, the point-wise method has been performed on drifter daily velocities to extract the collocated and contemporaneous velocities from the ADSCs. Four components of the velocities have been evaluated: the zonal (u), the meridional (v), the along-shelf, and the across-shelf velocities. These two latter velocity components, which are simply projections of the velocity vectors according to a geographic feature (here, the shelf break), tell us how the velocity field moves independently from the shape of the bathymetric gradients at the shelf break. Practically, this involves first defining a smoothed shelf break by using low-pass filtering of the 600 m isobath to remove the small-scale features such as submarine canyons and troughs. The cross-shelf vector (zcross, in a complex form) is obtained by connecting a straight line from the surface drifter to the closest point of the shelf break. The along-shelf angle is defined to be perpendicular to the across-shelf angle. Figure S2 provides a visual representation of this coordinate transformation. The velocity vector (zvelocity) at the position of the surface drifter is computed as zvelocity=u+iv. The angle θ between the cross-shelf vector and the drifter velocity vector is computed as

(1) θ = arg z cross z velocity v along = sin θ . | z velocity | v across = cos θ . | z velocity | .

Once the velocities have been rotated, we extract a time series of ADSC velocities along each drifter trajectory. For each drifter, we compare the standard deviation (normalized by drifter velocity standard deviation, σf^), the correlation coefficient r relative to the drifters' daily velocities, the root mean square error (RMSE), and the percent variance explained between the two time series (Eq. 2).

(2) Variance of X explained by Y (%) = 100 × 1 - σ 2 ( X - Y ) σ 2 ( X )

Taylor diagrams and skill scores S (Eq. 3) (Taylor2001) are used to concatenate and summarize the standard deviation, correlation coefficient, and RMSE. We compare these skill scores to the TOPAZ4 reanalysis product, which serves as a baseline from which to compare the ADSC values.

(3) S = 4 ( 1 + r ) ( σ f ^ + 1 / σ f ^ ) 2 ( 1 + r 0 ) ,

with r0 as the maximum expected correlation coefficient (here taken equal to 1).

Figure 3 Schematic of the averaging of velocities. This method is intended for Lagrangian instruments, but step nos. 4–6 can be applied to the Eulerian product using the velocity at the nearest grid cell and corresponding in time to the Lagrangian location.


2.4.2 Eulerian gridding of Lagrangian velocities

The Eulerian gridding consists of mapping the Lagrangian data into a Eulerian gridded field. The gridded field is obtained by averaging the Lagrangian velocities on the same 14 grid as the ADSCs.

Rather than the classical averaging of all Lagrangian devices available together in the same grid, we use an alternative method that corrects for the number of data points reported in a given area by slowly and quickly moving instruments. Otherwise, slow instruments have more position points in each bin and thus exert a stronger influence on the average.

To perform the mapping (Fig. 3), we proceed sequentially through each individual drifter and produce a unique gridded map based on data from a single drifter. For each Lagrangian trajectory s, all the cells are scanned to check if the drifter has passed inside. If the drifter has not passed through the cell, it is left empty. Otherwise, all velocities within the cell are averaged. This average is weighted to attribute more weight to the observations close to the center using a Gaussian distribution.

The grid U obtained for each Lagrangian instrument therefore contains empty cells (where it has not passed) and cells containing the average velocity measured inside (Eq. 4). A similar grid D is filled with the average distances from the center of the cells calculated using the same method (Eq. 5). This processing is repeated for each of the 34 drifters, and the velocities across the drifters are averaged together according to their weighted distances from the center of the cell (Eq. 6).

Mathematically, this process can be explained by taking u(s,x), the velocity of the instrument s at location x=(x,y). The gridded velocity field U(s,c) of the Lagrangian trajectory s at the grid cell c located at xc=(xc,yc) is computed as

(4) U ( s , c ) = x u ( s , x ) . q ( s , x ) x q ( s , x ) ,

and the grid of mean distance D is computed as

(5) D ( s , c ) = x d ( s , x ) . q ( s , x ) x q ( s , x ) ,

with x=(x,y) and only if xxc-Δx2;xc+Δx2yyc-Δy2;yc+Δy2 x means sum all locations verifying the previous condition (selecting only locations in the cell c) and q(s,x) is the weight of the location x of instrument s depending on the distance from the center with, as an example, a Gaussian weighting like

(6) q ( s , x ) = exp - d 2 0.5 ( Δ x ) 2 + ( Δ y ) 2 .

This processing is repeated for each of the drifters, and the 34 grids are then averaged. For each cell, the average of the 34 grids in that cell is calculated (Eq. 7) and weighted by the averaged distances to the center (Eq. 8). The final Eulerian velocity grid 𝓤 combines the 34 grids and thus the trajectories of all the drifters.


The velocity field is not entirely covered by Lagrangian instruments at all times. Our drifters were initially located on the southeast Greenland shelf and moved progressively around Cape Farewell to the west. So the drifters did not sample the southwest Greenland shelf initially and the southeast Greenland shelf near the end of the 2-week period. To properly account for this heterogeneous sampling of the drifters, we constructed a comparison data set by subsampling the ADSC fields in a fashion similar to the drifters sampling the real ocean. This procedure accounted for the temporal biasing that often occurs when comparing in situ observations with other more homogenously sampled data sets. To assess whether our results were sensitive to this procedure, we compared the results with this subsampled ADSC product with the mean of the entire ADSC fields over the same time period, and they were nearly identical. The gridded Lagrangian and Eulerian products are then compared to evaluate the differences relative to spatial patterns of velocity and directions and magnitudes of currents.

As pointed out by LaCasce (2008) the choice of bin size in this type of analysis is not trivial. Bins that are too large will smooth out the field, and bins that are too small will make the mean very sensitive to eddies and other fine-scale structures. In most of the cases, a convenient choice to facilitate the comparison is to use a grid identical to the original Eulerian gridded products. However, to evaluate the influence of gridding and the consistency of the 14 drifter gridding with the original drifter data, we compared the results at 14 resolution to 112. This higher resolution was chosen to test the sensitivity of the results to the resolution of the grid and ensure that an average of four drifters passed through each cell. A higher resolution would imply fewer drifters on average. We also tested another gridding methodology based on a k-means clustering as proposed by Koszalka and LaCasce (2010) that has the advantage of directly using the density of data to determine the location of grid nodes and does not require arbitrary choices regarding the grid resolution and bounds. However, due to the heterogeneous distribution of drifters on the region, the clustering method led to spurious results.

2.4.3 Observed and synthetic trajectories

The third part of the comparison framework consists of computing synthetic Lagrangian trajectories from the ADSC Eulerian field (Fig. 4). This corresponds to converting the Eulerian data into the Lagrangian frame of reference. The main idea is to compute a synthetic trajectory by advecting a Lagrangian particle through the Eulerian velocity field (Liu and Weisberg2011). We use the Parcels v2.0 (Delandmeter and van Sebille2019) particle tracking software to simulate the trajectories. The resulting synthetic trajectories must then be compared to the observed trajectories of drifters and evaluated using a metric. Various methods could be used such as the distance between the ending points of trajectories, the main direction of trajectories, or a cloud of points. We have found the methodology proposed by Liu and Weisberg (2011) particularly instructive: at each location of the Lagrangian data, a new synthetic particle is released and advected through the velocity field for 3 d. Then, each 3 d trajectory is compared to the actual trajectory over those 3 d, and a skill score is assessed for each particle launch. This skill score (Eqs. 9, 10) combines the cumulative distance traveled by observed particles (Lagrangian devices) dli and the cumulative separation distance between synthetic and observed particles di. It is important to note that this skill score is quite rigorous; to produce a perfect skill score (equal to 1), the simulated particle must not only end up at the same location after 3 d but also follow the exact pathway to the end point.

(9)c=i=13dii=13(j=1idlj),(10)s=1-c,if c10,otherwise

A disadvantage of this metric is its sensitivity to the duration of the advection of the particles. As proposed by Liu and Weisberg (2011), we tested the sensitivity of our results to the duration of the particles’ advection by repeating the experiment with various durations spanning from 2 to 14 d. We selected 3 d because it identified differences in trajectories over the synoptic wind-forcing timescale in the region. Another potential issue with this skill score derived in Liu and Weisberg (2011) is that it sets negative skill score values to 0 and alias regions of poor skill toward better scores. One example of an alternative skill score is introduced by Révelard et al. (2021); it is nearly identical to that in Liu and Weisberg (2011) but retains the negative skill scores. We compare both metrics and do not identify large differences.

Figure 4 Illustration of the method involving a comparison between observed and synthetic particles using the approach of Liu and Weisberg (2011). (a) Skill scores (colored circles) calculated along a Lagrangian instrument trajectory (black line). Blue lines represent synthetic particle trajectories run over 3 d. Dot color represents skill score at the daily observed location where the synthetic particle has been released. The orange star represents the location where the drifter was released. Gray contours are isobaths. (b) A schematic of the various variables used to calculate the skill score.

Figure 5 Taylor diagrams of the point-wise comparison between the surface drifter velocities and either ADSCs or TOPAZ4 for the (a) zonal, (b) meridional, (c) along-shelf, and (d) across-shelf components. The x and y axes depict the normalized standard deviation (ratio of the ADSC standard deviation to the surface drifter standard deviation). Polar coordinates represent the correlation coefficient between product and drifter time series. Contours represent Taylor skill scores. Each trajectory for each product is represented with a point. The stars represent the average across all 34 trajectories of normalized standard deviations and correlation coefficients for each product. Note the non-linearity of scales used for correlation and skill score contours in the figure. The Taylor diagrams are modified versions of an initial Python code by Copin (2012). (e) Root mean square error (RMSE) for each velocity component. (f) Percentage of variance explained for each velocity component.


3 Results

3.1 Applying the point-wise comparison to surface drifters and ADSCs around Cape Farewell

We present here the results of the point-wise comparison between the surface drifters and ADSCs around Cape Farewell. Figure 5 shows the Taylor diagrams of the evaluated products (ADSCs and TOPAZ4). Differences between high scores are visually exaggerated with non-linear axes in order to detail differences between the products. Meridional velocities (v) are well resolved by the altimetry-derived products (Fig. 5b) with an average skill score of 0.87 computed from the average correlation coefficients of 0.75 (Geostrophy) and 0.79 (Geostrophy + Ekman) and the average normalized standard deviations of 0.84 (Geostrophy + Ekman) and 0.92 (Geostrophy). The zonal component (u) has lower scores overall (Fig. 5a) with average skill scores of 0.63 (Geostrophy) and 0.69 (Geostrophy + Ekman) computed from the average correlation coefficients of 0.62 (Geostrophy) and 0.64 (Geostrophy + Ekman) and the average normalized standard deviations of 0.60 (Geostrophy) and 0.65 (Geostrophy + Ekman). This difference in scores between the u and v components can be explained by the shape of the coast, which makes meridional velocities easier to reproduce. As the drifters follow the shelf around Cape Farewell, the meridional velocities shift from high negative velocities (southward) to positive velocities (northward). If the gridded Eulerian product simulates this general behavior, i.e., the current follows the bathymetry, then the skill scores of the meridional velocities will be higher. Another potential explanation of this difference is that sea surface height gridded at 14 used to compute geostrophic currents does not have the same spatial resolution in latitude and in longitude. A degree of latitude always corresponds to the same distance, 111 km, but a degree of longitude at 60 N corresponds to half the latitudinal resolution, or 55 km. In a 14 product, the grid cells are 27 km in latitude and 13 km in longitude. The Rossby radius in this region is on the order of 10 km; thus the 14 product approaches the Rossby radius in longitudinal span. Meridional geostrophic velocities, computed as vg=gfηx, benefit from the better longitudinal resolution. The difference between altimetry-derived velocity and a reanalysis product is very important for meridional velocities, though it should be noted that the velocities in TOPAZ4 do not experience such improvement from u to v (average skill score for the zonal component is 0.59, and for the meridional component, it is 0.60).

The along- and across-shelf velocities present lower average Taylor skill scores (Fig. 5c and d), respectively, 0.51 (Geostrophy) and 0.60 (Geostrophy + Ekman) and 0.60 (Geostrophy) and 0.64 (Geostrophy + Ekman). This decrease from u and v to along-shelf and across-shelf velocities is due to the fact that the coordinate system follows the bathymetry and thus only considers variability distinct from the bathymetric contours. The across-shelf Taylor diagram shows that Geostrophy + Ekman performs better than other products, especially compared to Geostrophy only.

The RMSE calculated for each velocity component shows very good results for Geostrophy compared to TOPAZ4 (Fig. 5e). The RMSE was lower for the across-shelf component, which was not expected given that the geostrophic flow is generally along isobaths rather than across them. However, this lower RMSE could be explained by the smaller magnitude of the across-shelf velocities, which implies a lower signal-to-error ratio. This hypothesis is reinforced by the percent variance explained (Fig. 5f) that normalizes these errors by the total amount of variance in each direction. When the amount of total variance is accounted for, the along-shelf and across-shelf metrics are similar, and the difference between the u and v fields and the along- and across-shelf velocities is accentuated. In addition, the effect of the Ekman component is more apparent in the percent-variance-explained plot. In the along-shelf direction, the addition of the Ekman component nearly doubles the percent variance explained. Similarly, in the across-shelf direction, the Ekman component makes an important contribution, but, surprisingly, the change in percent variance explained when the Ekman component is added is larger in the along-shelf direction than the across-shelf direction. Given that the flow along bathymetric contours (along-shelf direction) is largely thought to be geostrophic, while deviations to it (the across-shelf direction) would be ageostrophic, it is surprising that the Ekman component is more important in the along-shelf velocities than the across-shelf velocities. The Ekman component can influence the along-shelf currents by winds that are misaligned with the currents. This misalignment occurs frequently in the vicinity of Cape Farewell as the dominant wind patterns are no longer constrained by the topography of Greenland and the winds across the shelf. Thus there is good reason for why the Ekman contribution is large in the along-shelf direction.

The Geostrophy + Ekman data set is also available at 6 h frequency due to the higher-frequency wind products, so we also investigated the difference between the 6 h data and the 6 h drifter data. Recall that we chose the daily Geostrophy + Ekman data for a more direct comparison to the fields for Geostrophy only, which are only available at daily resolution. The 6 h Eulerian data compared to 6 h drifter data have comparable correlation coefficients but weaker normalized standard deviations than the other altimetry-derived products. This data set also has a larger RMSE (Fig. 5e), but the percent variance explained is close to the daily data. Thus the higher temporal resolution does not directly lead to better results because the 6 h drifter data contain more variability. However, the higher temporal resolution is likely more representative of the total velocity fields and thus preferable for further studies.

Some drifter tracks scored very poorly in the Taylor diagrams, reflecting the fact that it is difficult for the ADSCs to reproduce exactly the same velocities as those experienced by drifters. Indeed, some of them may for example encounter very small-scale features near the coast or due to bathymetric changes that may lead to velocities different than the main current. However, we also observed some trajectories with very high scores and a good general representation of the velocities summarized through the mean scores.

We tested the sensitivity of our results to the interpolation methods. When comparing the “nearest-neighbor” method with linear interpolation, the correlation improved on average by 0.02, and the percent variance explained increased by 2 %, but the normalized standard deviation and the Taylor skill score both decreased by 0.03 on average. Thus interpolating reduced the variance in the gridded data while improving the correlations. The gridded altimetry is already a smoothed product from the along-track altimetry. We felt that retaining the natural variance is important. Given the similarity between the results, our conclusions from this first step are robust for this choice.

Figure 6 Gridded surface currents derived from drifter data (red), the Geostrophy product (cyan), and the Geostrophy + Ekman product (blue).

Figure 7 Differences in (a, b) angle and (c, d) magnitude between surface drifter velocities and (a, c) Geostrophy and (b, d) Geostrophy + Ekman velocities. In (a) and (b), color represents the angle and the size of dots represents the number of different drifters passing through a cell. Only velocities greater than 1 cm s−1 are considered to ignore the error due to random characteristics of vectors with small magnitudes. The mean corresponds to the average angle over the dots weighted with the number of drifters passing in a cell. In (c) and (d), color represents the ratio of ADSCs over surface drifter velocities and the size of dots represents the number of different drifters passing through a cell. Score written in the figure corresponds to the mean absolute error weighted with the number of drifters passing in the cells. Differences in (e) angle and (f) magnitude between the Geostrophy and Geostrophy + Ekman velocity vectors. Note the different color scales for these two panels than the top four.

3.2 Eulerian gridding of surface drifter velocities and ADSCs

The drifters and Geostrophy and Geostrophy + Ekman products have been gridded using the method prescribed above in Sect. 2.4.2. The vector map (Fig. 6) shows that drifters mainly sampled the shelf around Cape Farewell. The general pattern of the current around the southern tip of Greenland shows a good consistency concerning the shape of the current between the ADSC and drifter-derived field. From a spatial pattern point of view, the along-shore current and the detachment of the coastal current from the coast at Cape Farewell are clearly visible and represented in all fields. To investigate this consistency more specifically, two quantities have been computed at each grid point: (1) the angle between the ADSC and the drifter-derived velocities and (2) the ratio of the ADSC over the drifter-derived velocities.

On the maps presenting angles between drifter velocities and ADSC directions (Fig. 7a and b), the angles are very small over most of the domain (smaller than 20), especially close to the shelf break, where the direction of currents is a key point to resolve the shelf–basin exchanges. Some grid points very close to the coast present larger angles (larger than 60), but only few drifters passed inside these cells, and the statistical significance is therefore smaller. The average angle over the region, weighted by the number of drifter traveling the cells, is 14.5 for Geostrophy and 12.3 for Geostrophy + Ekman.

The maps showing velocity ratios (Fig. 7c and d) depict areas of strong under- and overestimation. The regions south of Cape Farewell and Cape Desolation show an underestimation of about 50 %, indicating that the ADSC fields have trouble resolving the zonal velocity in the coastal current. Again, this result is likely due to the large meridional separation between 14 grid cells that can be easily visualized in Fig. 7. Importantly, the magnitude of this error drops significantly at the shelf break, which is the key component for representing flux between the shelf and the ocean basin. So while the coastal current is likely not well resolved south of Cape Farewell, the shelf break jet remains present. The mean absolute errors observed are 37.2 % for Geostrophy and 35.6 % for Geostrophy + Ekman. These percentages capture the magnitude of the errors, both in underestimation and in overestimation.

The difference between Geostrophy and Geostrophy + Ekman is small in these averaged fields. Indeed, the mean angle between the vectors is around 5.2 (Fig. 7e), and the mean difference between velocities (Fig. 7f) is about 7.1 % of Geostrophy + Ekman velocities. The main area where directions are different is on the west Greenland shelf, around Cape Desolation, with differences up to 12–15. Alternatively, we observe two hot spots where the magnitudes of velocities are different: southeast of Cape Farewell, where the Ekman contribution increases the magnitudes, and northwest of Cape Farewell, around Julianehaab Bight, where the Ekman contribution decreases them.

The ADSCs seem able to reproduce the current’s direction with high accuracy, but the magnitude is underestimated. We suggest that this underestimation is due to the gridding of the along-track altimetry data that smoothes the velocities and likely underestimates their true variance. The Ekman contribution seems to be a small contribution to the total velocity but improves both the direction and magnitude consistently in both the along-shelf and across-shelf directions.

We compare 112 to 14 data to evaluate the differences due to a change in resolution (Fig. 8). The 112 velocity grid shows more details and small-scale features (Fig. 8a). This product identifies the very strong velocity regions on the shelf south of Cape Farewell and on the shelf break south of Cape Desolation. Current direction and magnitude seem consistent, and the mean kinetic energy (Fig. 8b, d) and eddy kinetic energy (Fig. 8c, e) fields show close values in general and in the energy hot spots as well. The 14 resolution therefore gets the right values at the right locations. In general, this exercise in testing various spatial resolutions of the gridded surface drifter velocity fields has shown us that future improvements to the gridded altimetry record from higher-resolution products and the SWOT mission will likely improve the comparisons to the surface drifter velocities and allow for more detailed studies of the circulation structure, but, surprisingly, the 14 gridded fields capture the majority of the circulation features important to the shelf–basin exchange in this region.

Figure 8 Comparison of gridded flows derived from drifter trajectories with varying resolution. (a) Comparison of 14 (blue) and 112 (black) vectors. (b, d) Mean kinetic energy (MKE) maps obtained on (b) 14 and (d) 112 grids. (c, e) Eddy kinetic energy (EKE) maps obtained on (c) 14 and (e) 112 grids. Values provided in (b–e) are the means over the color-coded boxes.

In this section, we presented a methodology to compute mean flows from ADSCs that closely simulates the progressive sampling of the domain by the surface drifters so that the Eulerian fields are averaged in the same spatiotemporal domain as the drifters. We also evaluated how sensitive this result is using a simple time average of the ADSCs. Figure S3 in the Supplement compares the Geostrophy and Geostrophy + Ekman velocity fields computed with these two different methods. They yield very similar fields, except in the region near the shelf break south of Cape Desolation. In this region, Geostrophy + Ekman computed with the “drifter-following” methodology clearly differs from the other fields. In our example, the methodology does not change much for the data set of Geostrophy only, and the time-averaged Geostrophy + Ekman field is very close to those Geostrophy fields. This could indicate that the Ekman contribution diminishes when averaged over a long enough period of time in this region where the synoptic storm timescale is around 3–5 d. The short-timescale wind events likely play a role in the cross-shelf exchange as highlighted by the vector field obtained with the drifter-following method.

3.3 Observed and synthetic trajectories around Cape Farewell

ADSC velocities are now evaluated in the Lagrangian frame by applying the methodology described in Sect. 2.3. The sensitivity of the skill scores to the duration of the particles' advection has been tested by repeating the experiment with various durations (from 2 to 14 d). We chose to use 3 d as proposed by Liu and Weisberg (2011) so that we could compare our values to theirs, as well as to avoid too short experiments which could be insignificant and too long ones where local information would be lost. Liu et al. (2014) applied this methodology to the evaluation of various remote-sensing products in the Gulf of Mexico and found a mean skill score of 0.50 in the open ocean and 0.41 on the shelf. The ADSC data they used were a combination AVISO (Archiving, Validation and Interpretation of Oceanographic Satellite data) (delayed time), CNES-CLS09 MDT (Rio et al.2011), and Ekman, which corresponds to the Geostrophy + Ekman product evaluated in our work. When the Ekman component is removed from their ADSC product, the open-ocean skill score does not change, but the shelf one decreases to 0.35. We can use their results as a reference to interpret the results obtained in the present work.

Here, the vast majority of our drifters traveled on the shelf, and we calculate a mean skill score of 0.47 for Geostrophy only and 0.50 for Geostrophy + Ekman (Fig. 9a, b). Our skill scores on the shelf with the Ekman component (0.50) exceed the Liu et al. (2014) skill scores on the shelf (0.41) and are comparable to those obtained in the open ocean. Without Ekman, the mean score obtained with data mainly located on the shelf is 0.47. There is also a clear spatial signal in the skill scores, with an area of lower skill scores (red hatch) around Julianehaab Bight. These seem to correspond to the area of slowly eddying shelf flow described by Duyck and De Jong (2021). The average associated skill score is 0.38 for Geostrophy and 0.39 for Geostrophy + Ekman. This turbulent area, which is supposed to be more difficult to be reproduce because of its particular dynamics, shows results almost the same as Liu et al. (2014) for the Ekman contribution on the shelf and even better results compared to Geostrophy only. The rest of the shelf shows particularly high mean skill scores of 0.55 for Geostrophy only and 0.58 for the Ekman contribution. The results obtained here with the remote-sensing products are thus particularly good, especially considering the Ekman contribution. The impact of this contribution on the skill score is investigated by computing the difference between the score obtained with both products (Fig. 9c). The Ekman contribution improves the main skill score by 0.03, and large improvements particularly appear on the east Greenland shelf for trajectories not located against the coast and south of Cape Desolation close to the coast.

Figure 9(a, b, e, f) Synthetic trajectories evaluation for (a, e) Geostrophy and (b, f) Geostrophy + Ekman. The dots represent the daily locations of drifters, and the colors represent the score obtained for the trajectory computed with the 3 following days. An area with lower scores is hatched in red, and the rest of the shelf with higher skill scores is hatched in green. Confidence intervals are computed using a bootstrap methodology. (a–d) Results computed with the original (Liu et al.2014) methodology, which sets negative values to 0. (e–h) Results based on an alternative metric proposed in Révelard et al. (2021), which conserves the negative values. (c, g) Skill score difference due to the Ekman contribution, for the (c) original and the (g) alternative metrics. Blue dots denote increases in the skill score when the Ekman velocities are included, and red dots indicate decreases. (c, g) Histograms showing the distribution of the skill score for the different subregions (represented by colors) and using Geostrophy (colored bars) or Geostrophy + Ekman (hatched bars), for the (d) original and the (h) alternative metrics. C.I.: confidence interval.

We extended the analysis by reproducing the results with the alternative metric (Fig. 9d, e) suggested by Révelard et al. (2021), which conserved the negative values instead of replacing them with 0. The mean results decreased from 0.47 and 0.5, respectively, for Geostrophy and Geostrophy + Ekman to 0.33 and 0.36, which remain good compared to skill scores obtained by Liu et al. (2014) on the shelf (0.35 and 0.41) without accounting for negative values. The skill scores in the large shelf area, highlighted by green hatches, remain very good and only decrease from 0.55 and 0.58 to 0.53 and 0.55 with this new metric. The results in the red shelf area present larger decreases from 0.38 and 0.39, respectively, for Geostrophy and Geostrophy + Ekman to 0.25 and 0.27, but it remains on the order of the results of Liu et al. (2014) without negative values. Values using this metric reported in Révelard et al. (2021) from the Ibiza channel were considerably lower, with some regions hitting −0.6. The proportion of negative scores, proposed in Révelard et al. (2021), is also a very informative metric to understand the full picture. In our case for Geostrophy (Geostrophy + Ekman), we obtained negative scores of 13.1 % (12.2 %) in the entire region, 7.3 % (7.1 %) in the red area, and 3.4 % (3.6 %) in the green one. This low proportion of negative values explained the relatively small change in the skill scores when accounting for negative values, especially in the area of good consistency highlighted in green. The impact of the Ekman contribution for this metric (Fig. 9f) remains very close to the result obtained using the first metric (Fig. 9c).

Finally, it is interesting to look at the distribution of scores for the different products and subregions and with the two metrics (Fig. 9g, h). Using the first metric (Fig. 9g) we observed a peak between 0 and 0.1 that is not present with the second metric (Fig. 9h) and simply represents the zeroing of all negative scores. In addition, we can observe that the scores in the red subregion are more spread compared to scores in the green subregion, especially considering the Révelard et al. (2021) metric. We can also observe that for the entire domain and the green subregion the number of scores smaller than 0.7 decreases and the number of scores between 0.7 and 0.8, and to a lesser extent between 0.8 and 0.9, strongly increases. The low average score in the red subregion is explained by the large variety of scores obtained with a combination of good and bad scores. Conversely, the good mean score in the green subregion is explained by the presence of almost only good scores.

From those results, we conclude that ADSCs seem able to reproduce the trajectories of surface drifters in the region of Cape Farewell and that they can reproduce trajectories in the open ocean in the Gulf of Mexico and Mediterranean. We expect the skill scores to be better for the open ocean than the shelf because the scales of motion are larger in the open ocean and the flow is more geostrophic. Thus, it is impressive that the skill scores found on the shelf in the current study are comparable to the open-ocean skill scores from Liu et al. (2014) and Révelard et al. (2021). Furthermore, our study is at higher latitudes, where the Rossby radius of deformation is smaller, and thus the comparison of skill scores would favor more southerly latitudes. However, the tracks of polar-orbiting altimeters converge at high latitudes, so there is better along-track coverage, and we are also using an updated MDT data set in our analysis (2018 version of the CNES CLS MDT compared to the 2009 version), both of which could contribute to this favorable comparison of our results.

4 Conclusions

The combination of all three steps demonstrates that the altimetry-derived surface currents are largely capable of recovering the spatial structure of the flow field on the south Greenland shelf and can mimic the Lagrangian nature of the flow as observed from surface drifters. This good agreement is especially strong for the meridional velocities, likely due to the strong bathymetric constraints and the meridional orientation of the shelf in the area, as well as the higher spatial resolution in the zonal dimension of the gridded altimetry product. The Taylor skill scores drop for the along-shelf and across-shelf velocities, but they remain high considering that the shelf coordinate system accounts for the bathymetric steering of the flow. So, the skill scores for these velocity components essentially estimate how well the flow field is characterized beyond the bathymetric control, and thus the skill scores are quite high given that assumption. ADSCs show a particular ability to reproduce the direction of currents around Cape Farewell with errors in directions around 12 to 14 in average, depending on the product used.

Overall, the addition of the Ekman velocities to the Geostrophy product improved the comparisons, though not by as much as originally hypothesized and not for specific velocity components. It is likely that during high wind events this Ekman component is more noteworthy. Though our drifter observations were limited to a relatively short period in August and September of 2021, when winds are at their climatological minima, we note that the conditions experienced during this period in August and September of 2021 were noteworthy for their strong winds for this time of year. The drifters directly sampled upwelling and downwelling on the east Greenland shelf, as well as upwelling on the west Greenland shelf. Finally, the simulations of particle trajectories highlight the particularly good ability of ADSCs to mimic the displacement of surface drifters in the region and by extension the displacement of water masses.

The main sources of error between the ADSC and the surface-derived velocities lie in the magnitude of velocities, especially the zonal velocities around the southernmost tip of Greenland. Here, the direction of the flow is well constrained, but the magnitude is about half of the observed velocity. This result is concerning, though the direction of the vectors remains high on the shelf, and the direction and magnitudes are well resolved at the shelf break. As the shelf break circulation is critical for evaluating shelf–basin fluxes, we remain confident that the ADSCs are doing a good job of tracking this exchange. Interestingly, the gridded products (step no. 2 in the methodology) seem to imply that the ADSCs would be better suited to tracking the shelf–basin exchange than the exchange from the coastal current to the shelf break currents. But the Lagrangian simulations (step no. 3 in the methodology) imply that the ADSCs are quite good on the southeast Greenland shelf, where we observe quite a bit of coastal current–shelf break exchange, especially under the initial upwelling-favorable winds. Thus, we conclude that the ADSCs are capable of tracking the shelf–basin exchange, and one should specify exactly what is of interest in the coastal current–shelf break exchange prior to using the ADSC fields. Though our initial “null hypothesis” was that the ADSCs would have trouble resolving the shelf circulation, the common result across all three steps of this proposed methodology was that the ADSCs capture critical components of the circulation.

A large caveat of our results here has been that these results are specific to a roughly 2-week period in August–September 2021 and may or may not be representative of the longer-term variability. To address this concern, we used a database of 34 drifters from the 6 h data set (Lumpkin and Centurioni2019) derived from Global Drifter Program (GDP) data that crossed onto the shelf from 1993–2021 (Supplement Sect. S5; Figs. S4, S5, S6). The results from extending the temporal scope show a good coherence between surface drifters' trajectories and altimetry-derived surface currents from 1996 to 2020 in all seasons investigated, leading us to believe that this good correspondence is not specific to our brief study period, albeit with limited data to test. How far back in time one can reliably reconstruct the shelf circulation with ADSCs, specifically as the number of altimeters decreases significantly prior to 2000, remains to be answered.

The results of this assessment pave the way for a long-term study of currents around the southern tip of Greenland based on satellite observations. They have the potential to improve our understanding of freshwater exchange between the shelf and the ocean basin by adding 30 years of observations to the results of modeling work and possibly eliminating some model disagreements.

Code availability

The Python codes to implement the presented methodology are available on the ADSC-SVP-Comparison repository (, Coquereau2023).

Data availability

The surface geostrophic-current data set Global Ocean Gridded L4 Sea Surface Heights And Derived Variables Nrt is publicly available on the CMEMS website at (Global Ocean Gridded L4 Sea Surface Heights And Derived Variables Nrt2022). The version including the Ekman contribution, named Global Total Surface and 15 m Current (COPERNICUS-GLOBCURRENT) from Altimetric Geostrophic Current and Modeled Ekman Current Processing, is available at (Global Total Surface and 15 m Current from Altimetric Geostrophic Current and Modeled Ekman Current Processing2022). Finally, the CMEMS website also provides a TOPAZ4 (Sakov et al.2012) open-access reanalysis at (Arctic Ocean Physics Analysis and Forecast2022). Global Drifter Program quality-controlled 6-hour interpolated data from Lumpkin and Centurioni (2019) are publicly available at (Lumpkin and Centurioni2019). The comparison with GDP data was performed on the reprocessed versions of altimetry-derived surface currents available on CMEMS as Global Ocean Gridded L4 Sea Surface Heights And Derived Variables Reprocessed 1993 Ongoing available at (Global Ocean Gridded L4 Sea Surface Heights And Derived Variables Reprocessed 1993 Ongoing2022) and Global Total Surface and 15 m Current (COPERNICUS-GLOBCURRENT) from Altimetric Geostrophic Current and Modeled Ekman Current Reprocessing available at (, Global Total Surface and 15 m Current from Altimetric Geostrophic Current and Modeled Ekman Current Reprocessing2022). The deployed drifters from GFWE and TERIFIC have been incorporated into the aforementioned GDP data set.


The supplement related to this article is available online at:

Author contributions

AC conducted the analysis and co-wrote the manuscript. NPF devised the study and co-wrote the manuscript.

Competing interests

The contact author has declared that neither of the authors has any competing interests.


Publisher’s note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.


This work was funded by the National Science Foundation (grant no. 2047952). TERIFIC was funded by the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation program (grant agreement no. 803140). The authors would like to thank Lucas Drumetz, Eleanor Frajka-Williams, Renske Gelderloos, Christophe Maes, Bob Pickart, Florian Sévellec, and Pierre Tandeo for their input and guidance during this project. We would also like to thank the editor, Erik van Sebille, and the anonymous reviewers who helped us improve our paper considerably.

Financial support

This research was funded by the National Science Foundation Division of Ocean Sciences grant 2047952.

Review statement

This paper was edited by Erik van Sebille and reviewed by two anonymous referees.


Arbic, B. K., Scott, R. B., Chelton, D. B., Richman, J. G., and Shriver, J. F.: Effects of Stencil Width on Surface Ocean Geostrophic Velocity and Vorticity Estimation from Gridded Satellite Altimeter Data, J. Geophys. Res.-Oceans, 117,, 2012. a

Arctic Ocean Physics Analysis and Forecast: E.U. Copernicus Marine Service Information (CMEMS), Marine Data Store (MDS), [data set],, last access: 26 April 2022. a

Böning, C. W., Behrens, E., Biastoch, A., Getzlaff, K., and Bamber, J. L.: Emerging Impact of Greenland Meltwater on Deepwater Formation in the North Atlantic Ocean, Nat. Geosci., 9, 523–527,, 2016. a, b

Castelao, R. M., Luo, H., Oliver, H., Rennermalm, A. K., Tedesco, M., Bracco, A., Yager, P. L., Mote, T. L., and Medeiros, P. M.: Controls on the Transport of Meltwater From the Southern Greenland Ice Sheet in the Labrador Sea, J. Geophys. Res.-Oceans, 124, 3551–3560,, 2019. a

Copin, Y.: Taylor Diagram for Python/Matplotlib, Zenodo,, 2012. a

Coquereau, A.: Coquereau/ADSC-SVP-Comparison: ADSC-SVP-Comparison v.1.0.0 (v.1.0.0), Zenodo [code],, 2023. a

Delandmeter, P. and van Sebille, E.: The Parcels v2.0 Lagrangian framework: new field interpolation schemes, Geosci. Model Dev., 12, 3571–3584,, 2019. a

Dukhovskoy, D. S., Myers, P. G., Platov, G., Timmermans, M.-L., Curry, B., Proshutinsky, A., Bamber, J. L., Chassignet, E., Hu, X., Lee, C. M., and Somavilla, R.: Greenland Freshwater Pathways in the Sub-Arctic Seas from Model Experiments with Passive Tracers, J. Geophys. Res.-Oceans, 121, 877–907,, 2016. a

Dukhovskoy, D. S., Yashayaev, I., Proshutinsky, A., Bamber, J. L., Bashmachnikov, I. L., Chassignet, E. P., Lee, C. M., and Tedstone, A. J.: Role of Greenland Freshwater Anomaly in the Recent Freshening of the Subpolar North Atlantic, J. Geophys. Res.-Oceans, 124, 3333–3360,, 2019. a, b

Dukhovskoy, D. S., Yashayaev, I., Chassignet, E. P., Myers, P. G., Platov, G., and Proshutinsky, A.: Time Scales of the Greenland Freshwater Anomaly in the Subpolar North Atlantic, J. Climate, 34, 8971–8987,, 2021. a

Duyck, E. and De Jong, M. F.: Circulation Over the South-East Greenland Shelf and Potential for Liquid Freshwater Export: A Drifter Study, Geophys. Res. Lett., 48, e2020JB020886,, 2021. a, b

Duyck, E., Gelderloos, R., and de Jong, M. F.: Wind-Driven Freshwater Export at Cape Farewell, J. Geophys. Res.-Oceans, 127, e2021JC018309,, 2022. a, b

Elipot, S., Lumpkin, R., Perez, R. C., Lilly, J. M., Early, J. J., and Sykulski, A. M.: A Global Surface Drifter Data Set at Hourly Resolution, J. Geophys. Res.-Oceans, 121, 2937–2966,, 2016. a

Faugère, Y., Taburet, G., Ballarotta, M., Pujol, I., Legeais, J. F., Maillard, G., Durand, C., Dagneau, Q., Lievin, M., Sanchez Roman, A., and Dibarboure, G.: DUACS DT2021: 28 years of reprocessed sea level altimetry products, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-7479,, 2022. a

Fichefet, T., Poncin, C., Goosse, H., Huybrechts, P., Janssens, I., and Le Treut, H.: Implications of Changes in Freshwater Flux from the Greenland Ice Sheet for the Climate of the 21st Century, Geophys. Res. Lett., 30,, 2003. a

Foukal, N. P., Gelderloos, R., and Pickart, R. S.: A Continuous Pathway for Fresh Water along the East Greenland Shelf, Sci. Adv., 6, eabc4254,, 2020. a

Garcia-Quintana, Y., Courtois, P., Hu, X., Pennelly, C., Kieke, D., and Myers, P. G.: Sensitivity of Labrador Sea Water Formation to Changes in Model Resolution, Atmospheric Forcing, and Freshwater Input, J. Geophys. Res.-Oceans, 124, 2126–2152,, 2019. a, b

Gillard, L. C., Hu, X., Myers, P. G., and Bamber, J. L.: Meltwater Pathways from Marine Terminating Glaciers of the Greenland Ice Sheet, Geophys. Res. Lett., 43, 10873–10882,, 2016. a, b

Glikson, A. Y.: North Atlantic and Sub-Antarctic Ocean Temperatures: Possible Onset of a Transient Stadial Cooling Stage, Clim. Change, 155, 311–321,, 2019. a

Global Ocean Gridded L4 Sea Surface Heights And Derived Variables Nrt: E.U. Copernicus Marine Service Information (CMEMS), Marine Data Store (MDS) [data set],, last access: 19 April 2022. a

Global Ocean Gridded L4 Sea Surface Heights And Derived Variables Reprocessed 1993 Ongoing: E.U. Copernicus Marine Service Information (CMEMS), Marine Data Store (MDS) [data set],, last access: 15 December 2022. a

Global Total Surface and 15 m Current (COPERNICUS-GLOBCURRENT) from Altimetric Geostrophic Current and Modeled Ekman Current Processing: E.U. Copernicus Marine Service Information (CMEMS), Marine Data Store (MDS) [data set],, last access: 19 April 2022. a

Global Total Surface and 15 m Current (COPERNICUS-GLOBCURRENT) from Altimetric Geostrophic Current and Modeled Ekman Current Reprocessing: E.U. Copernicus Marine Service Information (CMEMS), Marine Data Store (MDS) [data set],, last access: 15 December 2022. a

Gou, R., Feucher, C., Pennelly, C., and Myers, P. G.: Seasonal Cycle of the Coastal West Greenland Current System Between Cape Farewell and Cape Desolation From a Very High-Resolution Numerical Model, J. Geophys. Res.-Oceans, 126, e2020JC017017,, 2021. a

Gou, R., Pennelly, C., and Myers, P. G.: The Changing Behavior of the West Greenland Current System in a Very High-Resolution Model, J. Geophys. Res.-Oceans, 127, e2022JC018404,, 2022. a

Haine, T. W. N., Gelderloos, R., Jimenez-Urias, M. A., Siddiqui, A. H., Lemson, G., Medvedev, D., Szalay, A., Abernathey, R. P., Almansi, M., and Hill, C. N.: Is Computational Oceanography Coming of Age?, B. Am. Meteorol. Soc., 102, 1481–1493,, 2021. a

Hansen, J., Sato, M., Hearty, P., Ruedy, R., Kelley, M., Masson-Delmotte, V., Russell, G., Tselioudis, G., Cao, J., Rignot, E., Velicogna, I., Tormey, B., Donovan, B., Kandiano, E., von Schuckmann, K., Kharecha, P., Legrande, A. N., Bauer, M., and Lo, K.-W.: Ice melt, sea level rise and superstorms: evidence from paleoclimate data, climate modeling, and modern observations that 2 °C global warming could be dangerous, Atmos. Chem. Phys., 16, 3761–3812,, 2016. a

Håvik, L., Pickart, R. S., Våge, K., Torres, D., Thurnherr, A. M., Beszczynska-Möller, A., Walczowski, W., and von Appen, W.-J.: Evolution of the East Greenland Current from Fram Strait to Denmark Strait: Synoptic Measurements from Summer 2012, J. Geophys. Res.-Oceans, 122, 1974–1994,, 2017. a

Hersbach, H., Bell, B., Berrisford, P., Hirahara, S., Horányi, A., Muñoz-Sabater, J., Nicolas, J., Peubey, C., Radu, R., Schepers, D., Simmons, A., Soci, C., Abdalla, S., Abellan, X., Balsamo, G., Bechtold, P., Biavati, G., Bidlot, J., Bonavita, M., De Chiara, G., Dahlgren, P., Dee, D., Diamantakis, M., Dragani, R., Flemming, J., Forbes, R., Fuentes, M., Geer, A., Haimberger, L., Healy, S., Hogan, R. J., Hólm, E., Janisková, M., Keeley, S., Laloyaux, P., Lopez, P., Lupu, C., Radnoti, G., de Rosnay, P., Rozum, I., Vamborg, F., Villaume, S., and Thépaut, J.-N.: The ERA5 Global Reanalysis, Q. J. Roy. Meteor. Soc., 146, 1999–2049,, 2020. a, b

Johnson, G. C., Hosoda, S., Jayne, S. R., Oke, P. R., Riser, S. C., Roemmich, D., Suga, T., Thierry, V., Wijffels, S. E., and Xu, J.: Argo – Two Decades: Global Oceanography, Revolutionized, Annu. Rev. Mar. Sci., 14, 379–403,, 2022. a

Koszalka, I. M. and LaCasce, J. H.: Lagrangian Analysis by Clustering, Ocean Dynam., 60, 957–972,, 2010. a

Kundu, P. K.: Ekman Veering Observed near the Ocean Bottom, J. Phys. Ocean., 6, 238–242,<0238:EVONTO>2.0.CO;2, 1976. a

LaCasce, J. H.: Statistics from Lagrangian Observations, Prog. Oceanogr., 77, 1–29,, 2008. a

Liu, Y. and Weisberg, R. H.: Evaluation of Trajectory Modeling in Different Dynamic Regions Using Normalized Cumulative Lagrangian Separation, J. Geophys. Res.-Oceans, 116, C09013,, 2011. a, b, c, d, e, f, g, h

Liu, Y., Weisberg, R. H., Vignudelli, S., and Mitchum, G. T.: Evaluation of Altimetry-Derived Surface Current Products Using Lagrangian Drifter Trajectories in the Eastern Gulf of Mexico, J. Geophys. Res.-Oceans, 119, 2827–2842,, 2014. a, b, c, d, e, f, g, h

Lumpkin, R. and Centurioni, L.: Global Drifter Program qualitycontrolled 6-hour interpolated data from ocean surface drifting buoys. NOAA National Centers for Environmental Information, [data set],, 2019. a, b

Lumpkin, R. and Johnson, G. C.: Global Ocean Surface Velocities from Drifters: Mean, Variance, El Niño–Southern Oscillation Response, and Seasonal Cycle, J. Geophys. Res.-Oceans, 118, 2992–3006,, 2013. a

Lumpkin, R. and Centurioni, L.: Global Drifter Program quality-controlled 6-hour interpolated data from ocean surface drifting buoys. NOAA National Centers for Environmental Information, [data set],, 2019. a

Luo, H., Castelao, R. M., Rennermalm, A. K., Tedesco, M., Bracco, A., Yager, P. L., and Mote, T. L.: Oceanic Transport of Surface Meltwater from the Southern Greenland Ice Sheet, Nat. Geosci., 9, 528–532,, 2016. a, b

Marsh, R., Desbruyères, D., Bamber, J. L., de Cuevas, B. A., Coward, A. C., and Aksenov, Y.: Short-term impacts of enhanced Greenland freshwater fluxes in an eddy-permitting ocean model, Ocean Sci., 6, 749–760,, 2010. a

Morrison, T. J., Dukhovskoy, D. S., McClean, J. L., Gille, S. T., and Chassignet, E. P.: Mechanisms of Heat Flux Across the Southern Greenland Continental Shelf in 1/10 and 1/12Ocean/Sea Ice Simulations, J. Geophys. Res.-Oceans, 128, e2022JC019021,, 2023. a

Mulet, S., Rio, M.-H., Etienne, H., Artana, C., Cancet, M., Dibarboure, G., Feng, H., Husson, R., Picot, N., Provost, C., and Strub, P. T.: The new CNES-CLS18 global mean dynamic topography, Ocean Sci., 17, 789–808,, 2021. a, b

Oliver, H., Luo, H., Castelao, R. M., van Dijken, G. L., Mattingly, K. S., Rosen, J. J., Mote, T. L., Arrigo, K. R., Rennermalm, Å. K., Tedesco, M., and Yager, P. L.: Exploring the Potential Impact of Greenland Meltwater on Stratification, Photosynthetically Active Radiation, and Primary Production in the Labrador Sea, J. Geophys. Res.-Oceans, 123, 2570–2591,, 2018. a

Pacini, A. and Pickart, R. S.: Wind-Forced Upwelling Along the West Greenland Shelfbreak: Implications for Labrador Sea Water Formation, J. Geophys. Res.-Oceans, 128, e2022JC018952,, 2023. a

Pennelly, C., Hu, X., and Myers, P. G.: Cross-Isobath Freshwater Exchange Within the North Atlantic Subpolar Gyre, J. Geophys. Res.-Oceans, 124, 6831–6853,, 2019. a, b

Poulain, P.-M., Gerin, R., Mauri, E., and Pennel, R.: Wind Effects on Drogued and Undrogued Drifters in the Eastern Mediterranean, J. Atmos. Ocean. Technol., 26, 1144–1156,, 2009. a

Pujol, M.-I., Faugère, Y., Taburet, G., Dupuy, S., Pelloquin, C., Ablain, M., and Picot, N.: DUACS DT2014: the new multi-mission altimeter data set reprocessed over 20 years, Ocean Sci., 12, 1067–1090,, 2016. a, b

Pujol, M.-I., Dupuy, S., Vergara, O., Sánchez Román, A., Faugère, Y., Prandi, P., Dabat, M.-L., Dagneaux, Q., Lievin, M., Cadier, E., Dibarboure, G., and Picot, N.: Refining the Resolution of DUACS Along-Track Level-3 Sea Level Altimetry Products, Remote Sens., 15, 793,, 2023.  a

Rahmstorf, S., Box, J. E., Feulner, G., Mann, M. E., Robinson, A., Rutherford, S., and Schaffernicht, E. J.: Exceptional Twentieth-Century Slowdown in Atlantic Ocean Overturning Circulation, Nat. Clim. Change, 5, 475–480,, 2015. a

Révelard, A., Reyes, E., Mourre, B., Hernández-Carrasco, I., Rubio, A., Lorente, P., Fernández, C. D. L., Mader, J., Álvarez-Fanjul, E., and Tintoré, J.: Sensitivity of Skill Score Metric to Validate Lagrangian Simulations in Coastal Areas: Recommendations for Search and Rescue Applications, Front. Mar. Sci., 8, 2296–7745,, 2021. a, b, c, d, e, f, g

Rio, M.-H.: Use of Altimeter and Wind Data to Detect the Anomalous Loss of SVP-Type Drifter's Drogue, J. Atmos. Ocean. Tech., 29, 1663–1674,, 2012. a

Rio, M.-H. and Hernandez, F.: High-Frequency Response of Wind-Driven Currents Measured by Drifting Buoys and Altimetry over the World Ocean, J. Geophys. Res.-Oceans, 108,, 2003. a

Rio, M. H. and Santoleri, R.: Improved Global Surface Currents from the Merging of Altimetry and Sea Surface Temperature Data, Remote Sens. Environ., 216, 770–785,, 2018. a

Rio, M. H., Guinehut, S., and Larnicol, G.: New CNES-CLS09 Global Mean Dynamic Topography Computed from the Combination of GRACE Data, Altimetry, and in Situ Measurements, J. Geophys. Res.-Oceans, 116, C07018,, 2011. a, b, c, d, e

Rio, M.-H., Mulet, S., and Picot, N.: Beyond GOCE for the Ocean Circulation Estimate: Synergetic Use of Altimetry, Gravimetry, and in Situ Data Provides New Insight into Geostrophic and Ekman Currents, Geophys. Res. Lett., 41, 8918–8925,, 2014. a, b

Sakov, P., Counillon, F., Bertino, L., Lisæter, K. A., Oke, P. R., and Korablev, A.: TOPAZ4: an ocean-sea ice data assimilation system for the North Atlantic and Arctic, Ocean Sci., 8, 633–656,, 2012. a, b

Schulze Chretien, L. M. and Frajka-Williams, E.: Wind-driven transport of fresh shelf water into the upper 30 m of the Labrador Sea, Ocean Sci., 14, 1247–1264,, 2018. a, b

Tagklis, F., Bracco, A., Ito, T., and Castelao, R. M.: Submesoscale Modulation of Deep Water Formation in the Labrador Sea, Sci. Rep., 10, 17489,, 2020. a, b, c

Taylor, K. E.: Summarizing Multiple Aspects of Model Performance in a Single Diagram, J. Geophys. Res.-Atmos., 106, 7183–7192,, 2001. a

Short summary
Understanding meltwater circulation around Greenland is crucial as it could influence climate variability but difficult as data are scarce. Here, we use 34 surface drifters to evaluate satellite-derived surface currents and show that satellite data recover the general structure of the flow and can recreate the pathways of particles around the southern tip of Greenland. This result permits a wide range of work to proceed looking at long-term changes in the circulation of the region since 1993.