Research Article  Open Access
Reconstruction of Typhoon Structure Using 3Dimensional Doppler Radar Radial Velocity Data with the Multigrid Analysis: A Case Study in an Idealized Simulation Context
Abstract
Extracting multiplescale observational information is critical for accurately reconstructing the structure of mesoscale circulation systems such as typhoon. The Space and Time Mesoscale Analysis System (STMAS) with multigrid data assimilation developed in Earth System Research Laboratory (ESRL) in National Oceanic and Atmospheric Administration (NOAA) has addressed this issue. Previous studies have shown the capability of STMAS to retrieve multiscale information in 2dimensional Doppler radar radial velocity observations. This study explores the application of 3dimensional (3D) Doppler radar radial velocities with STMAS for reconstructing a 3D typhoon structure. As for the first step, here, we use an idealized simulation framework. A twoscale simulated “typhoon” field is constructed and referred to as “truth,” from which randomly distributed conventional wind data and 3D Doppler radar radial wind data are generated. These data are used to reconstruct the synthetic 3D “typhoon” structure by the STMAS and the traditional 3D variational (3DVar) analysis. The degree by which the “truth” 3D typhoon structure is recovered is an assessment of the impact of the data type or analysis scheme being evaluated. We also examine the effects of weak constraint and strong constraint on STMAS analyses. Results show that while the STMAS is superior to the traditional 3DVar for reconstructing the 3D typhoon structure, the strong constraint STMAS can produce better analyses on both horizontal and vertical velocities.
1. Introduction
Doppler radar has long been a valuable observational tool in meteorology. Threedimensional (3D) Doppler radar radial velocity data can provide an opportunity to estimate both horizontal and vertical velocities. Therefore, in recent years, Doppler radar data assimilation for shortterm numerical weather forecasting or called nowcasting has become a focal point of research [1–4]. Lots of techniques have been developed to retrieve wind field from Doppler radar radial velocity observations [2–24].
In Doppler radar radial velocity data assimilation used in the above literatures, in a threedimensional variational (3DVar) framework, a background error covariance matrix is always needed to determine the spatial spreading of observational information. It is well known that an analysis field at different locations may have different correlation scales [25], which are difficult to be estimated. Unfortunately, the traditional 3DVar always employs an empirical and static background error covariance matrix and therefore usually can only correct singlescale wavelength error. However, the errors in short wavelength scales cannot be sufficiently corrected until the long waves are corrected [25, 26].
To minimize the errors of long and short waves in turn, a sequential 3DVar approach has been proposed by Xie et al. [25, 26], implemented by either a recursive filter [27] or a multigrid technique [28] at Global Systems Division (GSD) of Earth System Research Laboratory (ESRL) in National Oceanic and Atmospheric Administration (NOAA) for a Federal Aviation Agency (FAA) project joined by the research team from the Lincoln Laboratory in Massachusetts Institute of Technology (MIT). Since this system also uses the temporal observation information, it is called a Space and Time Mesoscale Analysis System (STMAS, thereafter; see Xie et al. [26]). The STMAS has been applied to assimilating 2dimensional (2D) Doppler radar radial velocity data to improve the wind field analyses [29].
Here, we study the analysis of 3D Doppler radar radial velocities using the STMAS to reconstruct the 3D wind structure. As for the first step, this study is performed in a twin experiment framework. In the next section, we first briefly review the theory of the multigrid 3DVar data assimilation scheme in the STMAS. Some important aspects of the STMAS techniques such as smoothing, constraint, and Doppler radar radial wind operators used in the cost function of the STMAS multigrid 3DVar are described. Section 3 first introduces the twin experiment framework for 3D Doppler radar radial wind data assimilation with the STMAS and then gives the evaluation by comparing it to the traditional 3DVar. Section 3 also examines the performance of the STMAS in weak and strong constraints for 3D Doppler radar radial velocity analysis. Conclusions and discussions are given in Section 4.
2. Smoothing, Constraint, and Radar Radial Wind Operators in STMAS
In this study, the STMAS implemented by the multigrid 3DVar is applied to the analysis of 3D Doppler radar radial velocities. This method can extract long and short wavelength information in turn efficiently from observations and provide objective and accurate analysis. The basic idea of this multigrid implementation can be referred to Li et al. [28–30].
To assimilate 3D Doppler radar radial velocities, with the control variables being , where and represent zonal and meridional components of wind vector, the cost functional for the th level grid iswhere the subscript denotes the background term and the smooth term, the conventional observation data term, and the radar radial wind observation data term. The smooth matrixes and in the smooth term are derived from the Laplacian of control variables and , respectively, at grid points. Let represent the vertical component of wind vector. The details of the conventional observation data term and the 3D Doppler radar radial wind observation data term are as follows:where is the amount of radar radial wind observations, is the azimuth angle of the radar beam relative to north with positive clockwise, and is elevation angle of the radar beam. Of course, since radar scans at nonzero elevation angles, the fall speed of precipitation particles should be taken into account, and the algorithm of Sun and Crook [6] can be used to calculate terminal velocity. But for this study, we just neglect this terminal velocity, which does not lose its generality. The matrix is an error covariance matrix for radar radial wind observation; its superscript stands for the reverse matrix, and its subscript represents the radar radial wind observation.
During the procedure of sequential multiscale analyses, the operators , , , and remain the same when the full observation dataset is used through all multigrid levels; therefore, the superscript is omitted from these operators.
To make a strong constraint on these three components of wind vector, incompressible continuity equation is employed and discretized to calculate vertical velocity from the other two horizontal components. The discretized incompressible continuity equation is as follows:The adjoint codes are recursively developed for represented by and .
3. Simulation Methodology
3.1. Synthetic Typhoon Structure
The study domain covers a square region with 10 km thickness. The Doppler radar locates at the center (250 km, 250 km) of the study domain. A simulated typhoon field can be constructed by using the following function which consists of two subsection functions:This formula allows that the 2order derivatives of this function exist. Let and , where km^{2} s^{−1}, km, km, km^{2} s^{−1}, km, and km; then a stream function of a twoscale typhoon field can be constructed. The amplitude of one largescale ( km) is scaled by and the amplitude of one smallscale ( km) is scaled by . The locations of the maximum wind horizontal velocity of these two scales are km for largescale and km for smallscale, respectively, from the typhoon center.
The typhoon center is set at 300 km and 150 km. Then, the stream function and velocity potential function can be constructed as follows:where , km, and km. The horizontal components can be expressed in terms of and :In this study, incompressibility is assumed. And the true vertical velocity field can be obtained by integrating the continuity equation . That is,where the bottom boundary condition is . The wind speed field in this simulated typhoon field contains two different scale information. The first one is about 35 km and the other is about 10 km. The radial wind, component, and component pattern of middle level of this simulated typhoon field and a section wind field across the center of this typhoon are shown in Figure 1. This typhoon pattern is located at the southeast part of the study domain, so only the southeast square part is shown for the detailed structure. This simulated typhoon wind field is referred to as the “truth” typhoon field in this twin experiment.
(a)
(b)
(c)
(d)
(e)
(f)
3.2. Observations
Then, Doppler radar radial velocity data are generated from the “truth” typhoon field with onedegree azimuth angle increment and 2500 m gate spacing and 2degree elevation angle increment from 1 degree to 20 degrees by interpolating the “truth” velocity field to the radial velocity observations’ points and using the equation , and the number of radial velocity data is 114480. Spatially coarse (1000), moderate (10000), and dense (100000) random distributed conventional observations are also generated, respectively, from the “truth” typhoon field by interpolating the “truth” velocity field to the conventional observations’ points.
3.3. Twin Experiments Setup
In the following, the above observational data are used to retrieve the simulated “typhoon” structure by the STMAS analysis method with weak or strong constraint and the traditional 3dimensional variance (3DVar) analysis with different correlation scales, respectively, and by comparing the analyzed results with the “truth,” performances of different analysis methods are discussed.
The error variances of radial velocity observations and conventional observations can be determined by the measurement error of instruments. But, here for simplicity, the same error variance is set for each kind of data. However, because the amount of radial velocity data is much larger than that of conventional data, a scaling scheme is used to balance the weights of these two types of observations. Thus, the conventional observation can have the same weight as that of radial velocity observation, which may comprise these two types of observations to get to a reasonable wind analysis.
The limited memory BFGS (BroydenFletcherGoldfarbShanno) method [31] to solve the bound constrained optimization problem [32] is used as the minimization method in this study.
3.4. STMAS Results
Three level grids are employed ranging from about 31.25 km × 31.25 km × 2.5 km (, and the number of grids is 17 × 17 × 5) to 7.8125 km × 7.8125 km × 0.625 km (, and the number of grids is 65 × 65 × 17) with grid ratio being 0.5. The background is set to be zero for simplicity. The conventional data will be added gradually to investigate the impact of radar radial wind observations.
The STMAS analyses with the above 3D Doppler radar radial velocity data or conventional data which vary from coarse to dense are shown in Figures 2–4 and Table 1 gives the corresponding root mean square errors (RMSEs). Because the Doppler radar radial velocity data really provide some useful information of radial wind, the STMAS results by only using this type of data can make a good analysis on radial wind (see Figure 2(a)) and the twoscale information in radial wind has been captured by the STMAS. There is no information of tangential wind, and no statistical or empirical correlation information between and is used in this study. Therefore, only using radial velocity data, the STMAS analysis only can match the radial wind, and it cannot make a good analysis on and . This leads to large and RMSEs (15.7 m/s and 14.5 m/s, resp., see Table 1). With the 1000 coarse distributed conventional data only, the STMAS analysis only shows the large pattern of the typhoon wind (Figure 2(b)). With the combination of conventional data and radar radial wind data, the STMAS analysis significantly improves the typhoon structure (Figure 2(c)). The RMSE with the combined datasets is much smaller than the other one with an individual dataset (Table 1). Apparently, the radial velocity data provide additional information on radial direction to help the STMAS capture the detailed typhoon structure. If the conventional data change from the coarsest (totally 1000) to the densest (totally 100000), the detailed information on the typhoon is enhanced gradually (see Figures 2(b) and 2(c) corresponding to 1,000 conventional data and 1,000 conventional data combined with radar data, resp.; Figures 3(a) and 3(b) corresponding to 10,000 conventional data and 10,000 conventional data combined with radar data, resp.; Figures 4(a) and 4(b) corresponding to 100,000 conventional data and 100,000 conventional data combined with radar data, resp.). With the combination of radar radial velocity data and dense distributed conventional data, the STMAS can make almost perfect analysis for the typhoon structure (Figure 4(b)) and the RMSE is very small. And because dense distributed conventional data provide enough information, the RMSE improvements from the MG_DENSE to MG_RADAR_DENSE experiments are really negligible for the three wind components.

