With the aid of a global barotropic model, the role of the interaction of the synoptic-scale disturbance and the planetary flow in block onset is examined by a 4-dimensional variational approach. A cost function is de...With the aid of a global barotropic model, the role of the interaction of the synoptic-scale disturbance and the planetary flow in block onset is examined by a 4-dimensional variational approach. A cost function is defined to measure the squared errors of the forecasted stream functions during block onset period (day 4 and day 5 in this study) over a selected blocking domain. The sensitivity of block onset with respect to the initial synoptic-scale disturbance is studied by examining the gradient of the defined cost function with respect to the initial (during the first 24 hours) vorticity forcing, which is evaluated by the adjoint integration. Furthermore, the calculated cost function and gradient are connected with the limited-memory quasi-Newton optimization algorithm for solving the optimal initial vorticity forcing for block onset. For two studied cases of block onset (northern Atlantic and northern Pacific) introducing the optimal initial vorticity forcing, the nonlinear barotropic advection process mostly reconstructs these blocking onset processes. The results show that the formation of blocking can be correctly described by a barotropic nonlinear advection process, in which the wave- (synoptic-scale) flow (planetary-scale) interaction plays a very important role. On an appropriate planetary-scale flow, a certain synoptic-scale disturbance can cause the blocking onset by the interaction between the synoptic scale perturbations and the planetary scale basic flows. The extended forecasts show that the introduction of the optimal initial vorticity forcing can predict the blocking process up to the 7th or 8th day in this simple model case. The experimental results in this study show that the 4-dimensional variational approach has a good potential to be applied to study the dynamics of the medium-range weather processes. This simple model case study is only an initial trial. Applying the framework in this study to a complex model will further our understanding of the mechanism of the atmospheric/oceanic processes and improve their prediction.展开更多
The local convection initiation(CI)mechanisms of a convective case that occurred on5 August 2017 in Cangzhou,northern China,were studied using Doppler radar and automatic weather station observational analysis,along w...The local convection initiation(CI)mechanisms of a convective case that occurred on5 August 2017 in Cangzhou,northern China,were studied using Doppler radar and automatic weather station observational analysis,along with Variational Doppler Radar Analysis System assimilation analysis.During the convective process,a gust front appeared ahead of two existing convective systems,respectively.In the warm and moist environment ahead of the gust fronts in the south,there was a mesoscale air mass boundary.With the process of a gust front moving southward,approaching the mesoscale air mass boundary,the convergence intensified in the area between the gust front and the mesoscale air mass boundary.Finally,the strong convergent updraft exceeded the level of free convection and triggered the new convection.展开更多
The impact of diabatic processes on 4-dimensional variational data assimilation (4D-Var) was studied using the 1995 version of NCEP's global spectral model with and without full physics.The adjoint was coded manua...The impact of diabatic processes on 4-dimensional variational data assimilation (4D-Var) was studied using the 1995 version of NCEP's global spectral model with and without full physics.The adjoint was coded manually.A cost function measuring spectral errors of 6-hour forecasts to 'observation' (the NCEP reanalysis data) was minimized using the L-BFGS (the limited memory quasi-Newton algorithm developed by Broyden,Fletcher,Goldfard and Shanno) for optimizing parameters and initial conditions.Minimization of the cost function constrained by an adiabatic version of the NCEP global model converged to a minimum with a significant amount of decrease in the value of the cost function.Minimization of the cost function using the diabatic model, however,failed after a few iterations due to discontinuities introduced by physical parameterizations.Examination of the convergence of the cost function in different spectral domains reveals that the large-scale flow is adjusted during the first 10 iterations,in which discontinuous diabatic parameterizations play very little role.The adjustment produced by the minimization gradually moves to relatively smaller scales between 10-20th iterations.During this transition period,discontinuities in the cost function produced by 'on-off' switches in the physical parameterizations caused the cost function to stay in a shallow local minimum instead of continuously decreasing toward a deeper minimum. Next,a mixed 4D-Var scheme is tested in which large-scale flows are first adiabatically adjusted to a sufficient level,followed by a diabatic adjustment introduced after 10 to 20 iterations. The mixed 4D-Var produced a closer fit of analysis to observations,with 38% and 41% more decrease in the values of the cost function and the norm of gradient,respectively,than the standard diabatic 4D-Var,while the CPU time is reduced by 21%.The resulting optimal initial conditions improve the short-range forecast skills of 48-hour statistics.The detrimental effect of parameterization discontinuities on minimization was also reduced.展开更多
文摘With the aid of a global barotropic model, the role of the interaction of the synoptic-scale disturbance and the planetary flow in block onset is examined by a 4-dimensional variational approach. A cost function is defined to measure the squared errors of the forecasted stream functions during block onset period (day 4 and day 5 in this study) over a selected blocking domain. The sensitivity of block onset with respect to the initial synoptic-scale disturbance is studied by examining the gradient of the defined cost function with respect to the initial (during the first 24 hours) vorticity forcing, which is evaluated by the adjoint integration. Furthermore, the calculated cost function and gradient are connected with the limited-memory quasi-Newton optimization algorithm for solving the optimal initial vorticity forcing for block onset. For two studied cases of block onset (northern Atlantic and northern Pacific) introducing the optimal initial vorticity forcing, the nonlinear barotropic advection process mostly reconstructs these blocking onset processes. The results show that the formation of blocking can be correctly described by a barotropic nonlinear advection process, in which the wave- (synoptic-scale) flow (planetary-scale) interaction plays a very important role. On an appropriate planetary-scale flow, a certain synoptic-scale disturbance can cause the blocking onset by the interaction between the synoptic scale perturbations and the planetary scale basic flows. The extended forecasts show that the introduction of the optimal initial vorticity forcing can predict the blocking process up to the 7th or 8th day in this simple model case. The experimental results in this study show that the 4-dimensional variational approach has a good potential to be applied to study the dynamics of the medium-range weather processes. This simple model case study is only an initial trial. Applying the framework in this study to a complex model will further our understanding of the mechanism of the atmospheric/oceanic processes and improve their prediction.
基金supported by the Beijing Municipal Science and Technology Project [grant number 2171100004417008]the National Natural Science Foundation of China [grant numbers 41575050,41875049,and 41805034]
文摘The local convection initiation(CI)mechanisms of a convective case that occurred on5 August 2017 in Cangzhou,northern China,were studied using Doppler radar and automatic weather station observational analysis,along with Variational Doppler Radar Analysis System assimilation analysis.During the convective process,a gust front appeared ahead of two existing convective systems,respectively.In the warm and moist environment ahead of the gust fronts in the south,there was a mesoscale air mass boundary.With the process of a gust front moving southward,approaching the mesoscale air mass boundary,the convergence intensified in the area between the gust front and the mesoscale air mass boundary.Finally,the strong convergent updraft exceeded the level of free convection and triggered the new convection.
基金NSF grant ATM-9812729NOAA grant NA77WA0571Qiao is also supported by the Chinese National Key Basic Research Project under Contract G1999043809
文摘The impact of diabatic processes on 4-dimensional variational data assimilation (4D-Var) was studied using the 1995 version of NCEP's global spectral model with and without full physics.The adjoint was coded manually.A cost function measuring spectral errors of 6-hour forecasts to 'observation' (the NCEP reanalysis data) was minimized using the L-BFGS (the limited memory quasi-Newton algorithm developed by Broyden,Fletcher,Goldfard and Shanno) for optimizing parameters and initial conditions.Minimization of the cost function constrained by an adiabatic version of the NCEP global model converged to a minimum with a significant amount of decrease in the value of the cost function.Minimization of the cost function using the diabatic model, however,failed after a few iterations due to discontinuities introduced by physical parameterizations.Examination of the convergence of the cost function in different spectral domains reveals that the large-scale flow is adjusted during the first 10 iterations,in which discontinuous diabatic parameterizations play very little role.The adjustment produced by the minimization gradually moves to relatively smaller scales between 10-20th iterations.During this transition period,discontinuities in the cost function produced by 'on-off' switches in the physical parameterizations caused the cost function to stay in a shallow local minimum instead of continuously decreasing toward a deeper minimum. Next,a mixed 4D-Var scheme is tested in which large-scale flows are first adiabatically adjusted to a sufficient level,followed by a diabatic adjustment introduced after 10 to 20 iterations. The mixed 4D-Var produced a closer fit of analysis to observations,with 38% and 41% more decrease in the values of the cost function and the norm of gradient,respectively,than the standard diabatic 4D-Var,while the CPU time is reduced by 21%.The resulting optimal initial conditions improve the short-range forecast skills of 48-hour statistics.The detrimental effect of parameterization discontinuities on minimization was also reduced.