Next Article in Journal
Mapping Soil Organic Carbon for Airborne and Simulated EnMAP Imagery Using the LUCAS Soil Database and a Local PLSR
Previous Article in Journal
Changing Patterns of Lakes on The Southern Tibetan Plateau Based on Multi-Source Satellite Data
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Technical Note

Remote Cross-Calibration of Wave Buoys Based on Significant Wave Height Observations of Altimeters in the Northern Hemisphere

1
National Marine Environmental Forecasting Center, Beijing 100081, China
2
Department of Marine and Oceanography, Meteo France, 31057 Toulouse, France
*
Author to whom correspondence should be addressed.
Remote Sens. 2020, 12(20), 3447; https://0-doi-org.brum.beds.ac.uk/10.3390/rs12203447
Submission received: 23 September 2020 / Revised: 15 October 2020 / Accepted: 18 October 2020 / Published: 20 October 2020
(This article belongs to the Section Ocean Remote Sensing)

Abstract

:
Consistency between national wave buoy networks is extremely important for wave climate studies and verification of global operational wave forecasting systems; however, it is insufficiently investigated. The validation of altimeter significant wave heights (SWHs) with the wave buoy networks of China, Europe and the National Data Buoy Center (NDBC) show significant divergence in assessments. This reveals a negative bias and larger root mean square error and scatter index from the Chinese buoy network than from the European and NDBC buoy networks. A remote cross-calibration method is presented using the collocations between altimeters and buoys to match the buoy observations from different networks. The Chinese buoys are found to yield a negative bias of −0.127 m compared to European/NDBC buoy networks. The cross-calibration equation is achieved by regression of the SWHs between the Chinese and European/NDBC buoy networks. The use of this remote cross-calibration significantly reduces the inconsistency between the Chinese and European/NDBC buoys in the validation of SWH from altimeter HY2B.

Graphical Abstract

1. Introduction

Obtaining accurate wave conditions is a crucial issue for both marine activities and geographic science, such as the ocean climate. Significant wave height (SWH) is the most commonly used and effective parameter not only for describing sea surface dynamics but also for wave forecasting.
Providing in situ observations, wave buoys are the traditional means of observing waves. Wave buoys can measure diverse and comprehensive wave information, such as the SWH, mean period and even directional wave spectra. Because of their direct measurement of the ocean surface, buoy observations are always considered one of the most accurate ways to obtain SWH, which tends to be used as the ground truth for wave conditions in a wide variety of related applications. Wave observations are widely used in the validation/calibration phase of wave remote sensing missions to improve the retrieval algorithm [1], and buoy data are used as reliable references to validate and improve the physical processes, parametric schemes or model settings implemented in wave models [2,3,4]. As the traditional method of wave monitoring, buoys record long-term data that can indicate trends in climate change [5]. Therefore, wave buoy observations are one of the fundamentals of both oceanography and atmospheric sciences, and the accuracy and consistency of wave buoys are both extremely important for further usage.
However, a wave buoy is a complex system that is composed of multiple instruments and can be classified into various types or models. There are significant differences in the shapes and sizes of the hulls, measurement principles of the sensors, calculation methods, etc., and all these aspects have certain impacts on the accuracy of the buoy measurements [6]. It is known that a national wave buoy network is always a combination of various types of buoys, in addition to the diverse models used by different countries. The buoy network supported by the National Data Buoy Center (NDBC) under the National Oceanic and Atmospheric Administration (NOAA), which is widely used in operations and research, contains several different hull types, such as 3, 6 and 12 m diameters [7]. The European buoy network, which is operated by the Met Office and Meteo-France, mainly uses 3 m buoys in the Atlantic and some vessels in the English Channel [8]. Meanwhile, the Chinese operational wave buoy network, which belongs to the State Oceanic Administration (SOA) and Ministry of Natural Resources (MNR), is made up of 10m discus buoys that weight more than 50 tons each. It is necessary and valuable to calibrate the different buoy types or instrumental conditions before applying their data. Therefore, adequate works and studies focus on this objective. Both hardware and algorithm upgrades have been performed to improve the accuracy and maintain the consistency of NDBC buoys [9,10]. The errors in the NDBC buoy data are removed before these data are used in the analysis of long-term trends in wave height [11]. The wave measurements from light vessels were indicated to be less accurate than those from buoys based on the comparison to the SWH of the wave model, and vessels impart a significant negative bias due to their weights and sizes when compared to regular wave buoys [12]. Apart from the difference in buoy attributes, biological attachments also impact SWH observations if the buoy experiences biofouling [13].
Consistency between different national wave buoy systems is equally important as, if not more important than, the intercomparison within a single buoy network. This is because inconsistency between different networks prompts misleading directions of global wave forecast model tuning [14] or incorrect regional wave climate trends. Using parallel observations in field experiments is the most direct and effective way to obtain the consistency of different buoys by making simultaneous comparisons at close distances [15]. However, there would be many difficulties (e.g., excessive cost) in performing parallel comparisons between buoys from different national buoy networks, e.g., the NDBC and the Chinese operational buoy observation systems.
Therefore, a remote cross-calibration method dealing with widely separated wave buoys is essential for the subsequent usage of buoy wave data. The relative accuracy of buoys can be obtained by comparing both to another independent wave data source. The SWH results from a wave model are used as a reference to characterize the difference between 3 m buoys and light vessels [12]; however, the accuracy of the wave model is strongly related to regional metrological, climate or topographic features. Therefore, this method, which treats the wave model as a stable reference, may not work for buoys located excessively far apart.
After three decades of rapid development and improvement, space-borne altimeters such as the Jason [16,17,18] and HY2 series [19,20] now provide global wave observations under all weather conditions. The SWH from the altimeter nadir is accurate enough to be applied in validation or improvement of numerical wave models [21,22,23] and in wave climate studies [24]. Furthermore, the inversion principle of altimeters, which obtains the SWH from the slope of the leading edge of the returned signal, means that the accuracy of altimeters is not significantly related to regional metrological conditions. Therefore, the SWH from altimeters can be used as a reliable standalone reference in characterizing the buoy errors. Currently, Jason-3, Sentinel-3A and 3B, SARAL and HY-2B are routinely used in operational service with high and stable accuracy.
In this paper, we present how we use multi-mission altimeter SWHs to perform remote cross-calibration of buoys belonging to different networks. We first introduce the altimeter observations and the buoy data that are used. We then show the significant differences that result from the assessments of the same altimeter against different buoy networks. Then, we explain how to remotely “connect” buoys from different networks by collocation of altimeters and accomplish the remote cross-calibration of these buoys. The effectiveness of the remote calibration method is proven by assessing the HY2B altimeter SWH based on different buoy networks. Finally, the discussion and conclusions are presented.

