An Application of the Multi-Physics Ensemble Kalman Filter to Typhoon Forecast

This study examines the roles of the multi-physics approach in accounting for model errors for typhoon forecasts with the local ensemble transform Kalman filter (LETKF). Experiments with forecasts of Typhoon Conson (2010) using the weather research and forecasting (WRF) model show that use of the WRF’s multiple physical parameterization schemes to represent the model uncertainties can help the LETKF provide better forecasts of Typhoon Conson in terms of the forecast errors, the ensemble spread, the root mean square errors, the cross-correlation between mass and wind field as well as the coherent structure of the ensemble spread along the storm center. Sensitivity experiments with the WRF model show that the optimum number of the multi-physics ensemble is roughly equal to the number of combinations of different physics schemes assigned in the multi-physics ensemble. Additional idealized experiments with the Lorenz 40-variable model to isolate the dual roles of the multi-physics ensemble in correcting model errors and expanding the local ensemble space show that the multi-physics approach appears to be more essential in augmenting the local rank representation of the LETKF algorithm rather than directly accounting for model errors during the early cycles. The results in this study suggest that the multi-physics approach is a good option for short-range forecast applications with full physics models in which the spinup of the ensemble Kalman filter may take too long for the ensemble spread to capture efficiently model errors and cross-correlations among model variables.     

Conditional nonlinear optimal perturbation: Applications to stability,sensitivity, and predictability

Conditional nonlinear optimal perturbation (CNOP) is a nonlinear generalization of linear singular vector (LSV) and features the largest nonlinear evolution at prediction time for the initial perturbations in a given constraint. It was proposed initially for predicting the limitation of predictability of weather or climate. Then CNOP has been applied to the studies of the problems related to predictability for weather and climate. In this paper, we focus on reviewing the recent advances of CNOP’s applications, which involves the ones of CNOP in problems of ENSO amplitude asymmetry, block onset, and the sensitivity analysis of ecosystem and ocean’s circulations, etc. Especially, CNOP has been primarily used to construct the initial perturbation fields of ensemble forecasting, and to determine the sensitive area of target observation for precipitations. These works extend CNOP’s applications to investigating the nonlinear dynamical behaviors of atmospheric or oceanic systems, even a coupled system, and studying the problem of the transition between the equilibrium states. These contributions not only attack the particular physical problems, but also show the superiority of CNOP to LSV in revealing the effect of nonlinear physical processes. Consequently, CNOP represents the optimal precursors for a weather or climate event; in predictability studies, CNOP stands for the initial error that has the largest negative effect on prediction; and in sensitivity analysis, CNOP is the most unstable (sensitive) mode. In multi-equilibrium state regime, CNOP is the initial perturbation that induces the transition between equilibriums most probably. Furthermore, CNOP has been used to construct ensemble perturbation fields in ensemble forecast studies and to identify sensitive area of target observation. CNOP theory has become more and more substantial. It is expected that CNOP also serves to improve the predictability of the realistic predictions for weather and climate events plays an increasingly important role in exploring the nonlinear dynamics of atmospheric, oceanic and coupled atmosphere-ocean system. Supported by National Basic Research Program of China (Grant Nos. 2006CB403606, 2007CB411800), National Natural Science Foundation of China (Grant Nos. 40830955, 40675030, 40505013), Institute of Atmospheric Physics, Chinese Academy of Sciences (Grant No. IAP07202), and LASG State Key Laboratory Special Fund     

Ensemble prediction experiments using conditional nonlinear optimal perturbation

