赣江流域多种数据同化方案的径流模拟比较

为了探讨水文模型在不同水文数据同化方案下的径流模拟差异,本文采用集合卡尔曼滤波算法,以遥感蒸散发产品、实测径流为观测数据,构建了基于新安江模型的数据同化框架。基于此框架设计了4种不同同化方案(DA-ET、DAET(K)、DA-ET-Q、DA-ET-Q(K))以及1种对照方案OL,以赣江流域开展实例研究,评估了水文数据同化中遥感蒸散发产品的时间分辨率、模型蒸散发相关参数时变与否以及多源数据同化对径流模拟的影响。结果表明:在DA-ET方案下,同化两种不同时间分辨率的蒸散发产品均能提高模型整体的径流模拟精度,且时间分辨率更高的产品的同化效果更好;在DA-ET方案的基础上,考虑加入实测径流进行同化能够提升模型径流模拟精度,且DA-ET(K)与DA-ET-Q(K)方案所得径流相对误差的减幅均超过了20%,说明在蒸散发同化过程中同时考虑蒸散发参数动态变化的结果更优;相较于OL方案,4种同化方案均能不同程度地提高模型对径流高水部分的模拟能力,但DA-ET-Q(K)方案表现最差,而其余方案差异并不显著。本研究有助于进一步了解不同数据同化方案在径流模拟中的差异,从而为水资源高效利用与科学管理提供科学依据...     

基于热层电离层耦合数据同化的热层参量估计

本文采用高效集合卡尔曼滤波（EnKF）算法和背景场热层电离层理论模式NCAR-TIEGCM,开发了热层电离层数据同化系统.基于全球空地基GNSS电离层斜TEC观测、CHAMP和TIMED/GUVI热层参量观测构型设计了系列观测系统模拟实验,对热层参量进行估计.实验结果表明,（1）通过集合卡尔曼滤波算法同化电离层TEC观测能够较好地优化热层参量.（2）中性质量密度优化效果在整个同化阶段均有提升,提升百分比能达到40%.（3）积分氧氮比在同化阶段也能得到较好的优化,但在电子密度水平梯度变化剧烈区域效果较差.最后本文对中性质量密度进行了预报评估,结果表明,由于中性成分优化,在地磁平静条件下其预报时间尺度可长达24 h.

     

湖泊富营养化响应与流域优化调控决策的模型研究进展

湖泊富营养化是全球水环境领域面临的长期挑战,富营养化响应与流域优化决策模型是制定经济和高效调控方案的关键.然而已有的模型研究综述主要集中于模型开发、案例应用、敏感性分析、不确定性分析等单一方面,而缺少针对非线性响应、生态系统长期演变等最新湖泊治理挑战的研究总结.本文对数据驱动的统计模型、因果驱动的机理模型和决策导向的优化模型进行了综述.其中,统计模型包含经典统计、贝叶斯统计和机器学习模型,常用于建立响应关系、时间序列特征分析以及预报预警;机理模型包含流域的水文与污染物输移模拟以及湖泊的水文、水动力、水质、水生态等过程的模拟,用于不同时空尺度的变化过程模拟,其中复杂机理模型的敏感性分析、参数校验、模型不确定性等需要较高的计算成本;优化模型结合机理模型形成“模拟优化”体系,在不确定性条件下衍生出随机、区间优化等多种方法,通过并行计算、简化与替代模型可一定程度上解决计算时间成本的瓶颈.本文识别了湖泊治理面临的挑战,包括:①如何定量表征外源输入的非线性叠加和湖泊氮、磷、藻变化的非均匀性?②如何提高优化调控决策和水质目标的关联与精准性?③如何揭示湖泊生态系统的长期变化轨迹与驱动因素?最后,本文针对这些挑战提出研究展望,主要包括:①基于多源数据融合与机器学习算法以提升湖泊的短期水质预测精度;②以生物量为基础的机理模型与行为驱动的个体模型的升尺度或降尺度耦合以表达多种尺度的物质交互过程;③机器学习算法与机理模型的直接耦合或数据同化以降低模拟误差;④时空尺度各异的多介质模拟模型融合以实现精准和动态的优化调控.     

