Data assimilation methods in the Earth sciences

Although remote sensing data are often plentiful, they do not usually satisfy the users’ needs directly. Data assimilation is required to extract information about geophysical fields of interest from the remote sensing observations and to make the data more accessible to users. Remote sensing may provide, for example, measurements of surface soil moisture, snow water equivalent, snow cover, or land surface (skin) temperature. Data assimilation can then be used to estimate variables that are not directly observed from space but are needed for applications, for instance root zone soil moisture or land surface fluxes. The paper provides a brief introduction to modern data assimilation methods in the Earth sciences, their applications, and pertinent research questions. Our general overview is readily accessible to hydrologic remote sensing scientists. Within the general context of Earth science data assimilation, we point to examples of the assimilation of remotely sensed observations in land surface hydrology.     

Comparison of land surface scheme simulations with field observations versus atmospheric model output as forcing

ABSTRACT

The low density of meteorological stations in parts of Canada necessitates using numerical weather prediction (NWP)/assimilation output for hydrological modelling. In this study, comparisons are made of simulated land surface variables when using field observations versus NWP output as forcing for two well-instrumented sites: the mountainous and forested Marmot Creek Basin (MCRB) in the Canadian Rocky Mountains, and a prairie cropland/grassland site (Kenaston). The Canadian Land Surface Scheme 3.6 (CLASS) was used for modelling. The Global Environmental Multiscale (GEM) model with Canadian Precipitation Analysis (CaPA) was also used as forcing. There was good agreement between observed meteorology and GEM/CaPA, though some deficiencies in GEM/CaPA were identified: the effects of sub-grid topography on incoming radiation and wind speed were not accounted for at MCRB, and CaPA did not capture some convective rainfall events at Kenaston. CLASS simulations using both sets of forcing showed difficulties in simulating snow depth, soil moisture and evapotranspiration; certain difficulties were linked to GEM/CaPA deficiencies and/or CLASS. Both sets of forcing tended to overestimate the duration of snow cover at MCRB, but during different years. With GEM/CaPA as forcing, CLASS overestimated the duration of frozen soils. The GEM/CaPA precipitation difficulties at Kenaston degraded soil moisture simulations.

EDITOR A. Castellarin; ASSOCIATE EDITOR E. Volpi     

Assimilation and High Resolution Forecasts of Surface and Near Surface Conditions for the 2010 Vancouver Winter Olympic and Paralympic Games

A dynamical model was experimentally implemented to provide high resolution forecasts at points of interests in the 2010 Vancouver Olympics and Paralympics Region. In a first experiment, GEM-Surf, the near surface and land surface modeling system, is driven by operational atmospheric forecasts and used to refine the surface forecasts according to local surface conditions such as elevation and vegetation type. In this simple form, temperature and snow depth forecasts are improved mainly as a result of the better representation of real elevation. In a second experiment, screen level observations and operational atmospheric forecasts are blended to drive a continuous cycle of near surface and land surface hindcasts. Hindcasts of the previous day conditions are then regarded as today’s optimized initial conditions. Hence, in this experiment, given observations are available, observation driven hindcasts continuously ensure that daily forecasts are issued from improved initial conditions. GEM-Surf forecasts obtained from improved short-range hindcasts produced using these better conditions result in improved snow depth forecasts. In a third experiment, assimilation of snow depth data is applied to further optimize GEM-Surf’s initial conditions, in addition to the use of blended observations and forecasts for forcing. Results show that snow depth and summer temperature forecasts are further improved by the addition of snow depth data assimilation.     

Connecting Satellite Observations with Water Cycle Variables Through Land Data Assimilation: Examples Using the NASA GEOS-5 LDAS

A land data assimilation system (LDAS) can merge satellite observations (or retrievals) of land surface hydrological conditions, including soil moisture, snow, and terrestrial water storage (TWS), into a numerical model of land surface processes. In theory, the output from such a system is superior to estimates based on the observations or the model alone, thereby enhancing our ability to understand, monitor, and predict key elements of the terrestrial water cycle. In practice, however, satellite observations do not correspond directly to the water cycle variables of interest. The present paper addresses various aspects of this seeming mismatch using examples drawn from recent research with the ensemble-based NASA GEOS-5 LDAS. These aspects include (1) the assimilation of coarse-scale observations into higher-resolution land surface models, (2) the partitioning of satellite observations (such as TWS retrievals) into their constituent water cycle components, (3) the forward modeling of microwave brightness temperatures over land for radiance-based soil moisture and snow assimilation, and (4) the selection of the most relevant types of observations for the analysis of a specific water cycle variable that is not observed (such as root zone soil moisture). The solution to these challenges involves the careful construction of an observation operator that maps from the land surface model variables of interest to the space of the assimilated observations.     

