On regularized time varying gravity field models based on grace data and their comparison with hydrological models

Determination of spherical harmonic coefficients of the Earth’s gravity field is often an ill-posed problem and leads to solving an ill-conditioned system of equations. Inversion of such a system is critical, as small errors of data will yield large variations in the result. Regularization is a method to solve such an unstable system of equations. In this study, direct methods of Tikhonov, truncated and damped singular value decomposition and iterative methods of ν, algebraic reconstruction technique, range restricted generalized minimum residual and conjugate gradient are used to solve the normal equations constructed based on range rate data of the gravity field and climate experiment (GRACE) for specific periods. Numerical studies show that the Tikhonov regularization and damped singular value decomposition methods for which the regularization parameter is estimated using quasioptimal criterion deliver the smoothest solutions. Each regularized solution is compared to the global land data assimilation system (GLDAS) hydrological model. The Tikhonov regularization with L-curve delivers a solution with high correlation with this model and a relatively small standard deviation over oceans. Among iterative methods, conjugate gradient is the most suited one for the same reasons and it has the shortest computation time.     

Predation of the beetle Rhantus sikkimensis (Coleoptera: Dytiscidae) on the larvae of Chironomus Meigen (Diptera: Chironomidae) of the Darjeeling Himalayas of India

The dytiscid beetle Rhantus sikkimensis, Regimbart, 1899 (Coleoptera: Dytiscidae), a member of the freshwater insect communities of the Darjeeling Himalayas, were noted to predate on the coexisting larvae of Chironomus sp. Meigen. Evaluation of predation by R. sikkimensis on Chironomus sp. larvae, in the laboratory, revealed that a single adult morph of R. sikkimensis could kill and consume on an average 10–90 and 10–78 numbers of small and large Chironomus sp. larvae, respectively, per day, depending on the prey density. The attack rate ranged between 520 and 537, and the handling time ranged between 4.3 and 8.6 depending on the size of the preys. The predation varied with respect to predator density also, with a maximum of 151 larvae killed by three predators per day. Two indices of predation, ingestion rate (IR) ranging between 13.33 and 74.15 larvae/day/predator and clearance rate (CR) ranging between 19.67 and 39.99 L prey/day/predator, varied with the prey size and predator density, significantly, when the predation was observed for 9 consecutive days, at two predator densities. It was also noted that R. sikkimensis predated on an average 9.8 larvae of Chironomus sp. and 1 larva of Culex sp., when the larvae of both the species are present together as preys, showing a preference for the Chironomus sp. larvae.     

Exploring the effect of data assimilation by WRF‐3DVar for numerical rainfall prediction with different types of storm events

The mesoscale Numerical Weather Prediction (NWP) model is gaining popularity among the hydrometeorological community in providing high‐resolution rainfall forecasts at the catchment scale. Although the performance of the model has been verified in capturing the physical processes of severe storm events, the modelling accuracy is negatively affected by significant errors in the initial conditions used to drive the model. Several meteorological investigations have shown that the assimilation of real‐time observations, especially the radar data can help improve the accuracy of the rainfall predictions given by mesoscale NWP models. The aim of this study is to investigate the effect of data assimilation for hydrological applications at the catchment scale. Radar reflectivity together with surface and upper‐air meteorological observations is assimilated into the Weather Research and Forecasting (WRF) model using the three‐dimensional variational data‐assimilation technique. Improvement of the rainfall accumulation and its temporal variation after data assimilation is examined for four storm events in the Brue catchment (135.2 km²) located in southwest England. The storm events are selected with different rainfall distributions in space and time. It is found that the rainfall improvement is most obvious for the events with one‐dimensional evenness in either space or time. The effect of data assimilation is even more significant in the innermost domain which has the finest spatial resolution. However, for the events with two‐dimensional unevenness of rainfall, i.e. the rainfall is concentrated in a small area and in a short time period, the effect of data assimilation is not ideal. WRF fails in capturing the whole process of the highly convective storm with densely concentrated rainfall in a small area and a short time period. A shortened assimilation time interval together with more efficient utilisation of the weather radar data might help improve the effectiveness of data assimilation in such cases. Copyright © 2012 John Wiley & Sons, Ltd.     

Feed-forward vs recurrent neural network models for non-stationarity modelling using data assimilation and adaptivity

Abstract

