Groundwater contaminant source identification based on iterative local update ensemble smoother

Identification of the location and intensity of groundwater pollution source contributes to the effect of pollution remediation, and is called groundwater contaminant source identifcation. This is a kind of typical groundwater inverse problem, and the solution is usually ill-posed. Especially considering the spatial variability of hydraulic conductivity field, the identification process is more challenging. In this paper, the solution framework of groundwater contaminant source identification is composed with groundwater pollutant transport model (MT3DMS) and a data assimilation method (Iterative local update ensemble smoother, ILUES). In addition, Karhunen-Loève expansion technique is adopted as a PCA method to realize dimension reduction. In practical problems, the geostatistical method is usually used to characterize the hydraulic conductivity feld, and only the contaminant source information is inversely calculated in the identifcation process. In this study, the identification of contaminant source information under Kriging K-field is compared with simultaneous identification of source information and K-field. The results indicate that it is necessary to carry out simultaneous identification under heterogeneous site, and ILUES has good performance in solving high-dimensional parameter inversion problems.     

基于Whittaker滤波的陕西省植被物候特征

运用Whittaker滤波重构MODIS NDVI时序数列,利用地理探测器对比滤波前后影像信噪比,采用动态阈值法获取2000-2012年陕西省植被的3个关键物候参数(返青期、枯黄期和生长周期),在此基础上分析了该区植被物候参数空间分布特征。结果表明:(1)Whittaker滤波能够平滑原始NDVI曲线,有效减少原始影像的噪声,提高影像辨识度,并且参数设置简单;(2)陕西省植被物候地区分异明显,不同气候区划类植被物候表现出中温带半干旱区-暖温带半干旱区-暖温带半湿润区-北亚热带湿润区的递变规律:返青期逐步提前,枯黄期逐步推迟;(3)植被物候受高程和纬度影响,并且纬度影响更显著。海拔每升高200 m,返青期推迟1.3 d,枯黄期提前0.6 d;纬度每升高0.5°,返青期推迟3.6 d,枯黄期提前1.2 d。     

The representer method for two-phase flow in porous media

The representer method is applied to a one-dimensional two-phase flow model in porous media; capillary pressure and gravity are neglected. The Euler–Lagrange equations must be linearized, and one such linearization is presented here. The representer method is applied to the linear system iteratively until convergence, though a rigorous proof of convergence is out of reach. The linearization chosen is easy to calculate but does not converge for certain weights; however, a simple damping restores convergence at the cost of extra iterations. Numerical experiments are performed that illustrate the method, and quick comparison to the ensemble Kalman smoother is made. This research was supported by NSF grant EIA-0121523.     

A time series for monitoring vegetation activity and phenology at 10-daily time steps covering large parts of South America

Abstract

It is widely accepted that natural resources should only be sustainably exploited and utilized to effectively preserve our planet for future generations. To better manage the natural resources, and to better understand the closely linked Earth systems, the concept of Digital Earth has been strongly promoted since US Vice President Al Gore's speech in 1998. One core element of Digital Earth is the use and integration of remote sensing data. Only satellite imagery can cover the entire globe repeatedly at a sufficient high-spatial resolution to map changes in land cover and land use, but also to detect more subtle changes related for instance to climate change. To uncover global change effects on vegetation activity and phenology, it is important to establish high quality time series characterizing the past situation against which the current state can be compared. With the present study we describe a time series of vegetation activity at 10-daily time steps between 1998 and 2008 covering large parts of South America at 1 km spatial resolution. Particular emphasis was put on noise removal. Only carefully filtered time series of vegetation indices can be used as a benchmark and for studying vegetation dynamics at a continental scale. Without temporal smoothing, subtle spatio-temporal patterns in vegetation composition, density and phenology would be hidden by atmospheric noise and undetected clouds. Such noise is immanent in data that have undergone solely a maximum value compositing. Within the present study, the Whittaker smoother (WS) was applied to a SPOT VGT time series. The WS balances the fidelity to the observations with the roughness of the smoothed curve. The algorithm is extremely fast, gives continuous control over smoothness with only one parameter, and interpolates automatically. The filtering efficiently removed the negatively biased noise present in the original data, while preserving the overall shape of the curves showing vegetation growth and development. Geostatistical variogram analysis revealed a significantly increased signal-to-noise ratio compared to the raw data. Analysis of the data also revealed spatially consistent key phenological markers. Extracted seasonality parameters followed a clear meridional trend. Compared to the unfiltered data, the filtered time series increased the separability of various land cover classes. It is thus expected that the data set holds great potential for environmental and vegetation related studies within the frame of Digital Earth.     