A land surface data assimilation framework using the land information system: Description and applications

The Land Information System (LIS) is an established land surface modeling framework that integrates various community land surface models, ground measurements, satellite-based observations, high performance computing and data management tools. The use of advanced software engineering principles in LIS allows interoperability of individual system components and thus enables assessment and prediction of hydrologic conditions at various spatial and temporal scales. In this work, we describe a sequential data assimilation extension of LIS that incorporates multiple observational sources, land surface models and assimilation algorithms. These capabilities are demonstrated here in a suite of experiments that use the ensemble Kalman filter (EnKF) and assimilation through direct insertion. In a soil moisture experiment, we discuss the impact of differences in modeling approaches on assimilation performance. Provided careful choice of model error parameters, we find that two entirely different hydrological modeling approaches offer comparable assimilation results. In a snow assimilation experiment, we investigate the relative merits of assimilating different types of observations (snow cover area and snow water equivalent). The experiments show that data assimilation enhancements in LIS are uniquely suited to compare the assimilation of various data types into different land surface models within a single framework. The high performance infrastructure provides adequate support for efficient data assimilation integrations of high computational granularity.     

Numerical simulation and data assimilation of the water-energy cycle over semiarid northeastern China

The default fractional vegetation cover and terrain height were replaced by the estimated fractional vegetation cover, which was calculated by the Normalized Difference Vegetation Index(NDVI) of Earth Observing System Moderate-Resolution Imaging Spectroradiometer(EOS-MODIS) and the Digital Elevation Model of the Shuttle Radar Topography Mission(SRTM) system. The near-surface meteorological elements over northeastern China were assimilated into the three-dimensional variational data assimilation system(3DVar) module in the Weather Research and Forecasting(WRF) model. The structure and daily variations of air temperature, humidity, wind and energy fields over northeastern China were simulated using the WRF model. Four groups of numerical experiments were performed, and the simulation results were analyzed of latent heat flux, sensible heat flux, and their relationships with changes in the surface energy flux due to soil moisture and precipitation over different surfaces. The simulations were compared with observations of the stations Tongyu, Naiman, Jinzhou, and Miyun from June to August, 2009. The results showed that the WRF model achieves high-quality simulations of the diurnal characteristics of the surface layer temperature, wind direction, net radiation, sensible heat flux, and latent heat flux over semiarid northeastern China in the summer. The simulated near-surface temperature, relative humidity, and wind speed were improved in the data assimilation case(Case 2) compared with control case(Case 1). The simulated sensible heat fluxes and surface heat fluxes were improved by the land surface parameterization case(Case 3) and the combined case(Case 4). The simulated temporal variations in soil moisture over the northeastern arid areas agree well with observations in Case 4, but the simulated precipitation should be improved in the WRF model. This study could improve the land surface parameters by utilizing remote sensing data and could further improve atmospheric elements with a data assimilation system. This work provides an effective attempt at combining multi-source data with different spatial and temporal scales into numerical simulations. The assimilation datasets generated by this work can be applied to research on climate change and environmental monitoring of arid lands, as well as research on the formation and stability of climate over semiarid areas.     

青藏高原春季积雪在南海夏季风爆发过程中的作用

本文应用欧洲中期预报中心（ECMWF,European Centre for Medium\|Range Weather Forecasts—ERA\|40）资料和美国国家环境预测中心和国家大气研究中心(NCEP/NCAR, National Centers for Environmental Prediction/National Center for Atmospheric Research)资料,研究了青藏高原雪深变化对南海夏季风爆发的影响和ENSO对青藏高原降雪的影响.结果表明：（1）ECMWF的雪深资料是可信的,可以用来研究青藏高原雪深变化对南海夏季风爆发的影响;（2）青藏高原的积雪异常影响到500 hPa以上的温度异常和印度洋与大陆间的气温对比,一方面使上层的南亚高压移动速度发生变化,另一方面也影响到低层大气的运动和东西向风异常,在青藏高原少雪年,东印度洋产生西风异常和一个气旋对,而在青藏高原多雪年,东印度洋产生东风异常和一个反气旋对;（3）ENSO与青藏高原春季积雪关系密切.东太平洋SST正异常时,东印度洋和南海气压偏高,从而导致该区海陆经向压强梯度增强和西风异常.另外,此时青藏高原北部气压偏高,北风偏强,副热带锋面增强,同时,印度洋的SST偏高,为青藏高原降雪提供了水汽保障,这些都有利于青藏高原的降雪.     

