Indian Ocean Crossing Swells: New Insights from “Fireworks” Perspective Using Envisat Advanced Synthetic Aperture Radar

: Synthetic Aperture Radar (SAR) in wave mode is a powerful sensor for monitoring the swells propagating across ocean basins. Here, we investigate crossing swells in the Indian Ocean using 10-years Envisat SAR wave mode archive spanning from December 2003 to April 2012. Taking the beneﬁt of the unique “ﬁreworks” analysis on SAR observations, we reconstruct the origins and propagating routes that are associated with crossing swell pools in the Indian Ocean. Besides, three different crossing swell mechanisms are discriminated from space by the comparative analysis between results from “ﬁreworks” and original SAR data: (1) in the mid-ocean basin of the Indian Ocean, two remote southern swells form the crossing swell; (2) wave-current interaction; and, (3) co-existence of remote Southern swell and shamal swell contribute to the crossing swells in the Agulhas Current region and the Arabian Sea.


Introduction
It is estimated that more than 75% waves across the world seas are generated by remote storms that propagate as long-period ocean waves or swells [1]. In particular, crossing swells [2], defined as complex sea state of at least two swell systems with different directions, are of importance in navigation safety. Furthermore, an understanding of climatology of crossing swells would be helpful in interpreting the global wave climate change.
The strongly swell-dominated areas, termed "swell pools", in the ocean basins have been identified by the means of altimeters decades ago [3]. To date, Synthetic Aperture Radar (SAR), particularly equipped with the specific Wave Mode (WM), is the only spaceborne sensor that is able to provide the directional ocean wave spectra on a global scale in the long term [4,5]. In contrast to altimeters only providing wave height without directional information, wave spectra from SAR WM could readily distinguish swell from local wind sea and even separate different swell systems. Hence, based on 10-years Envisat Advanced Synthetic Aperture Radar (ASAR) WM observations, a global view on crossing swell pools has been revealed [6]. However, as reported by Jiang et al. [7], some limitations for crossing swell on SAR description still exist, mainly owing to the questionable ASAR retrievals (e.g., [8]). Besides, the in-depth understanding of the formation mechanism of crossing swell remains improved from the interpretation using SAR WM.
Taking advantage of space-time consistency of swells in open ocean discovered by [9], the so-called "fireworks" technique, in which the integral properties of swells are derived from SAR WM observations, and used to both monitor and predict their evolution across the globe, has been developed by Collard et al. [10]. Because the "fireworks" technique could be employed to enhance the capability of ocean swell monitoring [11][12][13], a question Furthermore, the quality-controlled ASAR swells are processed into a higher level (Level-3) using the "fireworks" algorithm detailed in [10,17]. The main steps are briefly, as follows. The observed swells are firstly retro-propagated along great circles until land/island in step of 3 h, neglecting wave-current interactions. Secondly, the swell origins (usually storms) are identified by seeking the spatiotemporal convergence of the retropropagated trajectories. As a result, the "fireworks" analysis separates ASAR Level-2 swells into: • Refocused swell observations: the converging swell partitions. From validation [18], this sub-dataset of swells has proven to be of better qualify due to space-time consistency. In addition to SAR observed swell as snapshot, the whole swell evolutions (i.e., generation date/area and three-hourly propagation from origin towards coast) are also provided. The percentage is around 30% in the Indian Ocean. • Non-refocused data: the inconsistent observed swell partitions, which failed to be converged in our "fireworks" analysis and, thus, irrelevant to any storm events.
It is worth noting that the relative low rate (30%) of retained in "fireworks" reported here is consistent with that (25%, globally) regarding Level-3 swell products of Sentinel-1A/B in the framework of Copernicus Marine Environment Monitoring Service (CMEMS) officially released by ESA [19]. The explanation is mainly that our "fireworks" qualitycontrols can filter most of the suspicious SAR retrievals (badly inversed swell partitions and 180 • ambiguity, for instance, see [7,18]). Additionally, it is also reasonable to raise the hypothesis that the majority of data (~70%) that are considered as non-refocused are not all noise, and we might mine geophysical meanings.

Auxiliary Data
Blended winds from Institut français de recherche pour l'exploitation de la mer (Ifremer) [20], i.e., six-hourly ocean winds from different scatterometers via the kriging technique, are used here. Regarding ocean current, we rely on the daily geostrophic currents from the ESA GlobCurrent project (http://www.globcurrent.org accessed on 10 February 2021), which is an optimal analysis on multiple altimeters. They are both available on global 0.25 • grid.