Two methods for initialization of ensemble forecasts are compared, namely, singular vector (SV) and conditional nonlinear optimal perturbation (CNOP). The comparison is done for forecast lengths of up to 10 days with a three-level quasi-geostrophic (QG) atmospheric model in a perfect model scenario. Ten cases are randomly selected from 1982/1983 winter to 1993/1994 winter (from December to the following February). Anomaly correlation coefficient (ACC) is adopted as a tool to measure the quality of the predicted ensembles on the Northern Hemisphere 500 hPa geopotential height. The results show that the forecast quality of ensemble samples in which the first SV is replaced by CNOP is higher than that of samples composed of only SVs in the medium range, based on the occurrence of weather regime transitions in Northern Hemisphere after about four days. Besides, the reliability of ensemble forecasts is evaluated by the Rank Histograms. The above conclusions confirm and extend those reached earlier by the authors, which stated that the introduction of CNOP improves the forecast skill under the condition that the analysis error belongs to a kind of fast-growing error by using a barotropic QG model. Supported by Knowledge Innovation Project of the Chinese Academy of Sciences (Grant No. KZCX3-SW-230), National Natural Science Foundation of China (Grant Nos. 40675030, 40633016)     

The application of the orthogonal conditional nonlinear optimal perturbations method to typhoon track ensemble forecasts

The orthogonal conditional nonlinear optimal perturbations (CNOPs) method, orthogonal singular vectors (SVs)method and CNOP+SVs method, which is similar to the orthogonal SVs method but replaces the leading SV (LSV) with the first CNOP, are adopted in both the Lorenz-96 model and Pennsylvania State University/National Center for Atmospheric Research (PSU/NCAR) Fifth-Generation Mesoscale Model (MM5) for ensemble forecasts. Using the MM5, typhoon track ensemble forecasting experiments are conducted for strong Typhoon Matsa in 2005. The results of the Lorenz-96 model show that the CNOP+SVs method has a higher ensemble forecast skill than the orthogonal SVs method, but ensemble forecasts using the orthogonal CNOPs method have the highest forecast skill. The results from the MM5 show that orthogonal CNOPs have a wider horizontal distribution and better describe the forecast uncertainties compared with SVs. When generating the ensemble mean forecast, equally averaging the ensemble members in addition to the anomalously perturbed forecast members may contribute to a higher forecast skill than equally averaging all of the ensemble members. Furthermore, for given initial perturbation amplitudes, the CNOP+SVs method may not have an ensemble forecast skill greater than that of the orthogonal SVs method, but the orthogonal CNOPs method is likely to have the highest forecast skill. Compared with SVs, orthogonal CNOPs fully consider the influence of nonlinear physical processes on the forecast results; therefore, considering the influence of nonlinearity may be important when generating fast-growing initial ensemble perturbations. All of the results show that the orthogonal CNOP method may be a potential new approach for ensemble forecasting.     

Short-term optimal operation of water systems using ensemble forecasts

Short-term water system operation can be realized using Model Predictive Control (MPC). MPC is a method for operational management of complex dynamic systems. Applied to open water systems, MPC provides integrated, optimal, and proactive management, when forecasts are available. Notwithstanding these properties, if forecast uncertainty is not properly taken into account, the system performance can critically deteriorate.Ensemble forecast is a way to represent short-term forecast uncertainty. An ensemble forecast is a set of possible future trajectories of a meteorological or hydrological system. The growing ensemble forecasts’ availability and accuracy raises the question on how to use them for operational management.The theoretical innovation presented here is the use of ensemble forecasts for optimal operation. Specifically, we introduce a tree based approach. We called the new method Tree-Based Model Predictive Control (TB-MPC). In TB-MPC, a tree is used to set up a Multistage Stochastic Programming, which finds a different optimal strategy for each branch and enhances the adaptivity to forecast uncertainty. Adaptivity reduces the sensitivity to wrong forecasts and improves the operational performance.TB-MPC is applied to the operational management of Salto Grande reservoir, located at the border between Argentina and Uruguay, and compared to other methods.     

High-resolution ensemble forecasting for the Gulf of Mexico eddies and fronts