Is the soil moisture precipitation feedback enhanced by heterogeneity and dry soils? A comparative study

The interaction between the land surface and the atmosphere is a crucial driver of atmospheric processes. Soil moisture and precipitation are key components in this feedback. Both variables are intertwined in a cycle, that is, the soil moisture – precipitation feedback for which involved processes and interactions are still discussed. In this study the soil moisture – precipitation feedback is compared for the sempiternal humid Ammer catchment in Southern Germany and for the semiarid to subhumid Sissili catchment in West Africa during the warm season, using precipitation datasets from the Climate Hazards Group InfraRed Precipitation with Station data (CHIRPS), from the German Weather Service (REGNIE) and simulation datasets from the Weather Research and Forecasting (WRF) model and the hydrologically enhanced WRF-Hydro model. WRF and WRF-Hydro differ by their representation of terrestrial water flow. With this setup we want to investigate the strength, sign and variables involved in the soil moisture – precipitation feedback for these two regions. The normalized model spread between the two simulation results shows linkages between precipitation variability and diagnostic variables surface fluxes, moisture flux convergence above the surface and convective available potential energy in both study regions. The soil moisture – precipitation feedback is evaluated with a classification of soil moisture spatial heterogeneity based on the strength of the soil moisture gradients. This allows us to assess the impact of soil moisture anomalies on surface fluxes, moisture flux convergence, convective available potential energy and precipitation. In both regions the amount of precipitation generally increases with soil moisture spatial heterogeneity. For the Ammer region the soil moisture – precipitation feedback has a weak negative sign with more rain near drier patches while it has a positive signal for the Sissili region with more rain over wetter patches. At least for the observed moderate soil moisture values and the spatial scale of the Ammer region, the spatial variability of soil moisture is more important for surface-atmosphere interactions than the actual soil moisture content. Overall, we found that soil moisture heterogeneity can greatly affect the soil moisture – precipitation feedback.     

GPS大气掩星技术在全球气候变化研究中的应用

人类活动引起全球变暖,衡量全球气候变化的指标有陆地、大气和海洋温度,水汽含量等等.研究对流层底层大气温度和水汽含量变化的传统方法是用数值天气预报模型和微波声纳,尚未实现用全球均匀覆盖的数据来做精确的定量研究.和GNSS系列卫星计划比较,最近发射的COSMIC卫星气象探测数据的空间、时间以及垂直分辨率都大大提高.采用COSMIC数据可以改进和量化南极洲的大气压力模型,并综合GNSS系列卫星测量的水汽和温度剖面研究全球气候变化.用一维协方差算法估计南极洲及附近海洋的大气压、温度和湿度剖面.把COSMIC卫星密集测量期间演算得到的大气折射率和GNSS系列卫星的结果进行比较.再和独立测量数据进行比较,包括南极洲自动气象观测站资料,数值天气预报模型资料,多种测高卫星水汽资料和海洋表面温度资料以及区域GPS水汽图.上述工作将改进发展中的气象遥感技术并应用于天气预报和空间天气预报及全球气候变化研究.     

Remote sensing in hydrology

Remote sensing provides a means of observing hydrological state variables over large areas. The ones which we will consider in this paper are land surface temperature from thermal infrared data, surface soil moisture from passive microwave data, snow cover using both visible and microwave data, water quality using visible and near-infrared data and estimating landscape surface roughness using lidar. Methods for estimating the hydrometeorlogical fluxes, evapotranspiration and snowmelt runoff, using these state variables are also described.     

One-dimensional soil moisture profile retrieval by assimilation of near-surface observations: a comparison of retrieval algorithms