Core Crossing Swell Pool and Its Properties
In this study, the criterion of defining a crossing swell event include (1) the angle between the two swell systems is larger than 40 • , (2) the wave height of the less energetic swell is larger than 1/2 of another swell, and both components >0.5 m. These thresholds are equivalent to [6,7].
In each location, we first calculated the probability of crossing swells using the number of crossing swell cases divided by the total number of records. Subsequently, we further normalized the probability through dividing by the peak value, making the comparative analysis between the results from "fireworks" and original SAR easier, as shown in Figure 1.
the ESA GlobCurrent project (http://www.globcurrent.org accessed on 10 February 2021), which is an optimal analysis on multiple altimeters. They are both available on global 0.25° grid.

Core Crossing Swell Pool and Its Properties
In this study, the criterion of defining a crossing swell event include (1) the angle between the two swell systems is larger than 40°, (2) the wave height of the less energetic swell is larger than 1/2 of another swell, and both components >0.5 m. These thresholds are equivalent to [6,7].
In each location, we first calculated the probability of crossing swells using the number of crossing swell cases divided by the total number of records. Subsequently, we further normalized the probability through dividing by the peak value, making the comparative analysis between the results from "fireworks" and original SAR easier, as shown in Figure 1.
As a result, statistics on 10-years Envisat wave mode ASAR yields the normalized occurrence map of the Indian Ocean crossing swells. Generally, both of the maps (using all ASAR observed swells and only refocused) highlight the core of crossing swell pool in the mid-ocean basin (hereafter, Region A: 15°S-5°N, 75°-90°E). Although consistent with previous studies [6,7] regarding geographical distribution, extra properties of the pool and its forming mechanism (e.g., origins and pathways) could be further investigated using "fireworks". We first illustrate this by the representative crossing swell case that occurred in Region A ( Figure 2). This event was captured by refocused ASAR "fireworks" data on 17 April 2008 at 12.36°S/83.03°E (black box in Figure  2c). The crossing swell, i.e., two systems propagating northwestward and northeastward, is clearly presented in the ASAR spectral retrieval, and js also visible from the imagette (blue and green lines indicating wave crests of two swells respectively overlaid on Figure  2a).
The generation and propagation of this event could be revealed from space, as depicted in Figure 2c. Based on "fireworks" analysis, the origins of the two swell systems were identified in the Southern Ocean, being represented as red disks. Accordingly, two extra-tropical storms (up to 35 knots) could be found from blended winds (Figure 2d) in Kerguelen Islands and west Australia, respectively, five days before the ASAR acquisition. Swells that were generated from these intense storms radiated and then crossed. Here, we predict the swell arrival dates and locations (colored lines) apart from As a result, statistics on 10-years Envisat wave mode ASAR yields the normalized occurrence map of the Indian Ocean crossing swells. Generally, both of the maps (using all ASAR observed swells and only refocused) highlight the core of crossing swell pool in the mid-ocean basin (hereafter, Region A: 15 • S-5 • N, 75 • -90 • E).
Although consistent with previous studies [6,7] regarding geographical distribution, extra properties of the pool and its forming mechanism (e.g., origins and pathways) could be further investigated using "fireworks". We first illustrate this by the representative crossing swell case that occurred in Region A ( Figure 2). This event was captured by refocused ASAR "fireworks" data on 17 April 2008 at 12.36 • S/83.03 • E (black box in Figure 2c). The crossing swell, i.e., two systems propagating northwestward and northeastward, is clearly presented in the ASAR spectral retrieval, and js also visible from the imagette (blue and green lines indicating wave crests of two swells respectively overlaid on Figure 2a).
The generation and propagation of this event could be revealed from space, as depicted in Figure 2c. Based on "fireworks" analysis, the origins of the two swell systems were identified in the Southern Ocean, being represented as red disks. Accordingly, two extra-tropical storms (up to 35 knots) could be found from blended winds ( Figure 2d) in Kerguelen Islands and west Australia, respectively, five days before the ASAR acquisition. Swells that were generated from these intense storms radiated and then crossed. Here, we predict the swell arrival dates and locations (colored lines) apart from storms, following linear wave propagation theory with the group velocity derived from ASAR swell periods (12 s and 14.5 s for two systems). As a result, the prediction lines from different swell generation area intersect close to ASAR acquisition in both space and time, indicating the Remote Sens. 2021, 13, 670 4 of 12 agreement between the satellite observation and the analysis from "fireworks". Therefore, in addition to the detection of crossing swell events, origins, and propagating routes that are associated with crossing swell pool could be reconstructed, providing a mechanism of these complex sea states from space, especially by our unique "fireworks" technique.
ASAR actual observations, see colored lines), we could further infer that swells from SW (SE) of the basin travel six to eight days before they cross each other in Region A, and then reach the coast in the Bay of Bengal (the Arabian Sea) two (4) days later.
Furthermore, the distribution of SAR observed swell wavelength and direction are investigated. Crossed swells in Region A generally propagate following swell incident regimes towards northeast (20-50°) and northwest (320-350°), coinciding with the two corridor-like pathways that are presented in colored map of Figure 3a, as shown by the spectral statistics presented in Figure 3b. Note that more swells are coming from southwest, and these swells have a longer wavelength (350 m approximately) than those from southeast (around 300 m) statistically. Both of these features indicate that Southern swells have more influence on the Bay of Bengal than the Arabian Sea, in line with previous finding from modeling [21,22]. Here, we statically analyze the Region A using 10-years "fireworks" data. Ocean waves that are responsible for crossing swells mainly originate from storms in the Southern Ocean westerly around 40 • S (black dots), as shown in Figure 3a. In addition, the normalized occurrences of swell traveling routes (i.e., dots backward and forward propagated from SAR observations in "fireworks" for every 3 h) are overlaid on the map. From the density, the propagation occurrences peak in the core crossing swell pool, decreasing towards the Southern Ocean (swell origins) and the north Indian Ocean coast (swell destinations). This characterizes the two propagating routes from the southwest (SW) and southeast (SE) Indian Ocean towards the Bay of Bengal and the Arabian Sea, respectively. From the time difference ("fireworks" three-hourly propagations minus the ASAR actual observations, see colored lines), we could further infer that swells from SW (SE) of the basin travel six to eight days before they cross each other in Region A, and then reach the coast in the Bay of Bengal (the Arabian Sea) two (4) days later.
Furthermore, the distribution of SAR observed swell wavelength and direction are investigated. Crossed swells in Region A generally propagate following swell incident regimes towards northeast (20-50 • ) and northwest (320-350 • ), coinciding with the two corridor-like pathways that are presented in colored map of Figure 3a, as shown by the spectral statistics presented in Figure 3b. Note that more swells are coming from southwest, and these swells have a longer wavelength (350 m approximately) than those from southeast (around 300 m) statistically. Both of these features indicate that Southern swells have more influence on the Bay of Bengal than the Arabian Sea, in line with previous finding from modeling [21,22]. Ifremer blended winds (12 April 2008 12:00). Red disks and white crosses indicate origins of crossing swell systems inferred from "fireworks" analysis in (c,d).