High-resolution models and realistic boundary conditions are necessary to reproduce the mesoscale dynamics of the Gulf of Mexico (GOM). In order to achieve this, we use a nested configuration of the Hybrid Coordinate Ocean Model (HYCOM), where the Atlantic TOPAZ system provides lateral boundary conditions to a high-resolution (5 km) model of the GOM . However, such models cannot provide accurate forecasts of mesoscale variability, such as eddy shedding event, without data assimilation. Eddy shedding events involve the rapid growth of nonlinear instabilities that are difficult to forecast. The known sources of error are the initial state, the atmospheric condition, and the lateral boundary condition. We present here the benefit of using a small ensemble forecast (10 members) for providing confidence indices for the prediction, while using a data assimilation scheme based on optimal interpolation. Our set of initial states is provided by using different values of a data assimilation parameter, while the atmospheric and lateral boundary conditions are perturbed randomly. Changes in the data assimilation parameter appear to control the main position of the large features of the GOM in the initial state, whereas changes in the boundary conditions (lateral and atmospheric) appears to control the propagation of cyclonic eddies at their boundary. The ensemble forecast is tested for the shedding of Eddy Yankee (2006). The Loop Current and eddy fronts observed from ocean color and altimetry are almost always within the estimated positions from the ensemble forecast. The ensemble spread is correlated both in space and time to the forecast error, which implies that confidence indices can be provided in addition to the forecast. Finally, the ensemble forecast permits the optimization of a data assimilation parameter for best performance at a given forecast horizon.     

A SVD-based ensemble projection algorithm for calculating the conditional nonlinear optimal perturbation

Conditional nonlinear optimal perturbation(CNOP) is an extension of the linear singular vector technique in the nonlinear regime.It represents the initial perturbation that is subjected to a given physical constraint,and results in the largest nonlinear evolution at the prediction time.CNOP-type errors play an important role in the predictability of weather and climate.Generally,when calculating CNOP in a complicated numerical model,we need the gradient of the objective function with respect to the initial perturbations to provide the descent direction for searching the phase space.The adjoint technique is widely used to calculate the gradient of the objective function.However,it is difficult and cumbersome to construct the adjoint model of a complicated numerical model,which imposes a limitation on the application of CNOP.Based on previous research,this study proposes a new ensemble projection algorithm based on singular vector decomposition(SVD).The new algorithm avoids the localization procedure of previous ensemble projection algorithms,and overcomes the uncertainty caused by choosing the localization radius empirically.The new algorithm is applied to calculate the CNOP in an intermediate forecasting model.The results show that the CNOP obtained by the new ensemble-based algorithm can effectively approximate that calculated by the adjoint algorithm,and retains the general spatial characteristics of the latter.Hence,the new SVD-based ensemble projection algorithm proposed in this study is an effective method of approximating the CNOP.     

MOS based forecast of 6-hourly area precipitation

A statistical post-processing methodology for application to numerical weather prediction (NWP) model outputs for precipitation forecast is proposed. The post-processing is based on the model output statistics approach. The statistical relationships are described by the multiple linear regression model, which is complemented by an iteration procedure to further correct the regression outputs. Prognostic fields of the ALADIN/LACE (Aire Limitée Adaptation Dynamique Développement InterNational/Limited Area Modelling in Central Europe) NWP model are used for the forecast of 6-hourly areal precipitation amounts at 15 river basins. The NWP model integration starts at 00UTC and forecasts are calculated for lead times of +12, +18, +24 and +30 hours. The post-processing models are developed separately for each lead time and for separate warm (April to September) and cool (October to March) seasons. The forecasts are focused on large precipitation amounts. Using all the combinations, data from four years (1999–2002) are divided into calibration data (3 years), where the models are developed, and verification data. The models are evaluated by examining the root-mean-square error (RMSE), bias, and correlation coefficient (CC) on the verification data samples. The results show that the additional iteration procedure increases the forecast accuracy for a given range of precipitation amounts and simultaneously does not deteriorate the bias, a situation which can arise when negative regression outputs are set to zero. The post-processing method improves the forecast of the NWP model in terms of RMSE and CC. For large precipitation amounts during the summer season, the decrease of RMSE reaches 10% to 20% depending upon the applied method of verification. For the cool season, the decrease is somewhat smaller (7% to 15%).     

Improved method to derive equivalent current systems from global MHD simulations

