scholarly journals An Evapotranspiration Assimilation Method Based on Ensemble Kalman Filter and À Trous Wavelet

2013 ◽  
Vol 2013 ◽  
pp. 1-13
Author(s):  
Shaohui Chen ◽  
Jianwei Qi ◽  
Xiaomin Sun ◽  
Xiangzheng Deng ◽  
Jing Tian

It is challenging to assimilate the evapotranspiration product (EP) retrieved from satellite data into land surface models (LSMs). In this paper, a perturbed ensemble Kalman filter (PEKF) and à trous wavelet transform (AWT) integrated method are proposed to implement the evapotranspiration assimilation. In this method, the AWT is used to decompose the EPs into multiple channels since it is very powerful in fusing high frequency spatial information of multisource data, and then the Kalman filter is performed in the AWT domain. The proposed method combines the advantages of the PEKF that is capable of accommodating model error and observation error, and the AWT can effectively perform multiresolution fusion. Assimilation experiment conducted with the Noah model and the EP retrieved from the MODIS data shows that the proposed method performs better than the traditional ensemble Kalman filter (EnKF) and PEKF methods. The analysis results fit well with the evapotranspiration observation at two field sites with different land surface conditions. These indicate that the proposed method is promising for assimilating regional scale satellite retrieved EP into LSMs.

2014 ◽  
Vol 7 (1) ◽  
pp. 283-302 ◽  
Author(s):  
B. Gaubert ◽  
A. Coman ◽  
G. Foret ◽  
F. Meleux ◽  
A. Ung ◽  
...  

Abstract. An ensemble Kalman filter (EnKF) has been coupled to the CHIMERE chemical transport model in order to assimilate ozone ground-based measurements on a regional scale. The number of ensembles is reduced to 20, which allows for future operational use of the system for air quality analysis and forecast. Observation sites of the European ozone monitoring network have been classified using criteria on ozone temporal variability, based on previous work by Flemming et al. (2005). This leads to the choice of specific subsets of suburban, rural and remote sites for data assimilation and for evaluation of the reference run and the assimilation system. For a 10-day experiment during an ozone pollution event over Western Europe, data assimilation allows for a significant improvement in ozone fields: the RMSE is reduced by about a third with respect to the reference run, and the hourly correlation coefficient is increased from 0.75 to 0.87. Several sensitivity tests focus on an a posteriori diagnostic estimation of errors associated with the background estimate and with the spatial representativeness of observations. A strong diurnal cycle of both these errors with an amplitude up to a factor of 2 is made evident. Therefore, the hourly ozone background error and the observation error variances are corrected online in separate assimilation experiments. These adjusted background and observational error variances provide a better uncertainty estimate, as verified by using statistics based on the reduced centered random variable. Over the studied 10-day period the overall EnKF performance over evaluation stations is found relatively unaffected by different formulations of observation and simulation errors, probably due to the large density of observation sites. From these sensitivity tests, an optimal configuration was chosen for an assimilation experiment extended over a three-month summer period. It shows a similarly good performance as the 10-day experiment.


2016 ◽  
Author(s):  
Hongjuan Zhang ◽  
Harrie-Jan Hendricks Franssen ◽  
Xujun Han ◽  
Jasper Vrugt ◽  
Harry Vereecken

Abstract. Land surface models (LSMs) contain a suite of different parameters and state variables to resolve the water and energy balance at the soil-atmosphere interface. Many of the parameters of these models cannot be measured directly in the field, and require calibration against flux and soil moisture data. In this paper, we use the Variable Infiltration Capacity Hydrologic Model (VIC) and the Community Land Model (CLM) to simulate temporal variations in soil moisture content at 5, 20 and 50 cm depth in the Rollesbroich experimental watershed in Germany. Four different data assimilation (DA) methods are used to jointly estimate the spatially distributed water content values, and hydraulic and/or thermal properties of the resolved soil domain. This includes the Ensemble Kalman Filter (EnKF) using state augmentation or dual estimation, the Residual Resampling Particle Filter (RRPF) and Markov chain Monte Carlo Particle Filter (MCMCPF). These four DA methods are tuned and calibrated for a five month data period, and subsequently evaluated for another five month period. Our results show that all the different DA methods improve the fit of the VIC and CLM model to the observed water content data, particularly if the maximum baseflow velocity (VIC), soil hydraulic (VIC) properties and/or soil texture (CLM) are jointly estimated along with the model states. In the evaluation period, the augmentation and dual estimation method performed slightly better than RRPF and MCMCPF. The differences in simulated soil moisture values between the CLM and VIC model were larger than variations among the data assimilation algorithms. The best performance for the Rollesbroich site was observed for the CLM model. The strong underestimation of the soil moisture values of the third VIC-layer are likely explained by an inadequate parameterization of groundwater drainage.


