Qiu Chong-jian

Development of a Chinese land data assimilation system: its progress and prospects

Abstract The Objective of land data assimilation is to merge multi-source observations into the dynamics of land surface model for improving the estimation of land surface states. We have developed a land data assimilation system for Chinas land territory. In this system, the Common Land Model and Simple Biosphere Model 2 are used to simulate land surface processes. The radiative transfer models of thawed and frozen soil, snow, lake, and vegetations are used as observation operators to transfer model predictions into estimated brightness temperatures. A Monte-Carlo sequential filter, the ensemble Kalman filter, is implemented as data assimilation method to integrate modeling and observation. The system is capable of assimilating passive microwave remotely sensed data such as special sensor microwave imager (SSM I), TRMM microwave imager (TMI), and advanced microwave scanning radiometer enhanced for EOS (AMSRE) and the conventional in situ measurements of soil and snow. A spatiotemporally consistent assim...

Estimating the soil moisture profile by assimilating near-surface observations with the ensemble Kalman filter (EnKF)

The paper investigates the ability to retrieve the true soil moisture profile by assimilating near-surface soil moisture into a soil moisture model with an ensemble Kaiman filter (EnKF) assimilation scheme, including the effect of ensemble size, update interval and nonlinearities in the profile retrieval, the required time for full retrieval of the soil moisture profiles, and the possible influence of the depth of the soil moisture observation. These questions are addressed by a desktop study using synthetic data. The “true” soil moisture profiles are generated from the soil moisture model under the boundary condition of 0.5 cm d−1 evaporation. To test the assimilation schemes, the model is initialized with a poor initial guess of the soil moisture profile, and different ensemble sizes are tested showing that an ensemble of 40 members is enough to represent the covariance of the model forecasts. Also compared are the results with those from the direct insertion assimilation scheme, showing that the EnKF is superior to the direct insertion assimilation scheme, for hourly observations, with retrieval of the soil moisture profile being achieved in 16 h as compared to 12 days or more. For daily observations, the true soil moisture profile is achieved in about 15 days with the EnKF, but it is impossible to approximate the true moisture within 18 days by using direct insertion. It is also found that observation depth does not have a significant effect on profile retrieval time for the EnKF. The nonlinearities have some negative influence on the optimal estimates of soil moisture profile but not very seriously.

An improved method for doppler wind and thermodynamic retrievals

A variational method is developed to retrieve winds in the first step and then thermodynamic fields in the second step from Doppler radar observations. In the first step, wind fields are retrieved at two time levels: the beginning and ending times of the data assimilation period, simultaneously from two successive volume scans by using the weak form constraints provided by the mass continuity and vorticity equations. As the retrieved wind fields are expressed by Legendre polynomial expansions at the beginning and ending times, the time tendency term in the vorticity equation can be conveniently formulated, and the retrieved winds can be compared with the radar observed radial winds in the cost function at the precise time and position of each radar beam. In the second step, the perturbation pressure and temperature fields at the middle time are then derived from the retrieved wind fields and the velocity time tendency by using the weak form constraints provided by the three momentum equations. The merits of the new method are demonstrated by numerical experiments with simulated radar observations and compared with the traditional least squares methods which consider neither the precise observation times and positions nor the velocity time tendency. The new method is also applied to real radar data for a heavy rainfall event during the 2001 Meiyu season in China.

Estimating heat fluxes by merging profile formulae and the energy budget with a variational technique

A variational technique (VT) is applied to estimate surface sensible and latent heat fluxes based on observations of air temperature, wind speed, and humidity, respectively, at three heights (1 m, 4 m, and 10 m), and the surface energy and radiation budgets by the surface energy and radiation system (SERBS). The method fully uses all information provided by the measurements of air temperature, wind, and humidity profiles, the surface energy budget, and the similarity profile formulae as well. Data collected at Feixi experiment station installed by the China Heavy Rain Experiment and Study (HeRES) Program are used to test the method. Results show that the proposed technique can overcome the well-known unstablility problem that occurs when the Bowen method becomes singular; in comparison with the profile method, it reduces both the sensitivities of latent heat fluxes to observational errors in humidity and those of sensible heat fluxes to observational errors in temperature, while the estimated heat fluxes approximately satisfy the surface energy budget. Therefore, the variational technique is more reliable and stable than the two conventional methods in estimating surface sensible and latent heat fluxes.

Kinematic structure of a heavy rain event from dual-Doppler radar observations

The detailed kinematic structure of a heavy rain event that occurred in the middle reaches of the Yangtze River was investigated using dual-Doppler radar observation. A variational analysis method was developed to obtain the three-dimensional wind fields. Before the analysis, a data preprocessing procedure was carried out, in which the temporal variation with the scanning time interval and the effect of the earth curvature on the data position were taken into account. The analysis shows that a shear line in the lower and middle levels played an important role in the rainfall event. The precipitation fell mainly on the south end of the shear line where southerly flow prevailed and convergence and updraft were obvious. With the movement and decay of the shear line, the precipitation moved and decayed correspondingly.