Derivation of equivalent current systems(ECS)from a global magnetospheric magnetohydrodynamics(MHD)model is very useful in studying magnetosphere-ionosphere coupling,ground induction effects,and space weather forecast.In this study we introduce an improved method to derive the ECS from a global MHD model,which takes account of the obliqueness of the magnetic field lines.By comparing the ECS derived from this improved method and the previous method,we find that the main characteristics of the ECS derived from the two methods are generally consistent with each other,but the eastward-westward component of the geomagnetic perturbation calculated from the ECS derived from the improved method is much stronger than that from the previous method.We then compare the geomagnetic perturbation as a function of the interplanetary magnetic field(IMF)clock angle calculated from the ECS derived from both methods with the observations.The comparison indicates that the improved method can improve the performance of the simulation.Furthermore,it is found that the incomplete counterbalance of the geomagnetic effect produced by the ionospheric poloidal current and field-aligned current(FAC)contributes to most of the eastward-westward component of geomagnetic perturbation.     

Ensemble rainfall forecasting with numerical weather prediction and radar‐based nowcasting models

The overall objective of this study is to improve the forecasting accuracy of the precipitation in the Singapore region by means of both rainfall forecasting and nowcasting. Numerical Weather Predication (NWP) and radar‐based rainfall nowcasting are two important sources for quantitative precipitation forecast. In this paper, an attempt to combine rainfall prediction from a high‐resolution mesoscale weather model and a radar‐based rainfall model was performed. Two rainfall forecasting methods were selected and examined: (i) the weather research and forecasting model (WRF); and (ii) a translation model (TM). The WRF model, at a high spatial resolution, was run over the domain of interest using the Global Forecast System data as initializing fields. Some heavy rainfall events were selected from data record and used to test the forecast capability of WRF and TM. Results obtained from TM and WRF were then combined together to form an ensemble rainfall forecasting model, by assigning weights of 0.7 and 0.3 weights to TM and WRF, respectively. This paper presented results from WRF and TM, and the resulting ensemble rainfall forecasting; comparisons with station data were conducted as well. It was shown that results from WRF are very useful as advisory of anticipated heavy rainfall events, whereas those from TM, which used information of rain cells already appearing on the radar screen, were more accurate for rainfall nowcasting as expected. The ensemble rainfall forecasting compares reasonably well with the station observation data. Copyright © 2012 John Wiley & Sons, Ltd.     

The Local Ensemble Transform Kalman Filter with the Weather Research and Forecasting Model: Experiments with Real Observations

The local ensemble transform Kalman filter (LETKF) is implemented with the Weather Research and Forecasting (WRF) model, and real observations are assimilated to assess the newly-developed WRF-LETKF system. The WRF model is a widely-used mesoscale numerical weather prediction model, and the LETKF is an ensemble Kalman filter (EnKF) algorithm particularly efficient in parallel computer architecture. This study aims to provide the basis of future research on mesoscale data assimilation using the WRF-LETKF system, an additional testbed to the existing EnKF systems with the WRF model used in the previous studies. The particular LETKF system adopted in this study is based on the system initially developed in 2004 and has been continuously improved through theoretical studies and wide applications to many kinds of dynamical models including realistic geophysical models. Most recent and important improvements include an adaptive covariance inflation scheme which considers the spatial and temporal inhomogeneity of inflation parameters. Experiments show that the LETKF successfully assimilates real observations and that adaptive inflation is advantageous. Additional experiments with various ensemble sizes show that using more ensemble members improves the analyses consistently.     

A One-dimensional Ensemble Forecast and Assimilation System for Fog Prediction

A probabilistic fog forecast system was designed based on two high resolution numerical 1-D models called COBEL and PAFOG. The 1-D models are coupled to several 3-D numerical weather prediction models and thus are able to consider the effects of advection. To deal with the large uncertainty inherent to fog forecasts, a whole ensemble of 1-D runs is computed using the two different numerical models and a set of different initial conditions in combination with distinct boundary conditions. Initial conditions are obtained from variational data assimilation, which optimally combines observations with a first guess taken from operational 3-D models. The design of the ensemble scheme computes members that should fairly well represent the uncertainty of the current meteorological regime. Verification for an entire fog season reveals the importance of advection in complex terrain. The skill of 1-D fog forecasts is significantly improved if advection is considered. Thus the probabilistic forecast system has the potential to support the forecaster and therefore to provide more accurate fog forecasts.     