Discussion
From Figure 1, the major discrepancy between the two maps is the missing of crossing swells detected from refocused data in Level-3 "fireworks". This is because some ASAR observed swell partitions consisting the crossing sea that were non-refocused in "fireworks" processing. In this study, we focus on three noticeable regions: In Figure 4a-f, we statistically analyze the swell spectral distribution from refocused and non-refocused ASAR swell partitions for these regions. All of the distributions from refocused data (Figure 4a-c) are different from the distinct bi-directional incident swell signature that is presented in Figure 3b. Moreover, non-refocused data distributions (Figure 4d-f) show various patterns, indicating different reasons for filtering out through "fireworks" analysis.
In addition, average altimeter-derived current speeds (Figure 4g) illustrate that Region B is well located in the main stream of the Agulhas current, and Region C is in the north of Agulhas Return Current (ARC). It implies the possibility of relation to this strong current, different from Region D in Arabian Sea (Figure 4h). Hence, we discuss Region B, C, and D separately below.

Discussion
From Figure 1, the major discrepancy between the two maps is the missing of crossing swells detected from refocused data in Level-3 "fireworks". This is because some ASAR observed swell partitions consisting the crossing sea that were non-refocused in "fireworks" processing. In this study, we focus on three noticeable regions: In Figure 4a-f, we statistically analyze the swell spectral distribution from refocused and non-refocused ASAR swell partitions for these regions. All of the distributions from refocused data (Figure 4a-c) are different from the distinct bi-directional incident swell signature that is presented in Figure 3b. Moreover, non-refocused data distributions (Figure 4d-f) show various patterns, indicating different reasons for filtering out through "fireworks" analysis.
In addition, average altimeter-derived current speeds (Figure 4g) illustrate that Region B is well located in the main stream of the Agulhas current, and Region C is in the north of Agulhas Return Current (ARC). It implies the possibility of relation to this strong current, different from Region D in Arabian Sea (Figure 4h). Hence, we discuss Region B, C, and D separately below.