2016 ◽  
Vol 35 (1) ◽  
pp. 30-37 ◽  
Author(s):  
Danian Liu ◽  
Ping Shi ◽  
Yeqiang Shu ◽  
Jinglong Yao ◽  
Dongxiao Wang ◽  
...  

2011 ◽  
Vol 139 (6) ◽  
pp. 2008-2024 ◽  
Author(s):  
Brian C. Ancell ◽  
Clifford F. Mass ◽  
Gregory J. Hakim

Abstract Previous research suggests that an ensemble Kalman filter (EnKF) data assimilation and modeling system can produce accurate atmospheric analyses and forecasts at 30–50-km grid spacing. This study examines the ability of a mesoscale EnKF system using multiscale (36/12 km) Weather Research and Forecasting (WRF) model simulations to produce high-resolution, accurate, regional surface analyses, and 6-h forecasts. This study takes place over the complex terrain of the Pacific Northwest, where the small-scale features of the near-surface flow field make the region particularly attractive for testing an EnKF and its flow-dependent background error covariances. A variety of EnKF experiments are performed over a 5-week period to test the impact of decreasing the grid spacing from 36 to 12 km and to evaluate new approaches for dealing with representativeness error, lack of surface background variance, and low-level bias. All verification in this study is performed with independent, unassimilated observations. Significant surface analysis and 6-h forecast improvements are found when EnKF grid spacing is reduced from 36 to 12 km. Forecast improvements appear to be a consequence of increased resolution during model integration, whereas analysis improvements also benefit from high-resolution ensemble covariances during data assimilation. On the 12-km domain, additional analysis improvements are found by reducing observation error variance in order to address representativeness error. Removing model surface biases prior to assimilation significantly enhances the analysis. Inflating surface wind and temperature background error variance has large impacts on analyses, but only produces small improvements in analysis RMS errors. Both surface and upper-air 6-h forecasts are nearly unchanged in the 12-km experiments. Last, 12-km WRF EnKF surface analyses and 6-h forecasts are shown to generally outperform those of the Global Forecast System (GFS), North American Model (NAM), and the Rapid Update Cycle (RUC) by about 10%–30%, although these improvements do not extend above the surface. Based on these results, future improvements in multiscale EnKF are suggested.


2011 ◽  
Vol 15 (8) ◽  
pp. 2437-2457 ◽  
Author(s):  
S. Nie ◽  
J. Zhu ◽  
Y. Luo

Abstract. The performance of the ensemble Kalman filter (EnKF) in soil moisture assimilation applications is investigated in the context of simultaneous state-parameter estimation in the presence of uncertainties from model parameters, soil moisture initial condition and atmospheric forcing. A physically based land surface model is used for this purpose. Using a series of identical twin experiments in two kinds of initial parameter distribution (IPD) scenarios, the narrow IPD (NIPD) scenario and the wide IPD (WIPD) scenario, model-generated near surface soil moisture observations are assimilated to estimate soil moisture state and three hydraulic parameters (the saturated hydraulic conductivity, the saturated soil moisture suction and a soil texture empirical parameter) in the model. The estimation of single imperfect parameter is successful with the ensemble mean value of all three estimated parameters converging to their true values respectively in both NIPD and WIPD scenarios. Increasing the number of imperfect parameters leads to a decline in the estimation performance. A wide initial distribution of estimated parameters can produce improved simultaneous multi-parameter estimation performances compared to that of the NIPD scenario. However, when the number of estimated parameters increased to three, not all parameters were estimated successfully for both NIPD and WIPD scenarios. By introducing constraints between estimated hydraulic parameters, the performance of the constrained three-parameter estimation was successful, even if temporally sparse observations were available for assimilation. The constrained estimation method can reduce RMSE much more in soil moisture forecasting compared to the non-constrained estimation method and traditional non-parameter-estimation assimilation method. The benefit of this method in estimating all imperfect parameters simultaneously can be fully demonstrated when the corresponding non-constrained estimation method displays a relatively poor parameter estimation performance. Because all these constraints between parameters were obtained in a statistical sense, this constrained state-parameter estimation scheme is likely suitable for other land surface models even with more imperfect parameters estimated in soil moisture assimilation applications.


