Tropical cyclone(TC) rainfall forecast has remained a challenge. To create initial conditions with high quality for simulation, the present study implemented a data assimilation scheme based on the EnKF method to inge...Tropical cyclone(TC) rainfall forecast has remained a challenge. To create initial conditions with high quality for simulation, the present study implemented a data assimilation scheme based on the EnKF method to ingest the satellite-retrieved cloud water path(C_(w)) and tested it in WRF. The scheme uses the vertical integration of forecasted cloud water content to transform control variables to the observation space, and creates the correlations between C_(w) and control variables in the flow-dependent background error covariance based on all the ensemble members, so that the observed cloud information can affect the background temperature and humidity. For two typhoons in 2018(Yagi and Rumiba), assimilating C_(w) significantly increases the simulated rainfalls and TC intensities. In terms of the average equitable threat score of daily moderate to heavy rainfall(5-120 mm), the improvements are over 130%, and the dry biases are cut by about 30%. Such improvements are traced down to the fact that C_(w) assimilation increases the moisture content, especially that further away from the TC center, which provides more precipitable water for the rainfall,strengthens the TC and broadens the TC size via latent heat release and internal wind field adjustment.展开更多
A 3D compressible nonhydrostatic dynamic core based on a three-point multi-moment constrained finite-volume (MCV) method is developed by extending the previous 2D nonhydrostatic atmospheric dynamics to 3D on a terrain...A 3D compressible nonhydrostatic dynamic core based on a three-point multi-moment constrained finite-volume (MCV) method is developed by extending the previous 2D nonhydrostatic atmospheric dynamics to 3D on a terrainfollowing grid. The MCV algorithm defines two types of moments: the point-wise value (PV) and the volume-integrated average (VIA). The unknowns (PV values) are defined at the solution points within each cell and are updated through the time evolution formulations derived from the governing equations. Rigorous numerical conservation is ensured by a constraint on the VIA moment through the flux form formulation. The 3D atmospheric dynamic core reported in this paper is based on a three-point MCV method and has some advantages in comparison with other existing methods, such as uniform third-order accuracy, a compact stencil, and algorithmic simplicity. To check the performance of the 3D nonhydrostatic dynamic core, various benchmark test cases are performed. All the numerical results show that the present dynamic core is very competitive when compared to other existing advanced models, and thus lays the foundation for further developing global atmospheric models in the near future.展开更多
A positivity-preserving conservative semi-Lagrangian transport model by multi-moment finite volume method has been developed on the cubed-sphere grid.Two kinds of moments(i.e.,point values(PV moment) at cell interface...A positivity-preserving conservative semi-Lagrangian transport model by multi-moment finite volume method has been developed on the cubed-sphere grid.Two kinds of moments(i.e.,point values(PV moment) at cell interfaces and volume integrated average(VIA moment) value) are defined within a single cell.The PV moment is updated by a conventional semi-Lagrangian method,while the VIA moment is cast by the flux form formulation to assure the exact numerical conservation.Different from the spatial approximation used in the CSL2(conservative semi-Lagrangian scheme with second order polynomial function) scheme,a monotonic rational function which can effectively remove non-physical oscillations is reconstructed within a single cell by the PV moments and VIA moment.To achieve exactly positive-definite preserving,two kinds of corrections are made on the original conservative semi-Lagrangian with rational function(CSLR)scheme.The resulting scheme is inherently conservative,non-negative,and allows a Courant number larger than one.Moreover,the spatial reconstruction can be performed within a single cell,which is very efficient and economical for practical implementation.In addition,a dimension-splitting approach coupled with multi-moment finite volume scheme is adopted on cubed-sphere geometry,which benefitsthe implementation of the 1 D CSLR solver with large Courant number.The proposed model is evaluated by several widely used benchmark tests on cubed-sphere geometry.Numerical results show that the proposed transport model can effectively remove nonphysical oscillations and preserve the numerical nonnegativity,and it has the potential to transport the tracers accurately in a real atmospheric model.展开更多
基金National Key R&D Project of China(2018YFC1507001)。
文摘Tropical cyclone(TC) rainfall forecast has remained a challenge. To create initial conditions with high quality for simulation, the present study implemented a data assimilation scheme based on the EnKF method to ingest the satellite-retrieved cloud water path(C_(w)) and tested it in WRF. The scheme uses the vertical integration of forecasted cloud water content to transform control variables to the observation space, and creates the correlations between C_(w) and control variables in the flow-dependent background error covariance based on all the ensemble members, so that the observed cloud information can affect the background temperature and humidity. For two typhoons in 2018(Yagi and Rumiba), assimilating C_(w) significantly increases the simulated rainfalls and TC intensities. In terms of the average equitable threat score of daily moderate to heavy rainfall(5-120 mm), the improvements are over 130%, and the dry biases are cut by about 30%. Such improvements are traced down to the fact that C_(w) assimilation increases the moisture content, especially that further away from the TC center, which provides more precipitable water for the rainfall,strengthens the TC and broadens the TC size via latent heat release and internal wind field adjustment.
基金supported by the National Key Research and Development Program of China (Grant Nos. 2017YFC1501901 and 2017YFA0603901)the Beijing Natural Science Foundation (Grant No. JQ18001)
文摘A 3D compressible nonhydrostatic dynamic core based on a three-point multi-moment constrained finite-volume (MCV) method is developed by extending the previous 2D nonhydrostatic atmospheric dynamics to 3D on a terrainfollowing grid. The MCV algorithm defines two types of moments: the point-wise value (PV) and the volume-integrated average (VIA). The unknowns (PV values) are defined at the solution points within each cell and are updated through the time evolution formulations derived from the governing equations. Rigorous numerical conservation is ensured by a constraint on the VIA moment through the flux form formulation. The 3D atmospheric dynamic core reported in this paper is based on a three-point MCV method and has some advantages in comparison with other existing methods, such as uniform third-order accuracy, a compact stencil, and algorithmic simplicity. To check the performance of the 3D nonhydrostatic dynamic core, various benchmark test cases are performed. All the numerical results show that the present dynamic core is very competitive when compared to other existing advanced models, and thus lays the foundation for further developing global atmospheric models in the near future.
基金supported by the National Key Research and Development Program of China (Grant Nos.2017YFC1501901 and 2017YFA0603901)the Beijing Natural Science Foundation (Grant No.JQ18001)。
文摘A positivity-preserving conservative semi-Lagrangian transport model by multi-moment finite volume method has been developed on the cubed-sphere grid.Two kinds of moments(i.e.,point values(PV moment) at cell interfaces and volume integrated average(VIA moment) value) are defined within a single cell.The PV moment is updated by a conventional semi-Lagrangian method,while the VIA moment is cast by the flux form formulation to assure the exact numerical conservation.Different from the spatial approximation used in the CSL2(conservative semi-Lagrangian scheme with second order polynomial function) scheme,a monotonic rational function which can effectively remove non-physical oscillations is reconstructed within a single cell by the PV moments and VIA moment.To achieve exactly positive-definite preserving,two kinds of corrections are made on the original conservative semi-Lagrangian with rational function(CSLR)scheme.The resulting scheme is inherently conservative,non-negative,and allows a Courant number larger than one.Moreover,the spatial reconstruction can be performed within a single cell,which is very efficient and economical for practical implementation.In addition,a dimension-splitting approach coupled with multi-moment finite volume scheme is adopted on cubed-sphere geometry,which benefitsthe implementation of the 1 D CSLR solver with large Courant number.The proposed model is evaluated by several widely used benchmark tests on cubed-sphere geometry.Numerical results show that the proposed transport model can effectively remove nonphysical oscillations and preserve the numerical nonnegativity,and it has the potential to transport the tracers accurately in a real atmospheric model.