Crossing Swells in Agulhas Current Regions
From the spectral statistics presented in Figure 4a,d, the ASAR observed crossing swells in Region B have two sorts: (1) refocused ASAR swells mainly originated from Southern Ocean and propagated towards northeastward and (2) non-refocused swell with a broad directional pattern scattering from north to east.
The fact that most of the observed crossing swells here did not converge in our "fireworks" procedure could be explained by the wave-current interaction. This effect includes wave spectrum modulation, wave energy enhancement, or attenuation by the co-existing surface current [23,24]. The equatorial currents could refract swells up to 10° according to buoy records, as documented by [25]. Here, Region B is characterized as the Agulhas Current main stream, the strongest western boundary in Southern Hemisphere [26]. The average current speeds are >1.2 m/s during the 10-years Envisat mission, as illustrated in the GlobCurrent current map shown in Figure 4g. Because of the reflection of the strong Agulhas, incoming swells from the Southern Ocean would deviate from their incident directions and, thus, severely interfere with the refocusing algorithm of retropropagated great circle trajectories. Recall that ignoring ocean current is one of the assumptions for our "fireworks" technique.
The complex crossing sea states that were caused by wave-current interaction were also documented by other remotely sensed data. For instance, investigation using Sentinel-2 optical imagery has reported the spectra splitting of incident swell caused by Agulhas Current [27]. Here, this formation mechanism of crossing swell is highlighted through inter-comparison between the remaining and rejected data from the SAR "fireworks".
This could also be demonstrated by the typical case presented in Figure 5. Inferred from "fireworks", the extra-tropical storm (>50 knots wind speed from Ifremer blended winds) generated swells that were captured by 12 ASAR WM imagettes (numbered from 0 to 11) along descending satellite orbit on 18 April 2007 in the Agulhas region. Figure 6 presents the corresponding ocean spectra retrieved from ASAR WM imagettes. ASAR observations #3 and 6 are located in the stronger Agulhas current areas, when compared