2. Data and Validation of Altimeters against Buoy Networks

In this paper, the altimeter L2 geophysical data records (GDRs) for Jason-3, SARAL, and Sentinel-3A from 1 January 2018, to 31 December 2019, were used to perform both the validation and cross-calibration of buoy networks. The L2 products from April 2019 to July 2020 of HY2B were used in the assessment to show the difference in validation results before and after buoy cross-calibration. HY2B was launched in October 2018 and is now the world’s newest altimeter and the calibration/validation phase has indicated good accuracy of wave observations [25,26].
The wave observations of the buoys from the NDBC, European and Chinese networks are collocated with altimeters. The locations of all buoys mentioned are shown by red circles marked in Figure 1. We used 45 buoys from the NDBC, 4 buoys from the European network and 16 buoys from the Chinese network in this paper. To remove the land effect on the altimeter signal, only buoys with distances offshore greater than 50 km were considered. These buoys provided hourly SWH and mean wave period observations.
The validations of the SWHs from the altimeters were performed against the buoys from the NDBC, European and Chinese networks to verify whether the three buoy networks showed consistent results. The observations from the altimeters and buoys were matched according to the following rules: spatially, we retained only the altimeter footprint nearest to the buoy location and the distance between the crossovers was considered to be less than 50 km; temporally, the gap in the SWH observation times between the altimeter and the collocated buoy had to be less than 1 h (±30 min). During the validation of altimeter SWH data, statistical analyses based on the bias, normalized root mean square error (NRMSE), root mean square error (RMSE) and scatter index (SI) were implemented. These statistical parameters are defined as follows:
b i a s = i = 1 N ( T i T N i ) N ;
N R M S E = 100 % × i = 1 N ( T i T N i ) 2 / N T N i ¯ ;
R M S E = i = 1 N ( T i T N i ) 2 N ;
S I = 100 % × i = 1 N [ ( T i T ¯ ) ( T N i T ¯ N ) 2 / N T N ¯ ;
where T i represents the SWH from altimeters and T N i is the matched buoy observation.
The assessments of Jason-3, Sentinel-3A, SARAL and HY2B SWHs were performed against the three buoy networks. Figure 2 first shows the assessment of the Jason-3 SWH from the L2 GDR products. From the scatter diagrams in Figure 2a–c, in general, the Jason-3 SWH shows a reasonable scatter against the three buoy systems. However, one can find a positive shift in the scatter against the Chinese network compared to the other two buoy networks. The statistical parameters also indicate this shift. Limited biases of 0.020 m and 0.047 m were obtained from the comparison against the European and NDBC buoy networks, respectively, while a 0.168 m positive bias was obtained from the Chinese buoy network. Significant differences in NRMSE and SI were also found in the assessment results between the Chinese buoys and NDBC/European buoys, which were 25.1% to 11.5%/12.2% and 20.6% to 11.5%/12.0%, respectively. We also investigated the variation of the bias and NRMSE with the SWH ranges, as illustrated in Figure 2c–e. Six SWH ranges were considered. We clearly see that the positive bias from the validation against the Chinese network was observed for all SWH ranges, while the bias distributions from the European and NDBC buoys were close to each other. This shows that the buoys from the European and NDBC networks give similar bias estimation with respect to Jason-3, but the Chinese buoys give a much more positive bias. Figure 2e shows the percentage of the bias obtained from dividing by the mean SWH in each SWH range. The Chinese buoys again showed a significant deviation when compared to the European and NDBC networks. In the section from 1 m to 1.5 m, the Chinese buoys showed that the proportion of the bias accounts for 30% of the average SWH observed, which is two times larger than those from the European and NDBC buoys. A similar situation was also found in the NRMSE distribution. The NRMSE from Chinese buoys again stands out when the SWH was lower than 3 m, after which all the curves became centralized. From the assessment of Jason-3 SWH against the three buoy networks, the European and NDBC buoys show similar validation results; however, an obvious difference was found from the Chinese buoys which mean an inconsistency among the buoy networks can be inferred.
If the Jason-3 SWH is used as the reference, the relative bias distributions with the mean wave period of the three buoy networks are indicated in Figure 3. As shown in Figure 3a, the Chinese buoys gave the largest negative bias throughout the mean wave period of 4s to 9s. The European network and NDBC buoys yield limited bias within −0.08 m to 0.07 m, and their bias was distributed similarly throughout the whole mean wave period region. A more significant discrepancy between Chinese buoys and European/NDBC is shown in Figure 3b. The negative bias reaches −14% to −23% of the mean SWH observed when the mean period is less than 6 s, which is likely due to wind sea conditions. The bias percentages of the three buoy networks were centralized when the mean period increased to more than 7 s, which indicates a more swell-dominated condition. The related NRMSEs of the three buoy networks showed the same trend as the mean wave period increases, as shown in Figure 3c. That is, a large discrepancy was found between Chinese buoys and the other two networks and the curves merged together during the likely swell cases.
In addition to Jason-3, validations of Sentinel-3A, SARAL and HY2B SWH against three buoy networks were performed. The statistics for these altimeters, including Jason-3, are listed in Table 1. From the comparison of the biases of all the altimeters involved, Chinese buoys showed the highest positive bias, while European and NDBC buoys gave relatively limited bias values, and the biases of these two buoy networks were similar to each other. The same situation occurred in the statistics of NRMSE and SI. The assessments based on the collocated Chinese buoys achieved the largest NRMSE and SI, approximately 15% to 25%, but they were approximately 10% for the European and NDBC buoy networks. A possible explanation for the large difference in the altimeter validations between the Chinese and European/NDBC buoy networks is the different sizes and weights of the buoys. The Chinese network is formed by 10 m discus buoys that weigh more than 50 tons, while the European/NDBC networks mainly use 3 m buoys that weigh much less. Like the heavy vessels in the English Chanel [12], the wave-following property of Chinese buoys is quite limited, making the large Chinese buoys less easily shaken than other buoys under light wave conditions, especially in the wind sea. The insufficient calibration and parallel observation experiment made the situation worse because the calibration is not always carried out when the payload is installed on or under the buoy, and the SWH corrections due to the poor wave-following property are mainly achieved by very limted experiments. The effects of their corrections now seem to be inadequate. Therefore, the Chinese buoys always give an underestimated SWH compared to the European/NDBC buoys, and the negative bias decreases with an increase in wave height.

3. Remote Cross-Calibration Method

Since a significant inconsistency was found between the Chinese buoy network and European/NDBC buoys, cross calibration was seriously required. During the collocated assessment of altimeter SWH, the following equations hold:
S W H b = T + b b ( T ) + σ b ( T ) ;
S W H a l t = T + b a l t ( T ) + σ a l t ( T ) ;
where, S W H b and S W H a l t are the SWH observations of the collocated buoy and altimeter, respectively; T is the true wave SWH; b b ( T ) and b a l t ( T ) are the biases of the buoy and altimeter, respectively, when the SWH is T ; and σ b ( T ) and σ a l t ( T ) are the random errors of the buoy and altimeter, respectively, when the SWH is T . We assumed that the bias and randomness were functions that only related to the true SWH.
Combining (5) and (6), we can eliminate the true SWH:
S W H a l t = S W H b + b a l t ( T ) b b ( T ) + σ a l t ( T ) σ b ( T ) .
Then, taking mathematical expectations on both sides of the equation, the random error should be 0, which leads to:
E ( S W H a l t ) = E ( S W H b ) + E ( b a l t ) E ( b b ) ,
where E ( b a l t ) and E ( b b ) are actually the biases of the altimeter and buoy, respectively.
Similarly, when this altimeter collocated with another buoy, we can also have:
E ( S W H a l t ) = E ( S W H b 2 ) + E ( b a l t ) E ( b b 2 ) ,
where the subscript “b2” indicates the parameters of another buoy.
Now, it should be noted that for a fixed SWH range (for example, 1, 2 m), if there are “enough” and “random” samples of altimeter observations distributed in this range, the expectation of altimeter observations collocated with these two buoys should be identical, and the biases of altimeters in this range should also be equal. Then, from Equations (8) and (9), we can infer that:
E ( S W H b ) E ( S W H b 2 ) = E ( b b ) E ( b b 2 ) .
Now, we make a connection of the two buoy observations by collocating these two buoys with altimeter observations. E ( b b ) E ( b b 2 ) indicates the relative difference between these two buoys. There is no dependency upon altimeter attributes in Equation (10).
The following section describes how we implemented the remote matchup of buoys through altimeter collocation. Considering that the relative difference between buoys may relate to the SWH and that the numerical resolution of buoy is always 0.1 m, several ranges of the collocated SWHs from the altimeter were divided. They are: (0.35, 0.45 m), (0.45, 0.55 m), (0.55, 0.65 m), … (3.95, 4.05 m), of which the expectations are 0.4 m, 0.5 m, 0.6 m, …and 4.0 m, respectively, under a fully random distribution. One example of the range (0.65, 0.75 m) is shown in Figure 4, indicating the distribution of Chinese buoys and NDBC buoys when Jason-3 gives the expectation of 0.7 m SWH.
From Figure 4, we note that when Jason-3 gives an expected SWH of 0.7 m, the collocated Chinese buoys and NDBC buoys show different SWH distributions. Shown as the blue bars in Figure 4, the Chinese buoys tend to give a lower centered distribution than the NDBC buoys, which again proves that the Chinese buoys produce a significant negative bias compared to the NDBC buoys. Then, the expectations of Chinese buoys and NDBC buoys can be calculated, and the values are 0.50 m and 0.65 m, respectively. This means that when the NDBC buoys give a SWH of 0.65 m, the Chinese buoys would give a SWH of 0.50 m.
Then, all SWH ranges were processed as described above, and because the method does not depend on the type and accuracy characteristics of altimeters, the collocations with Jason-3, Sentinel-3A and SARAL were all used to avoid unsatisfied statistical assumptions caused by an insufficient number of samples. Since high consistency is found between the European buoys and NDBC buoys, in the following cross-matching and calibrations, we considerthe buoys from the European and NDBC networks together.
Based on all the collocations of buoys and Jason-3, Sentinel-3A and SARAL, the remote matches of the Chinese buoys and European/NDBC buoys are shown as scatters in Figure 5. It can be easily noted that all the matchups from the three altimeters show the same pattern in which the Chinese buoys underestimate the SWH compared to European/NDBC buoys. The negative bias of Chinese buoys is −0.127 m. The least squares regression method is applied to obtain the cross-calibration equation to correct the Chinese buoys:
S W H c a l = 1.117 × S W H c 0.931 ,
where S W H c a l and S W H c are the calibrated and original SWH observations from Chinese buoys, respectively.

4. Validation of HY2B with Calibrated and Original SWHs

As shown in Table 1, obvious discrepancies can be found in the bias, RMSE, and scatter index between the Chinese buoys and the European/NDBC buoys. New validation of the HY2B SWH was performed based on the calibrated Chinese buoys using Equation (11) and compared to the original Chinese buoys and European/NDBC buoys. The results are shown in Figure 6.
From Figure 6a–c, one of the most significant improvements after the calibration of the Chinese buoy SWH was the reduction of the bias difference between the original Chinese buoy SWH and the European/NDBC network value. Because of the decreases in the bias, the RMSE also decreased to 0.230 m from the original 0.301 m which led to a decrease in the NRMSE from 24.1% to 16.9%, which is closer to the value of 12.3% from the European/NDBC buoys. The scatter index also decreased to 14.3% from 15.4% but was still greater than that from European/NDBC buoys. This is mainly because the majority of the collocated SWHs from Chinese buoys are less than 2.5 m, which results in a large scatter index, as shown in Figure 6f. There was also significant improvement when looking at the bias distribution with the SWH. As indicated in Figure 6d, the bias distribution curve from the calibrated Chinese buoy achieves much more consistency than the original data. The same conclusion for the bias percentage is also revealed in Figure 6e. The discrepancy above 3.5 m is mainly due to the limited number of collocated SWHs above 3.5 m. Significant consistency improvement in NRMSE is also shown in Figure 6f, although the calibration does not remove the inconsistency completely. In order to quantitatively measure the correlations between European/NDBC buoys and Chinese buoys with/without calibration, the RMSEs of the error parameters at each SWH ranges were introduced by using the error estimations of European/NDBC buoy SWH as a reference. A lower RMSE means a decrease in the difference between the Chinese and European/NDBC buoy. These RMSEs are indicated in Figure 6d–f. One can note that all the RMSEs of bias, bias percentage and NRMSE after the calibration significantly decreased compared to those without calibration, which are 0.097 m to 0.166 m, 3.1% to 8.1% and 4.1% to 8.4%. Therefore, the calibration of Chinese buoys does significantly improve consistency with the European/NDBC buoys.
The improvement in the consistency distribution with the mean wave period was also investigated, and the results are shown in Figure 7. The calibration method works well for the bias distribution with the mean wave period. The bias of HY2B against calibrated Chinese buoys showed good consistency with European/NDBC buoys in both short and long mean wave periods, which indicate dominant windsea and swell wave regimes, respectively.

5. Conclusions and Discussion

5.1. Conclusions

Wave buoy data are the most commonly used observations to verify the performance of wave models in different ocean regions. For global wave forecasting, there is a need to obtain consistent SWHs provided by different national buoy networks. Consistency among buoys is extremely important because the insitu SWH from a buoy is always seen as the observed truth. Therefore, the internal calibrations of the buoys are performed to maintain the consistency of the whole buoy network, which usually contains different types of buoys. However, because of the limitations of cost and implementation, comparisons or cross-calibrations are rarely performed between different national buoy networks. Such comparisons generally require simultaneous parallel observation experiments using buoys from different countries. The cross-calibration is critical because SWH inconsistency between buoy networks would lead to serious misinterpretation in the conclusions of regional wave studies.
Unfortunately, significant inconsistency between the buoy networks is found during the validation of operational altimeter SWH observations. The altimeters obtain different bias, RMSE and SI against the Chinese buoy network compared to those against the NDBC and European buoy networks. The large weight and inadequate calibration of Chinese buoys may lead to negative bias with respect to the other buoy networks, especially under low sea state or wind sea dominant conditions.
A remote cross-calibration method has been implemented based on the collocations between altimeters and buoys from different networks. The collocated altimeter data were used to remotely connect the buoy networks. The method has the advantage of matching the buoys independent of the altimeter accuracy. The cross-calibration of Chinese buoys was achieved by least squares regression based on scatter diagrams from the buoy matchups with Jason-3, Sentinel-3A and SARAL. The cross-calibration was tested for the validation of SWH from altimeter HY2B. This test shows that the negative SWH bias of Chinese buoys is effectively removed, indicating a significant improvement in the consistency among the European/NDBC and Chinese buoy networks.

5.2. Discussion

The method presented in this study will be useful for applications analyzing the regional wave climate in the China Sea and establishing a benchmark for indicators of sea state forecasts. For operational use, such calibration is relevant to favor the assimilation of SWHs from the Chinese network into coastal wave models and thus improve the reliability of marine warning bulletins during the typhoon season.
Finally, we would like to emphasize that the cross-calibration of different national buoy networks still requires further experiments and studies.The remote cross-calibration method in this paper is cost effective but must be updated regularly by adequate collocations between buoys and altimeters. This method is suitable for estimating the consistency between buoy systems. Most importantly, the routine calibration or maintenance of buoys is critical to preserve their consistencies and to avoid serious issues occurring in Chinese buoys due to inadequate studies or calibrations. Among the methods of calibration, the parallel observation experiment is still the most direct way to obtain detailed cross-calibration because a buoy network always contains buoys of multiple types; of course, massive support would be needed.

Author Contributions

Conceptualization, J.W. (Jiuke Wang) and L.A.; methodology and formal analysis, J.W. (Jiuke Wang); data acquisition and preprocessing, X.W. and J.W. (Juanjuan Wang); writing—original draft preparation, J.W. (Jiuke Wang); writing—review and editing, L.A.; project administration, B.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Key Research and Development Program of China (grant No. 2016YFC1401701 and No.2017YFC1403306) and the Laboratory for Regional Oceanography and Numerical Modeling, Pilot National Laboratory for Marine Science and Technology (grant No.2019B03).

Acknowledgments

The authors would like to thank the National Data Buoy Center for providing the buoy data. The authors would like to thank Christophe Bataille for providing the European buoy data from the GTS system.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Jiang, H.; Zheng, H.; Mu, L. Improving Altimeter Wind Speed Retrievals Using Ocean Wave Parameters. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2020, 13, 1917–1924. [Google Scholar] [CrossRef]
  2. Wang, X.; Yao, C.; Gao, G.; Jiang, H.; Xu, D.; Chen, G.; Zhang, Z. Simulating tropical cyclone waves in the East China Sea with an event-based, parametric-adjusted model. J. Oceanogr. 2020, 1–19. [Google Scholar] [CrossRef]
  3. Crosby, S.C.; O’Reilly, W.C.; Guza, R.T. Modeling long-period swell in Southern California: Practical boundary conditions from buoy observations and global wave model predictions. J. Atmos. Ocean. Technol. 2016, 33, 1673–1690. [Google Scholar] [CrossRef]
  4. Doble, M.J.; Bidlot, J.-R. Wave buoy measurements at the Antarctic sea ice edge compared with an enhanced ECMWF WAM: Progress towards global waves-in-ice modelling. Ocean Model. 2013, 70, 166–173. [Google Scholar] [CrossRef]
  5. Gemmrich, J.; Thomas, B.; Bouchard, R. Observational changes and trends in northeast Pacific wave records. Geophys. Res. Lett. 2011, 38. [Google Scholar] [CrossRef]
  6. Meindl, A. Guide to Moored Buoys and Other Ocean Data Acquisition Systems DBCP Technical Document No. 8. Intergov. Oceanogr. Comm. World Meteorol. Organ. 1996, 3. [Google Scholar]
  7. Mettlach, T.; Teng, C.-C. Concepts for an ideal ocean wave-measuring buoy. In Proceedings of the OCEANS 2010 MTS/IEEE SEATTLE, Seattle, WA, USA, 20–23 September 2010; pp. 1–8. [Google Scholar]
  8. Turton, J.; Fenna, P. Observations of extreme wave conditions in the north-east Atlantic during December 2007. Weather 2008, 63, 352–355. [Google Scholar] [CrossRef]
  9. Hall, C.; Bouchard, R.H.; Riley, R.; Stewart, R.; Wang, D.; DiNapoli, S. Emerging National Data Buoy Center (NDBC) Wave Systems. In Proceedings of the Technical Report Presentation at DBCP-34 Meeting, Cape Town, South Africa, 22 October 2018. [Google Scholar]
  10. Teng, C.-C.; Bernard, L.; Taft, B.; Burdette, M. A compact wave and ocean data buoy system. In Proceedings of the Proceedings of OCEANS 2005 MTS/IEEE, Washington, DC, USA, 17–23 September 2005; pp. 1249–1254. [Google Scholar]
  11. Mettlach, T.; Teng, C.-C.; Weir, R.; LaRue, K. Quality of archived NDBC data as climate records. In Proceedings of the OCEANS, Quebec City, QC, Canada, 15–18 September 2008; pp. 1–10. [Google Scholar]
  12. Anderson, G.; Carse, F.; Turton, J.; Saulter, A. Quantification of bias of wave measurements from lightvessels. J. Oper. Oceanogr. 2016, 9, 93–102. [Google Scholar] [CrossRef]
  13. Islam, H.; Campos, R.M.; Ferreira, T.R.S.; Soares, C.G. Hydrodynamic assessment of a biofouled wave buoy in coastal zone. In Proceedings of the 39th International Conference on Ocean, Offshore & Arctic Engineering—OMAE 2020, Virtual, Online, 3–7 August 2020. [Google Scholar]
  14. Wang, J.; Li, B.; Gao, Z.; Wang, J. Comparison of ECMWF significant wave height forecasts in the China sea with buoy data. Weather Forecast. 2019, 34, 1693–1704. [Google Scholar] [CrossRef]
  15. Bender III, L.C.; Guinasso, N.L., Jr.; Walpert, J.N.; Howden, S.D. A comparison of methods for determining significant wave heights—Applied to a 3-m discus buoy during Hurricane Katrina. J. Atmos. Ocean. Technol. 2010, 27, 1012–1028. [Google Scholar] [CrossRef]
  16. Nerem, R.S.; Chambers, D.P.; Choe, C.; Mitchum, G.T. Estimating mean sea level change from the TOPEX and Jason altimeter missions. Mar. Geod. 2010, 33, 435–446. [Google Scholar] [CrossRef]
  17. Lee, H.; Shum, C.K.; Emery, W.; Calmant, S.; Deng, X.; Kuo, C.-Y.; Roesler, C.; Yi, Y. Validation of Jason-2 altimeter data by waveform retracking over California coastal ocean. Mar. Geod. 2010, 33, 304–316. [Google Scholar] [CrossRef]
  18. Yang, J.; Zhang, J.; Jia, Y.; Fan, C.; Cui, W. Validation of Sentinel-3A/3B and Jason-3 Altimeter Wind Speeds and Significant Wave Heights Using Buoy and ASCAT Data. Remote Sens. 2020, 12, 2079. [Google Scholar] [CrossRef]
  19. Liu, Q.; Babanin, A.V.; Guan, C.; Zieger, S.; Sun, J.; Jia, Y. Calibration and validation of HY-2 altimeter wave height. J. Atmos. Ocean. Technol. 2016, 33, 919–936. [Google Scholar] [CrossRef] [Green Version]
  20. Jiang, M.; Xu, K.; Liu, Y. Global statistical assessment and cross-calibration with Jason-2 for reprocessed HY-2 A altimeter data. Mar. Geod. 2018, 41, 289–312. [Google Scholar] [CrossRef]
  21. Lionello, P.; Günther, H.; Janssen, P.A.E.M. Assimilation of altimeter data in a global third-generation wave model. J. Geophys. Res. Ocean. 1992, 97, 14453–14474. [Google Scholar] [CrossRef]
  22. Bauer, E.; Hasselmann, S.; Hasselmann, K.; Graber, H.C. Validation and assimilation of Seasat altimeter wave heights using the WAM wave model. J. Geophys. Res. Ocean. 1992, 97, 12671–12682. [Google Scholar] [CrossRef]
  23. Breivik, L.-A.; Reistad, M. Assimilation of ERS-1 altimeter wave heights in an operational numerical wave model. Weather Forecast. 1994, 9, 440–451. [Google Scholar] [CrossRef] [Green Version]
  24. Ribal, A.; Young, I.R. 33 years of globally calibrated wave height and wind speed data based on altimeter observations. Sci. Data 2019, 6, 1–15. [Google Scholar]
  25. Wang, J.; Aouf, L.; Jia, Y.; Zhang, Y. Validation and Calibration of Significant Wave Height and Wind Speed Retrievals from HY2B Altimeter Based on Deep Learning. Remote Sens. 2020, 12, 2858. [Google Scholar] [CrossRef]
  26. Jia, Y.; Yang, J.; Lin, M.; Zhang, Y.; Ma, C.; Fan, C. Global Assessments of the HY-2B Measurements and Cross-Calibrations with Jason-3. Remote Sens. 2020, 12, 2470. [Google Scholar] [CrossRef]
Figure 1. The locations of buoys from the three networks used in this paper. The locations of the NDBC, European and Chinese buoys are indicated in subfigure (ac), respectively.
Figure 1. The locations of buoys from the three networks used in this paper. The locations of the NDBC, European and Chinese buoys are indicated in subfigure (ac), respectively.
Remotesensing 12 03447 g001
Figure 2. The validation of Jason-3 SWH against buoys from the three buoy networks. (ac) indicate the scatters and the statistical parameters; (d) indicates the bias distribution with the SWH value; (e) indicates the percentage of bias in each SWH level; (f) shows the NRMSE distribution with the significant wave height (SWH) value. Among (df), the legends C, E and N represent the Chinese network, European network and theNDBC network, respectively.
Figure 2. The validation of Jason-3 SWH against buoys from the three buoy networks. (ac) indicate the scatters and the statistical parameters; (d) indicates the bias distribution with the SWH value; (e) indicates the percentage of bias in each SWH level; (f) shows the NRMSE distribution with the significant wave height (SWH) value. Among (df), the legends C, E and N represent the Chinese network, European network and theNDBC network, respectively.
Remotesensing 12 03447 g002
Figure 3. Validation of buoy SWH observations against Jason-3 SWH. (a) indicates the biases of the three buoy networks distributed with the mean wave period; (b) indicates the bias percentage distribution with the mean wave period; (c) indicates the NRMSE distribution with the mean wave period. The legends C, E and N represent the Chinese network, European network and NDBC, respectively.
Figure 3. Validation of buoy SWH observations against Jason-3 SWH. (a) indicates the biases of the three buoy networks distributed with the mean wave period; (b) indicates the bias percentage distribution with the mean wave period; (c) indicates the NRMSE distribution with the mean wave period. The legends C, E and N represent the Chinese network, European network and NDBC, respectively.
Remotesensing 12 03447 g003
Figure 4. The SWH distribution of Chinese buoys and NDBC buoys when Jason-3 SWH is 0.7 m.
Figure 4. The SWH distribution of Chinese buoys and NDBC buoys when Jason-3 SWH is 0.7 m.
Remotesensing 12 03447 g004
Figure 5. Scatters of the Chinese buoys and European/NDBC buoys. Circles, crosses and stars indicate the matchups from Jason-3, Sentinel-3A and SARAL, respectively. The dashed line represents the least squares method fitting from the scatters from all three altimeters.
Figure 5. Scatters of the Chinese buoys and European/NDBC buoys. Circles, crosses and stars indicate the matchups from Jason-3, Sentinel-3A and SARAL, respectively. The dashed line represents the least squares method fitting from the scatters from all three altimeters.
Remotesensing 12 03447 g005
Figure 6. The validation of HY2B SWH against buoys from the calibrated and original Chinese buoys and European/NDBC buoys. (ac) indicate the scatters and the statistical parameters; (d) indicates the bias distribution with the SWH value; (e) indicates the percentage of bias in each SWH range; (f) shows the NRMSE distribution with the SWH value. Among (df), the legends C, CC and ALL represent Original Chinese buoy SWH, Calibrated Chinese buoy SWH and European/NDBC buoy SWH, respectively.
Figure 6. The validation of HY2B SWH against buoys from the calibrated and original Chinese buoys and European/NDBC buoys. (ac) indicate the scatters and the statistical parameters; (d) indicates the bias distribution with the SWH value; (e) indicates the percentage of bias in each SWH range; (f) shows the NRMSE distribution with the SWH value. Among (df), the legends C, CC and ALL represent Original Chinese buoy SWH, Calibrated Chinese buoy SWH and European/NDBC buoy SWH, respectively.
Remotesensing 12 03447 g006
Figure 7. Validation of buoy SWH observations against HY2B SWH. (a) indicates the biases of the original and calibrated Chinese and European/NDBC buoy networks distributed with mean wave period; (b) indicates the bias percentage distributions with the mean wave period; (c) indicates the NRMSE distributions with mean wave period. The legends C, CC and ALL represent original Chinese buoy SWH, Calibrated Chinese buoy SWH and European/NDBC buoy SWH, respectively.
Figure 7. Validation of buoy SWH observations against HY2B SWH. (a) indicates the biases of the original and calibrated Chinese and European/NDBC buoy networks distributed with mean wave period; (b) indicates the bias percentage distributions with the mean wave period; (c) indicates the NRMSE distributions with mean wave period. The legends C, CC and ALL represent original Chinese buoy SWH, Calibrated Chinese buoy SWH and European/NDBC buoy SWH, respectively.
Remotesensing 12 03447 g007
Table 1. Validations of altimeter SWHs against the buoy networks.
Table 1. Validations of altimeter SWHs against the buoy networks.
AltimeterJason-3Sentinel-3ASARALHY2B
Bias-NDBC0.0470.1170.0800.135
Bias-European0.0200.1390.1160.126
Bias-Chinese0.1680.1900.1710.232
NRMSE-NDBC12.2%13.6%11.9%12.3%
NRMSE-European11.5%13.0%13.2%10.9%
NRMSE-Chinese25.1%23.9%21.5%24.1%
SI-NDBC12.0%12.0%11.1%9.4%
SI-European11.5%11.8%12.3%10.8%
SI-Chinese20.6%18.2%16.1%15.4%
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wang, J.; Aouf, L.; Wang, X.; Li, B.; Wang, J. Remote Cross-Calibration of Wave Buoys Based on Significant Wave Height Observations of Altimeters in the Northern Hemisphere. Remote Sens. 2020, 12, 3447. https://0-doi-org.brum.beds.ac.uk/10.3390/rs12203447

AMA Style

Wang J, Aouf L, Wang X, Li B, Wang J. Remote Cross-Calibration of Wave Buoys Based on Significant Wave Height Observations of Altimeters in the Northern Hemisphere. Remote Sensing. 2020; 12(20):3447. https://0-doi-org.brum.beds.ac.uk/10.3390/rs12203447

Chicago/Turabian Style

Wang, Jiuke, Lotfi Aouf, Xianqiao Wang, Benxia Li, and Juanjuan Wang. 2020. "Remote Cross-Calibration of Wave Buoys Based on Significant Wave Height Observations of Altimeters in the Northern Hemisphere" Remote Sensing 12, no. 20: 3447. https://0-doi-org.brum.beds.ac.uk/10.3390/rs12203447

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop