Assimilation of ASCAT Sea Surface Wind Retrievals with Correlated Observation Errors

考虑观测误差相关的ASCAT海表风场反演资料同化

+ Author Affiliations + Find other works by these authors
  • Corresponding author: Weimin ZHANG, wmzhang104@139.com
  • Funds:

    Supported by the National Natural Science Foundation of China (41675097 and 41375113) and Key Research and Development Program of Hainan Province (ZDYF2017167)

  • doi: 10.1007/s13351-021-1007-0

PDF

  • Data assimilation systems usually assume that the observation errors of wind components, i.e., u (the longitudinal component) and v (the latitudinal component), are uncorrelated. However, since wind components are derived from observations in the form of wind speed and direction (spd and dir), the observation errors of u and v are correlated. In this paper, an explicit expression of the observation errors and correlation for each pair of wind components are derived based on the law of error propagation. The new data assimilation scheme considering the correlated error of wind components is implemented in the Weather Research and Forecasting Data Assimilation (WRFDA) system. Besides, adaptive quality control (QC) is introduced to retain the information of high wind-speed observations. Results from real data experiments assimilating the Advanced Scatterometer (ASCAT) sea surface winds suggest that analyses from the new data assimilation scheme are more reasonable compared to those from the conventional one, and could improve the forecasting of Typhoon Noru.

    资料同化系统通常假定纬向风分量u与经向风分量v的观测误差不相关。然而,由于风分量是从观测的风速和风向推导出来的,所以风分量的观测误差实际上是相关的。本文根据误差传播定律,导出了已知风速和风向误差的前提下,风分量的观测误差和相关系数的显示表达式,并在WRF模式同化系统(WRFDA)中实现了考虑风分量观测误差相关的直接同化。此外,为了更好的利用高风条件下的风场信息,引入了自适应的质量控制方案。基于ASCAT海表风场反演的初步实验结果表明,与常规方法相比,本文的方法更加合理,并且改进了台风“奥鹿”的预报效果。

  • 加载中
  • Fig. 1.  Region where the wind “truth” might be located when given the wind observation with uncertainties. The red arrow is the wind observation, while the blue ones are the corresponding wind components. The red shadow region is the place where the wind “truth” might be located given the uncertainties of the wind observation. Panel (a) is from the wind-vector perspective, in which the uncertainties of ${\rm{spd}}$ and ${\rm{dir}}$ (${\delta _{{\rm{spd}}}}$ and ${\delta _{{\rm{dir}}}}$) are known. Panel (b) is from the wind-components perspective, in which the uncertainties of $u$ and $v$ (${\delta _u}$ and ${\delta _v}$) are known.

    Fig. 2.  Standard errors and corresponding correlation of $u$ and $v$ within a certain range of the observed values: (a) standard errors of $u$, denoted as ${\sigma _u}$; (b) standard errors of $v$, denoted as ${\sigma _v}$; and (c) error correlation of $u$ and $v$, denoted as ${\rho _{u,v}}$.

    Fig. 3.  As in Fig. 2, but after variance inflation.

    Fig. 4.  ASCAT wind field of Typhoon Noru with a grid size of 25 km: (a) wind speed (spd) field, (b) standard deviation of $u$ component, (c) standard deviation of $v$ component, and (d) error correlation of wind components.

    Fig. 5.  Ensemble mean of the wind fields of forecasts and corresponding flow-dependent error of the $u$ and $v$ components: (a) ensemble average of the wind fields of forecasts, (b) standard deviation of the $u$ component, and (c) standard deviation of the $v$ component.

    Fig. 6.  QC using different assimilation schemes, in which the blue dots represent the wind observations that are successfully assimilated into the data assimilation system: (a) ($u$, $v$) with independent errors, (b) (spd, dir) with independent errors, and (c) ($u$, $v$) with correlated errors.

    Fig. 7.  Wind speed (spd) analysis errors of different experiments: (a) background, (b) ($u$, $v$) components with independent error, (c) (spd, dir) with independent error, and (d) ($u$, $v$) components with correlated error.

    Fig. 8.  Mean error and root-mean-square (RMS) of spd analysis errors of different experiments: (a) background, (b) ($u$, $v$) components with independent error, (c) (spd, dir) with independent error, and (d) ($u$, $v$) components with correlated error.

    Fig. 9.  As in Fig. 7, but for pressure analysis errors.

    Fig. 10.  As in Fig. 8, but for pressure analysis errors.

    Fig. 11.  The forecast (a) tracks and (b) track errors of Typhoon Noru (2017).

    Fig. 12.  (a) Minimum pressure and (b) maximum spd forecast errors of the typhoon eye.

    Table 1.  Different assimilation schemes for the ASCAT wind observation

    No.Assimilation schemeQC schemeObservation type
    1Control experimentNo QCNone
    2(u, v) with independent errorsAdaptive QCASCAT wind
    3(spd, dir) with independent errorsAdaptive QCASCAT wind
    4(u, v) with correlated errorsAdaptive QCASCAT wind
    Download: Download as CSV

    Table 2.  QC and number of iterations of minimization in different assimilation schemes

    No.Assimilation schemeObservation number before QCObservation number after QCIteration
    1(u, v) with independent errors1681159511
    2(spd, dir) with independent errors1681168111
    3(u, v) with correlated errors16811658 7
    Download: Download as CSV

    Table 3.  Mean error and RMS of analysis errors in different experiments

    No.Assimilation scheme${ {e} }_{{\rm{speed}}}^{{\rm{mean}}}$${ {e} }_{{\rm{speed}}}^{{\rm{RMS}}}$${ {e} }_{{\rm{pressure}}}^{{\rm{mean}}}$${ {e} }_{{\rm{pressure}}}^{{\rm{RMS}}}$
    1Control experiment4.77−6.31 −7.11 7.80
    2(u, v) with independent errors3.05 4.13 −6.42 7.15
    3(spd, dir) with independent errors5.33 7.32−10.1711.12
    4(u, v) with correlated errors1.60 3.74 −5.87 6.66
    Download: Download as CSV

    Table 4.  Mean error of the typhoon forecast in different experiments

    No.Assimilation scheme${ {\overline{e}_{ {\rm{track} } } } }$${ {\overline{e}_{ {\rm{speed} } } } }$${ {\overline{e}_{ {\rm{pressure} } } } }$
    1Control experiment78.44−15.0022.66
    2(u, v) with independent errors40.96 −7.3211.52
    3(spd, dir) with independent errors82.90−14.8120.91
    4(u, v) with correlated errors41.85 −6.33 9.53
    Download: Download as CSV
  • [1]

    Anderson, J., T. Hoar, K. Raeder, et al., 2009: The Data Assimilation Research Testbed: A community facility. Bull. Amer. Meteor. Soc., 90, 1283–1296. doi: 10.1175/2009BAMS2618.1.
    [2]

    Andersson, E., J. Haseler, P. Undén, et al., 1998: The ECMWF implementation of three-dimensional variational assimilation (3D-Var). III: Experimental results. Quart. J. Roy. Meteor. Soc., 124, 1831–1860. doi: 10.1002/qj.49712455004.
    [3]

    Atlas, R., and R. N. Hoffman, 2000: The use of satellite surface wind data to improve weather analysis and forecasting at the NASA Data Assimilation Office. Elsevier Oceanography Series, 63, 57–78. doi: 10.1016/S0422-9894(00)80005-7.
    [4]

    Atlas, R., R. N. Hoffman, S. M. Leidner, et al., 2001: The effects of marine winds from scatterometer data on weather analysis and forecasting. Bull. Amer. Meteor. Soc., 82, 1965–1990. doi: 10.1175/1520-0477(2001)082<1965:TEOMWF>2.3.CO;2.
    [5]

    Barker, D., X.-Y. Huang, Z. Q. Liu, et al., 2012: The Weather Research and Forecasting Model’s Community Variational/Ensemble Data Assimilation System: WRFDA. Bull. Amer. Meteor. Soc., 93, 831–843. doi: 10.1175/BAMS-D-11-00167.1.
    [6]

    Barker, D. M., W. Huang, Y.-R. Guo, et al., 2004: A three-dimensional variational data assimilation system for MM5: Implementation and initial results. Mon. Wea. Rev., 132, 897–914. doi: 10.1175/1520-0493(2004)132<0897:ATVDAS>2.0.CO;2.
    [7]

    Bi, L., J. A. Jung, M. C. Morgan, et al., 2011: Assessment of assimilating ASCAT surface wind retrievals in the NCEP Global Data Assimilation System. Mon. Wea. Rev., 139, 3405–3421. doi: 10.1175/2011MWR3391.1.
    [8]

    Bonavita, M., M. Dahoui, P. Lopez, et al., 2017: On the Initialization of Tropical Cyclones. ECMWF Technical Memorandum No. 810, ECMWF, Shinfield Park, Reading, 39 pp.
    [9]

    Candy, B., 2001: The Assimilation of Ambiguous Scatterometer Winds Using a Variational Technique: Method and Forecast Impact. Forecasting Research Technical Report No. 349, Met Office, Exeter, UK, 23 pp. Available online at https://library.metoffice.gov.uk/Portal/Default/en-GB/RecordView/Index/46757. Accessed on 24 May 2021.
    [10]

    Dee, D. P., S. M. Uppala, A. J. Simmons, et al., 2011: The ERA-Interim reanalysis: configuration and performance of the data assimilation system. Quart. J. Roy. Meteor. Soc., 137, 553–597. doi: 10.1002/qj.828.
    [11]

    Dobson, F., L. Hasse, and R. Davis, 1980: Air–Sea Interaction: Instruments and Methods. Springer, Boston, 1–4, doi: 10.1007/978-1-4615-9182-5.
    [12]

    Ebuchi, N., H. C. Graber, and M. J. Caruso, 2002: Evaluation of wind vectors observed by QuikSCAT/SeaWinds using ocean buoy data. J. Atmos. Oceanic Technol., 19, 2049–2062. doi: 10.1175/1520-0426(2002)019<2049:EOWVOB>2.0.CO;2.
    [13]

    Gao, F., X. Y. Zhang, N. A. Jacobs, et al., 2012: Estimation of TAMDAR observational error and assimilation experiments. Wea. Forecasting, 27, 856–877. doi: 10.1175/WAF-D-11-00120.1.
    [14]

    Hilton, F., A. Collard, V. Guidard, et al., 2009: Assimilation of IASI radiances at European NWP centres. Proc. Workshop on the Assimilation of IASI in NWP, ECMWF/EUMETSAT Numerical Weather Prediction-Satellite Application Facility (NWP-SAF), Reading, UK, 73 pp.
    [15]

    Huang, X.-Y., F. Gao, N. A. Jacobs, et al., 2013: Assimilation of wind speed and direction observations: A new formulation and results from idealised experiments. Tellus A, 65, 19936. doi: 10.3402/tellusa.v65i0.19936.
    [16]

    Isaksen, L., and A. Stoffelen, 2000: ERS scatterometer wind data impact on ECMWF’s tropical cyclone forecasts. IEEE Trans. Geosci. Remote Sens., 38, 1885–1892. doi: 10.1109/36.851771.
    [17]

    Lin, W. M., M. Portabella, A. Stoffelen, et al., 2017: On the improvement of ASCAT wind data assimilation in global NWP. Proc. 2017 IEEE International Geoscience and Remote Sensing Symposium, IEEE, Fort Worth, TX, 402–405, doi: 10.1109/IGARSS.2017.8126980.
    [18]

    Ochoa, B., and S. Belongie, 2006: Covariance propagation for guided matching. Proc. Workshop on Statistical Methods in Multi-Image and Video Processing, University of California, San Diego, USA, 12 pp.
    [19]

    Plant, W. J., 2000: Effects of wind variability on scatterometry at low wind speeds. J. Geophys. Res. Oceans, 105, 16,899–16,910. doi: 10.1029/2000JC900043.
    [20]

    Ricciardulli, L., and F. Wentz, 2012: Development of consistent geophysical model functions for different scatterometer missions: Ku and C-band. Proc. 2012 NASA International Ocean Vector Wind Science Team Meeting, NASA, Utrecht, Netherlands, 28 pp.
    [21]

    Ricciardulli, L., and F. J. Wentz, 2016: Remote Sensing Systems ASCAT C-2015 Daily Ocean Vector Winds on 0.25 Deg Grid, Version 02.1. Remote Sensing Systems, Santa Rosa. Available online at www.remss.com/missions/ascat. Accessed on 12 May 2021.
    [22]

    Schwartz, B., and S. G. Benjamin, 1995: A comparison of temperature and wind measurements from ACARS-equipped aircraft and rawinsondes. Wea. Forecasting, 10, 528–544. doi: 10.1175/1520-0434(1995)010<0528:ACOTAW>2.0.CO;2.
    [23]

    Stoffelen, A., and P. van Beukering, 1997: Implementation of Improved ERS Scatterometer Data Processing and Its Impact on HIRLAM Short Range Weather Forecasts. HIRLAM Technical Report 31, IMET, Dublin, Ireland, 77 pp. Available online at http://www.hirlam.org/index.php/hirlam-documentation/doc_download/1322-hirlam-technical-report-no-31. Accessed on 24 May 2021.
    [24]

    Stoffelen, A. C. M., and G. J. Cats, 1991: The impact of Seasat-A scatterometer data on high-resolution analyses and forecasts: The development of the QE II storm. Mon. Wea. Rev., 119, 2794–2802. doi: 10.1175/1520-0493(1991)119<2794:TIOASD>2.0.CO;2.
    [25]

    Wang, X. G., D. M. Barker, C. Snyder, et al., 2008: A hybrid ETKF–3DVAR data assimilation scheme for the WRF model. Part I: Observing system simulation experiment. Mon. Wea. Rev., 136, 5116–5131. doi: 10.1175/2008MWR2444.1.
    [26]

    Ying, M., W. Zhang, H. Yu, et al., 2014: An overview of the China Meteorological Administration tropical cyclone database. J. Atmos. Oceanic Technol., 31, 287–301. doi: 10.1175/JTECH-D-12-00119.1.
    [27]

    Zhang, F. Q., Y. H. Weng, J. A. Sippel, et al., 2009: Cloud-resolving hurricane initialization and prediction through assimilation of Doppler radar observations with an ensemble Kalman filter. Mon. Wea. Rev., 137, 2105–2125. doi: 10.1175/2009MWR2645.1.
  • Weimin ZHANG and Boheng DUAN.pdf

  • 加载中
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Assimilation of ASCAT Sea Surface Wind Retrievals with Correlated Observation Errors

    Corresponding author: Weimin ZHANG, wmzhang104@139.com
  • 1. College of Meteorology and Oceanology, National University of Defense Technology, Changsha 410073
  • 2. Laboratory of Software Engineering for Complex Systems, Changsha 410073
  • 3. State Key Laboratory of Remote Sensing Science, Chinese Academy of Sciences, Beijing 100101
  • 4. Hainan Key Laboratory of Earth Observation, Sanya 572029
  • 5. Beijing Institute of Applied Meteorology, Beijing 100029
Funds: Supported by the National Natural Science Foundation of China (41675097 and 41375113) and Key Research and Development Program of Hainan Province (ZDYF2017167)

Abstract: 

Data assimilation systems usually assume that the observation errors of wind components, i.e., u (the longitudinal component) and v (the latitudinal component), are uncorrelated. However, since wind components are derived from observations in the form of wind speed and direction (spd and dir), the observation errors of u and v are correlated. In this paper, an explicit expression of the observation errors and correlation for each pair of wind components are derived based on the law of error propagation. The new data assimilation scheme considering the correlated error of wind components is implemented in the Weather Research and Forecasting Data Assimilation (WRFDA) system. Besides, adaptive quality control (QC) is introduced to retain the information of high wind-speed observations. Results from real data experiments assimilating the Advanced Scatterometer (ASCAT) sea surface winds suggest that analyses from the new data assimilation scheme are more reasonable compared to those from the conventional one, and could improve the forecasting of Typhoon Noru.

考虑观测误差相关的ASCAT海表风场反演资料同化

资料同化系统通常假定纬向风分量u与经向风分量v的观测误差不相关。然而,由于风分量是从观测的风速和风向推导出来的,所以风分量的观测误差实际上是相关的。本文根据误差传播定律,导出了已知风速和风向误差的前提下,风分量的观测误差和相关系数的显示表达式,并在WRF模式同化系统(WRFDA)中实现了考虑风分量观测误差相关的直接同化。此外,为了更好的利用高风条件下的风场信息,引入了自适应的质量控制方案。基于ASCAT海表风场反演的初步实验结果表明,与常规方法相比,本文的方法更加合理,并且改进了台风“奥鹿”的预报效果。

    • Wind observations, especially those derived from satellites, are of great importance for improving the numerical prediction of tropical cyclones. Observation products from scatterometers constitute one of the most important sources of sea surface wind data. The ECMWF was the first numerical weather forecasting center in the world to successfully incorporate the first European Remote Sensing Satellite (ERS-1) scatterometer wind data into a global three-dimensional variational system (Andersson et al., 1998). Previous studies have shown that wind data derived from scatterometers are of great importance for weather forecasting and climate monitoring (e.g., Stoffelen and Cats, 1991; Stoffelen and van Beukering, 1997; Atlas and Hoffman, 2000; Isaksen and Stoffelen, 2000; Atlas et al., 2001; Candy, 2001; Lin et al., 2017). In particular, the usefulness of such data has been demonstrated in the prediction of tropical cyclones (Isaksen and Stoffelen, 2000) and extratropical cyclones (Stoffelen and van Beukering, 1997).

      Before assimilation, wind observations are usually first transformed into longitudinal and latitudinal components, $u$ and $v$, although most wind observations are observed or derived in the form of wind speed (${\rm{spd}}$) and wind direction (${\rm{dir}}$). Since the $u$ and $v$ components are derived from observations of ${\rm{spd}}$ and ${\rm{dir}}$, the observation errors of the $u$ and $v$ components are actually correlated and directly impacted by the errors of ${\rm{spd}}$ and ${\rm{dir}}$. However, most current data assimilation systems—for example, the Weather Research and Forecasting Data Assimilation (WRFDA) system (Barker et al., 2012)—assume for the sake of simplicity that the observation errors of wind components are uncorrelated (Huang et al., 2013). This simplification impairs the results of analysis since correlated observations bring in redundant information to the data assimilation. Huang et al. (2013) proposed a method for assimilating wind observations in their observed form, and thus the observation errors of spd and dir could be directly taken into account in the assimilation process. Although this method can help avoid the introduction of correlation into the wind components in the process of transformation, it may also introduce propagation error due to the observation operator that transforms background wind components to ${\rm{spd}}$ and ${\rm{dir}}$. This also impacts the process of quality control (QC), which is based on the departure between the observation and background. In light of these problems, it is necessary to find a way to estimate the errors and correlation of the wind components derived from ${\rm{spd}}$ and ${\rm{dir}}$, so as to make it possible to assimilate the correlated wind components directly into the data assimilation system.

      The rest of this paper is organized as follows. In Section 2, we discuss the error characteristics of wind observation from different viewpoints, and give a definite expression of the observation errors and correlation for each pair of wind components based on the law of error propagation. In Section 3, we describe the possibility of directly assimilating the correlated wind components. The results from real data experiments with Advanced Scatterometer (ASCAT) wind observations are given in Section 4. Finally, conclusions are drawn in Section 5.

    2.   Error characteristics of wind components
    • In most data assimilation systems, wind observations are assimilated in the form of $u$ and $v$ components. For one wind vector collected in the form of ${\rm{spd}}$ and ${\rm{dir}}$, the wind components can be acquired by using the equation

      $$ \begin{array}{l} u = - {\rm{spd}} \cdot \sin({\rm{dir}}), \\ v = - {\rm{spd}} \cdot \cos({\rm{dir}}) . \end{array} $$ (1)

      It is certain from Eq. (1) that both the ${\rm{spd}}$ error and the ${\rm{dir}}$ error can impact the uncertainties of $u$ and $v$. However, data assimilation systems (Barker et al., 2012) usually use the observation errors of $u$ and $v$ estimated only from the spd error. It is obvious that the errors of $u$ and $v$ are correlated to each other, but they are assumed to be independent in data assimilation systems.

      Next, we discuss the impact of considering the difference in the wind observation errors between the wind-vector (${\rm{spd}},\;{\rm{dir}}$) form and the wind-components ($u,v$) form. Figure 1 shows the region where the wind “truth” might be located when given the wind observation with uncertainties. Figure 1a is from the wind-vector perspective where the uncertainties (${\delta _{{\rm{spd}}}}$, ${\delta _{{\rm{dir}}}}$) of one wind-vector pair (${\rm{spd}}^{\rm{o}}$, ${\rm{dir}}^{\rm{o}}$) are given. Figure 1b, meanwhile, is from the wind-components perspective, which assumes that the uncertainties (${\delta _u}$, ${\delta _v}$) of the (${u^{\rm{o}}},{v^{\rm{o}}}$) are fixed and uncorrelated. For most observing systems acquiring wind information in the form of ${\rm{spd}}$ and ${\rm{dir}}$, Fig. 1a is more reasonable for the real world. Under this circumstance, the uncertainties of wind components are no longer fixed values, but vary along with the values of ${\rm{spd}}$ and ${\rm{dir}}$ when their uncertainties are given. As shown in Fig. 1a, the possible region where the wind “truth” might be located will change in different (${\rm{spd}}^{\rm{o}},\;{\rm{dir}}^{\rm{o}}$) pairs, as will the range of ${u^{\rm{t}}}$ and ${v^{\rm{t}}}$ (where the superscript ${\rm{t}}$ means “truth”). In other words, the real errors of wind components do not only depend on the values of ${\rm{spd}}$ error and ${\rm{dir}}$ error, but also on the observed values of ${\rm{spd}}$ and ${\rm{dir}}$, and the errors of each (${u^{\rm{o}}},{v^{\rm{o}}}$) pair should be correlated.

      Figure 1.  Region where the wind “truth” might be located when given the wind observation with uncertainties. The red arrow is the wind observation, while the blue ones are the corresponding wind components. The red shadow region is the place where the wind “truth” might be located given the uncertainties of the wind observation. Panel (a) is from the wind-vector perspective, in which the uncertainties of ${\rm{spd}}$ and ${\rm{dir}}$ (${\delta _{{\rm{spd}}}}$ and ${\delta _{{\rm{dir}}}}$) are known. Panel (b) is from the wind-components perspective, in which the uncertainties of $u$ and $v$ (${\delta _u}$ and ${\delta _v}$) are known.

    • In this subsection, we calculate the correlated errors of observed wind components based on the law of error propagation (Ochoa and Belongie, 2006). The law of error propagation is widely used in the area of mapping and surveying. It maps the error (or errors) from one variable (or variables) to another. A brief introduction to this theory is given as follows.

      Let ${{x}} \in {{\bf{R}}^n}$ be a random vector with mean ${\mu _x}$ and covariance matrix ${\Delta _x}$, and let $f:{{\bf{R}}^n} \to {{\bf{R}}^m}$ be a nonlinear function. Assume that ${{y}} = f({{x}})$, and ${{J}} \in {{\bf{R}}^{m \times n}}$ is the Jacobian matrix $\partial f/\partial {{x}}$ evaluated at ${\mu _x}$. Then, the random vector ${{y}} \in {{\bf{R}}^m}$ has mean ${\mu _y} \approx f({\mu _x})$ and covariance ${\Delta _y} \approx {{J}}{\Delta _x}{{{J}}^{\rm{T}}}$ [see Ochoa and Belongie (2006) for further details], where the superscript ${\rm{T}}$ denotes the transpose.

      For the wind vector transformation equation from (${\rm{spd}},\;{\rm{dir}}$) to ($u,v$) [see Eq. (1)], we have ${{x}} = {({\rm{spd}},\;{\rm{dir}})^{\rm{T}}}$ and ${{y}} = {(u,v)^{\rm{T}}}$, and the Jacobian matrix is

      $$\begin{aligned} {{J}} = \partial {{y}}/\partial {{x}} = & \left({\begin{array}{*{20}{c}} {\partial u/\partial {\rm{spd}}}&{\partial u/\partial {\rm{dir}}} \\ {\partial v/\partial {\rm{spd}}}&{\partial v/\partial {\rm{dir}}} \end{array}} \right) \\ = & \left({\begin{array}{*{20}{c}} { - \sin({\rm{dir}})}&{ - {\rm{spd}} \cdot \cos({\rm{dir}})} \\ { - \cos({\rm{dir}})}&{{\,\,\,\,\rm{spd}} \cdot \sin({\rm{dir}})} \end{array}} \right). \end{aligned} \hspace{38pt} $$ (2)

      Assume that $\sigma _{{\rm{spd}}}^2$ and $\sigma _{{\rm{dir}}}^2$ are the error variances of ${\rm{spd}}$ and ${\rm{dir}}$, respectively, and their correlation can be ignored. Then, the error covariance matrix of ${{x}}$ is

      $$ {\Delta _x} = \left({\begin{array}{*{20}{c}} {\sigma _{{\rm{spd}}}^2}&{} \\ {}&{\sigma _{{\rm{dir}}}^2} \end{array}} \right). \hspace{135pt}$$ (3)

      The covariance matrix of ${{y}}$ can thus be written as

      $$\begin{aligned} {\Delta _y} = & \left({\begin{array}{*{20}{c}} {\sigma _u^2}&{{\sigma _{v,u}}} \\ {{\sigma _{u,v}}}&{\sigma _v^2} \end{array}} \right) \approx \left({\begin{array}{*{20}{c}} { - \sin({\rm{dir}})}&{ - {\rm{spd}} \cdot \cos({\rm{dir}})} \\ { - \cos({\rm{dir}})}&{{\,\,\,\,\rm{spd}} \cdot \sin({\rm{dir}})} \end{array}} \right)\\ & \left({\begin{array}{*{20}{c}} {\sigma _{{\rm{spd}}}^2}&{} \\ {}&{\sigma _{{\rm{dir}}}^2} \end{array}} \right)\left({\begin{array}{*{20}{c}} { - \sin({\rm{dir}})}&{ - \cos({\rm{dir}})} \\ { - {\rm{spd}} \cdot \cos({\rm{dir}})}&{{\rm{spd}} \cdot \sin({\rm{dir}})} \end{array}} \right), \end{aligned}$$ (4)

      where $\sigma _u^2$ and $\sigma _v^2$ are the error variances of $u$ and $v$, respectively; and ${\sigma _{u,v}}$ equals ${\sigma _{v,u}}$, which is the error covariance. Given the error variances and the observed values of (${\rm{spd}},\;{\rm{dir}}$), the error variances and covariance of the corresponding ($u,v$) can be derived by using the equation above, which are

      $$\begin{array}{*{20}{c}} {\sigma _u^2 = {\sin^2}({\rm{dir}}) \cdot \sigma _{{\rm{spd}}}^2 + {\rm{sp}}{{\rm{d}}^2} \cdot {\cos^2}({\rm{dir}}) \cdot \sigma _{{\rm{dir}}}^2,} \\ {\sigma _v^2 = {\cos^2}({\rm{dir}}) \cdot \sigma _{{\rm{spd}}}^2 + {\rm{sp}}{{\rm{d}}^2} \cdot {\sin^2}({\rm{dir}}) \cdot \sigma _{{\rm{dir}}}^2,} \\ {{\,\,\,\sigma _{u,v}} = \sin({\rm{dir}}) \cdot \cos({\rm{dir}}) \cdot (\sigma _{{\rm{spd}}}^2 - {\rm{sp}}{{\rm{d}}^2} \cdot \sigma _{{\rm{dir}}}^2).} \end{array}$$ (5)

      The correlation coefficient ${\rho _{u,v}}$ can also be calculated by using the equation

      $${\,\,\,\,\,\,\rho _{u,v}} = \frac{{{\sigma _{u,v}}}}{{\sqrt {\sigma _u^2\sigma _v^2} }}. \hspace{110pt}$$ (6)

      To better demonstrate how the error characteristics of wind components would vary along with the observed values of wind, we calculate the standard deviations of the wind components within a certain range of the observed values. Furthermore, in this study, an assumption is made that the standard deviations of ${\rm{spd}}$ and ${\rm{dir}}$ are unchanged and are set to 2 m s−1 and 20°, respectively. The standard deviations of $u$ and $v$ along with their correlation are shown in Fig. 2; the observed wind vectors vary within the range of −25 to 25 m s−1 with an interval of 0.05 m s−1. Each point in the $(u,v)$ coordinate space corresponds to a $({\rm{spd}},\;{\rm{dir}})$ vector. We can see from Fig. 2a that the standard error of $u$ does not remain constant but grows rapidly in its vertical (latitudinal) direction. This is because the magnitude of the standard deviation is mainly determined by its second term, which changes rapidly with the spd if the value of $v$ is nonzero. The same is true for $v$ (see Fig. 2b). Furthermore, the errors of the wind components are highly correlated, especially in the diagonal direction, as shown in Fig. 2c. Errors of the $(u,v)$ pair are negatively correlated; while if the signs are opposite, the correlation is positive. The correlation is zero only when the value of one component is zero. However, the absolute value of the correlation would always be less than one.

      Figure 2.  Standard errors and corresponding correlation of $u$ and $v$ within a certain range of the observed values: (a) standard errors of $u$, denoted as ${\sigma _u}$; (b) standard errors of $v$, denoted as ${\sigma _v}$; and (c) error correlation of $u$ and $v$, denoted as ${\rho _{u,v}}$.

      Based on the above analysis, we can see that the new scheme improves the wind assimilation process from two aspects. On the one hand, it changes the fixed weight assigned to the wind observation based on the observation error in the assimilation process. The observation error of the $u$ wind component can reach 7 m s−1 when the spd exceeds 20 m s−1, as can be seen from Fig. 2a. Under strong convective weather conditions, such as in the area of the typhoon center, the spd is much higher, and so too is the uncertainty of the observed wind, meaning that reducing the weight of such wind observations is appropriate. In addition, consideration of the correlation between the observation errors of wind components can be seen as another correction to the observation weight. On the other hand, the increase in the magnitude of the observation error is equivalent to appropriately relaxing the threshold condition of QC. This ensures that more observations can be incorporated into the assimilation system, especially for strong convective weather conditions.

    3.   Assimilation of correlated wind components
    • After the new observation errors and corresponding correlation of the pair of $(u,v)$ components are obtained, the QC scheme must first be adjusted. The threshold for the criterion $y - h(x) < \alpha \cdot {(\sigma _{\rm{o}}^2 + \sigma _{\rm{b}}^2)^{1/2}}$ is no longer a fixed value, but varies with the observation error. Another important factor is the accuracy of the background error covariance, since it plays an important role in the QC procedure and assigning weight to the minimization of the cost function. For extreme weather conditions, background variances are typically more than 10 times larger than the local climatology (Bonavita et al., 2017). Thus, we use a flow-dependent background error variance based on the hybrid ensemble transform Kalman filter–three-dimensional variational data assimilation (ETKF–3DVAR) scheme developed for WRF (Wang et al., 2008), which will also make a difference to the QC. The cost function of the assimilation also needs to be adjusted since the error correlation between the pair of $(u,v)$ components is considered.

    • The mathematical formulation of the hybrid method in WRFDA is described in Wang et al. (2008). Therefore, we only give a brief introduction here. The final analysis increment $\delta {{x}}$ of the hybrid is a sum of two terms,

      $$\delta {{x}} = \delta {{{x}}_1} + \sum\limits_{k = 1}^K {({{{a}}_k} \text{◦} {{x}}_k^{\rm{e}})} , \hspace{80pt}$$ (7)

      where the first term $\delta {{{x}}_1}$ is the analysis increment associated with the 3DVAR static background covariance and the second term is the increment associated with the flow-dependent ensemble covariance. The vector ${{{a}}_k}$ denotes the extended control variable for the kth ensemble member. The symbol ◦ denotes the Schur product of the vectors ${{{a}}_k}$ and ${{x}}_k^{\rm{e}}$, and ${{x}}_k^{\rm{e}}$ is the kth perturbation normalized by $\sqrt {K - 1} $, where $K$ is the ensemble size:

      $$ \hspace{-52pt} {{x}}_k^{\rm{e}} = ({{{x}}_k} - \bar {{x}})/\sqrt {K - 1} . \hspace{100pt}$$ (8)

      Here, ${{{x}}_k}$ is the kth ensemble forecast and $\bar {{x}}$ is the mean of the K-member ensemble forecasts.

      The analysis increment $\delta {{x}}$ is obtained by minimizing the following cost function:

      $$\begin{aligned}[b] \!\!\!\!\!\! J(\delta {{x}}) = \; & {\beta _1} \cdot {J_{\rm{s}}}(\delta {{{x}}_1}) + {\beta _2} \cdot {J_{\rm{e}}}({{a}}) + {J_{\rm{o}}}(\delta {{x}}) \\ = \; & {\beta _1} \cdot \frac{1}{2}\delta {{x}}_1^{\rm{T}}{{{B}}^{ - 1}}\delta {{{x}}_1} + {\beta _2} \cdot \frac{1}{2}{{{a}}^{\rm{T}}}{{{A}}^{ - 1}}{{a}} \\ & + \frac{1}{2}{({{H}}\delta {{x}} - {{d}})^{\rm{T}}}{{{R_{\rm o}^{ - 1}}}}({{H}}\delta {{x}} - {{d}}). \end{aligned} \hspace{60pt} $$ (9)

      Compared to the normal 3DVAR cost function, the weighted sum of the terms ${J_{\rm{s}}}$ and ${J_{\rm{e}}}$ replaces the usual background term. Here, ${J_{\rm{s}}}$ is the traditional 3DVAR background term associated with the static covariance ${{B}}$, and ${J_{\rm{e}}}$ is the background term associated with the flow-dependent ensemble covariance. Coefficients ${\beta _1}$ and ${\beta _2}$ are the weights assigned to the static background-error covariance and the ensemble covariance and are constrained by

      $${\beta _1} + {\beta _2} = 1. \hspace{171pt}$$ (10)

      Matrix ${{B}}$ is the background error covariance matrix, and the superscripts $ - 1$ and ${\rm{T}}$ denote the inverse and adjoint, respectively, of a matrix or linear operator. Here, ${{a}}$ is formed by $K$ vectors ${{{a}}_k}$, in which ${{{a}}^{\rm{T}}} = ({{a}}_1^{\rm{T}},\;{{a}}_2^{\rm{T}},\; \cdots, $$\;{{a}}_k^{\rm{T}})$; ${{A}}$ is a block-diagonal matrix that defines the spatial covariance of ${{a}}$.

      Here, we mainly concentrate on the observation term ${J_{\rm{o}}}$, of which ${{d}} = {{y}} - H({{{x}}^{\rm{b}}})$ is the innovation vector and ${{y}}$ is the observation vector. In a hybrid data assimilation scheme, the background ${{{x}}^{\rm{b}}}$ is the ensemble mean forecast, $H$ is the nonlinear observation operator, ${{H}}$ is the linearized $H$, and ${{R_{\rm o}}}$ is the observation error covariance matrix.

      We only consider the correlation between the pair of wind components within one wind vector observation, and assume that errors between different wind vector observations are independent. Therefore, the minimization of the cost function of the observation can be calculated independently. For wind components, the cost function of one observation can be written as

      $${J_{\rm{o}}}(\delta {{x}}) = \frac{1}{2}(\delta u - {d_{{u^{\rm{o}}}}},\delta v - {d_{{v^{\rm{o}}}}}){\left({\begin{array}{*{20}{c}} {\sigma _u^2}&{{\sigma _{v,u}}} \\ {{\sigma _{u,v}}}&{\sigma _v^2} \end{array}} \right)^{ - 1}}\left({\begin{array}{*{20}{c}} {\delta u - {d_{{u^{\rm{o}}}}}} \\ {\delta v - {d_{{v^{\rm{o}}}}}} \end{array}} \right),$$ (11)

      where $\delta u = {u^{\rm{a}}} - {u^{\rm{b}}}$ and $\delta v = {v^{\rm{a}}} - {v^{\rm{b}}}$ are the analysis increments of $u$ and $v$, respectively; and ${d_{{u^{\rm{o}}}}} = {u^{\rm{o}}} - {u^{\rm{b}}}$ and ${d_{{v^{\rm{o}}}}} = $${v^{\rm{o}}} - {v^{\rm{b}}}$ are the innovations. Assuming that ${{X}} = {({x_1},{x_2})^{\rm{T}}}$, with ${x_1} = \delta u - {d_{{u^{\rm{o}}}}}$ and ${x_2} = \delta v - {d_{{v^{\rm{o}}}}}$, then the cost function ${J_{\rm{o}}}(\delta {{x}})$ is a quadratic function of vector ${{X}}$, which is

      $$\hspace{-135pt} {J_{\rm{o}}}(\delta {{x}}) = f({{X}}) = \frac{1}{2}{{{X}}^{\rm{T}}}{{A}}{{X}}, $$ (12)

      where ${{A}}$ is the inverse of the error covariance matrix. The condition that the quadratic function $f({{X}})$ has a minimum is that the matrix ${{A}}$ is a positive definite matrix, which is the same as saying that the inverse of ${{A}}$ is a positive definite matrix. From the equation above, we have

      $$\hspace{-135pt}{\rho _{u,v}} = - 1 < \frac{{{\sigma _{u,v}}}}{{\sqrt {\sigma _u^2\sigma _v^2} }} < 1, $$ (13)

      such that

      $$\hspace{-184pt} \sigma _{v,u}^2 < \sigma _u^2\sigma _v^2. $$ (14)

      Then, for the error covariance matrix, we have

      $$\hspace{-49pt} \sigma _u^2 > 0,{\kern 1pt} {\kern 1pt} {\kern 1pt} {\rm{and}}{\kern 1pt} {\kern 1pt} {\kern 1pt} \left| {\begin{array}{*{20}{c}} {\sigma _u^2}&{{\sigma _{v,u}}} \\ {{\sigma _{u,v}}}&{\sigma _v^2} \end{array}} \right| = \sigma _u^2\sigma _v^2 - \sigma _{v,u}^2 > 0. $$ (15)

      Therefore, the error covariance matrix is a positive definite matrix; that is, the quadratic function $f({{X}})$ has its minimum.

    • The observation error matrix ${{R}}_{\rm o}$ is a diagonal matrix under the independent error assumption, so the minimization of the cost function of wind components can be calculated separately. Unlike the traditional scheme, the minimization for the new scheme should be considered simultaneously since the wind components are correlated. The same is true for the gradient of the cost function, in which the gradient for the new scheme is a vector rather than a scalar, as shown by

      $$\hspace{-40pt} \nabla f({{X}}) = {{AX}} = \left({\begin{array}{*{20}{c}} {{g_1}} \\ {{g_2}} \end{array}} \right) = {\left({\begin{array}{*{20}{c}} {\sigma _u^2}&{{\sigma _{v,u}}} \\ {{\sigma _{u,v}}}&{\sigma _v^2} \end{array}} \right)^{ - 1}}\left({\begin{array}{*{20}{c}} {{x_1}} \\ {{x_2}} \end{array}} \right), $$ (16)

      where ${g_1}$ and ${g_2}$ denote the gradients of $u$ and $v$, respectively. The inverse of the error covariance matrix can be derived as follows:

      $$\hspace{-1pt} \begin{aligned}[b] \left({\begin{array}{*{20}{c}} {{g_1}} \\ {{g_2}} \end{array}} \right) = & \left({\begin{array}{*{20}{c}} {\dfrac{1}{{\sigma _u^2(1 - {\rho ^2})}}}&{ - \dfrac{\rho }{{{\sigma _u}{\sigma _v}(1 - {\rho ^2})}}} \\ { - \dfrac{\rho }{{{\sigma _u}{\sigma _v}(1 - {\rho ^2})}}}&{\dfrac{1}{{\sigma _v^2(1 - {\rho ^2})}}} \end{array}} \right)\left({\begin{array}{*{20}{c}} {{x_1}} \\ {{x_2}} \end{array}} \right)\\ =& \left({\begin{array}{*{20}{c}} {\dfrac{{\sigma _v^2{x_1} - \rho {\sigma _u}{\sigma _v}{x_2}}}{{\sigma _u^2\sigma _v^2(1 - {\rho ^2})}}} \\ {\dfrac{{\sigma _u^2{x_2} - \rho {\sigma _u}{\sigma _v}{x_1}}}{{\sigma _u^2\sigma _v^2(1 - {\rho ^2})}}} \end{array}} \right). \end{aligned} $$ (17)

      Here, $\rho $ is short for ${\rho _{u,v}}$. We can see from the equation above that when $\rho = 0$, the gradient vector becomes

      $$\left({\begin{array}{*{20}{c}} {{g_1}} \\ {{g_2}} \end{array}} \right) = \left({\begin{array}{*{20}{c}} {\dfrac{{{x_1}}}{{\sigma _u^2}}} \\ {\dfrac{{{x_2}}}{{\sigma _v^2}}} \end{array}} \right), \hspace{140pt} $$ (18)

      which is equal to the scheme with independent error assumption. However, when the absolute value of $\rho $ is close to one, the magnitude of the gradient value will be extremely large and slows down the convergence of the minimization process. To avoid this problem, we must limit the value of $\rho $ to a certain range. In this study, we adopt the following constraint:

      $$\hspace{-130pt} \bar \rho = \left\{ {\begin{array}{*{20}{c}} {\rho,{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\rm{if}}{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} |\rho | < 0.9} , \\ {0.9,{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\rm{if}}{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} |\rho | \geqslant 0.9} , \end{array}} \right. $$ (19)

      where $\bar \rho $ is the correlation we used in WRFDA. For the wind components with $|\rho | \geqslant 0.9$, we take a measure called variance inflation, which is widely used in satellite observation assimilation (Hilton et al., 2009). Then, the new variance ${{\overline\sigma} _x ^2}$ we used in WRFDA becomes

      $${ {\overline\sigma} _x ^2} = \left\{ {\begin{array}{*{20}{c}} {\sigma _x^2,{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt}{\kern 1pt} {\kern 1pt} {\rm{if}}\;|\rho | < 0.9}, \\ {\dfrac{{|\rho |}}{{0.9}} \cdot \sigma _x^2,\;\;{\rm{if}}\;|\rho | \geqslant 0.9} , \end{array}} \right. \hspace{80pt}$$ (20)

      where the subscript $x$ denotes $u$ or $v$. The reason for inflation is to reduce the weight of the observation in the analysis as a compensation for the reduction of correlation. The variance inflation guarantees that the matrix ${{A}}$ will still be a positive definite matrix.

      The new variances and correlation are shown in Fig. 3. The inflation has no impact on the low spd observations, because wind components with absolute values of correlation larger than 0.9 are mainly distributed in the region where both the absolute values of $u$ and $v$ are larger than 13 m s−1. Therefore, this strategy will make a difference to the wind observations in high-wind conditions, especially for tropical cyclones.

      Figure 3.  As in Fig. 2, but after variance inflation.

    • Conventional QC procedures based on selective blacklisting and first-guess rejections are too blunt to deal with the task of whether to retain an observation with a large initial departure from the model background (Bonavita et al., 2017). Observations with a large departure from the background do not mean that they have large errors. Sometimes, the background also has large errors, especially in extreme weather conditions. The strategy of QC should be considered carefully to assign a reasonable weight to the observation, since the observation error is varying with the observed value.

      For traditional QC using an error cutoff method, an observation that satisfies $y - h(x) < \alpha \cdot {(\sigma _{\rm{o}}^2 + \sigma _{\rm{b}}^2)^{1/2}}$ can be retained by the QC step, where $y$ is the observation and $h$ is the observation (or forward model) operator linking the model state $x$ to the observation; ${\sigma _{\rm{o}}}$ and ${\sigma _{\rm{b}}}$ represent the standard deviation of the observation error and background equivalent error, respectively; and $\alpha $ is the cutoff factor. For the WRFDA system, $\alpha $ is usually set to 5, which means an observation with departure more than five times the value of ${(\sigma _{\rm{o}}^2 + \sigma _{\rm{b}}^2)^{1/2}}$ will be rejected by QC. However, in extreme weather conditions like those of a typhoon, it is often the case that high spd observations near the typhoon center are screened out by such a QC scheme. This is usually because of the incorrect expression of background error covariance in the meteorological environment of the typhoon.

      In situations in which we can retain the high-wind-speed observational information, as well as the incremental 3DVAR within the bounds of validity of the tangent linear hypothesis (see Section 3.1), we use an adaptive QC scheme that was proposed by ECMWF (Bonavita et al., 2017). Here, we give the new QC scheme directly as

      $$ {{\overline\sigma} _{\rm{o}} ^2} = \left\{ {\begin{array}{*{20}{l}} {\sigma _{\rm{o}}^2}, &{{\rm{if}}\;{{[y - h(x)]}^2} \leqslant \sigma _{\rm{o}}^2 + \sigma _{\rm{b}}^2},\\ {{{[y - h(x)]}^2} - \sigma _{\rm{b}}^2}, &{{\rm{if}}\;{{[y - h(x)]}^2} > \sigma _{\rm{o}}^2 + \sigma _{\rm{b}}^2}, \end{array}} \right.$$ (21)

      where ${\overline {\sigma} _{\rm{o}} ^2}$ is the new observation error variance and $\sigma _{\rm{b}}^2$ is the flow-dependent variance of background equivalent $h(x)$. As can be seen from these equations, for small departures, the adjustment is not used. However, for extreme weather conditions where large observation departures occur, the observation errors will increase with the amplitude of the departures. In such a circumstance, regardless of whether the background errors are overestimated or underestimated, the new observation error variance will reduce the weight of the observation and make the analysis more robust.

    4.   Real data experiments
    • The wind observation data used in this study were produced by Remote Sensing Systems (Ricciardulli and Wentz, 2016), which were derived from ASCAT, which is one of the instruments carried on board the MetOp polar satellites operated by the European Organisation for the Exploitation of Meteorological Satellites (EUMETSAT) (Bi et al., 2011). The data files contain the measurements for one orbit of the satellite around the earth and are organized by observation cells that are perpendicular to the direction of the satellite travels.

      Scatterometer wind products have been used in many operations, such as those of ECMWF (Isaksen and Stoffelen, 2000), the UK Met Office (Candy, 2001), and the NCEP (Bi et al., 2011). ASCAT obtains three independent backscatter measurements using three different viewing directions separated by a short time delay. Then, the surface spd and dir are calculated using these “triplets” within the C-2015 geophysical model function (GMF; Ricciardulli and Wentz, 2012).

      The ASCAT wind product used in this study is from the MetOp-A satellite, with a grid size of 25 km. The spd error and dir error are assumed uncorrelated in this study and are set to 2 m s−1 and 20° for experimental use, respectively. Figure 4a shows the ASCAT wind field of the center of Typhoon Noru (2017) at 2100 UTC 25 July 2017.

      Figure 4.  ASCAT wind field of Typhoon Noru with a grid size of 25 km: (a) wind speed (spd) field, (b) standard deviation of $u$ component, (c) standard deviation of $v$ component, and (d) error correlation of wind components.

      Figure 4 also shows the new error characteristics of the $u$ and $v$ components, in which Figs. 4b, c present the standard deviations of the observation errors of $u$ and $v$, respectively, which are derived from the standard deviations of the ${\rm{spd}}$ and ${\rm{dir}}$; and Fig. 4d shows the corresponding correlation. We can see from Fig. 4 that the error standard variances of wind components are affected by both spd and dir, which is consistent with the theory presented in Section 2. However, the error correlation is highly direction-related, and it reaches its maximum when the $u$ and $v$ have the same amplitude.

    • This study uses Typhoon Noru (2017) as a numerical example, which generated in the Northwest Pacific Ocean on 20 July 2017 and further intensified into a typhoon on 23 July 2017. Here, we implement the new data assimilation scheme that considers the correlation of the u and v components and the adaptive QC in the WRFDA system. The WRFDA was developed by NCAR (Barker et al., 2004) and is a widely used operational system. It can produce a multivariate incremental analysis in the WRF model space (Zhang et al., 2009). The setup of the assimilation region is $80\; \times \;80$ with a resolution of 25 km, and the vertical discretization is 51 layers. The time of assimilation is 2100 UTC 25 July 2017. We use the FNL (final) global analysis data provided by NCEP as the initial field and boundary conditions, and take the ensemble mean of 40 forecasts [generated using the Data Assimilation Research Testbed (Anderson et al., 2009)] adjustment from 1800 UTC 23 to 2100 UTC 25 July 2017 as the background field of the assimilation system. After the assimilation, a 24-h forecast is made, which is a forecast to 2100 UTC 26 July 2017. A set of assimilation and comparison experiments are carried out, as shown in Table 1.

      No.Assimilation schemeQC schemeObservation type
      1Control experimentNo QCNone
      2(u, v) with independent errorsAdaptive QCASCAT wind
      3(spd, dir) with independent errorsAdaptive QCASCAT wind
      4(u, v) with correlated errorsAdaptive QCASCAT wind

      Table 1.  Different assimilation schemes for the ASCAT wind observation

      The control experiment is the natural run of the forecast without assimilation of any observations. Accordingly, we have three groups of assimilation experiments, including two wind component ($u,v$) schemes and one wind vector (${\rm{spd}},\;{\rm{dir}}$) scheme.

    5.   Results
    • Figure 5 gives the ensemble mean of the wind fields of forecasts and corresponding flow-dependent error of the $u$ and $v$ components. The background error of the wind field shows very similar characteristics to that of the observation in the central area of the typhoon. As shown in Figs. 4b, c; 5b, c, the error distribution is highly correlated to the values of spd and dir.

      Figure 5.  Ensemble mean of the wind fields of forecasts and corresponding flow-dependent error of the $u$ and $v$ components: (a) ensemble average of the wind fields of forecasts, (b) standard deviation of the $u$ component, and (c) standard deviation of the $v$ component.

      The results of QC by using different assimilation schemes are shown in Table 2 and Fig. 6. The reason for the observation gap in the center of the typhoon area for the wind-component schemes is mainly because of the position of the typhoon center in the background wind field (27.5°N, 157.5°E) being different from that of the observed wind field (28°N, 158°E). The misalignment of the typhoon vortex structure will result in a large deviation between the wind components of the background and observation. However, the use of the adaptive QC loosens the observation error in the high-wind region, which is equivalent to relaxing the QC threshold, and thus more observations can be reserved in the QC step (as shown in Fig. 6). For the wind vector scheme, the adaptive QC in the (${\rm{spd}},\;{\rm{dir}}$) form results in a much larger QC threshold than that of the wind-component schemes. Reasonable observation error also means better weight assignment for the cost function of assimilation, which is of great significance for the minimization process of assimilation. As listed in Table 2, the number of iterations of minimization for the new scheme is much smaller than that of the traditional schemes.

      No.Assimilation schemeObservation number before QCObservation number after QCIteration
      1(u, v) with independent errors1681159511
      2(spd, dir) with independent errors1681168111
      3(u, v) with correlated errors16811658 7

      Table 2.  QC and number of iterations of minimization in different assimilation schemes

      Figure 6.  QC using different assimilation schemes, in which the blue dots represent the wind observations that are successfully assimilated into the data assimilation system: (a) ($u$, $v$) with independent errors, (b) (spd, dir) with independent errors, and (c) ($u$, $v$) with correlated errors.

      For verification, we use the reanalysis data of ERA-Interim as the “truth.” ERA-Interim is the latest global atmospheric reanalysis, available from 1979, produced by ECMWF, and is continuously updated in real time (Dee et al., 2011). Figure 7 presents the analysis errors (analysis minus “truth”) of the spd field, and Fig. 8 gives the corresponding mean error and root-mean-square (RMS) value of the different experiments. Compared to the background, the wind-component assimilation schemes improve the spd field and the new scheme has the best fit to the reanalysis. As shown in Fig. 7, the assimilation of ASCAT wind components has improved the wind field structure in the typhoon area [the typhoon center is located at (28°N, 158°E) at the analysis time], which is critical to the track forecast. Although more observations are assimilated in the wind vector scheme, the improper observation error assigned for the high-wind observations may cause disharmony between the model states and the observations.

      Figure 7.  Wind speed (spd) analysis errors of different experiments: (a) background, (b) ($u$, $v$) components with independent error, (c) (spd, dir) with independent error, and (d) ($u$, $v$) components with correlated error.

      Figure 8.  Mean error and root-mean-square (RMS) of spd analysis errors of different experiments: (a) background, (b) ($u$, $v$) components with independent error, (c) (spd, dir) with independent error, and (d) ($u$, $v$) components with correlated error.

      The sea surface pressure errors of different experiments are also compared in Figs. 9, 10. Assimilation of ASCAT wind observations using wind-component schemes improves the sea surface pressure field compared to the background. The scheme using wind components with correlated error shows the best result for the pressure field (see Fig. 10d for details). However, they all have a bias compared to the reanalysis.

      Figure 9.  As in Fig. 7, but for pressure analysis errors.

      Figure 10.  As in Fig. 8, but for pressure analysis errors.

      Table 3 gives the analysis errors of the different assimilation schemes. It shows that in all cases the analysis errors of the spd field and pressure field improve when using the new assimilation scheme, which suggests that the new method is very promising for data assimilation.

      No.Assimilation scheme${ {e} }_{{\rm{speed}}}^{{\rm{mean}}}$${ {e} }_{{\rm{speed}}}^{{\rm{RMS}}}$${ {e} }_{{\rm{pressure}}}^{{\rm{mean}}}$${ {e} }_{{\rm{pressure}}}^{{\rm{RMS}}}$
      1Control experiment4.77−6.31 −7.11 7.80
      2(u, v) with independent errors3.05 4.13 −6.42 7.15
      3(spd, dir) with independent errors5.33 7.32−10.1711.12
      4(u, v) with correlated errors1.60 3.74 −5.87 6.66

      Table 3.  Mean error and RMS of analysis errors in different experiments

      It is also instructive to examine how the forecast skill changes with time. Figure 11a shows the best track and forecast typhoon paths of the different experiments, and Fig. 11b gives the errors of the forecast typhoon paths. The best-track data are from the China Meteorological Administration’s tropical cyclone database (tcdata.typhoon.org.cn; Ying et al., 2014). The different assimilation schemes each show some improvement in terms of the forecast path compared to that of the control experiment. The new scheme shows very promising results in the initial stage of the forecast compared to the other schemes. The track forecast of the new scheme begins to diverge at 18 h; however, it is still effective for the 24-h forecast.

      Figure 11.  The forecast (a) tracks and (b) track errors of Typhoon Noru (2017).

      The intensity of the typhoon based on the different schemes is also compared (Fig. 12), from which we can see that the new scheme shows strong improvement in intensity forecasts compared to the traditional scheme. To better demonstrate the results of the different experiments, Table 4 gives the average forecast errors of the typhoon track (${{\overline{e}_{{\rm{track}}}}}$), minimum pressure (${{\overline{e}_{{\rm{pressure}}}}}$), and maximum spd (${{\overline{e}_{{\rm{speed}}}}}$) of the typhoon eye. In each case the new scheme shows a promising result.

      Figure 12.  (a) Minimum pressure and (b) maximum spd forecast errors of the typhoon eye.

      No.Assimilation scheme${ {\overline{e}_{ {\rm{track} } } } }$${ {\overline{e}_{ {\rm{speed} } } } }$${ {\overline{e}_{ {\rm{pressure} } } } }$
      1Control experiment78.44−15.0022.66
      2(u, v) with independent errors40.96 −7.3211.52
      3(spd, dir) with independent errors82.90−14.8120.91
      4(u, v) with correlated errors41.85 −6.33 9.53

      Table 4.  Mean error of the typhoon forecast in different experiments

    6.   Summary
    • This paper presents a new method for assimilating remotely sensed sea surface wind components by considering their correlation. We have theoretically demonstrated how different it can be when considering the error characteristics of wind observations in their original form and the wind-components form. Given the error variances of wind observations in the form of spd and dir, the error variances and correlation of corresponding wind components can be calculated by the law of error propagation. Results from real data experiments compared to the traditional ($u,v$) assimilation scheme and the (${\rm{spd}},\;{\rm{dir}}$) scheme verify the ability of this new method.

      Since the main focus of this study is on the correlation of the wind components, other observation error sources, such as representativity error, have not been addressed, which may also affect the configuration of the observation error variance matrix. We also assume in this study that the observation errors of spd and dir are not correlated; however, large direction errors in wind observation are correlated to lighter spd for in situ observations (e.g., Dobson et al., 1980; Schwartz and Benjamin, 1995; Gao et al., 2012), as well as remotely sensed wind observations (e.g., Plant, 2000; Ebuchi et al., 2002). We will deal with these aspects in future studies.

Reference (27)

Catalog

    /

    DownLoad:  Full-Size Img  PowerPoint
    Return
    Return