Data assimilation (DA) has been widely used in land surface models (LSM) to improve model state estimates. Among various DA methods, the particle filter (PF) with Markov chain Monte Carlo (MCMC) has become increasingly popular for estimating the states of the nonlinear and non-Gaussian LSMs. However, the standard PF always suffers from the particle impoverishment problem, characterized by loss of particle diversity. To solve this problem, an immune evolution particle filter with MCMC simulation inspired by the biological immune system, entitled IEPFM, is proposed for DA in this paper. The merit of this approach is in imitating the antibody diversity preservation mechanism to further improve particle diversity, thus increasing the accuracy of estimates. Furthermore, the immune memory function refers to promise particle evolution process towards optimal estimates. Effectiveness of the proposed approach is demonstrated by the numerical simulation experiment using a highly nonlinear atmospheric model. Finally, IEPFM is applied to a soil moisture (SM) assimilation experiment, which assimilates in situ observations into the Variable Infiltration Capacity (VIC) model to estimate SM in the MaQu network region of the Tibetan Plateau. Both synthetic and real case experiments demonstrate that IEPFM mitigates particle impoverishment and provides more accurate assimilation results compared with other popular DA algorithms.

Hydrological models play an essential role in data assimilation (DA) systems. However, it is a challenging task to acquire the distributed hydrological model parameters that affect the accuracy of the simulations at a grid scale. Remote sensing data provide an ideal observation for DA to estimate parameters and state variables. In this study, a special assimilation scheme was proposed to jointly estimate parameters and soil moisture (SM) by assimilating brightness temperature (TB) from the Soil Moisture and Ocean Salinity (SMOS) mission. Variable infiltration capacity (VIC) hydrological model and L-band microwave emission of the biosphere model (L-MEB) are coupled as model and observation operators, respectively. The scheme combines two stages of estimators, one for the static model parameters and the other for the dynamic state variables. The estimators approximate the posterior probability distribution of an unknown target through sequential Monte Carlo (SMC) sampling. Markov chain Monte Carlo (MCMC) and immune evolution strategy are embedded in both stages to solve particle impoverishment problems. To evaluate the effectiveness of the scheme, the estimated SM sets are compared with in-situ observations and SMOS products in Maqu on the Tibetan Plateau. Specifically, the root mean square error decreased from 0.126 to 0.087 m3m−3 for surface SM, with a slight impact on the root zone. The temporal correlation between DA results and in-situ measurements increased to 0.808 and 0.755 for surface SM (+0.057) and root zone SM (+0.040), respectively. The results demonstrate that assimilating TB has tremendous potential as an approach to improve the estimation of distributed model parameters and SMs of surface and root zone at a grid scale, and the immune evolution strategy is effective for increasing the accuracy of approximation in sampling.

Remote sensing images are widely used in earth observation. However, the existence of clouds seriously affects the interpretation of remote sensing images. In order to improve the accuracy of cloud detection, it’s usually necessary to complete the segmentation of the cloud boundary before cloud detection. Based on the Simple Linear Iterative Clustering (SLIC) algorithm, an improved SLIC superpixels segmentation method based on density feature is proposed to realize the segmentation of cloud remote sensing images. First, to generate the initial clustering center , the density peak clustering method is used instead of the uniform setting method. Then, in the calculation of distance measurement, we added the local density distance term. Finally, we get the superpixels segmented image by iteration. Four remote sensing images with different underlying surfaces were selected as the test data. The comparison experiments with other two algorithms show that the algorithm promoted in this paper shows superior performance in boundary recall (BR) rate and the error rate is lower than other algorithms in cloud detection.