Artificial neural networks (ANN) are nonlinear models widely investigated in hydrology due to their properties of universal approximation and parsimony. Their performance during the training phase is very good, and their ability to generalize can be improved by using regularization methods such as early stopping and cross-validation. In our research, two kinds of generic models are implemented: the feed-forward model and the recurrent model. At first glance, the feed-forward model would seem to be more effective than the recurrent one on non-stationary datasets, because measured information on the state of the system (measured discharge) is used as input, thereby implementing a kind of data assimilation. This study investigates the feasibility and effectiveness of data assimilation and adaptivity when implemented in both feed-forward and recurrent neural networks. Based on the IAHS Workshop held in Göteborg, Sweden (July 2013), the hydrological behaviour of two watersheds of different sizes and different kind of non-stationarity will be modelled: (a) the Fernow watershed (0.2 km²) in the USA, affected by significant modifications in land cover during the study period, and (b) the Durance watershed (2170 km²) in France, affected by an increase in temperature that is causing a decrease in the extent of glaciers. Two methods were applied to evaluate the ability of ANN to adapt on the test set: (i) adaptivity using observed data to adapt parameter values in real time; and (ii) data assimilation using observed data to modify inaccurate inputs in real time. The goal of the study is thus re-analysis and not forecasting. This study highlights how effective the feed-forward model is compared to the recurrent model for dealing with non-stationarity. It also shows that adaptivity and data assimilation improve the recurrent model considerably, whereas improvement is marginal for the feed-forward model in the same conditions. Finally, this study suggests that adaptivity is effective in the case of changing conditions of the watershed, whereas data assimilation is better in the case of climate change (inputs modification).     

Snow water equivalent prediction using Bayesian data assimilation methods

Using the U.S. National Weather Service’s SNOW-17 model, this study compares common sequential data assimilation methods, the ensemble Kalman filter (EnKF), the ensemble square root filter (EnSRF), and four variants of the particle filter (PF), to predict seasonal snow water equivalent (SWE) within a small watershed near Lake Tahoe, California. In addition to SWE estimation, the various data assimilation methods are used to estimate five of the most sensitive parameters of SNOW-17 by allowing them to evolve with the dynamical system. Unlike Kalman filters, particle filters do not require Gaussian assumptions for the posterior distribution of the state variables. However, the likelihood function used to scale particle weights is often assumed to be Gaussian. This study evaluates the use of an empirical cumulative distribution function (ECDF) based on the Kaplan–Meier survival probability method to compute particle weights. These weights are then used in different particle filter resampling schemes. Detailed analyses are conducted for synthetic and real data assimilation and an assessment of the procedures is made. The results suggest that the particle filter, especially the empirical likelihood variant, is superior to the ensemble Kalman filter based methods for predicting model states, as well as model parameters.     

Neural networks in auroral data assimilation

Data assimilation is an essential step for improving space weather forecasting by means of a weighted combination between observational data and data from a mathematical model. In the present work data assimilation methods based on Kalman filter (KF) and artificial neural networks are applied to a three-wave model of auroral radio emissions. A novel data assimilation method is presented, whereby a multilayer perceptron neural network is trained to emulate a KF for data assimilation by using cross-validation. The results obtained render support for the use of neural networks as an assimilation technique for space weather prediction.     

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

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

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

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

Assimilation of sea surface temperature into CECOM by flux correction

Sea surface temperature (SST) from a near real-time data set produced from satellites data has been assimilated into a coupled ice–ocean forecasting model (Canadian East Coast Ocean Model) using an efficient data assimilation method. The method is based on an optimal interpolation scheme by which SST is melded into the model through the adjustment of surface heat flux. The magnitude and space–time variation of the adjustment depend on the depth of heat diffusion into the water column in response to changes in surface flux, the correlation time scale of the data, and model and data errors. The diffusion depth is scaled by the eddy diffusivity for temperature. The ratio of the model and data errors is treated as an adjustable parameter. To evaluate the quality of the assimilation, the results from the model with and without assimilation are compared to independent ship data from the Atlantic Zone Monitoring Program and the World Ocean Circulation Experiment. It is shown that the assimilation has a significant impact on the modeled SST, reducing the root mean square difference (RMSD) between the model SST and the ship SST by 0.63°C or 37%. The RMSD of the assimilated SST is smaller than that of the satellite SST by 0.23°C. This suggests that model simulations or predictions with data assimilation can provide the best estimate of the true SST. A sensitivity study is performed to examine the change of the model RMSD with the adjustable parameter in the assimilation equation. The results show that there is an optimal value of the parameter and the model SST is not very sensitive to the parameter.     

Improving hydrodynamic modeling of river networks by incorporating data assimilation using a particle filter