Towards improved post‐processing of hydrologic forecast ensembles

Forecast ensembles of hydrological and hydrometeorologial variables are prone to various uncertainties arising from climatology, model structure and parameters, and initial conditions at the forecast date. Post‐processing methods are usually applied to adjust the mean and variance of the ensemble without any knowledge about the uncertainty sources. This study initially addresses the drawbacks of a commonly used statistical technique, quantile mapping (QM), in bias correction of hydrologic forecasts. Then, an auxiliary variable, the failure index (γ), is proposed to estimate the ineffectiveness of the post‐processing method based on the agreement of adjusted forecasts with corresponding observations during an analysis period prior to the forecast date. An alternative post‐processor based on copula functions is then introduced such that marginal distributions of observations and model simulations are combined to create a multivariate joint distribution. A set of 2500 hypothetical forecast ensembles with parametric marginal distributions of simulated and observed variables are post‐processed with both QM and the proposed multivariate post‐processor. Deterministic forecast skills show that the proposed copula‐based post‐processing is more effective than the QM method in improving the forecasts. It is found that the performance of QM is highly correlated with the failure index, unlike the multivariate post‐processor. In probabilistic metrics, the proposed multivariate post‐processor generally outperforms QM. Further evaluation of techniques is conducted for river flow forecast of Sprague River basin in southern Oregon. Results show that the multivariate post‐processor performs better than the QM technique; it reduces the ensemble spread and is a more reliable approach for improving the forecast. Copyright © 2012 John Wiley & Sons, Ltd.     

The comparison of ensemble or deterministic dispersion modeling on global dispersion during Fukushima Dai-ichi nuclear accident

Ensemble forcasting,originally developed for weather prediction,is lately being extended to atmospheric dispersion applications,which is a new,effective methodology for improving the atmospheric dispersion numerical modeling.In March 2011,due to the massive 9.0 earthquakes and ensuing tsunami that struck off the northern coast of the island of Honshu,the Fukushima Nuclear Plant I had the substantial leak of radioactive materials into surrounding environment and atmosphere.To aim at the global dispersion modeling of atmospheric radionuclides from Fukushima Nuclear Accident,this paper presents two approaches of atmospheric dispersion forecasting:ensemble dispersion modeling(EDM) and deterministic dispersion modeling(DDM),conducts the globally dispersion modeling cases for Fukushima nuclear accident,and analyzes and evaluates the simulation results using observation data.In this paper,EDM includes three different perturbation methods:meteorological perturbation method,turbulence perturbation method,and physical parameterization ensemble forecasting method.The simulation results show that the trajectories from EDM have a better performance,which is in better agreement with the atmospheric circulation and observation data; the spread from DDM is slower and not as far as EDM.Additionally,the results from EDM display a better performance in the modeling of transport from Japan to China East Sea on April 4.The reasons for these results are:the techniques of MET and TUR are performed by adding perturbations on mean wind and turbulent velocity,respectively; the various different flow fields will result in far spreading in horizontal and the simulation results closer to observation; PHY is performed by using different diffusion physical parameterizations and produces the perturbations on vertical wind,which results the spreading in smaller range and discontinuous in horizontal.Finally,the comparative analysis between modeling results and observation data shows that all cases results are in good agreement with trends of observed radionuclides surface concentration; however,the modeling surface concentration is smaller than observation,especially in DDM and PHY.Furthermore,the EDM results show that MET and TUR are of more evolutionary advantage than PHY in modeling of average and maximum concentration.Therefore,this study can serve as a reference to atmospheric dispersion and environmental emergency response(EER).     

Seasonal Sensitivity on COBEL-ISBA Local Forecast System for Fog and Low Clouds