This paper investigates the ability to retrieve the true soil moisture and temperature profiles by assimilating near-surface soil moisture and surface temperature data into a soil moisture and heat transfer model. The direct insertion and Kalman filter assimilation schemes have been used most frequently in assimilation studies, but no comparisons of these schemes have been made. This study investigates which of these approaches is able to retrieve the soil moisture and temperature profiles the fastest, over what depth soil moisture observations are required, and the effect of update interval on profile retrieval. These questions are addressed by a desktop study using synthetic data. The study shows that the Kalman filter assimilation scheme is superior to the direct insertion assimilation scheme, with retrieval of the soil moisture profile being achieved in 12 h as compared to 8 days or more, depending on observation depth, for hourly observations. It was also found that profile retrieval could not be realised for direct insertion of the surface node alone, and that observation depth does not have a significant effect on profile retrieval time for the Kalman filter. The observation interval was found to be unimportant for profile retrieval with the Kalman filter when the forcing data is accurate, whilst for direct insertion the continuous Dirichlet boundary condition was required for an increasingly longer period of time. It was also found that the Kalman filter assimilation scheme was less susceptible to unstable updates if volumetric soil moisture was modelled as the dependent state rather than matric head, because the volumetric soil moisture state is more linear in the forecasting model.     

积雪在El Niño影响东亚夏季气候异常中的作用

本文利用1948-2010年Global Land Data Assimilation System（GLDAS）NOAH陆面模式资料、GPCC月平均降水资料和NCAR/NCEP全球月平均再分析资料,采用滤波、距平合成和线性相关等方法,分析了El Niño成熟位相冬季欧亚大陆积雪异常的分布特征,研究了关键区积雪融化对后期春、夏季土壤湿度、土壤温度以及大气环流与降水的影响,揭示了El Niño事件通过关键区积雪储存其强迫信号并影响东亚夏季气候异常的机制和过程.主要结论如下：El Niño成熟阶段冬季伊朗高原、巴尔喀什湖东北部和青藏高原南麓区域是雪深异常的三个关键区,这些区域的雪深、雪融和土壤湿度有明显的正相关;这三个关键区雪深异常通过春季融雪将冬季El Niño信号传递给春、夏季局地土壤湿度,通过减少感热通量和增加潜热通量对大气环流产生影响;春末夏初伊朗高原土壤湿度异常对东亚夏季气候异常的影响最大,其引起的降水异常与El Niño次年夏季降水异常分布基本一致,春夏季青藏高原南麓和巴尔喀什湖附近土壤湿度也都明显增加,均会对中国华北降水增加有显著正贡献.总之,在利用El Niño事件研究和预测东亚夏季气候异常时,还应考虑关键区雪深异常对El Niño信号的存储和调制作用.     

Improving the estimation of hydrological states in the SWAT model via the ensemble Kalman smoother: Synthetic experiments for the Heihe River Basin in northwest China

Data assimilation as a method to predict variables, reduce uncertainties and explicitly handle various sources of uncertainties has recently received widespread attention and has been utilized to combine in situ and remotely sensed measurements with hydrological models. However, factors that significantly influence the capability of data assimilation still need testing and verifying. In this paper, synthetic surface soil moisture data are assimilated into the Soil and Water Assessment Tool (SWAT) model to evaluate their impact on other hydrological variables via the ensemble Kalman smoother (EnKS), using data from the Heihe River Basin, northwest China. The results show that the assimilation of surface soil moisture can moderately improve estimates of deep layer soil moisture, surface runoff and lateral flow, which reduces the negative influences of erroneous forcing and inaccurate parameters. The effects of the spatially heterogeneous input data (land cover and soil type) on the performance of the data assimilation technique are noteworthy. Moreover, the approaches including inflation and localization are specifically diagnosed to further extend the capability of the EnKS.     

A New Integrated Weighted Model in SNOW-V10: Verification of Continuous Variables