Crossing Swells in Agulhas Current Regions
From the spectral statistics presented in Figure 4a,d, the ASAR observed crossing swells in Region B have two sorts: (1) refocused ASAR swells mainly originated from Southern Ocean and propagated towards northeastward and (2) non-refocused swell with a broad directional pattern scattering from north to east.
The fact that most of the observed crossing swells here did not converge in our "fireworks" procedure could be explained by the wave-current interaction. This effect includes wave spectrum modulation, wave energy enhancement, or attenuation by the co-existing surface current [23,24]. The equatorial currents could refract swells up to ±10 • according to buoy records, as documented by [25]. Here, Region B is characterized as the Agulhas Current main stream, the strongest western boundary in Southern Hemisphere [26]. The average current speeds are >1.2 m/s during the 10-years Envisat mission, as illustrated in the GlobCurrent current map shown in Figure 4g. Because of the reflection of the strong Agulhas, incoming swells from the Southern Ocean would deviate from their incident directions and, thus, severely interfere with the refocusing algorithm of retro-propagated great circle trajectories. Recall that ignoring ocean current is one of the assumptions for our "fireworks" technique.
The complex crossing sea states that were caused by wave-current interaction were also documented by other remotely sensed data. For instance, investigation using Sentinel-2 optical imagery has reported the spectra splitting of incident swell caused by Agulhas Current [27]. Here, this formation mechanism of crossing swell is highlighted through inter-comparison between the remaining and rejected data from the SAR "fireworks".
This could also be demonstrated by the typical case presented in Figure 5. Inferred from "fireworks", the extra-tropical storm (>50 knots wind speed from Ifremer blended winds) generated swells that were captured by 12 ASAR WM imagettes (numbered from 0 to 11) along descending satellite orbit on 18 April 2007 in the Agulhas region. Figure 6 presents the corresponding ocean spectra retrieved from ASAR WM imagettes. ASAR observations #3 and 6 are located in the stronger Agulhas current areas, when compared to the acquisitions of #1 and 4, as shown in GlobCurrent surface velocity colored in Figure 5 (right). Correspondingly, from Figure 6, in contrast to the unimodal swell spectra (#1 and #4) observed from ASAR in the absence of strong Agulhas current, the ocean swell spectra (#3 and #6) exhibit bimodality due to the spectra splitting by the wave-current interaction and, thus, present the crossing swell states. Although the surface current is not taken to consideration in our "fireworks" analysis in present version, the crossing swell mechanism of wave-current interaction could be highlighted by the careful examination of these rejected observations.
Here, we also discuss the possibility of including wave-current interactions in an improved "fireworks" algorithm. Even if great-circles back-tracing is most often in good agreement with the location of large storms, there are indeed some cases in which there are some inconsistencies (as mentioned by [25], and also highlighted here). These arise from wave refraction by surface currents. One feasible solution is based on the spherical ray tracing using reconstructed surface currents that are inferred from satellite altimetry, rather than the present simpler great-circle propagation. Such a technique has been tested for wave spectral observations from buoy, and it could also be implemented in an improved version of "fireworks" algorithm for SAR derived swells. However, this implementation requires a careful analysis, since, as it has experienced [25], directional fluctuations of the order of 10 • can be observed at the buoy location, which result from the arrival of wave packets that have travelled along the different rays. This can bias the directional average at the buoy, not to mention the significant wave height.    With respect to Region C, statistical analysis on non-refocused swell observations (in Figure 4e) shows a very scattered pattern, and the cluster containing the highest density of points representing the swell propagating towards directions that range 30-70 • (represented by black dashed lines). Examining the Globcurrent map shown in Figure 4g, one could see that the eastward flowing ARC (along 40 • S approximately, ranging from 20 • to 35 • E) is just located in the regime of non-refocused swell propagating pathway that corresponds to the spectral statistical analysis in Figure 4e. Thus, the relatively strong ARC (though weaker than core Agulhas current) could interfere the current-ignored retro-propagation in firework, and lead to the non-refocused, and thereby discrepancy, of crossing swell in Figure 1.

Crossing Swells in the Arabian Sea
From statistics of ASAR WM spectra that are shown in Figure 4c, only swells originating from the Southern Ocean and northward propagating towards the Arabian Sea could be refocused in "fireworks" analysis. With respect to the non-refocused (Figure 4f), the pattern of directional statistics indicates two kinds of waves that contribute the crossing swells that were identified from Level-2 observations, but they were rejected during "fireworks": (1) swells propagating from SE and (2) waves from northwest (NW). We examine these two cases of ASAR left-overs spectra and discuss, as follows.
The island blocking effect could explain the non-refocusing of ASAR observed swells from SE in the Arabian Sea. When the swells from directions of 140-160 • (regime defined by dash lines in Figure 4f,h) propagate across the Maldives islands, although attenuation and diffraction will occur, some swell energy will continuously travel and arrive at the Arabian Sea, thus crossing with other swells. However, because of the island obstructions, ASAR observed swells from these directions could not be retro-propagated and then associated with their origins successfully and, therefore, be regarded as left-overs in the "fireworks".
Regarding waves from NW, the reason of missing crossing swell may be attributed to shamal swell events. The northwest propagating shamal swells [28] in the Arabian Sea is in agreement with the ASAR derived spectra towards SE regarding the propagation direction ( Figure 4h). In fact, because shamal swells are generated by the strong northwesterly winds in the Arabian Peninsula [29], it is difficult to associate these swells with storms (the retro-propagation was stopped near the coast in our "fireworks" analysis). Therefore, another crossing swell formation mechanism, i.e., a combination of Southern Ocean swells and shamal swells could be revealed.
Using the buoy data along the west coast of India, Aboobacker et al. [28] reported a representative shamal swell event in the Arabian Sea throughout 2-4 February 2008. Here, we revisit this event using SAR data from space. This shamal swell case was captured by wave mode ASAR (Figure 7a) on 3 February 2008 at 16.26 • N/61.80 • E, as shown in Figure 7 (red box in Figure 7b). The SW propagating swell system from ASAR derived spectrum that is shown in Figure 7a (2nd partition) could be associated with that shamal event. From three-days successive Ifremer blended wind maps in Figure 7c-e, the northwesterly shamal wind event began on 1 February in the Persian Gulf, blew offshore the Arabian Peninsula on 2 February, and then vanished over the Arabian Sea on 3 February at the time of SAR observation. Additionally, one could see the topography effect of the Arabic Peninsula on the wind flow, i.e., the winds coming from the Persian Gulf were obstructed at the Hormuz Strait and separated into two flows in the Arabian Sea and then created the two swell systems (one from N and one from NW) that are presented in Figure 7a.