浅水湖泊湖沼学与太湖富营养化控制研究

自2007年无锡暴发饮用水危机事件以来,太湖经历了前所未有的高强度、大规模治理,各种治理措施累计投资已经超过千亿元.监测显示,在治理初期太湖的氮、磷浓度下降明显,水质有所好转,但最近几年关键水质指标总磷与浮游植物叶绿素α浓度出现了波动,蓝藻水华有所反弹.研究表明,太湖的外源负荷并没有减少,这与城镇用水量增加、污水排放标准偏低、面源污染削减不足有很大的关系;同时,内源负荷也因为蓝藻水华的持续而加重,浅水湖泊水深浅、扰动强的特点强化了磷的循环利用效率,加剧了内源负荷对湖泊富营养化和蓝藻水华的影响.气候变暖叠加营养盐富集的复合效应、流域风速下降以及暴雨事件频次和强度增加等气象水文条件变化,都促进了太湖蓝藻水华的暴发;蓝藻水华的时空分布特征则受湖泊水动力的决定性影响.太湖治理的曲折过程,凸显了大型浅水湖泊湖沼学研究的不断深入与发展,未来需要继续加强多学科交叉研究,特别是基于湖泊-流域系统的气象水文、生物地球化学和生物生态学的学科交叉.对于太湖生态环境的综合治理和管理,既要注重湖泊与流域相结合,更需要重视自然科学和人文科学的有机融合,才能真正达到控制太湖富营养化、维护流域水环境安全与社会经济可持...     

基于微波亮温及集合Kalman滤波的土壤湿度同化方案

基于集合Kalman滤波及SCE-UA(shuffled complex evolution)算法发展了能够直接同化微波亮温的土壤湿度同化方案. 该方案以陆面过程模式CLM 3.0中的土壤水模型作为预报算子, 以辐射传输模型作为观测算子. 整个同化过程分为参数优化和土壤湿度同化两个阶段, 利用SCE-UA算法优化辐射传输模型中难以确定的植被光学厚度参数和地表粗糙度参数, 并利用优化参数作为观测算子的模型参数进行同化. 通过人工理想试验表明该同化方案可以明显改善表层土壤湿度的模拟精度, 并且对深层土壤湿度的模拟也有一定程度的改善; 利用AMSR-E亮温(10.65 GHz垂直极化)所进行的实际同化试验表明顶层(0~10 cm)土壤湿度同化结果与观测的均方根误差(RMSE)由模拟的0.05052减小到0.03355, 相对减小了33.6%, 而较深层(10~50 cm)平均减小了20.9%. 这些同化试验显示该同化方案的合理性.     

富营养化湖泊中硫酸盐对蓝藻衰亡产甲烷过程的影响

随着外源性硫酸盐（SO₄^2-）的持续性输入,富营养化湖泊水体的SO₄^2-浓度持续升高.野外长期监测结果表明,近几十年太湖水体的SO₄^2-浓度逐渐升高,达到了96 mg/L的水平.此外,富营养化湖泊中蓝藻水华衰亡会产生并释放大量的甲烷（CH₄）,湖泊水体的SO₄^2-浓度升高是否会对沉积物产CH₄过程造成影响仍缺乏相关研究.本实验构建了蓝藻水体沉积物微宇宙系统,通过添加30、60、90、120和150 mg/L五组浓度的硫酸盐,探究不同SO₄^2-浓度下蓝藻衰亡过程中水体的SO₄^2-、还原性硫化物（∑S^2-）和CH₄的变化规律.结果表明,蓝藻聚积衰亡的第6~9天硫酸盐还原作用最为强烈,此时水体中的SO₄^2-浓度快速下降到最低值,依次为7.65、8.87、21.21、41.14和56.54 mg/L.伴随着硫酸盐还原过程的进行,水柱中∑S^2-的浓度不断上升,并达到最高值,依次为4.77、6.98、7.49、7.49和7.43 mg/L.蓝藻聚积衰亡的第10~21天水体中的SO₄^2-浓度维持在较低水平,∑S^2-浓度逐渐下降,并趋近于0.培养开始时,CH₄增长缓慢,SO₄^2-浓度下降之后,CH₄浓度逐渐上升,并在第6~9天迅速上升,培养结束时,CH₄的最终浓度随着水体初始SO₄^2-浓度的增加而降低,依次为546.39、207.24、79.61、37.25和5.56 μmol/L,CH₄的浓度与初始水体SO₄^2-浓度呈指数型负相关关系.因此,对于精准评估富营养化湖泊的产甲烷过程,需要考虑不断上升的SO₄^2-浓度所带来的影响.     

耦合敏感参数实时识别的新型数据同化算法研究——以湖泊藻类模拟为例

数据同化是提升复杂机理过程模型精度的关键技术之一,而湖泊藻类模型的敏感参数具有随时间动态变化的特征,导致数据同化过程中无法精准更新某一时段的敏感参数,影响数据同化的模型精度提升效果.针对上述问题,本研究耦合了参数敏感性分析与集合卡尔曼滤波,研发了一种能够实时识别模型敏感参数的新型数据同化算法;为验证研发算法的效率,依托巢湖的高频水质自动监测数据,测试算法对藻类动态模型的精度提升效果.测试结果表明：研发算法能够精准跟踪模型敏感参数的动态变化,并根据监测数据实时更新模型敏感参数,实现了水质高频自动监测数据与藻类动态模型的深度融合,藻类生物量模拟精度提升了55%,即纳什系数(NSE)从0.49提升到0.76,模拟精度提升效果也显著优于传统数据同化算法(NSE=0.63).研发算法可应用于其它水生态环境模型的数据同化,为水生态环境相关要素的精准模拟预测提供关键技术支撑.     

基于Gauss-Markov卡尔曼滤波的电离层数值同化现报预报系统的构建—以中国及周边地区为例的观测系统模拟试验

本文给出了一个基于Gauss-Markov卡尔曼滤波的电离层数据同化系统的初步构建和试验结果.我们选择中国及周边地区部分涉及电离层观测的台站(包括子午工程台站、中国地壳形变网和部分IGS台站)作为观测系统进行模拟试验,背景场利用IRI模式,观测值则由NeQuick模式计算得到.我们的同化结果表明,采用Kalman滤波算法,把部分斜TEC同化到背景模式当中,能够获得较好的同化结果,说明我们设计的算法可行、所选择的各种参数比较合理,采用Gauss-Markov假设进行短期预报也取得了较合理的结果.本项研究经过进一步的改进和完善,可以用来对中国地区的电离层进行现报和短期预报,一方面满足相关空间工程应用,另一方面可以提升现有观测系统的科学意义.     

太湖叶绿素a同化系统敏感性分析

太湖叶绿素a同化系统对于不同参数的敏感性将直接影响到该系统能否精确的估算太湖叶绿素a的浓度分布.利用2009年4月21日环境一号卫星(HJ-1B CCD2)影像数据反演太湖叶绿素a浓度场信息.以此作为背景场信息,结合基于集合均方根滤波的太湖叶绿素a同化系统,分析和评价了样本数目、同化时长、背景场误差、观测误差和模型误差对于同化系统性能的影响.结果表明:从计算成本、系统运行时间和同化效果等方面分析,当集合样本数目达到30~40左右时同化系统取得了较好的结果;同化系统对于背景场误差的估计变化不是很敏感,即初始场的估计是否准确对于同化系统的性能影响不是很大;同化系统对于模型误差和观测误差的变化较为敏感,不同的测试点位由于水体动力学性质不一,其敏感性的表现形式有所差异;利用数据同化方法可以有效地估算太湖叶绿素a浓度.     

基于集合四维变分方法PODEn4DVar的GRACE陆地水储量同化：方法与验证

地球重力场季节和年际变化主要来源于地球表层大气、海洋和陆地各系统间水的质量交换,由GRACE重力卫星探测地球重力场变化所反映陆地水储量的改变主要来自降雨、土壤蒸发蒸腾、河流输运以及向地下深层的渗透等过程．本研究利用陆面过程模式CLM3．5以及基于本征正交分解的集合四维变分同化方法PODEn4DVar,构建能够同化GRACE卫星重力场的陆面水文同化系统LDAS．G,实现对地球重力场所反映的大尺度陆地水储量变化在时间及垂直方向各分量的分解,并对垂直方向的水文变量进行同化,从而更好地估计陆面水循环要素变化并实现其监测．利用LDAS．G同化系统进行理想试验以及针对中国区域所进行的同化试验表明该同化系统能够改善对陆面水文要素变化的模拟,对大尺度陆面水文循环监测研究具有重要意义．     

Hydrological data assimilation with the Ensemble Square-Root-Filter: Use of streamflow observations to update model states for real-time flash flood forecasting

The objective of the study is to evaluate the potential of a data assimilation system for real-time flash flood forecasting over small watersheds by updating model states. To this end, the Ensemble Square-Root-Filter (EnSRF) based on the Ensemble Kalman Filter (EnKF) technique was coupled to a widely used conceptual rainfall-runoff model called HyMOD. Two small watersheds susceptible to flash flooding from America and China were selected in this study. The modeling and observational errors were considered in the framework of data assimilation, followed by an ensemble size sensitivity experiment. Once the appropriate model error and ensemble size was determined, a simulation study focused on the performance of a data assimilation system, based on the correlation between streamflow observation and model states, was conducted. The EnSRF method was implemented within HyMOD and results for flash flood forecasting were analyzed, where the calibrated streamflow simulation without state updating was treated as the benchmark or nature run. Results for twenty-four flash-flood events in total from the two watersheds indicated that the data assimilation approach effectively improved the predictions of peak flows and the hydrographs in general. This study demonstrated the benefit and efficiency of implementing data assimilation into a hydrological model to improve flash flood forecasting over small, instrumented basins with potential application to real-time alert systems.     

Data assimilation for distributed hydrological catchment modeling via ensemble Kalman filter

Catchment scale hydrological models are critical decision support tools for water resources management and environment remediation. However, the reliability of hydrological models is inevitably affected by limited measurements and imperfect models. Data assimilation techniques combine complementary information from measurements and models to enhance the model reliability and reduce predictive uncertainties. As a sequential data assimilation technique, the ensemble Kalman filter (EnKF) has been extensively studied in the earth sciences for assimilating in-situ measurements and remote sensing data. Although the EnKF has been demonstrated in land surface data assimilations, there are no systematic studies to investigate its performance in distributed modeling with high dimensional states and parameters. In this paper, we present an assessment on the EnKF with state augmentation for combined state-parameter estimation on the basis of a physical-based hydrological model, Soil and Water Assessment Tool (SWAT). Through synthetic simulation experiments, the capability of the EnKF is demonstrated by assimilating the runoff and other measurements, and its sensitivities are analyzed with respect to the error specification, the initial realization and the ensemble size. It is found that the EnKF provides an efficient approach for obtaining a set of acceptable model parameters and satisfactory runoff, soil water content and evapotranspiration estimations. The EnKF performance could be improved after augmenting with other complementary data, such as soil water content and evapotranspiration from remote sensing retrieval. Sensitivity studies demonstrate the importance of consistent error specification and the potential with small ensemble size in the data assimilation system.     

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.     

Review of the Kalman-type hydrological data assimilation

ABSTRACT

There is great potential in Data Assimilation (DA) for the purposes of uncertainty identification, reduction and real-time correction of hydrological models. This paper reviews the latest developments in Kalman filters (KFs), particularly the Extended KF (EKF) and the Ensemble KF (EnKF) in hydrological DA. The hydrological DA targets, methodologies and their applicability are examined. The recent applications of the EKF and EnKF in hydrological DA are summarized and assessed critically. Furthermore, this review highlights the existing challenges in the implementation of the EKF and EnKF, especially error determination and joint parameter estimation. A detailed review of these issues would benefit not only the Kalman-type DA but also provide an important reference to other hydrological DA types.

Editor D. Koutsoyiannis; Associate editor F. Pappenberger     

湖泊水环境模糊数学评价模型及其在富营养化排序中的应用

本文根据模糊水文学关于水体“清洁”与“污染”之间无明确界面,在识别过程中具有模糊性的论点,提出湖泊水环境评价的模糊数学模型。应用此模型对我国12个湖泊富营养化污染状况进行排序,结果与实际相符,文中提出的评价模型,原则上也适用于其他环境领域的评价工作。     

On ensemble representation of the observation-error covariance in the Ensemble Kalman Filter

Evensen (2003) presents a modification of the Ensemble Kalman Filter (EnKF), in which the observation-error and background-error covariance matrices are both represented by ensembles, in contrast to the usual practice, where only the background error is so represented. It is shown that this modification can cause the ensemble to collapse to a single member, in the common situation where the number of observations is more than twice the number of ensemble members, and to be rank-deficient when the number of observations is greater than or equal to the ensemble size. It is also shown that some further modifications to the scheme, presented by Evensen as offering numerical efficiencies, can prevent this collapse. However, these latter modifications are shown in some simple numerical examples to require tuning to produce acceptable results, which are nevertheless inferior to those of the standard EnKF.Acknowledgements The author acknowledges useful discussions with Peter Steinle, and other participants at the EnKF workshop held in BMRC in November, 2003.     

Flow relevant covariance localization during dynamic data assimilation using EnKF

Multiphase dynamic data integration into high resolution subsurface models is an integral aspect of reservoir and groundwater management strategies and uncertainty assessment. Over the past two decades, advances in computing and the development and implementation of robust algorithms for automatic history matching have considerably reduced the time and effort associated with subsurface characterization and reduced the subjectivity associated with manual model calibration. However, reliable and accurate subsurface characterization continues to be challenging due to the large number of model unknowns to be estimated using a relatively smaller set of measurements. For ensemble-based methods in particular, the difficulties are compounded by the need for a large number of model replicates to estimate sample-based statistical measures, specifically the covariances and cross-covariances that directly impact the spread of information from the measurement locations to the model parameters. Statistical noise resulting from modest ensemble sizes can overwhelm and degrade the model updates leading to geologically inconsistent subsurface models. In this work we propose to address the difficulties in the implementation of the ensemble Kalman filter (EnKF) for operational data integration problems. The methods described here use streamline-derived information to identify regions within the reservoir that will have a maximum impact on the dynamic response. This is achieved through spatial localization of the sample-based cross-covariance estimates between the measurements and the model unknowns using streamline trajectories. We illustrate the approach with a synthetic example and a large field-study that demonstrate the difficulties with the traditional EnKF implementation. In both the numerical experiments, it is shown that these challenges are addressed using flow relevant conditioning of the cross-covariance matrix. By mitigating sampling error in the cross-covariance estimates, the proposed approach provides significant computational savings through the use of modest ensemble sizes, and consequently offers the opportunity for use with large field-scale groundwater and reservoir characterization studies.     

Comparison of Ensemble Kalman Filter groundwater-data assimilation methods based on stochastic moment equations and Monte Carlo simulation

Traditional Ensemble Kalman Filter (EnKF) data assimilation requires computationally intensive Monte Carlo (MC) sampling, which suffers from filter inbreeding unless the number of simulations is large. Recently we proposed an alternative EnKF groundwater-data assimilation method that obviates the need for sampling and is free of inbreeding issues. In our new approach, theoretical ensemble moments are approximated directly by solving a system of corresponding stochastic groundwater flow equations. Like MC-based EnKF, our moment equations (ME) approach allows Bayesian updating of system states and parameters in real-time as new data become available. Here we compare the performances and accuracies of the two approaches on two-dimensional transient groundwater flow toward a well pumping water in a synthetic, randomly heterogeneous confined aquifer subject to prescribed head and flux boundary conditions.