This paper presents the verification results of nowcasts of four continuous variables generated from an integrated weighted model and underlying Numerical Weather Prediction (NWP) models. Real-time monitoring of fast changing weather conditions and the provision of short term forecasts, or nowcasts, in complex terrain within coastal regions is challenging to do with sufficient accuracy. A recently developed weighting, evaluation, bias correction and integration system was used in the Science of Nowcasting Olympic Weather for Vancouver 2010 project to generate integrated weighted forecasts (INTW) out to 6 h. INTW forecasts were generated with in situ observation data and background gridded forecasting data from Canadian high-resolution deterministic NWP system with three nested grids at 15-, 2.5- and 1-km horizontal grid-spacing configurations. In this paper, the four variables of temperature, relative humidity, wind speed and wind gust are treated as continuous variables for verifying the INTW forecasts. Fifteen sites were selected for the comparison of the model performances. The results of the study show that integrating surface observation data with the NWP forecasts produce better statistical scores than using either the NWP forecasts or an objective analysis of observed data alone. Overall, integrated observation and NWP forecasts improved forecast accuracy for the four continuous variables. The mean absolute errors from the INTW forecasts for the entire test period (12 February to 21 March 2010) are smaller than those from NWP forecasts with three configurations. The INTW is the best and most consistent performer among all models regardless of location and variable analyzed.     

Dual state-parameter estimation of root zone soil moisture by optimal parameter estimation and extended Kalman filter data assimilation

With well-determined hydraulic parameters in a hydrologic model, a traditional data assimilation method (such as the Kalman filter and its extensions) can be used to retrieve root zone soil moisture under uncertain initial state variables (e.g., initial soil moisture content) and good simulated results can be achieved. However, when the key soil hydraulic parameters are incorrect, the error is non-Gaussian, as the Kalman filter will produce a persistent bias in its predictions. In this paper, we propose a method coupling optimal parameters and extended Kalman filter data assimilation (OP-EKF) by combining optimal parameter estimation, the extended Kalman filter (EKF) assimilation method, a particle swarm optimization (PSO) algorithm, and Richards’ equation. We examine the accuracy of estimating root zone soil moisture through the optimal parameters and extended Kalman filter data assimilation method by using observed in situ data at the Meiling experimental station, China. Results indicate that merely using EKF for assimilating surface soil moisture content to obtain soil moisture content in the root zone will produce a persistent bias between simulated and observed values. Using the OP-EKF assimilation method, estimates were clearly improved. If the soil profile is heterogeneous, soil moisture retrieval is accurate in the 0-50 cm soil profile and is inaccurate at 100 cm depth. Results indicate that the method is useful for retrieving root zone soil moisture over large areas and long timescales even when available soil moisture data are limited to the surface layer, and soil moisture content are uncertain and soil hydraulic parameters are incorrect.     

A land surface soil moisture data assimilation framework in consideration of the model subgrid-scale heterogeneity and soil water thawing and freezing

The Ensemble Kalman Filter (EnKF) is well known and widely used in land data assimilation for its high precision and simple operation. The land surface models used as the forecast operator in a land data assimilation system are usually designed to consider the model subgrid-heterogeneity and soil water thawing and freezing. To neglect their effects could lead to some errors in soil moisture assimilation. The dual EnKF method is employed in soil moisture data assimilation to build a soil moisture data as- similation framework based on the NCAR Community Land Model version 2.0 (CLM 2.0) in considera- tion of the effects of the model subgrid-heterogeneity and soil water thawing and freezing: Liquid volumetric soil moisture content in a given fraction is assimilated through the state filter process, while solid volumetric soil moisture content in the same fraction and solid/liquid volumetric soil moisture in the other fractions are optimized by the parameter filter. Preliminary experiments show that this dual EnKF-based assimilation framework can assimilate soil moisture more effectively and precisely than the usual EnKF-based assimilation framework without considering the model subgrid-scale heteroge- neity and soil water thawing and freezing. With the improvement of soil moisture simulation, the soil temperature-simulated precision can be also improved to some extent.     

The subsurface–land surface–atmosphere connection under convective conditions

The dynamics of the free groundwater table influence land surface soil moisture and energy balance components, and are therefore also linked to atmospheric processes. In this study, the sensitivity of the atmosphere to groundwater table dynamics induced heterogeneity in land surface processes is examined under convective conditions. A fully coupled subsurface–land surface–atmosphere model is applied over a 150 km × 150 km study area located in Western Germany and ensemble simulations are performed over two convective precipitation events considering two separate model configurations based on groundwater table dynamics. Ensembles are generated by varying the model atmospheric initial conditions following the prescribed ensemble generation method by the German Weather Service in order to account for the intrinsic, internal atmospheric variability. The results demonstrate that especially under strong convective conditions, groundwater table dynamics affect atmospheric boundary layer height, convective available potential energy, and precipitation via the coupling with land surface soil moisture and energy fluxes. Thus, this study suggests that systematic uncertainties may be introduced to atmospheric simulations if groundwater table dynamics are neglected in the model.     