LOESS四维客观分析在中国近海的应用

利用全国海洋综合调查资料和世界大洋数据库(WOD)资料,运用局部加权回归法(LOESS)四维客观分析方案对中国近海的温度场进行了客观分析,将散布在空间中的温度资料点插值到均一的网格中。在分析中充分考虑了中国海曲折复杂的海岸特点,将底深调整系统(TAR)和障碍调整系统(BAR)整合到四维插值系统中,使得资料点空间权重的计算更为合理。计算结果显示,LOESS客观分析方案充分考虑了空间因素(经度、纬度、深度)和时间因素,能客观反映出中国近海一年四季的温度时间-空间结构。通过与Kriging插值方法的对比,发现LOESS方法能更精确地把中国近海温度结构反映出来。     

集合卡尔曼平滑和集合卡尔曼滤波在污染源反演中的应用

此文目的是讨论污染源反演问题的统计方法.基于Bayes估计理论,该文将资料同化中的集合平滑、集合卡尔曼平滑和集合卡尔曼滤波应用在污染源反演问题中.在详细给出污染源反演的集合平滑、集合卡尔曼平滑和集合卡尔曼滤波的严格数学表达后,用一个简单的模型演示了集合卡尔曼平滑和集合卡尔曼滤波在污染源反演中的可行性,并且通过对比理想试验结果比较了集合卡尔曼平滑和集合卡尔曼滤波方法在反演污染源排放的效果,讨论了观测误差和污染源先验误差估计对反演结果的影响.试验结果表明在观测间隔小和观测误差小的情况下,集合卡尔曼滤波和集合卡尔曼平滑都可以有效地反演出随时间变化的污染源排放.当观测误差增大时,集合卡尔曼滤波和集合卡尔曼平滑的反演效果都有一定降低,但是反演误差的增加少于观测误差的增加,同时集合卡尔曼平滑(Ensemble Kalman smoother,简称EnKS)对观测误差比集合卡尔曼滤波(Ensemble Kalman filter,简称EnKF)更为敏感.当观测时间间隔较大时,EnKF不能对没有观测时的污染源排放进行估计,仅能对有观测时的污染源排放进行较好的反演.而EnKS可以利用观测对观测时刻前的污染源排放进行反演,因此其效果明显好于EnKF,并且在观测时间间隔较大的情况下依然可以较好地反演出污染源排放.试验结果还显示污染源排放的先验误差估计对反演的结果有较大影响.     

Data assimilation in groundwater modelling: ensemble Kalman filter versus ensemble smoothers

Groundwater modelling calls for an effective and robust data integrating method to fill the gap between the model and observation data. The ensemble Kalman filter (EnKF), a real‐time data assimilation method, has been increasingly applied in multiple disciplines such as petroleum engineering and hydrogeology. In this approach, a groundwater model is updated sequentially with measured data such as hydraulic head and concentration. As an alternative to the EnKF, the ensemble smoother (ES) has been proposed for updating groundwater models using all the data together, with much less computational cost. To further improve the performance of the ES, an iterative ES has been proposed for continuously updating the model by assimilating measurements together. In this work, we compare the performance of the EnKF, the ES, and the iterative ES using a synthetic example in groundwater modelling. Hydraulic head data modelled on the basis of the reference conductivity field are used to inversely estimate conductivities at unsampled locations. Results are evaluated in terms of the characterization of conductivity and groundwater flow predictions. It is concluded that (a) the iterative ES works better than the standard ES because of its continuous updating and (b) the iterative ES could achieve results comparable with those of the EnKF, with less computational cost. These findings show that the iterative ES should be paid much more attention for data assimilation in groundwater modelling.