2019 ◽  
Vol 11 (7) ◽  
pp. 753 ◽  
Author(s):  
Guodong Zhang ◽  
Hongmin Zhou ◽  
Changjing Wang ◽  
Huazhu Xue ◽  
Jindi Wang ◽  
...  

Continuous, long-term sequence, land surface albedo data have crucial significance for climate simulations and land surface process research. Sensors such as the Moderate-Resolution Imaging Spectroradiometer (MODIS) and Visible Infrared Imaging Radiometer (VIIRS) provide global albedo product data sets with a spatial resolution of 500 m over long time periods. There is demand for new high-resolution albedo data for regional applications. High-resolution observations are often unavailable due to cloud contamination, which makes it difficult to obtain time series albedo estimations. This paper proposes an “amalgamation albedo“ approach to generate daily land surface shortwave albedo with 30 m spatial resolution using Landsat data and the MODIS Bidirectional Reflectance Distribution Functions (BRDF)/Albedo product MCD43A3 (V006). Historical MODIS land surface albedo products were averaged to obtain an albedo estimation background, which was used to construct the albedo dynamic model . The Thematic Mapper (TM) albedo derived via direct estimation approach was then introduced to generate high spatial-temporal resolution albedo data based on the Ensemble Kalman Filter algorithm (EnKF). Estimation results were compared to field observations for cropland, deciduous broadleaf forest, evergreen needleleaf forest, grassland, and evergreen broadleaf forest domains. The results indicated that for all land cover types, the estimated albedos coincided with ground measurements at a root mean squared error (RMSE) of 0.0085–0.0152. The proposed algorithm was then applied to regional time series albedo estimation; the results indicated that it captured spatial and temporal variation patterns for each site. Taken together, our results suggest that the amalgamation albedo approach is a feasible solution to generate albedo data sets with high spatio-temporal resolution.


2011 ◽  
Vol 139 (11) ◽  
pp. 3389-3404 ◽  
Author(s):  
Thomas Milewski ◽  
Michel S. Bourqui

Abstract A new stratospheric chemical–dynamical data assimilation system was developed, based upon an ensemble Kalman filter coupled with a Chemistry–Climate Model [i.e., the intermediate-complexity general circulation model Fast Stratospheric Ozone Chemistry (IGCM-FASTOC)], with the aim to explore the potential of chemical–dynamical coupling in stratospheric data assimilation. The system is introduced here in a context of a perfect-model, Observing System Simulation Experiment. The system is found to be sensitive to localization parameters, and in the case of temperature (ozone), assimilation yields its best performance with horizontal and vertical decorrelation lengths of 14 000 km (5600 km) and 70 km (14 km). With these localization parameters, the observation space background-error covariance matrix is underinflated by only 5.9% (overinflated by 2.1%) and the observation-error covariance matrix by only 1.6% (0.5%), which makes artificial inflation unnecessary. Using optimal localization parameters, the skills of the system in constraining the ensemble-average analysis error with respect to the true state is tested when assimilating synthetic Michelson Interferometer for Passive Atmospheric Sounding (MIPAS) retrievals of temperature alone and ozone alone. It is found that in most cases background-error covariances produced from ensemble statistics are able to usefully propagate information from the observed variable to other ones. Chemical–dynamical covariances, and in particular ozone–wind covariances, are essential in constraining the dynamical fields when assimilating ozone only, as the radiation in the stratosphere is too slow to transfer ozone analysis increments to the temperature field over the 24-h forecast window. Conversely, when assimilating temperature, the chemical–dynamical covariances are also found to help constrain the ozone field, though to a much lower extent. The uncertainty in forecast/analysis, as defined by the variability in the ensemble, is large compared to the analysis error, which likely indicates some amount of noise in the covariance terms, while also reducing the risk of filter divergence.


Sign in / Sign up

Export Citation Format

Share Document