Skillful low visibility forecasts are essential for air-traffic managers to effectively regulate traffic and to optimize air-traffic control at international airports. For this purpose, the COBEL-ISBA local numerical forecast system has been implemented at Paris CDG international airport. This local approach is robust owing to the assimilation of detailed local observations. However, even with dedicated observations and initialization, uncertainties remain in both initial conditions and mesoscale forcings. The goal of the research presented here is to address the sensitivity of COBEL-ISBA forecast to initial conditions and mesoscale forcings during the winter season 2002–2003. The main sources of uncertainty of COBEL-ISBA input parameters have been estimated and the evaluation of parameter uncertainty on the forecasts has been studied. A budget strategy is applied during the winter season to quantify COBEL-ISBA sensitivity. This study is the first step toward building a local ensemble prediction system based on COBEL-ISBA. The conclusions of this work point out the potential for COBEL-ISBA ensemble forecasting and quantify sources of uncertainty that lead to dispersion.     

Prediction and predictability of a catastrophic local extreme precipitation event through cloud-resolving ensemble analysis and forecasting with Doppler radar observations

Local extreme rain usually resulted in disasters such as flash floods and landslides. Upon today, it is still one of the most difficult tasks for operational weather forecast centers to predict those events accurately. In this paper, we simulate an extreme precipitation event with ensemble Kalman filter (EnKF) assimilation of Doppler radial-velocity observations, and analyze the uncertainties of the assimilation. The results demonstrate that, without assimilation radar data, neither a single initialization of deterministic forecast nor an ensemble forecast with adding perturbations or multiple physical parameterizations can predict the location of strong precipitation. However, forecast was significantly improved with assimilation of radar data, especially the location of the precipitation. The direct cause of the improvement is the buildup of a deep mesoscale convection system with EnKF assimilation of radar data. Under a large scale background favorable for mesoscale convection, efficient perturbations of upstream mid-low level meridional wind and moisture are key factors for the assimilation and forecast. Uncertainty still exists for the forecast of this case due to its limited predictability. Both the difference of large scale initial fields and the difference of analysis obtained from EnKF assimilation due to small amplitude of initial perturbations could have critical influences to the event's prediction. Forecast could be improved through more cycles of EnKF assimilation. Sensitivity tests also support that more accurate forecasts are expected through improving numerical models and observations.     

Development of a short-term river flood forecasting method for snowmelt driven floods based on wavelet and cross-wavelet analysis

In this study, a new method of stand-alone short-term spring snowmelt river flood forecasting was developed based on wavelet and cross-wavelet analysis. Wavelet and cross-wavelet analysis were used to decompose flow and meteorological time series data and to develop wavelet based constituent components which were then used to forecast floods 1, 2, and 6 days ahead. The newly developed wavelet forecasting method (WT) was compared to multiple linear regression analysis (MLR), autoregressive integrated moving average analysis (ARIMA), and artificial neural network analysis (ANN) for forecasting daily stream flows with lead-times equal to 1, 2, and 6 days. This comparison was done using data from the Rideau River watershed in Ontario, Canada. Numerical analysis was performed on daily maximum stream flow data from the Rideau River station and on meteorological data (rainfall, snowfall, and snow on ground) from the Ottawa Airport weather station. Data from 1970 to 1997 were used to train the models while data from 1998 to 2001 were used to test the models. The most significant finding of this research was that it was demonstrated that the proposed wavelet based forecasting method can be used with great accuracy as a stand-alone forecasting method for 1 and 2 days lead-time river flood forecasting, assuming that there are no significant trends in the amplitude for the same Julian day year-to-year, and that there is a relatively stable phase shift between the flow and meteorological time series. The best forecasting model for 1 day lead-time was a wavelet analysis model. In testing, it had the lowest RMSE value (13.8229), the highest R² value (0.9753), and the highest EI value (0.9744). The best forecasting model for 2 days lead-time was also a wavelet analysis model. In testing, it had the lowest RMSE value (31.7985), the highest R² value (0.8461), and the second highest EI value (0.8410). It was also shown that the proposed wavelet based forecasting method is not particularly accurate for longer lead-time forecasting such as 6 days, with the ANN method providing more accurate results. The best forecasting model for 6 days lead-time was an ANN model, with the wavelet model not performing as well. In testing, the wavelet model had an RMSE of 57.6917, an R² of 0.4835, and an EI of 0.4366.     