Numerical modeling is a well-recognized method for studying the hydrodynamic processes in river networks. Multi-source measurements also offer abundant information on the patterns and mechanisms within the processes. Therefore, improving hydrodynamic modeling of river networks through the use of data assimilation techniques has become a hot research topic in recent years. The particle filter (PF) is a commonly used data assimilation method and has been proven to be applicable to various nonlinear and non-Gaussian models. In the current study, an improved numerical hydrodynamic model for large-scale river networks is established by incorporating the advanced PF algorithm. Furthermore, the PF method based on the Gaussian likelihood function (GLF) and the method based on the Cauchy likelihood function (CLF) are compared for a complex river network scenario. The feasibility of the PF-based methods was evaluated through application to the Yangtze-Dongting River-lake Network (YDRN) by assimilating water stage data collected at six hydrometric stations during the entire hydrodynamic process in 2003. Additionally, the parameters used in the likelihood function, which affect the assimilation performance, also were explored in the current study. The study results found that the accuracy of the model-derived water stage data was improved when the PF-based methods are utilized, with improvement not only at the data assimilation (calibration) sites but also at three hydrometric stations not used in the data assimilation (i.e., verification sites). The highest average Nash-Sutcliffe Efficiency result for the six assimilation sites were 0.98 while the lowest summed root-mean-square-error result was 1.801 m. The comparison results also indicated that the CLF-based PF outperformed the GLF-based PF when high-accuracy observed data are available. Specifically, the CLF can effectively resolve the filtering failure problem and the dispersion problem of PFs, and further improve the accuracy of the filtering results for a river network scenario. In summary, the CLF-based PF method along with high-accuracy observation data shows promise to provide reliable reference and technical support for hydrodynamic modeling of large-scale river networks.     

Estimation of model errors generated by atmospheric forcings for ocean data assimilation: experiments in a regional model of the Bay of Biscay

The characterization of model errors is an essential step for effective data assimilation into open-ocean and shelf-seas models. In this paper, we propose an experimental protocol to properly estimate the error statistics generated by imperfect atmospheric forcings in a regional model of the Bay of Biscay, nested in a basin-scale North Atlantic configuration. The model used is the Hybrid Coordinate Ocean Model (HYCOM), and the experimental protocol involves Monte Carlo (or ensemble) simulations. The spatial structure of the model error is analyzed using the representer technique, which allows us to anticipate the subsequent impact in data assimilation systems. The results show that the error is essentially anisotropic and inhomogeneous, affecting mainly the model layers close to the surface. Even when the forcings errors are centered around zero, a divergence is observed between the central forecast and the mean forecast of the Monte Carlo simulations as a result of nonlinearities. The 3D structure of the representers characterizes the capacity of different types of measurement (sea level, sea surface temperature, surface velocities, subsurface temperature, and salinity) to control the circulation. Finally, data assimilation experiments demonstrate the superiority of the proposed methodology for the implementation of reduced-order Kalman filters.     

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.     

The role of prior model calibration on predictions with ensemble Kalman filter

This paper, based on a real world case study (Limmat aquifer, Switzerland), compares inverse groundwater flow models calibrated with specified numbers of monitoring head locations. These models are updated in real time with the ensemble Kalman filter (EnKF) and the prediction improvement is assessed in relation to the amount of monitoring locations used for calibration and updating. The prediction errors of the models calibrated in transient state are smaller if the amount of monitoring locations used for the calibration is larger. For highly dynamic groundwater flow systems a transient calibration is recommended as a model calibrated in steady state can lead to worse results than a noncalibrated model with a well-chosen uniform conductivity. The model predictions can be improved further with the assimilation of new measurement data from on-line sensors with the EnKF. Within all the studied models the reduction of 1-day hydraulic head prediction error (in terms of mean absolute error [MAE]) with EnKF lies between 31% (assimilation of head data from 5 locations) and 72% (assimilation of head data from 85 locations). The largest prediction improvements are expected for models that were calibrated with only a limited amount of historical information. It is worthwhile to update the model even with few monitoring locations as it seems that the error reduction with EnKF decreases exponentially with the amount of monitoring locations used. These results prove the feasibility of data assimilation with EnKF also for a real world case and show that improved predictions of groundwater levels can be obtained.     

An effective method based on dynamic sampling for data assimilation in a global wave model

