Twenty-Seven Years of Scatterometer Surface Wind Analysis over Eastern Boundary Upwelling Systems

: More than twelve satellite scatterometers have operated since 1992 through the present, providing the main source of surface wind vector observations over global oceans. In this study, these scatterometer winds are used in combination with radiometers and synthetic aperture radars (SAR) for the better determination and characterization of high spatial and temporal resolution of regional surface wind parameters, including wind speed and direction, wind stress components, wind stress curl, and divergence. In this paper, a 27-year-long (1992–2018) 6-h satellite wind analysis with a spatial resolution of 0.125 ◦ in latitude and longitude is calculated using spatial structure functions derived from high-resolution SAR data. The main objective is to improve regional winds over three major upwelling regions (the Canary, Benguela, and California regions) through the use of accurate and homogenized wind observations and region-speciﬁc spatial and temporal wind variation structure functions derived from buoy and SAR data. The long time series of satellite wind analysis over the California upwelling, where a signiﬁcant number of moorings is available, are used for assessing the accuracy of the analysis. The latter is close to scatterometer wind retrieval accuracy. This assessment shows that the root mean square difference between collocated 6-h satellite wind analysis and buoys is lower than 1.50 and 1.80 m s − 1 for offshore and nearshore locations, respectively. The temporal correlation between buoy and satellite analysis winds exceeds 0.90. The analysis accuracy is lower for 1992–1999 when satellite winds were mostly retrieved from ERS-1 and/or ERS-2 scatterometers. To further assess the improvement brought by this new wind analysis, its data and data from three independent products (ERA5, CMEMS, and CCMP) are compared with purely scatterometer winds over the Canary and Benguela regions. Even though the four products are generally similar, the new satellite analysis shows signiﬁcant improvements, particularly in the upwelling areas.