Root-zone soil moisture estimation from assimilation of downscaled Soil Moisture and Ocean Salinity data

The crucial role of root-zone soil moisture is widely recognized in land–atmosphere interaction, with direct practical use in hydrology, agriculture and meteorology. But it is difficult to estimate the root-zone soil moisture accurately because of its space-time variability and its nonlinear relationship with surface soil moisture. Typically, direct satellite observations at the surface are extended to estimate the root-zone soil moisture through data assimilation. But the results suffer from low spatial resolution of the satellite observation. While advances have been made recently to downscale the satellite soil moisture from Soil Moisture and Ocean Salinity (SMOS) mission using methods such as the Disaggregation based on Physical And Theoretical scale Change (DisPATCh), the assimilation of such data into high spatial resolution land surface models has not been examined to estimate the root-zone soil moisture. Consequently, this study assimilates the 1-km DisPATCh surface soil moisture into the Joint UK Land Environment Simulator (JULES) to better estimate the root-zone soil moisture. The assimilation is demonstrated using the advanced Evolutionary Data Assimilation (EDA) procedure for the Yanco area in south eastern Australia. When evaluated using in-situ OzNet soil moisture, the open loop was found to be 95% as accurate as the updated output, with the updated estimate improving the DisPATCh data by 14%, all based on the root mean square error (RMSE). Evaluation of the root-zone soil moisture with in-situ OzNet data found the updated output to improve the open loop estimate by 34% for the 0–30 cm soil depth, 59% for the 30–60 cm soil depth, and 63% for the 60–90 cm soil depth, based on RMSE. The increased performance of the updated output over the open loop estimate is associated with (i) consistent estimation accuracy across the three soil depths for the updated output, and (ii) the deterioration of the open loop output for deeper soil depths. Thus, the findings point to a combined positive impact from the DisPATCh data and the EDA procedure, which together provide an improved soil moisture with consistent accuracy both at the surface and at the root-zone.     

Utilization of satellite data in land surface hydrology: sensitivity and assimilation

This paper investigates the sensitivity of potential evapotranspiration to input meteorological variables, i.e. surface air temperature and surface vapor pressure. The sensitivity studies have been carried out for a wide range of land surface variables such as wind speed, leaf area index and surface temperatures. Errors in the surface air temperature and surface vapor pressure result in errors of different signs in the computed potential evapotranspiration. This result has implications for use of estimated values from satellite data or analysis of surface air temperature and surface vapor pressure in large‐scale hydrological modeling. The comparison of cumulative potential evapotranspiration estimates using ground observations and satellite observations over Manhattan, Kansas for a period of several months shows a variable difference between the two estimates. The use of satellite estimates of surface skin temperature in hydrological modeling to update the soil moisture using a physical adjustment concept is studied in detail, including the extent of changes in soil moisture resulting from the assimilation of surface skin temperature. The soil moisture of the 1 cm surface layer was adjusted by 0·9 mm over a 10‐day period as a result of a 3 K difference between the predicted and the observed surface temperature. This is a considerable amount given the fact that the top layer can hold only 5 mm of moisture. Copyright © 2001 John Wiley & Sons, Ltd.     

A dual-pass data assimilation scheme for estimating surface fluxes with FY3A-VIRR land surface temperature

In this work, a dual-pass data assimilation scheme is developed to improve predictions of surface flux. Pass 1 of the dual-pass data assimilation scheme optimizes the model vegetation parameters at the weekly temporal scale, and Pass 2 optimizes the soil moisture at the daily temporal scale. Based on ensemble Kalman filter(EnKF), the land surface temperature(LST) data derived from the new generation of Chinese meteorology satellite(FY3A-VIRR) are assimilated into common land model(CoLM) for the first time. Six sites, Daman, Guantao, Arou, BJ, Miyun and Jiyuan, are selected for the data assimilation experiments and include different climatological conditions. The results are compared with those from a dataset generated by a multi-scale surface flux observation system that includes an automatic weather station(AWS), eddy covariance(EC) and large aperture scintillometer(LAS). The results indicate that the dual-pass data assimilation scheme is able to reduce model uncertainties and improve predictions of surface flux with the assimilation of FY3A-VIRR LST data.