The ensemble Kalman filter (EnKF) performs well because that the covariance of background error is varying along time. It provides a dynamic estimate of background error and represents the reasonable statistic characters of background error. However, high computational cost due to model ensemble in EnKF is employed. In this study, two methods referred as static and dynamic sampling methods are proposed to obtain a good performance and reduce the computation cost. Ensemble adjustment Kalman filter (EAKF) method is used in a global surface wave model to examine the performance of EnKF. The 24-h interval difference of simulated significant wave height (SWH) within 1 year is used to compose the static samples for ensemble errors, and these errors are used to construct the ensemble states at each time the observations are available. And then, the same method of updating the model states in the EAKF is applied for the ensemble states constructed by a static sampling method. The dynamic sampling method employs a similar method to construct the ensemble states, but the period of the simulated SWH is changing with time. Here, 7 days before and after the observation time is used as this period. To examine the performance of three schemes, EAKF, static, or dynamic sampling method, observations from satellite Jason-2 in 2014 are assimilated into a global wave model, and observations from satellite Saral are used for validation. The results indicate that the EAKF performs best, while the static sampling method is relatively worse. The dynamic sampling method improves an assimilation effect dramatically compared to the static sampling method, and its overall performance is closed to the EAKF. In low latitudes, the dynamic sampling method has a slight advantage over the EAKF. In the dynamic or static sampling methods, only one wave model is required to run and their computational cost is reduced sharply. According to the performance of these three methods, the dynamic sampling method can treated as an effective alternative of EnKF, which could reduce the computational cost and provide a good performance of data assimilation.     

Environmental toxicology: population modeling of cod larvae shows high sensitivity to loss of zooplankton prey

Two factors determine whether pollution is likely to affect a population indirectly through loss of prey: firstly, the sensitivity of the prey to the pollutants, and secondly, the sensitivity of the predator population to loss of prey at the given life stage. We here apply a statistical recruitment model for Northeast Arctic cod to evaluate the sensitivity of cod cohorts to loss of zooplankton prey, for example following an oil spill. The calculations show that cod cohorts are highly sensitive to possible zooplankton biomass reductions in the distribution area of the cod larvae, and point to a need for more knowledge about oil-effects on zooplankton. Our study illustrates how knowledge about population dynamics may guide which indirect effects to consider in environmental impact studies.     

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.     

A finite element model for simulating surface run‐off and unsaturated seepage flow in the shallow subsurface

This work introduces water–air two‐phase flow into integrated surface–subsurface flow by simulating rainfall infiltration and run‐off production on a soil slope with the finite element method. The numerical model is formulated by partial differential equations for hydrostatic shallow flow and water–air two‐phase flow in the shallow subsurface. Finite element computing formats and solution strategies are presented to obtain a numerical solution for the coupled model. An unsaturated seepage flow process is first simulated by water–air two‐phase flow under the atmospheric pressure boundary condition to obtain the rainfall infiltration rate. Then, the rainfall infiltration rate is used as an input parameter to solve the surface run‐off equations and determine the value of the surface run‐off depth. In the next iteration, the pressure boundary condition of unsaturated seepage flow is adjusted by the surface run‐off depth. The coupling process is achieved by updating the rainfall infiltration rate and surface run‐off depth sequentially until the convergence criteria are reached in a time step. A well‐conducted surface run‐off experiment and traditional surface–subsurface model are used to validate the new model. Comparisons with the traditional surface–subsurface model show that the initiation time of surface run‐off calculated by the proposed model is earlier and that the water depth is larger, thus providing values that are closer to the experimental results.     

Evolutionary assimilation of streamflow in distributed hydrologic modeling using in-situ soil moisture data

This study has applied evolutionary algorithm to address the data assimilation problem in a distributed hydrological model. The evolutionary data assimilation (EDA) method uses multi-objective evolutionary strategy to continuously evolve ensemble of model states and parameter sets where it adaptively determines the model error and the penalty function for different assimilation time steps. The assimilation was determined by applying the penalty function to merge background information (i.e., model forecast) with perturbed observation data. The assimilation was based on updated estimates of the model state and its parameterizations, and was complemented by a continuous evolution of competitive solutions.The EDA was illustrated in an integrated assimilation approach to estimate model state using soil moisture, which in turn was incorporated into the soil and water assessment tool (SWAT) to assimilate streamflow. Soil moisture was independently assimilated to allow estimation of its model error, where the estimated model state was integrated into SWAT to determine background streamflow information before they are merged with perturbed observation data. Application of the EDA in Spencer Creek watershed in southern Ontario, Canada generates a time series of soil moisture and streamflow. Evaluation of soil moisture and streamflow assimilation results demonstrates the capability of the EDA to simultaneously estimate model state and parameterizations for real-time forecasting operations. The results show improvement in both streamflow and soil moisture estimates when compared to open-loop simulation, and a close matching between the background and the assimilation illustrates the forecasting performance of the EDA approach.     

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

为了探讨水文模型在不同水文数据同化方案下的径流模拟差异,本文采用集合卡尔曼滤波算法,以遥感蒸散发产品、实测径流为观测数据,构建了基于新安江模型的数据同化框架。基于此框架设计了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)方案表现最差,而其余方案差异并不显著。本研究有助于进一步了解不同数据同化方案在径流模拟中的差异,从而为水资源高效利用与科学管理提供科学依据...