Introduction
Some of the most productive marine areas are regions of eastern boundary upwelling systems where equatorward alongshore winds produce coastal upwelling due to offshore Ekman transport. Several previous studies have characterized their high variable space and time dynamics in attribution to the mean and variability characteristics of forcing winds. The wind forcing data that are needed to drive the coastal upwelling can be derived from multiple sources including blended satellite-model analyses (e.g., [1,2], atmospheric re-analyses (e.g., [3,4]), or high-resolution regional atmospheric models (such as Weather Research Forecasting, WRF [5]). Ground truth data for the above spatially resolving winds are normally provided by ocean mooring time series (e.g., [6]), ship observations (International Comprehensive Ocean-Atmosphere Data Set, ICOADS [7]), and/or from purely remote sensing winds (e.g., [8]).
The right representation of mean and transient wind characteristics is crucial for the correct upwelling representation. In particular, focusing on Bay St Hélène regional upwelling, [9] have emphasized the impact of proper spatial and especially temporal resolution of forcing winds for the correct representation of anomalous sea surface temperature (SST). Numerical simulations have shown the importance of the combined coastal and open ocean upwelling response [10]. The latter is driven by the surface wind stress curl, which requires high spatial resolution winds to capture the spatial wind inhomogeneity. Some important upwelling regions are located in regions of complex wind-coastal-orography interactions, including the formation of coastal wind jets. As shown by [11][12][13], the correct representation of the spatial structure and magnitude of these coastal wind jets is crucial for reliable simulations of upwelling and related eastern boundary current transport. In turn, strong SST contrasts in upwelling regions lead to important air-sea interactions whose proper representation again requires high spatial resolution. In particular, [1] have shown that the use of 0.25 • satellite winds improves the expected relationship between crosswind SST and wind stress curl (e.g., [14,15]) in comparison with coarser winds. Moreover, through twin numerical simulations, [2] have shown the need for high spatial and temporal surface wind resolution for better characterization of transient upwelling features known for occurring in the Canary and Benguela regions.
Research focusing on upwelling systems requires an improved wind forcing with better accuracy and regional representation of spatial and temporal variability. This study aims at meeting the above scientific requirements through the determination of accurate surface wind analyses at regional scales. It focuses on developing high-resolution winds based on the use of various satellite observations including measurements from radars (scatterometers, synthetic aperture radar (SAR)) and radiometers. This suite of remote sensing data is used for developing a new surface wind analysis that provides the mean structure and accompanying variability of wind vector, stress, curl, and divergence, and it is specifically designed for coastal upwelling regions.
This study represents our continuing efforts in improving surface wind analyses. It focuses on three upselling system regions: Canary, Benguela, and California ( Figure A1). It follows the previous global blended wind analysis of [16]. Although this previous analysis exhibits good accuracy, it needs improvements at regional scales and especially in high variable upwelling regions such as the Canary and Benguela eastern boundary regions. Consequently, this study aims at the improvement of surface wind analysis in terms of spatial and temporal resolutions, as well as its accuracy. The resulting wind analysis is expected to be suitable for the investigation of upwelling system dynamics and its modeling. To our best knowledge, such regional remote sensing-based surface wind analysis was not available before.
This paper is organized as follows. The remote sensing and numerical model data used in this study are described in Section 2. The method of determination of regular space and time surface wind analysis and the assessment of its accuracy characteristics are provided in Section 3. Section 4 focuses on the assessing of improvements of this new regional analysis through comparisons with purely scatterometer wind retrievals.

In Situ Data
In situ wind measurements from National Data Buoy Center (NDBC) buoys are used as a reference for assessing the accuracy of remote sensing data and analyses. This buoy Remote Sens. 2021, 13, 940 3 of 28 network, among others, provides the longest time series of air-sea interface parameters in various atmospheric and oceanic conditions (https://www.ndbc.noaa.gov/ (accessed on 2 March 2021); see Figure 1 for the buoy geographical locations). All available buoys ( Figure 1a) are used for assessing the long-term accuracy of radar and radiometer wind retrievals, while buoys from California upwelling (Figure 1b) are also used for assessing the accuracy of this wind analysis (see Section 3).

In Situ Data
In situ wind measurements from National Data Buoy Center (NDBC) buoys are used as a reference for assessing the accuracy of remote sensing data and analyses. This buoy network, among others, provides the longest time series of air-sea interface parameters in various atmospheric and oceanic conditions (https://www.ndbc.noaa.gov/ (accessed on 2 March 2021); see Figure 1 for the buoy geographical locations). All available buoys ( Figure 1a) are used for assessing the long-term accuracy of radar and radiometer wind retrievals, while buoys from California upwelling (Figure 1b) are also used for assessing the accuracy of this wind analysis (see Section 3, 2 m ≤ ℎ ≤ = ℎ = 10). decimation, while earlier buoy data are available only hourly. In addition to wind speed and direction measurements, several other atmospheric and oceanic parameters (e.g., SST, near-surface air temperature, Ta, relative humidity, Qa, and sea state) required for the bulk atmospheric boundary layer characterization are also available. These data are used for converting buoy wind measurements at z = h to equivalent neutral wind at z = 10 m (ENW) based on the COARE3.0 algorithm [17].

Remote Sensing Data
Since 1992, several satellite radars and radiometers have provided measurements of ocean surface parameters and near-surface atmosphere parameters (Table 1). Calibrated and validated geophysical parameter retrievals from each instrument (e.g., [18][19][20][21]) are homogenized and combined into long time series (e.g., [8,22]). Calculation and accuracy determination of surface wind analyses (gridded wind in space and time) are discussed in e.g., [16,18,23,24] Bentamy et al., 1999; Tomita and Kubota. The methods of determination of retrieval accuracy of various oceanic and atmospheric parameters can be found in e.g., [25][26][27]. As shown in Table 1, most scatterometers and radiometers operated within the study period  are included in the regional surface wind analysis. The details of radar and radiometer measurements such as their measurement physics, calibration and validation, and accuracy determination are not discussed in this study given the substantial body of publications that provide that information. In this study, eight scatterometers (ERS-1, ERS-2, QuikSCAT (QSCAT), ASCAT-A, ASCAT-B, RapidScat (RSCAT), HY-2A, ScatSat (SSCAT)), one polarimetric radiometer (WindSat), two synthetic aperture radars (SAR onboard of Sentinel-1a and 1b), and 10 radiometers (SSM/I F10 through F15, SSMI F16 through F18, AMSRE, and AMSR2) are used. Scatterometers and WindSat provide wind speed (W 10 ) and direction (ϕ) as 10 m equivalent neutral wind (ENW) at swath wind vector cells (WVC) of 25 km 2 or 12.5 km 2 over the global ocean. W 10 and ϕ are also available from SAR measurements over selected geographical regions but with a greatly improved spatial resolution of 1 km 2 . In this study, only SAR wind speed and direction retrieved from backscatter coefficients acquired in the interferometric wide (SAR IW) swath mode are used. Passive microwave radiometer measurements provide only 10 m wind speed. One may notice that only remotely sensed data occurring within the study period (January 1992 to December 2018) are used. Wind speed and direction units used throughout this study are m s −1 (or m/s) and degree (or • ), respectively.
The statistics characterizing the consistency between in situ NDBC and collocated remotely sensed wind data (Table 2) is based only on NDBC offshore buoys (50 km off coastlines) and spatial and temporal collocation criteria of 25 km and 1 h, respectively. It includes the difference between time means and the root mean square difference (bias and RMSD) of buoy and satellite data (in this order), scalar (ρ) and vector (ρ 2 ) correlation coefficients, symmetrical regression coefficient slope (bs), and intercept (as). One should notice that vector correlation varies between −2 and 2. The equations allowing the calculations of the aforementioned statistical parameters are reminded hereafter.
For any two ground-truth and satellite-derived scalar variables, X and Y, the following statistical measures are defined as: where the overbar stands for arithmetical mean.
Scalar correlation, The statistical parameters used for characterizing the comparisons of two wind directions, θ and ϕ, are estimated as follows: STD of wind direction difference = sin −1 (ε)(1 + 0.1547ε 3 ), where [28]. The wind direction correlation (also named vector correction) is defined based on the method described in [29].
Σ ij is the cross-covariance two-dimensional matrix associated with wind directions θ and ϕ, being dimensions 1 and 2, respectively.
Tr indicates matrix trace. ERS-1 and ERS-2 scatterometers exhibit higher RMSD and lower correlation that result from a relatively low matchup length and coarse original spatial resolution of 50 km 2 (resampled to 25 km 2 ). The rest of the satellite winds exhibit very good agreement with buoy data without any significant bias, RMSD below 1 m s −1 and 20 • , and scalar (vector) correlation higher than 0.97 (1.80). Interestingly, the RMSD is higher for SAR wind retrievals than that for scatterometers. Even though the high spatial resolution of SAR winds (1 km 2 ) allows for almost perfect spatial collocation, the bulk of SAR data were acquired close to coasts (50 to 100 km off coasts) where wind variability is higher.

Copernicus/Marine Environment Monitoring Service (CMEMS) L4 Wind Analyses
CMEMS (https://marine.copernicus.eu/ (accessed on 2 March 2021)) provides two wind analyses estimated from remotely sensed data. They are referenced as WIND_GLO_ WIND_L4_NRT_OBSERVATIONS_012_004 (http://marine.copernicus.eu/documents/ PUM/CMEMS-WIND-PUM-012-004.pdf (accessed on 2 March 2021)) and WIND_GLO_ WIND_L4_REP_OBSERVATIONS_012_006 (http://marine.copernicus.eu/documents/PUM/ CMEMS-WIND-PUM-012-006.pdf (accessed on 2 March 2021)). The first is calculated as a near-real-time product from available scatterometer and radiometer wind retrievals in combination with the European Centre for Medium-Range Weather Forecasts (ECMWF) surface wind forecasts. The second product is calculated from reprocessed and inter-calibrated scatterometer and radiometer wind retrievals [16]. It is an offline product named hereafter as CMEMS wind analysis and used in this study for assessing the improvement of satellite wind analyses at regional scales in Canary and Benguela zones. The CMEMS wind analysis is primarily based on scatterometer winds (Section 2.2), with ancillary data from radiometer winds. In combination with ECMWF ERA-Interim atmospheric reanalysis [30], these data are used to produce wind analysis on a global regular 0.25 • × 0.25 • space/time grid at synoptic times 00 h:00, 06 h:00, 12 h:00, and 18 h:00 UTC [16]. Its accuracy determined against 6-hourly buoy winds is similar to that of purely scatterometer wind retrievals.
It should be noted that the CMEMS wind analysis is quite similar to the satellite wind analysis (Section 3) developed in this study (also referred hereafter as the IFREMER satellite wind analysis). Indeed, both are determined from almost the same remote sensing data and share the objective analysis method. The anticipated differences between the two arise from different spatial and temporal structure functions and atmospheric reanalysis used as a background.

Cross-Calibrated Multi-Platform (CCMP) Wind Analysis
In this study, version 2 of cross-calibrated multi-platform (CCMP) gridded winds is used (http://www.remss.com/measurements/ccmp/ (accessed on 2 March 2021)). They are calculated as blended winds combining satellite observations, in situ measurements, and NWP model estimates. More specifically, CCMP wind analysis is computed from QuikSCAT and ASCAT wind velocity, version 7 RSS radiometer wind speed, and the former ECMWF re-analysis (ERA-Interim) wind velocity using the variational analysis method [24]. It is available over the global ocean at synoptic times (00 h:00, 06 h:00, 12 h:00, and 18 h:00 UTC) on a 0.25 • × 0.25 • grid.

Atmospheric Reanalysis
The latest, fifth-generation ECMWF global atmospheric reanalysis named ERA5 [31] used in this study is available at https://climate.copernicus.eu/climate-reanalysis (accessed on 2 March 2021). Details of ERA5 physics and data can be found in (https: //www.ecmwf.int/en/forecasts/datasets/reanalysis-datasets/era5 (accessed on 2 March 2021)). ERA5 is a data assimilation system that combines model data and observations. Most of the available buoy (including NDBC) and satellite (almost all shown in Table 1) data are assimilated by the ERA5 reanalysis. ERA5 fields used here are provided by Copernicus services (https://climate.copernicus.eu/climate-reanalysis (accessed on 2 March 2021)). For this study, we use hourly ERA5 data on a 0.25 • grid. Note that ERA5 10m winds are provided as real winds (not ENW). The ERA5 has several improvements in comparison with the previous ECMWF reanalysis that enable the better characterization of surface wind patterns. Although the RMSD between ERA5 and collocated scatterometer winds is 20% lower than that for the ERA-Interim, it is still lacking transient wind variability crucial for regional analyses [32].

Analysis Method
The space/time gridding procedure for satellite-derived surface winds is detailed in [33]. It utilizes an objective interpolation based on the kriging technique with the external drift method. In this study, the external drift is taken from the ERA5. Briefly, the scatterometer and radiometer winds are used to estimate 6-h averaged wind speed and direction on a regular 0.125 • . The contribution of each measurement is weighted depending on its spatial and temporal separations. Weights are determined from the regional spatial and temporal structure functions (variograms) for wind speed, zonal, and meridional wind components.

Spatial Structure Functions
Spatial structure functions are derived from high-resolution SAR IW images over each upwelling region using 2017-2019 data.
Specifically, for each hour of the day and each point on a 0.125 • × 0.125 • grid, SARbased wind covariances are estimated as a function of distance δh for 1 km ≤ δh ≤ 300 km at 5 km steps. The hourly statistics is estimated only if the sampling length of SAR retrievals is significant (≥30 ). Space distribution of the spatial scales for each region reflects the nature of local air-sea-land interactions and generally aligns with local topography (Figure 2). The presence of islands (such as in the Canary region) introduces apparent inhomogeneities to spatial scale maps ( Figure 2, middle row). Local air-sea-land interactions are also reflected in the spatial distribution of wind scale patterns that tend to be aligned with regional coastal configuration and topography. In coastal areas (100 km off coastlines), spatial scales of wind speed vary between 10 and 30 km, while they generally exceed 50 km in offshore areas. Although zonal and meridional wind scales exhibit similar spatial patterns as those for wind speed, they are larger and generally exceed 30 km, except at a few locations. March 2021)). For this study, we use hourly ERA5 data on a 0.25° grid. Note that ERA5 10m winds are provided as real winds (not ENW). The ERA5 has several improvements in comparison with the previous ECMWF reanalysis that enable the better characterization of surface wind patterns. Although the RMSD between ERA5 and collocated scatterometer winds is 20% lower than that for the ERA-Interim, it is still lacking transient wind variability crucial for regional analyses [32].

Analysis Method
The space/time gridding procedure for satellite-derived surface winds is detailed in [33]. It utilizes an objective interpolation based on the kriging technique with the external drift method. In this study, the external drift is taken from the ERA5. Briefly, the scatterometer and radiometer winds are used to estimate 6-h averaged wind speed and direction on a regular 0.125°. The contribution of each measurement is weighted depending on its spatial and temporal separations. Weights are determined from the regional spatial and temporal structure functions (variograms) for wind speed, zonal, and meridional wind components.

Spatial Structure Functions
Spatial structure functions are derived from high-resolution SAR IW images over each upwelling region using 2017-2019 data.
Specifically, for each hour of the day and each point on a 0.125° × 0.125° grid, SARbased wind covariances are estimated as a function of distance δh for 1 km ≤ ℎ ≤ 300 km at 5 km steps. The hourly statistics is estimated only if the sampling length of SAR retrievals is significant (≥ 30). Space distribution of the spatial scales for each region reflects the nature of local air-sea-land interactions and generally aligns with local topography ( Figure 2). The presence of islands (such as in the Canary region) introduces apparent inhomogeneities to spatial scale maps ( Figure 2, middle row). Local air-sea-land interactions are also reflected in the spatial distribution of wind scale patterns that tend to be aligned with regional coastal configuration and topography. In coastal areas (100 km off coastlines), spatial scales of wind speed vary between 10 and 30 km, while they generally exceed 50 km in offshore areas. Although zonal and meridional wind scales exhibit similar spatial patterns as those for wind speed, they are larger and generally exceed 30 km, except at a few locations.

Temporal Structure Functions
Due to satellite orbit repeat cycles normally exceeding 4 days, the estimation of temporal scales from a particular satellite is challenging. To circumvent this limitation, we use all satellite measurements assuming that they are properly homogenized [8,16]. In this study, the temporal scales are estimated from lagged temporal correlations ( computed from all wind retrievals over each region during the Northern Hemisphere winter and summer (January and July of 2000). Figure 3-left shows that wind speed correlation remains high for ≤ 2 h and starts decreasing at ≥ 3 h (middle panels). The correlation significantly decreases in all three upwelling areas at ≥ 6 h, which is expected from the diurnal cycle of coastal winds. The sampling length associated with 5 to 6h-long lags is poor and does not allow us to assess the correlation at some coastal locations (200 km off the coast). Inferring lagged temporal correlation from polar-orbiting satellite data requires a validation that is performed against the corresponding temporal scales estimated from NDBC buoy 10-min averages. Buoy-based lagged correlation is estimated from time series collected during the same period of the year 2000 and compared with satellite-based estimates at grid points located within 12.5 km of buoy locations. This validation indicates that satellite-based temporal scale estimates are consistent with buoy-based estimates ( Figure A2).

Temporal Structure Functions
Due to satellite orbit repeat cycles normally exceeding 4 days, the estimation of temporal scales from a particular satellite is challenging. To circumvent this limitation, we use all satellite measurements assuming that they are properly homogenized [8,16]. In this study, the temporal scales are estimated from lagged temporal correlations (1 h ≤ δt ≤ h) computed from all wind retrievals over each region during the Northern Hemisphere winter and summer (January and July of 2000). Figure 3-left shows that wind speed correlation remains high for δt ≤ 2 h and starts decreasing at δt ≥ 3 h (middle panels). The correlation significantly decreases in all three upwelling areas at δt ≥ 6 h, which is expected from the diurnal cycle of coastal winds. The sampling length associated with 5 to 6h-long lags is poor and does not allow us to assess the correlation at some coastal locations (200 km off the coast).
Inferring lagged temporal correlation from polar-orbiting satellite data requires a validation that is performed against the corresponding temporal scales estimated from NDBC buoy 10-min averages. Buoy-based lagged correlation is estimated from time series collected during the same period of the year 2000 and compared with satellite-based estimates at grid points located within 12.5 km of buoy locations. This validation indicates that satellite-based temporal scale estimates are consistent with buoy-based estimates ( Figure A2).

IFREMER satellite wind analyses
As mentioned above, the objective method used for determining the gridded IFREMER satellite wind analysis is quite similar to the method used in the CMEMS analysis (Section 2.3). The main difference is in the use of spatial and temporal structure functions (variograms) specifically estimated from SAR data for each upwelling region (Sections 3.1 and 3.2) that allows characterizing smaller spatial scales. The regional IFREMER satellite wind analysis is calculated on a 0.125° grid at synoptic times (00 h:00, 06 h:00, 12 h:00, and 18 h:00 UTC) from homogenized scatterometer and radiometer winds retrievals during January 1st 1992 through December 31st 2018. ERA5 reanalysis winds are used for the external drift correction [16]. One should notice that SAR IW data are only used for the determination of spatial structure functions but not in the calculation of IFREMER wind analysis.

Accuracy of Satellite Wind Analyses
Determination of the regional satellite wind analysis accuracy requires the availability of a "ground truth" reference. It is assessed only over the California upwelling region because this region has better NDBC buoy coverage than the two other regions. Indeed, to the best of our knowledge, no nearshore or offshore moorings are available in the Benguela and Canary upwelling zones. For comparison purposes, 6-h averages are calculated from each buoy's raw wind measurements and collocated in space and time with

IFREMER Satellite Wind Analyses
As mentioned above, the objective method used for determining the gridded IFRE-MER satellite wind analysis is quite similar to the method used in the CMEMS analysis (Section 2.3). The main difference is in the use of spatial and temporal structure functions (variograms) specifically estimated from SAR data for each upwelling region (Sections 3.1 and 3.2) that allows characterizing smaller spatial scales. The regional IFRE-MER satellite wind analysis is calculated on a 0.125 • grid at synoptic times (00 h:00, 06 h:00, 12 h:00, and 18 h:00 UTC) from homogenized scatterometer and radiometer winds retrievals during January 1st 1992 through December 31st 2018. ERA5 reanalysis winds are used for the external drift correction [16]. One should notice that SAR IW data are only used for the determination of spatial structure functions but not in the calculation of IFREMER wind analysis.

Accuracy of Satellite Wind Analyses
Determination of the regional satellite wind analysis accuracy requires the availability of a "ground truth" reference. It is assessed only over the California upwelling region because this region has better NDBC buoy coverage than the two other regions. Indeed, to the best of our knowledge, no nearshore or offshore moorings are available in the Benguela and Canary upwelling zones. For comparison purposes, 6-h averages are calculated from each buoy's raw wind measurements and collocated in space and time with IFREMER satellite wind analysis and ERA5 atmospheric reanalysis synoptic rime winds. Figures 4 and 5 show scatter diagrams for 10-m wind speed (W 10 ), zonal and meridional components (U and V) for offshore (≥50 km off coastlines) and nearshore buoys. The scatter of ERA5 versus NDBC winds serves as a reference for assessing the improvement introduced by the IFREMER satellite wind analysis. The IFREMER analysis winds are highly correlated with buoy winds with somewhat larger temporal correlation coefficients for W 10 , U, and V at offshore (0.95, 0.84, and 0.93) than at nearshore buoy locations (0.92, 0.84, and 0.90). These correlation coefficients are similar to those for the ERA5. The main improvement of the IFREMER satellite wind analysis is reflected in absolute wind values. While both IFREMER and ERA5 analyses underestimate strong winds, this underestimation is lower for the IFREMER satellite wind analysis. For all winds and buoys, the bias of IFREMER and ERA5 W 10 is 0.41 and 0.91 m s −1 , respectively. For high winds (W 10 > 15 m s −1 ), these values increase to 1.66 and 3.05 m s −1 W 10 > 15. This high wind bias presented in both analyses reflects the known lack of transient variability associated with high winds in the ERA5 reanalysis (Belmont and Stoffelen, 2019) and the effect of high wind underestimation by scatterometers.

Wind Vector Issues
There are two major approaches to producing regularly gridded wind fields. They may be obtained, among many other atmospheric parameters, from atmospheric re-analyses. In this approach, the analyzed gridded fields involve a balance between atmospheric model dynamics and assimilated data. Alternative no-model analyses, such as the current satellite wind analysis, rely on measurements only, but they also heavily rely on a priori knowledge of variograms. In this section, the IFREMER satellite analysis is compared with the ERA5 reanalysis, CMEMS satellite wind analysis, and CCMP analysis (a merge of observed and model winds) against scatterometer wind observations. The above wind data are collocated in space and time with scatterometer swath wind observations from ERS-2, QuikSCAT, and ASCAT for January and July of 1996, 2006, and 2016, respectively over the Canary ( Figure 6) and Benguela (Figure 7) upwelling regions. Figures 6 and 7 show spatial distributions of wind speed RMSD, while those for each wind component are shown in the Appendix A ( Figures A4-A6). RMSD magnitude and spatial distribution The accuracy of satellite wind analysis is also investigated separately for each synoptic time (00 h:00, 06 h:00, 12 h:00, and 18 h:00 UTC). No significant changes in the correlation coefficients as well as in the RMSD are found as a function of synoptic time for offshore as well as nearshore buoys. The satellite analysis captures the diurnal cycle of winds. Due to orographic impacts, it is particularly strong (about 2 m s −1 ) in the California region within the 34 • N-35 • N, 121 • W-120 • W box (Figure 1b). This enhanced diurnal variability is confirmed by comparisons with measurements from the WMO 46011 buoy (34.88 • N, 120.87 • W) located 38 km from the coast ( Figure A3).

Wind Vector Issues
There are two major approaches to producing regularly gridded wind fields. They may be obtained, among many other atmospheric parameters, from atmospheric re-analyses. In this approach, the analyzed gridded fields involve a balance between atmospheric model dynamics and assimilated data. Alternative no-model analyses, such as the current satellite wind analysis, rely on measurements only, but they also heavily rely on a priori knowledge of variograms. In this section, the IFREMER satellite analysis is compared with the ERA5 reanalysis, CMEMS satellite wind analysis, and CCMP analysis (a merge of observed and model winds) against scatterometer wind observations. The above wind data are collocated in space and time with scatterometer swath wind observations from ERS-2, QuikSCAT, and ASCAT for January and July of 1996, 2006, and 2016, respectively over the Canary ( Figure 6) and Benguela (Figure 7) upwelling regions. Figures 6 and 7 show spatial distributions of wind speed RMSD, while those for each wind component are shown in the Appendix A ( Figures A4-A6). RMSD magnitude and spatial distribution vary among the four wind products, with the IFREMER satellite analysis having the lowest deviation from collocated scatterometer data. For each product, the highest RMSD from scatterometer data occurs in 1996 and is related to the ERS-2 scatterometer sampling scheme, which was poorer than that for QuikSCAT (in 2006), and ASCAT (in 2016). Excluding ERS-2 comparisons, the highest IFREMER wind speed RMSD values in the Canary region are located in the vicinity of the Canary Islands and result from a combination of poorer scatterometer sampling and higher wind velocity variability in space and time due to orographic effects ( Figure 6). Impacts of orographic effects are also present in the Benguela region (Figure 7) where higher RMSD values (exceeding 2 m s −1 for the meridional wind velocity, Figure A7

Wind Stress Issues
In addition to the previous comparisons mostly focusing on wind speed, the quality of wind direction is investigated from wind stress curl that is crucial for the correct representation of the Ekman pumping. If not provided by a particular analysis, wind stress is calculated using the COARE3.0 [17], with atmospheric boundary layer parameters taken from the ERA5. Figures 8 and 9 show monthly mean wind stress curl patterns estimated from the QuickSCAT scatterometer and wind analyses for the summer of 2006. Corresponding patterns for the ASCAT period (2016) are illustrated in the Appendix A (Figures A8 and A9). The strongest wind stress curl features are present along continental coasts and around the Canary and Madeira Islands (Figure 8). The latter features are associated with island shadowing effects and wind acceleration in inter-island gaps (known as the Venturi effect) that produce downstream bands of accelerated and decelerated winds behind islands. These island shadowing and acceleration wind effects are properly resolved by all analyses but the CCMP (Figure 8). The positive wind vorticity pattern along the African coast is captured by all analyses with a somewhat weaker magnitude in the CCMP. Wind stress Even though the IFREMER and CMEMS analyses are estimated based on the use of the same remote sensing data, some important differences between the two are depicted. The main improvement of the IFREMER analysis (Figures 6b,c and 7b,c) in comparison with the CMEMS analysis (Figures 6h,i and 7h,i) is found in near coastal areas. This improvement is related first to the use of the updated spatial and temporal structure functions inferred from SAR data that better account for characteristics of fine-scale regional variability (Sections 3.1 and 3.2). Furthermore, they allow for an improved selection of remotely sensed observations for wind analysis that relies on the local spatial and temporal scales.
IFREMER and CCMP wind analyses (Figures 6b-k and 7b-k) exhibit quite similar RMSD from scatterometer observations in 2006 when both analyses use QuikSCAT wind retrievals. However, the deviation from ASCAT scatterometer winds (Figures 6c-l and 7c-l) is higher for the CCMP in high wind areas south of 30 • S (Figure 7l) and coastal areas. IFREMER satellite wind speed also has a lower (by a factor of 30%) deviation from scatterometer observations than the ERA5. This is expected, because the ERA5 atmospheric model is not perfect in capturing transient wind events [32].

Wind Stress Issues
In addition to the previous comparisons mostly focusing on wind speed, the quality of wind direction is investigated from wind stress curl that is crucial for the correct representation of the Ekman pumping. If not provided by a particular analysis, wind stress is calculated using the COARE3.0 [17], with atmospheric boundary layer parameters taken from the ERA5. Figures 8 and 9 show monthly mean wind stress curl patterns estimated from the QuickSCAT scatterometer and wind analyses for the summer of 2006. Corresponding patterns for the ASCAT period (2016) are illustrated in the Appendix A (Figures A8 and A9). The strongest wind stress curl features are present along continental coasts and around the Canary and Madeira Islands (Figure 8). The latter features are associated with island shadowing effects and wind acceleration in inter-island gaps (known as the Venturi effect) that produce downstream bands of accelerated and decelerated winds behind islands. These island shadowing and acceleration wind effects are properly resolved by all analyses but the CCMP (Figure 8). The positive wind vorticity pattern along the African coast is captured by all analyses with a somewhat weaker magnitude in the CCMP. Wind stress curl along the African coast is not homogeneous and is spatially modulated by orographic interactions over major capes along the Morocco-Mauritania coast (e.g., Cape Ghir, Sim, Juby [34]). All the above coastal features amplify in summer when the northeasterly trade winds strengthen. curl along the African coast is not homogeneous and is spatially modulated by orographic interactions over major capes along the Morocco-Mauritania coast (e.g., Cape Ghir, Sim, Juby [34]). All the above coastal features amplify in summer when the northeasterly trade winds strengthen.  Southern Hemisphere trade winds affecting the Benguela upwelling region also intensify in local summer (see January 2006 in Figure 9). As in the Canary region (Figure 8), the highest wind stress curl in the Benguela region is present along coastlines downstream of main capes (Cape Columbine and Cape of Good Hope). More details on wind stress curl patterns in the Canary and Benguela can be found in . Although the four wind sources exhibit qualitatively similar wind stress curl patterns in the Benguela region, all the analyses tend to smooth and underestimate the wind stress curl magnitude in comparison with purely scatterometer-based winds, with a stronger discrepancy in the CCMP. The CCMP has the weakest wind vorticity magnitude among the products. This suggests that Benguela upwelling simulations would be biased if forced by CCMP wind stress.

Assessment of the Local Wind Patterns
IFREMER satellite wind analysis captures significant diurnal variability of winds (Figures 10 and 11) present in both the Canary and Benguela upwelling regions. Southern Hemisphere trade winds affecting the Benguela upwelling region also intensify in local summer (see January 2006 in Figure 9). As in the Canary region (Figure 8), the highest wind stress curl in the Benguela region is present along coastlines downstream of main capes (Cape Columbine and Cape of Good Hope). More details on wind stress curl patterns in the Canary and Benguela can be found in . Although the four wind sources exhibit qualitatively similar wind stress curl patterns in the Benguela region, all the analyses tend to smooth and underestimate the wind stress curl magnitude in comparison with purely scatterometer-based winds, with a stronger discrepancy in the CCMP. The CCMP has the weakest wind vorticity magnitude among the products. This suggests that Benguela upwelling simulations would be biased if forced by CCMP wind stress.

Assessment of the Local Wind Patterns
IFREMER satellite wind analysis captures significant diurnal variability of winds (Figures 10 and 11) present in both the Canary and Benguela upwelling regions. In the Canary region, the main diurnal pattern associated with wind acceleration at 6:00 UTC and deceleration at 18:00 UTC develops along the Moroccan coast. It resembles nocturnal jets developing due to diurnal changes in atmospheric boundary layer stability and height. The strongest variability in the low wind jet velocity occurs by the end of nocturnal cooling (6:00 UTC, sunrise) and day-time heating (18:00 UTC, sunset). Its diurnal magnitude is modulated by the coastal orography and amplifies to about 2 m s −1 in the vicinity of Cape Ghir and Cape Blanc (21°N). Similar diurnal variability is present in the Benguela atmospheric low jet ( Figure 11). In contrast to the Canary region (Figure 10), the diurnal wind variability extends beyond the coastal boundary in the Benguela region, which is an observation for which we do not have any immediate explanation. As in the Canary region, the diurnal magnitude of Benguela winds is spatially modulated along the coastline, reaching a local maximum of 1.5 m s −1 in St. Helena Bay (32°S, South Africa). This regional amplification is probably linked to summer season radiative heating and its impact on changes of Ta, SST, land temperature, and land-sea interactions. In the Canary region, the main diurnal pattern associated with wind acceleration at 6:00 UTC and deceleration at 18:00 UTC develops along the Moroccan coast. It resembles nocturnal jets developing due to diurnal changes in atmospheric boundary layer stability and height. The strongest variability in the low wind jet velocity occurs by the end of nocturnal cooling (6:00 UTC, sunrise) and day-time heating (18:00 UTC, sunset). Its diurnal magnitude is modulated by the coastal orography and amplifies to about 2 m s −1 in the vicinity of Cape Ghir and Cape Blanc (21 • N). Similar diurnal variability is present in the Benguela atmospheric low jet ( Figure 11). In contrast to the Canary region (Figure 10), the diurnal wind variability extends beyond the coastal boundary in the Benguela region, which is an observation for which we do not have any immediate explanation. As in the Canary region, the diurnal magnitude of Benguela winds is spatially modulated along the coastline, reaching a local maximum of 1.5 m s −1 in St. Helena Bay (32 • S, South Africa). This regional amplification is probably linked to summer season radiative heating and its impact on changes of Ta, SST, land temperature, and land-sea interactions. Remote Sens. 2021, 13, x FOR PEER REVIEW 20 of 31 Figure 11. The same as in Figure 10 but for the Benguela region in July.

Summary
This study presents the description and quality assessment of a regional coastal wind analysis with enhanced spatial and temporal resolution. This analysis referred to as IFREMER satellite wind analysis is currently available only over two major eastern boundary upwelling systems, the Canary and Benguela. For accuracy estimation purposes, the IFREMER analysis is also calculated over the California upwelling system, where several ground-truth mooring time series used are available. To the best of our knowledge, the IFREMER analysis provides a uniquely long time series of regional surface winds estimated from well-calibrated and homogenized scatterometer and radiometer data. It provides 27 years (1992 to 2018) of 10 m wind speed, zonal, and meridional velocity components, wind stress vector, along with its curl and divergence. Analysis uncertainty data are also provided at each grid cell. Data are available at synoptic times (00 h:00, 06 h:00, 12 h:00, and 18 h:00 UTC) on a regular 0.125° × 0.125° grid.
Even though the approach used for the determination of IFREMER 6-h winds is pretty similar to that used in the Copernicus/CMEMS (L4) wind reprocessing [16], some significant changes should be underlined. In this study, the determination of spatial structure functions is based on high spatial resolution winds from Sentinel-1A/B SAR, which allow for better characterization of fine spatial scale variability at each grid cell, including Figure 11. The same as in Figure 10 but for the Benguela region in July.

Summary
This study presents the description and quality assessment of a regional coastal wind analysis with enhanced spatial and temporal resolution. This analysis referred to as IFREMER satellite wind analysis is currently available only over two major eastern boundary upwelling systems, the Canary and Benguela. For accuracy estimation purposes, the IFREMER analysis is also calculated over the California upwelling system, where several ground-truth mooring time series used are available. To the best of our knowledge, the IFREMER analysis provides a uniquely long time series of regional surface winds estimated from well-calibrated and homogenized scatterometer and radiometer data. It provides 27 years (1992 to 2018) of 10 m wind speed, zonal, and meridional velocity components, wind stress vector, along with its curl and divergence. Analysis uncertainty data are also provided at each grid cell. Data are available at synoptic times (00 h:00, 06 h:00, 12 h:00, and 18 h:00 UTC) on a regular 0.125 • × 0.125 • grid.
Even though the approach used for the determination of IFREMER 6-h winds is pretty similar to that used in the Copernicus/CMEMS (L4) wind reprocessing [16], some significant changes should be underlined. In this study, the determination of spatial structure functions is based on high spatial resolution winds from Sentinel-1A/B SAR, which allow for better characterization of fine spatial scale variability at each grid cell, including near coast cells. Temporal structure functions are also improved based on the use of consistent wind observations from all available scatterometers and radiometers ( Table 1).
The accuracy of the methodology as resulting 6-h wind analysis is assessed from comparisons with ground-truth 6h-averaged NDBC buoy wind measurements. IFREMER wind analysis exhibits high accuracy in representing the wind speed and the related components temporal variability, with high correlation coefficients and quite low biases and RMSD values.
At regional scales, the new IFREMER satellite wind analysis and the three wind analyses (CMEMS, CMMP, and ERA5) are compared with purely scatterometer winds in the Canary and Benguela upwelling regions. Since the four analyses use mostly the same scatterometer retrievals, their comparison indicates the ability of each particular surface wind analysis methodology to resolve the spatial and temporal variability present in scatterometer observations. The four analyses exhibit comparable deviations from pure scatterometer winds for wind speed and components, wind stress and components, and wind stress curl. The IFREMER wind analysis has the lowest RMSD from the scatterometry in offshore and nearshore areas. Furthermore, excluding comparisons based on the use of poor scatterometer samplings (e.g., ERS-2), the IFREMER wind analysis does not show any time dependence in its deviation from the pure scatterometry.
We believe that these unique data have the potential to resolve various spatial and temporal wind characteristics at local scales (e.g., in cape areas) of both Canary and Benguela regions, including diurnal, seasonal, and annual variability. Further investigations of the IFREMER wind analysis content will be performed through the forcing ocean circulation model process. The latter is expected to better address the role of the high spatial and temporal surface wind variability on the coastal circulation. In the future, the surface wind analyses should address the fast, highly variable, and high wind conditions. The latter tend to be smoother in all surface wind analyses. Further improvements, as required by oceanic and atmospheric modelers, include the enhancement of the temporal resolution of the scatterometer wind analyses. The latter would be achieved through the use of a new approach including numerical modeling and data learning methods. Additional scatterometer missions such as ASCAT-C, HY-2b, HY-2c, and CFOSAT will be considered in a reprocessed long time series of IFREMER wind analyses.
The long time series datasets will be available from IFREMER data center CERSAT (http://cersat.ifremer.fr/ (accessed on 2 March 2021)). In the meantime, the data can be downloaded from (ftp://ftp.ifremer.fr/ifremer/cersat/projects/upwesatwind (accessed on 2 March 2021)). For the Canary or Benguela regions, each 6-h analysis file contains wind speed and velocity components, wind stress and its components, wind stress divergence and curl on a 0.125 • × 0.125 • grid. Supplementary information such as error estimate is also provided. Data are in NetCDF 4 format.