A variational method for correcting non-systematic errors in numerical weather prediction

A variational method based on previous numerical forecasts is developed to estimate and correct non-systematic component of numerical weather forecast error. In the method, it is assumed that the error is linearly dependent on some combination of the forecast fields, and three types of forecast combination are applied to identifying the forecasting error: 1) the forecasts at the ending time, 2) the combination of initial fields and the forecasts at the ending time, and 3) the combination of the forecasts at...     

Assimilating Doppler radar observations with an ensemble Kalman filter for convection-permitting prediction of convective development in a heavy rainfall event during the pre-summer rainy season of South China

This study examines the effectiveness of an ensemble Kalman filter based on the weather research and forecasting model to assimilate Doppler-radar radial-velocity observations for convection-permitting prediction of convection evolution in a high-impact heavy-rainfall event over coastal areas of South China during the pre-summer rainy season. An ensemble of 40 deterministic forecast experiments(40 DADF) with data assimilation(DA) is conducted, in which the DA starts at the same time but lasts for different time spans(up to 2 h) and with different time intervals of 6, 12, 24, and 30 min. The reference experiment is conducted without DA(NODA).To show more clearly the impact of radar DA on mesoscale convective system(MCS)forecasts, two sets of 60-member ensemble experiments(NODA EF and exp37 EF) are performed using the same 60-member perturbed-ensemble initial fields but with the radar DA being conducted every 6 min in the exp37 EF experiments from 0200 to0400 BST. It is found that the DA experiments generally improve the convection prediction. The 40 DADF experiments can forecast a heavy-rain-producing MCS over land and an MCS over the ocean with high probability, despite slight displacement errors. The exp37 EF improves the probability forecast of inland and offshore MCSs more than does NODA EF. Compared with the experiments using the longer DA time intervals, assimilating the radial-velocity observations at 6-min intervals tends to produce better forecasts. The experiment with the longest DA time span and shortest time interval shows the best performance.However, a shorter DA time interval(e.g., 12 min) or a longer DA time span does not always help. The experiment with the shortest DA time interval and maximum DA window shows the best performance, as it corrects errors in the simulated convection evolution over both the inland and offshore areas. An improved representation of the initial state leads to dynamic and thermodynamic conditions that are more conducive to earlier initiation of the inland MCS and longer maintenance of the offshore MCS.     

A soil moisture assimilation scheme based on the ensemble Kalman filter using microwave brightness temperature

This study presents a soil moisture assimilation scheme, which could assimilate microwave brightness temperature directly, based on the ensemble Kalman filter and the shuffled complex evolution method (SCE-UA). It uses the soil water model of the land surface model CLM3.0 as the forecast operator, and a radiative transfer model (RTM) as the observation operator in the assimilation system. The assimilation scheme is implemented in two phases: the parameter calibration phase and the pure soil moisture assimilation phase. The vegetation optical thickness and surface roughness parameters in the RTM are calibrated by SCE-UA method and the optimal parameters are used as the final model parameters of the observation operator in the assimilation phase. The ideal experiments with synthetic data indicate that this scheme could significantly improve the simulation of soil moisture at the surface layer. Furthermore, the estimation of soil moisture in the deeper layers could also be improved to a certain extent. The real assimilation experiments with AMSR-E brightness temperature at 10.65 GHz (vertical polarization) show that the root mean square error (RMSE) of soil moisture in the top layer (0–10 cm) by assimilation is 0.03355 m³ · m⁻³, which is reduced by 33.6% compared with that by simulation (0.05052 m³ · m⁻³). The mean RMSE by assimilation for the deeper layers (10–50 cm) is also reduced by 20.9%. All these experiments demonstrate the reasonability of the assimilation scheme developed in this study.