(a)
(b)
(c)
(a)
(b)
(a)
(b)
(c)
3.5. Strong Constraint versus Weak Constraint
To compare the performance of the STMAS with a strong constraint or weak constraint, the continuity equation is added as penalty term to make a weak constraint case (otherwise, the STMAS analysis is a strong constraint of the three components of wind vector). Then, the control variables become :And background term is changed to The smooth term is Here, is smooth matrix for component of wind vector. The weak constraint term is where is the penalty coefficient for weak constraint. The other terms, and , are the same as (2) and (3), respectively.
Assimilating the whole simulated Doppler radar radial velocity data and dense conventional data, the STMAS analyses with continuity equation as strong constraint or weak constraint are shown in Figures 4(b) and 4(c), respectively. And the RMSE is shown in Table 1. The coefficient of penalty term for weak constraint is set to 100.
As shown in Figure 4(b), the strong constraint forces these three components of wind vector to exactly satisfy the continuity equation. Then, while the horizontal components can be analyzed very well, the vertical component is reconstructed very well too. The weak constraint (Figure 4(c)) can only constrain the horizontal and vertical components of wind vector to satisfy the continuity equation to some degree. Whereas the order of vertical velocities is much smaller than that of horizontal ones and no direct observation of vertical velocities is available, the minimization of cost function primarily focuses on the horizontal parts. This causes little improvement on the vertical velocities even if the horizontal components can be well analyzed. Therefore, the vertical component RMSE score of MG_RADAR_DENSE_WEAK experiment is worse than that of MG_RADAR_DENSE experiment. When a large penalty coefficient is used, although vertical velocity analysis can be improved a little bit suffering a lower accuracy of horizontal velocity, the accuracy of the analysis vertical velocities is still not high.
3.6. STMAS versus Traditional 3DVar
For a traditional 3DVar analysis, the cost function takes the formFor simplicity, here, we assume , and no covariance between the two horizontal velocity components is considered. The background error covariance matrix takes the following form [33]:where , , and are characteristic length scales which reflect the extent of spatial correlation; and are coordinates, and and are indexes of grid point; and is the firstguess error variance. and are similar to formulas (2) and (3), respectively, but with an omitted superscript . Strong constraint is also imposed on , , and component by using discretized continuity equation. To make a comparison with STMAS, the dense conventional data and Doppler radar radial velocity data are all used. The RMSE is shown in Table 1.
Based on Gaussian distribution, the traditional 3DVar using correlation scale usually constructs the background error covariance matrix by an empirical correlation scale. Therefore, the traditional 3DVar with a certain correlation scale only can analyze this kind of scale information. However, the “truth” typhoon field in this study contains two different spatial scales wind speed information (~35 km and ~10 km). The traditional 3DVar with 50 km horizontal correlation scale can only capture the main pattern of this typhoon field (i.e., the long wave information) but lose the smallscale information and produce a smooth analysis (Figure 5(a)). The traditional 3DVar with 25 km or 12.5 km horizontal correlation scale can capture some detailed information of the long wave as well as major short wave information, but detailed short wave information cannot be well analyzed (Figures 5(b) and 5(c)). On the contrary, the traditional 3DVar with 5 km correlation scale can analyze short wave features but incorrectly treats long wave information. This leads to an erroneous analysis (Figure 5(d)). The two scales of this true typhoon field are 35 km and 10 km, respectively, and the largescale (35 km) component covers most of study domain. In the experiment T25_RADAR_DENSE, the horizontal correlation scale is 25 km which is the closest to the largescale (35 km) of the true typhoon field among all these traditional experiments; therefore, the horizontal wind velocity analysis is the best.
(a)
(b)
(c)
(d)
From the vertical velocity distribution as well as wind vector shown in Figure 6, we can further verify the above argument. Near the typhoon center, there are two branches of vertical circulations on each side: one is narrow and the other is broad, corresponding to four upwelling zones in this section. Using the dense conventional data and Doppler radar radial velocity data, the STMAS can analyze these sets of vertical circulations and corresponding upwelling zones (Figure 6(a)). However, the traditional 3DVar with horizontal correlation scales being 50 km, 25 km, or 12.5 km can only analyze two (Figure 6(b)) or three (Figures 6(c) and 6(d)) upwelling zones, since these kinds of correlation scales filter out the smallscale information and merge these upwelling zones. Although the traditional 3DVar with 6.25 km horizontal correlation scale can distinguish the middle two upwelling zones, it incorrectly treats the largescale information (Figure 6(e)).
(a)
(b)
(c)
(d)
(e)
4. Conclusions and Discussions
Within an idealized simulation framework, the role of 3D Doppler radar radial velocity data for reconstructing 3D typhoon structures has been examined using the Space and Time Mesoscale Analysis System (STMAS). A twoscale simulated “typhoon” field is constructed and referred to as “truth,” from which randomly distributed conventional wind data and 3D Doppler radar radial wind data are generated. These data are used to reconstruct the synthetic 3D “typhoon” structure by the STMAS or the traditional 3D variational (3DVar) analysis. The degree by which the “truth” 3D typhoon structure is recovered is an assessment of the impact of data type or analysis scheme being evaluated. The effects of weak or strong constraint on STMAS analysis have also been examined. We found that () the STMAS is superior to traditional 3DVar for reconstructing the 3D typhoon structure, since the STMAS can retrieve multiscale information from observational network. () The radial velocity data provide additional useful information for the STMAS to reconstruct the detailed structure of 3D typhoon field. () Compared to a weak constraint, the strong constraint STMAS can produce better analyses on both horizontal and vertical velocities of the 3D typhoon structure.
This study gives us promising results. Challenges still remain when 3D radar radial velocity data are assimilated for the reconstruction and initialization of real typhoon structures in the future. First, given the fact that the real atmosphere is compressive, the wind vector may not satisfy the nondiffusivity continuity equation used in this study. Therefore, a full continuity equation should be used to make a more general strong constraint in future study. Second, the model error has not been taken into account in this study. The influence of model errors on typhoon reconstruction and initialization has to be addressed and how to deal with model errors could be an important research topic in the followup studies.
Competing Interests
The authors declare that there are no competing interests regarding the publication of this paper.
Acknowledgments
This research was jointly supported by grants of National Basic Research Program (2013CB430304), National Natural Science Foundation (41376013, 41376015, 41306006, 41541041, and 41506039), National HighTech R&D Program (2013AA09A505), and National Programme on Global Change and AirSea Interaction (GASI010112 and GASIIPOVAI04) of China.
References
 J. Sun and N. A. Crook, “Realtime lowlevel wind and temperature analysis using single WSR88D data,” Weather and Forecasting, vol. 16, no. 1, pp. 117–132, 2001. View at: Publisher Site  Google Scholar
 S. S. Weygandt, A. Shapiro, and K. K. Droegemeier, “Retrieval of model initial fields from singleDoppler observations of a supercell thunderstorm. Part I: singleDoppler velocity retrieval,” Monthly Weather Review, vol. 130, no. 3, pp. 433–453, 2002. View at: Publisher Site  Google Scholar
 S. S. Weygandt, A. Shapiro, and K. K. Droegemeier, “Retrieval of model initial fields from singleDoppler observations of a supercell thunderstorm. Part II: thermodynamic retrieval and numerical prediction,” Monthly Weather Review, vol. 130, no. 3, pp. 454–476, 2002. View at: Publisher Site  Google Scholar
 Y.C. Liou, J.L. Chiou, W.H. Chen, and H.Y. Yu, “Improving the model convective storm quantitative precipitation nowcasting by assimilating state variables retrieved from multipleDoppler radar observations,” Monthly Weather Review, vol. 142, no. 11, pp. 4017–4035, 2014. View at: Publisher Site  Google Scholar
 J. Sun and N. A. Crook, “Dynamical and microphysical retrieval from Doppler radar observations using a cloud model and its adjoint—part I: model development and simulated data experiments,” Journal of the Atmospheric Sciences, vol. 54, no. 12, pp. 1642–1661, 1997. View at: Publisher Site  Google Scholar
 J. Sun and N. A. Crook, “Dynamical and microphysical retrieval from Doppler radar observations using a cloud model and its adjoint. Part II: retrieval experiments of an observed Florida convective storm,” Journal of the Atmospheric Sciences, vol. 55, no. 5, pp. 835–852, 1998. View at: Publisher Site  Google Scholar
 J. Gao, M. Xue, A. Shapiro, and K. K. Droegemeier, “A variational method for the analysis of threedimensional wind fields from two Doppler radars,” Monthly Weather Review, vol. 127, no. 9, pp. 2128–2142, 1999. View at: Publisher Site  Google Scholar
 M. Lindskog, H. Järvinen, and D. B. Michelson, “Development of radar wind data assimilation for the HIRLAM 3DVar,” HIRLAM Technical Report 52, 2002. View at: Google Scholar
 J.D. Gao, M. Xue, K. Brewster, and K. K. Droegemeier, “A threedimensional variational data analysis method with recursive filter for Doppler radars,” Journal of Atmospheric and Oceanic Technology, vol. 21, no. 3, pp. 457–469, 2004. View at: Publisher Site  Google Scholar
 Q. Xiao, Y.H. Kuo, J. Sun et al., “Assimilation of Doppler radar observations with a regional 3DVAR system: impact of Doppler velocities on forecasts of a heavy rainfall case,” Journal of Applied Meteorology, vol. 44, no. 6, pp. 768–788, 2005. View at: Publisher Site  Google Scholar
 J. Gao, M. Xue, S.Y. Lee, A. Shapiro, Q. Xu, and K. K. Droegemeier, “A threedimensional variational singleDoppler velocity retrieval method with simple conservation equation constraint,” Meteorology and Atmospheric Physics, vol. 94, no. 1–4, pp. 11–26, 2006. View at: Publisher Site  Google Scholar
 M. Hu, M. Xue, and K. Brewster, “3DVAR and cloud analysis with WSR88D levelII data for the prediction of the Fort Worth, Texas, tornadic thunderstorms. Part I: cloud analysis and its impact,” Monthly Weather Review, vol. 134, no. 2, pp. 675–698, 2006. View at: Publisher Site  Google Scholar
 M. Hu, M. Xue, J. Gao, and K. Brewster, “3DVAR and cloud analysis with WSR88D levelII data for the prediction of the Fort Worth, Texas, tornadic thunderstorms—part II: impact of radial velocity analysis via 3DVAR,” Monthly Weather Review, vol. 134, no. 2, pp. 699–721, 2006. View at: Publisher Site  Google Scholar
 Q. Xiao and J. Sun, “Multipleradar data assimilation and shortrange quantitative precipitation forecasting of a squall line observed during IHOP_2002,” Monthly Weather Review, vol. 135, no. 10, pp. 3381–3404, 2007. View at: Publisher Site  Google Scholar
 S. Sugimoto, N. A. Crook, J. Sun, Q. Xiao, and D. Barker, “Assimilation of Doppler radar data with WRF 3DVAR: evaluation of its potential benefits to quantitative precipitation forecasting through Observing System Simulation Experiments,” Monthly Weather Review, vol. 137, pp. 4011–4029, 2009. View at: Google Scholar
 F. Q. Zhang, Y. Weng, J. A. Sippel, Z. Meng, and C. H. Bishop, “Cloudresolving hurricane initialization and prediction through assimilation of doppler radar observations with an ensemble Kalman filter,” Monthly Weather Review, vol. 137, no. 7, pp. 2105–2125, 2009. View at: Publisher Site  Google Scholar
 K. Zhao and M. Xue, “Assimilation of coastal Doppler radar data with the ARPS 3DVAR and cloud analysis for the prediction of Hurricane Ike (2008),” Geophysical Research Letters, vol. 36, no. 12, Article ID L12803, 2009. View at: Publisher Site  Google Scholar
 A. D. Schenkman, M. Xue, A. Shapiro, K. Brewster, and J. Gao, “The analysis and prediction of the 89 May 2007 Oklahoma tornadic mesoscale convective system by assimilating WSR88D and CASA radar data using 3DVAR,” Monthly Weather Review, vol. 139, no. 1, pp. 224–246, 2011. View at: Publisher Site  Google Scholar
 Y. Z. Li, X. G. Wang, and M. Xue, “Assimilation of radar radial velocity data with the wrf hybrid ensemble–3DVAR system for the prediction of hurricane ike (2008),” Monthly Weather Review, vol. 140, no. 11, pp. 3507–3524, 2012. View at: Publisher Site  Google Scholar
 N. Du, M. Xue, K. Zhao, and J. Min, “Impact of assimilating airborne Doppler radar velocity data using the ARPS 3DVAR on the analysis and prediction of Hurricane Ike (2008),” Journal of Geophysical Research Atmospheres, vol. 117, no. 17, Article ID D18113, 2012. View at: Publisher Site  Google Scholar
 S. Abhilash, A. K. Sahai, K. Mohankumar, J. P. George, and S. Das, “Assimilation of Doppler weather radar radial velocity and reflectivity observations in WRF3DVAR system for shortrange forecasting of convective storms,” Pure and Applied Geophysics, vol. 169, no. 11, pp. 2047–2070, 2012. View at: Publisher Site  Google Scholar
 J. Dong and M. Xue, “Assimilation of radial velocity and reflectivity data from coastal WSR88D radars using an ensemble Kalman filter for the analysis and forecast of landfalling hurricane Ike (2008),” Quarterly Journal of the Royal Meteorological Society, vol. 139, no. 671, pp. 467–487, 2013. View at: Publisher Site  Google Scholar
 L. Zhu, Q. Wan, X. Shen et al., “Prediction and predictability of highimpact Western Pacific landfalling Tropical Cyclone Vicente (2012) through convectionpermitting ensemble assimilation of doppler radar velocity,” Monthly Weather Review, vol. 144, no. 1, pp. 21–43, 2016. View at: Publisher Site  Google Scholar
 F. Shen, J. Min, and D. Xu, “Assimilation of radar radial velocity data with the WRF Hybrid ETKF3DVAR system for the prediction of Hurricane Ike (2008),” Atmospheric Research, vol. 169, pp. 127–138, 2016. View at: Publisher Site  Google Scholar
 Y. F. Xie, S. E. Koch, J. A. McGinley, S. Albers, and N. Wang, “A sequential variational analysis approach for mesoscale data assimilation,” in Proceedings of the 21st Conference on Weather Analysis and Forecasting/17th Conference on Numerical Weather Prediction, 15B.7, American Meteorological Society, Washington, DC, USA, 2005, http://ams.confex.com/ams/pdfpapers/93468.pdf. View at: Google Scholar
 Y. Xie, S. Koch, J. McGinley et al., “A spacetime multiscale analysis system: a sequential variational analysis approach,” Monthly Weather Review, vol. 139, no. 4, pp. 1224–1240, 2011. View at: Publisher Site  Google Scholar
 Z. J. He, Y. F. Xie, W. Li et al., “Application of the sequential threedimensional variational method to assimilating SST in a global ocean model,” Journal of Atmospheric and Oceanic Technology, vol. 25, no. 6, pp. 1018–1033, 2008. View at: Publisher Site  Google Scholar
 W. Li, Y. Xie, Z. He et al., “Application of the multigrid data assimilation scheme to the China seas' temperature forecast,” Journal of Atmospheric and Oceanic Technology, vol. 25, no. 11, pp. 2106–2116, 2008. View at: Publisher Site  Google Scholar
 W. Li, Y. F. Xie, S.M. Deng, and Q. Wang, “Application of the multigrid method to the twodimensional Doppler radar radial velocity data assimilation,” Journal of Atmospheric and Oceanic Technology, vol. 27, no. 2, pp. 319–332, 2010. View at: Publisher Site  Google Scholar
 W. Li, Y. F. Xie, and G. J. Han, “A theoretical study on the multigrid data assimilation scheme using a simple bilinear interpolation scheme,” Acta Oceanologica Sinica, vol. 32, no. 3, pp. 80–87, 2013. View at: Google Scholar
 R. Fletcher, Practical Methods of Optimization, John Wiley & Sons, Chichester, UK, 2nd edition, 1987. View at: MathSciNet
 R. H. Byrd, P. Lu, J. Nocedal, and C. Y. Zhu, “A limited memory algorithm for bound constrained optimization,” SIAM Journal on Scientific Computing, vol. 16, no. 5, pp. 1190–1208, 1995. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 J. Derber and A. Rosati, “A global oceanic data assimilation system,” Journal of Physical Oceanography, vol. 19, no. 9, pp. 1333–1347, 1989. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2016 Hongli Fu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.