Conclusions and Perspectives
In this communication we present the new insights into Indian Ocean crossing swells from 10 years Envisat ASAR data taking advantage of "fireworks" analysis. Our results present the three crossing-swell dominated areas, i.e., the mid-Indian Ocean basin, the Arabian Sea, and Agulhas region.
The "fireworks" analysis on SAR WM observations improves our understanding of crossing swells from space. On the one hand, knowledge on swell origins, propagating pathways for crossing swell pool (the mid-Indian Ocean basin), are revealed from the "fireworks" product. On the other hand, some of the rejected observations in the "fireworks analysis" indeed make geophysical sense, revealing the formation of crossing swell. Thus, the discrimination of crossing swell mechanisms could be revealed by comparative analysis taking opportunity of the SAR WM fireworks technique. In the Indian Ocean, the wave-current interaction and the combination of remote southern swell Besides, we could further identify the 1st swell partition in Figure 7a as a long-distance traveling swell system (the pathway depicted as a dashed line in Figure 7b) originating from the extra-tropical storm on 25 January in the Southern Ocean (>50 knots from Ifremer blended winds in bottom of Figure 7b) using ASAR WM "fireworks" analysis. Hence, in this case, the complex crossing swell state resulted from the co-existence of Southern Ocean swell and shamal swell could be interpreted from the SAR observation using "fireworks" analysis, complementary to the in situ that only documented the shamal swell [28].
Another interesting pattern shown in Figure 4f is that some points (close to North) representing swells from the Southern Ocean are still retained in the non-refocused subset.
This portion data failed to be refocused in "fireworks", owing to the similar wave-current interaction mechanism that was mentioned in the last section. Not Agulhas, but equatorial currents and eddies, in the midway of propagating route, refract some of these long-distance traveling swells. Therefore, in the present "fireworks", this wave-current interaction bends the back-propagated great-circle rays and leads to a focus on non-storm regions (over sea ice or even on Antarctica for instance). The examples of buoy measurements and related theory are detailed in [25].

Conclusions and Perspectives
In this communication we present the new insights into Indian Ocean crossing swells from 10 years Envisat ASAR data taking advantage of "fireworks" analysis. Our results present the three crossing-swell dominated areas, i.e., the mid-Indian Ocean basin, the Arabian Sea, and Agulhas region.
The "fireworks" analysis on SAR WM observations improves our understanding of crossing swells from space. On the one hand, knowledge on swell origins, propagating pathways for crossing swell pool (the mid-Indian Ocean basin), are revealed from the "fireworks" product. On the other hand, some of the rejected observations in the "fireworks analysis" indeed make geophysical sense, revealing the formation of crossing swell. Thus, the discrimination of crossing swell mechanisms could be revealed by comparative analysis taking opportunity of the SAR WM fireworks technique. In the Indian Ocean, the wavecurrent interaction and the combination of remote southern swell and local shamal swell are highlighted in the Agulhas current region and the Arabian Sea, respectively, as compared to the mechanism of two remote southern swells cross-over in the mid-Indian Ocean basin.
Based on this 10-year remotely sensed swell dataset processed by the unique "fireworks" technique, studies could be further performed regarding the spatio-temporal distributions of crossing swells that resulted from different mechanisms. This long-term analysis on a global scale could deepen our understanding of nature of ocean. Future work will also be dedicated to the improvement of our "fireworks" algorithm by considering wave-current interactions.
Although the Envisat mission ended in 2012, C-band SARs that are equipped with wave mode continue providing swell spectra nowadays, such as European Copernicus missions Sentinel-1A/B [19] and Chinese Gaofen-3 [30]. In addition, Chinese-French joint mission CFOSAT launched in 2018, which is also capable of providing 2D ocean wave spectra globally [31]. With the accumulation of more data in near future, it is encouraging to enlarge the temporal series of remotely sensed swell spectra archive and merging multisensor data, therefore improving the significance of our results. Data Availability Statement: Envisat ASAR data were provided by the ESA. Blended winds and ocean current data were obtained from Ifremer and ESA's GlobCurrent project.