Sea surface wind speed retrieval from Sentinel-1 HH polarization data using conventional and neural network methods
-
Abstract: Conventional retrieval and neural network methods are used simultaneously to retrieve sea surface wind speed (SSWS) from HH-polarized Sentinel-1 (S1) SAR images. The Polarization Ratio (PR) models combined with the CMOD5.N Geophysical Model Function (GMF) is used for SSWS retrieval from the HH-polarized SAR data. We compared different PR models developed based on previous C-band SAR data in HH-polarization for their applications to the S1 SAR data. The recently proposed CMODH, i.e., retrieving SSWS directly from the HH-polarized S1 data is also validated. The results indicate that the CMODH model performs better than results achieved using the PR models. We proposed a neural network method based on the backward propagation (BP) neural network to retrieve SSWS from the S1 HH-polarized data. The SSWS retrieved using the BP neural network model agrees better with the buoy measurements and ASCAT dataset than the results achieved using the conventional methods. Compared to the buoy measurements, the bias, root mean square error (RMSE) and scatter index (SI) of wind speed retrieved by the BP neural network model are 0.10 m/s, 1.38 m/s and 19.85%, respectively, while compared to the ASCAT dataset the three parameters of training set are –0.01 m/s, 1.33 m/s and 15.10%, respectively. It is suggested that the BP neural network model has a potential application in retrieving SSWS from Sentinel-1 images acquired at HH-polarization.
-
Key words:
- Sentinel-1 /
- HH-polarization /
- sea surface wind speed /
- retrieval methods
-
Figure 7. An example of SSWS retrieval for Sentinel-1 HH polarization data using Mouche-PR1 model (a), Mouche-PR2 model (b), Zhang-PR model (c), Liu-PR model (d), CMODH model (e), and BP neural network model (f). The collocated ERA5 wind vectors within the coverage of the Sentinel-1 image are overlaid. The red pentagram marks the location of the R/V Xuelong.
Table 1. Technical specifications of Sentinel-1 IW and EW modes
Mode Incidence angle/(°) Nominal resolution Swath width/km Polarization IW 29–46 5 m×20 m 250 HH+HV, HH, VH+VV, VV EW 19–47 25 m×40 m 400 HH+HV, HH, VH+VV, VV Table 2. Information table of buoy data used in this study
Station Latitude Longitude Height/m Station Latitude Longitude Height/m 46022 40.720°N 124.531°W 4 46047 32.398°N 119.498°W 5 46014 39.233°N 123.967°W 4 46029 46.143°N 124.485°W 5 46013 38.238°N 123.307°W 4 46041 47.353°N 124.742°W 5 46015 42.779°N 124.874°W 4 51003 19.289°N 160.569°W 5 46050 44.677°N 124.515°W 4 51000 23.535°N 153.781°W 5 46092 36.751°N 122.029°W 4 44258 44.500°N 63.400°W 5 46025 33.749°N 119.053°W 4 45138 49.540°N 65.710°W 5 51004 17.602°N 152.395°W 4 SMKF1 24.628°N 81.109°W 6.1 32303 5.000°N 95.000°W 4 VCAF1 24.711°N 81.107°W 6.5 43301 8.000°N 95.000°W 4 NPSF1 26.132°N 81.807°W 6.5 51002 17.037°N 157.696°W 4 MTBF1 27.661°N 82.594°W 6.7 46035 57.026°N 177.738°W 5 LONF1 24.844°N 80.864°W 7 46027 41.852°N 124.382°W 5 VCVA2 57.125°N 170.285°W 8.5 46089 45.925°N 125.771°W 5 CDRF1 29.136°N 83.029°W 10 46011 34.956°N 121.019°W 5 VENF1 27.072°N 82.453°W 11.6 46028 35.712°N 121.858°W 5 SANF1 24.456°N 81.877°W 14.6 46042 36.785°N 122.398°W 5 MLRF1 25.012°N 80.376°W 15.8 46069 33.674°N 120.212°W 5 PLSF1 24.693°N 82.773°W 17.7 46086 32.491°N 118.035°W 5 Table 3. Coefficients of Mouche-PR1 model
Coefficient Value A0 0.006 507 04 B0 0.128 983 00 C0 0.992 839 00 Aπ/2 0.007 821 94 Bπ/2 0.121 405 00 Cπ/2 0.992 839 00 Aπ 0.005 984 16 Bπ 0.140 952 00 Cπ 0.992 885 00 Table 4. Coefficients of PR model with incidence angle dependence
PR model A B C Mouche-PR2 0.00799793 0.125465 0.997379 Zhang-PR 0.282 8 0.045 1 0.289 1 Liu-PR 0.453 041 0 0.0324573 0.5243030 Table 5. Statistical parameters of SSWS retrieval using different methods
Models Bias/(m·s–1) RMSE/(m·s–1) SI/% Compared with buoy Mouche-PR1 0.11 1.60 22.89 Mouche-PR2 0.13 1.63 23.42 Zhang-PR –0.33 1.50 21.09 Liu-PR 0.25 1.78 25.36 CMODH –0.31 1.48 20.89 BPNN 0.10 1.38 19.85 Compared with ASCAT Mouche-PR1 0.27 1.57 17.73 Mouche-PR2 0.41 1.70 18.85 Zhang-PR –0.09 1.40 16.26 Liu-PR 0.91 2.52 26.86 CMODH –0.35 1.39 15.82 BPNN (training set) 0.03 1.33 15.39 BPNN (test set) –0.01 1.33 15.10 -
[1] Bentamy A, Croize-Fillon D, Perigaud C. 2008. Characterization of ASCAT measurements based on buoy and QuikSCAT wind vector observations. Ocean Science, 4(4): 265–274. doi: 10.5194/os-4-265-2008 [2] Elfouhaily T. 1996. Physical modeling of electromagnetic backscatter from the ocean surface; application to retrieval of wind fields and wind stress by remote sensing of the marine atmospheric boundary layer [dissertation]. Paris: University Paris VII [3] ESA (European Space Agency). 2016. Sentinel-1 product specification. https://sentinel.esa.int/web/sentinel/document-library/content/-/article/sentinel-1-product-specification, [2020-2-17] [4] Hersbach H. 2010. Comparison of C-band scatterometer CMOD5. N equivalent neutral winds with ECMWF. Journal of Atmospheric and Oceanic Technology, 27(4): 721–736. doi: 10.1175/2009JTECHO698.1 [5] Hersbach H, Stoffelen A, de Haan S. 2007. An improved C‐band scatterometer ocean geophysical model function: CMOD5. Journal of Geophysical Research, 112(C3): C03006 [6] Horstmann J, Schiller H, Schulz-Stellenfleth J, et al. 2003. Global wind speed retrieval from SAR. IEEE Transactions on Geoscience and Remote Sensing, 41(10): 2277–2286. doi: 10.1109/TGRS.2003.814658 [7] Komarov A S, Buehner M. 2017. Automated detection of ice and open water from dual-polarization RADARSAT-2 images for data assimilation. IEEE Transactions on Geoscience and Remote Sensing, 55(10): 5755–5769. doi: 10.1109/TGRS.2017.2713987 [8] Liu Guihong, Yang Xiaofeng, Li Xiaofeng, et al. 2013. A systematic comparison of the effect of polarization ratio models on sea surface wind retrieval from C-band synthetic aperture radar. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 6(3): 1100–1108. doi: 10.1109/JSTARS.2013.2242848 [9] Monaldo F M, Thompson D R, Pichel W G, et al. 2004. A systematic comparison of QuikSCAT and SAR ocean surface wind speeds. IEEE Transactions on Geoscience and Remote Sensing, 42(2): 283–291. doi: 10.1109/TGRS.2003.817213 [10] Moreira A. 1991. Improved multi-look techniques applied to SAR and Scan SAR imagery. IEEE Transactions on Geoscience and Remote Sensing, 29(4): 529–534. doi: 10.1109/36.135814 [11] Mouche A A, Hauser D, Daloze J F, et al. 2005. Dual-polarization measurements at C-band over the ocean: results from airborne radar observations and comparison with ENVISAT ASAR data. IEEE Transactions on Geoscience and Remote Sensing, 43(4): 753–769. doi: 10.1109/TGRS.2005.843951 [12] Olivier P, Vidal-Madjar D. 1994. Empirical estimation of the ERS-1 SAR radiometric resolution. International Journal of Remote Sensing, 15(5): 1109–1114. doi: 10.1080/01431169408954144 [13] Peixoto J P, Oort A H. 1992. Physics of Climate. New York: American Institute of Physics, 67 [14] Pierson Jr W J. 1990. Examples of, reasons for, and consequences of the poor quality of wind data from ships for the marine boundary layer: Implications for remote sensing. Journal of Geophysical Research, 95(C8): 13313–13340. doi: 10.1029/JC095iC08p13313 [15] Pond S, Pickard G L. 1983. Currents with friction; wind-driven circulation. In: Pond S, Pickard G, eds. Introductory Dynamical Oceanography. 2nd ed. Amsterdam: Elsevier, 100–162 [16] Quilfen Y, Chapron B, Elfouhaily T, et al. 1998. Observation of tropical cyclones by high-resolution scatterometry. Journal of Geophysical Research, 103(C4): 7767–7786. doi: 10.1029/97JC01911 [17] Richaume P, Badran F, Crepon M, et al. 2000. Neural network wind retrieval from ERS-1 scatterometer data. Journal of Geophysical Research, 105(C4): 8737–8751. doi: 10.1029/1999JC900225 [18] Schwerdt M, Schmidt K, Tous Ramon N, et al. 2014. Independent verification of the Sentinel-1A system calibration. In: Proceedings of 2014 Geoscience and Remote Sensing Symposium. Quebec City, QC, Canada: IEEE [19] Stoffelen A, Anderson D. 1997. Scatterometer data interpretation: estimation and validation of the transfer function CMOD4. Journal of Geophysical Research, 102(C3): 5767–5780. doi: 10.1029/96JC02860 [20] Thiria S, Mejia C, Badran F, et al. 1993. A neural network approach for modeling nonlinear transfer functions: application for wind retrieval from spaceborne scatterometer data. Journal of Geophysical Research, 98(C12): 22827–22841. doi: 10.1029/93JC01815 [21] Wang Lihua, Lu Peng, Ma Jiapei. 2017. Deriving sea surface wind from synthetic aperture radar based on Fourier transform and neural network. In: Proceedings of the 10th International Congress on Image and Signal Processing. Shanghai: IEEE, 1–6 [22] Yang Yonghong. 2009. Introduction to Synthetic Aperture Radar Ocean Remote Sensing (in Chinese). Beijing: China Ocean Press [23] Zhang Biao, Mouche A, Lu Yiru, et al. 2019. A geophysical model function for wind speed retrieval from C-band HH-polarized synthetic aperture radar. IEEE Geoscience and Remote Sensing Letters, 16(10): 1521–1525. doi: 10.1109/LGRS.2019.2905578 [24] Zhang Biao, Perrie W, He Yijun. 2011. Wind speed retrieval from RADARSAT-2 quad-polarization images using a new polarization ratio model. Journal of Geophysical Research, 116(C8): C08008 [25] Zhang Biao, Perrie W, Hwang P A, et al. 2010. A new polarization ratio model from C-band RADARSAT-2 fine Quad-Pol imagery. In: Proceedings of 2010 IEEE International Geoscience and Remote Sensing Symposium. Honolulu, HI, USA: IEEE, 1948–1951