Evaluating and comparing methods of sinkhole susceptibility mapping in the Ebro Valley evaporite karst (NE Spain)

Multiple sinkhole susceptibility models have been generated in three study areas of the Ebro Valley evaporite karst (NE Spain) applying different methods (nearest neighbour distance, sinkhole density, heuristic scoring system and probabilistic analysis) for each sinkhole type separately (cover collapse sinkholes, cover and bedrock collapse sinkholes and cover and bedrock sagging sinkholes). The quantitative and independent evaluation of the predictive capability of the models reveals that: (1) The most reliable susceptibility models are those derived from the nearest neighbour distance and sinkhole density. These models can be generated in a simple and rapid way from detailed geomorphological maps. (2) The reliability of the nearest neighbour distance and density models is conditioned by the degree of clustering of the sinkholes. Consequently, the karst areas in which sinkholes show a higher clustering are a priori more favourable for predicting new occurrences. (3) The predictive capability of the best models obtained in this research is significantly higher (12.5–82.5%) than that of the heuristic sinkhole susceptibility model incorporated into the General Urban Plan for the municipality of Zaragoza. Although the probabilistic approach provides lower quality results than the methods based on sinkhole proximity and density, it helps to identify the most significant factors and select the most effective mitigation strategies and may be applied to model susceptibility in different future scenarios.     

Simulated multi-scale watershed runoff and sediment production based on GeoWEPP model

The runoff and sediment yield data from the Qiaozidonggou, Qiaozixigou, and Lu＇ergou watersheds, in the Loess Plateau of China are used to calibrate and validate the runoff and sediment yield simulated by GeoWEPP model of the WEPP Model at watershed scale. The indices of relative error, R, correlation coefficient, Re, and Nash-Suttcliffe efficiency coefficient Ens are used to evaluate the model fit. The eco-hydrological responses in the Luoyugou and Lu＇ergou watersheds are also forecast based on the WEPP Model. Meanwhile, the relation between vegetation pattern changes and sediment yield in the watershed is discussed, and the responses of runoff and sediment yield in the watersheds concerning forest growth stages are studied. The results show that the relative errors of simulated values of runoff and sediment yield are below 30%, the correlation coefficients axe above 0.90, and the Nash-Suttcliffe efficiency coefficients axe above 0.80. The simulation results present satisfactory performance, thus, the model could be used to simulate the runoff and sediment yield in these small watersheds. It is also observed that soil erosion tended to become severe as precipitation increased in the watershed, while soil erosion has a decreasing trend as forest cover increases and vegetation pattern is optimized. When the watershed is fully covered by forest, erosion and sediment yield are minimized. When the forest cover is about 30% and evenly distributed in the watershed, the erosion intensity is lower than if the forest cover is collectively distributed in the watershed. Erosion varies with different forest growth stages in the watershed; it is more serious at the young and near planting stage and is the smallest at the mature forest stage.     

A comparative modeling study of a dual tracer experiment in a large lysimeter under atmospheric conditions

In this paper, five model approaches with different physical and mathematical concepts varying in their model complexity and requirements were applied to identify the transport processes in the unsaturated zone. The applicability of these model approaches were compared and evaluated investigating two tracer breakthrough curves (bromide, deuterium) in a cropped, free-draining lysimeter experiment under natural atmospheric boundary conditions. The data set consisted of time series of water balance, depth resolved water contents, pressure heads and resident concentrations measured during 800 days. The tracer transport parameters were determined using a simple stochastic (stream tube model), three lumped parameter (constant water content model, multi-flow dispersion model, variable flow dispersion model) and a transient model approach. All of them were able to fit the tracer breakthrough curves. The identified transport parameters of each model approach were compared. Despite the differing physical and mathematical concepts the resulting parameters (mean water contents, mean water flux, dispersivities) of the five model approaches were all in the same range. The results indicate that the flow processes are also describable assuming steady state conditions. Homogeneous matrix flow is dominant and a small pore volume with enhanced flow velocities near saturation was identified with variable saturation flow and transport approach. The multi-flow dispersion model also identified preferential flow and additionally suggested a third less mobile flow component. Due to high fitting accuracy and parameter similarity all model approaches indicated reliable results.     

Parameter-independent model reduction of transient groundwater flow models: Application to inverse problems

A new methodology is proposed for the development of parameter-independent reduced models for transient groundwater flow models. The model reduction technique is based on Galerkin projection of a highly discretized model onto a subspace spanned by a small number of optimally chosen basis functions. We propose two greedy algorithms that iteratively select optimal parameter sets and snapshot times between the parameter space and the time domain in order to generate snapshots. The snapshots are used to build the Galerkin projection matrix, which covers the entire parameter space in the full model. We then apply the reduced subspace model to solve two inverse problems: a deterministic inverse problem and a Bayesian inverse problem with a Markov Chain Monte Carlo (MCMC) method. The proposed methodology is validated with a conceptual one-dimensional groundwater flow model. We then apply the methodology to a basin-scale, conceptual aquifer in the Oristano plain of Sardinia, Italy. Using the methodology, the full model governed by 29,197 ordinary differential equations is reduced by two to three orders of magnitude, resulting in a drastic reduction in computational requirements.     

Intercomparison of the JULES and CABLE land surface models through assimilation of remotely sensed soil moisture in southeast Australia

Numerous land surface models exist for predicting water and energy fluxes in the terrestrial environment. These land surface models have different conceptualizations (i.e., process or physics based), together with structural differences in representing spatial variability, alternate empirical methods, mathematical formulations and computational approach. These inherent differences in modeling approach, and associated variations in outputs make it difficult to compare and contrast land surface models in a straight-forward manner. While model intercomparison studies have been undertaken in the past, leading to significant progress on the improvement of land surface models, additional framework towards identification of model weakness is needed. Given that land surface models are increasingly being integrated with satellite based estimates to improve their prediction skill, it is practical to undertake model intercomparison on the basis of soil moisture data assimilation. Consequently, this study compares two land surface models: the Joint UK Land Environment Simulator (JULES) and the Community Atmosphere Biosphere Land Exchange (CABLE) for soil moisture estimation and associated assessment of model uncertainty. A retrieved soil moisture data set from the Soil Moisture and Ocean Salinity (SMOS) mission was assimilated into both models, with their updated estimates validated against in-situ soil moisture in the Yanco area, Australia. The findings show that the updated estimates from both models generally provided a more accurate estimate of soil moisture than the open loop estimate based on calibration alone. Moreover, the JULES output was found to provide a slightly better estimate of soil moisture than the CABLE output at both near-surface and deeper soil layers. An assessment of the updated membership in decision space also showed that the JULES model had a relatively stable, less sensitive, and more highly convergent internal dynamics than the CABLE model.     

Short-term optimal operation of water systems using ensemble forecasts

Short-term water system operation can be realized using Model Predictive Control (MPC). MPC is a method for operational management of complex dynamic systems. Applied to open water systems, MPC provides integrated, optimal, and proactive management, when forecasts are available. Notwithstanding these properties, if forecast uncertainty is not properly taken into account, the system performance can critically deteriorate.Ensemble forecast is a way to represent short-term forecast uncertainty. An ensemble forecast is a set of possible future trajectories of a meteorological or hydrological system. The growing ensemble forecasts’ availability and accuracy raises the question on how to use them for operational management.The theoretical innovation presented here is the use of ensemble forecasts for optimal operation. Specifically, we introduce a tree based approach. We called the new method Tree-Based Model Predictive Control (TB-MPC). In TB-MPC, a tree is used to set up a Multistage Stochastic Programming, which finds a different optimal strategy for each branch and enhances the adaptivity to forecast uncertainty. Adaptivity reduces the sensitivity to wrong forecasts and improves the operational performance.TB-MPC is applied to the operational management of Salto Grande reservoir, located at the border between Argentina and Uruguay, and compared to other methods.     

紫坪铺水库造成孔隙弹性耦合变化及其对2008年汶川地震触发作用

位于川西地区龙门山断裂带附近的紫坪铺水库于2005年9月开始蓄水.约2.7年后,2008年5月12日,M_w7.9级汶川地震在龙门山断裂带上发生,两事件在时间和空间上的接近,揭示其可能相互关联,但前人的诸多研究给出了不同甚至是相反的结果.本研究基于完全耦合孔隙弹性理论,利用二维有限元模型(FEM),模拟水库蓄水造成的区域孔隙压力场和应力场的演化过程,基于库仑应力演化探讨其对龙门山断裂带活动的影响.模拟结果显示紫坪铺水库蓄水打破了原来的区域孔隙压力平衡,形成孔隙压力梯度源,向周围地壳传播;进而造成龙门山断裂带上库仑应力正值范围不断扩大,由浅入深影响到整条断层,尤其对浅层范围的加载作用明显,达上百千帕,为整个断层面的失稳提供了基础.震源区域库仑应力呈持续增长趋势,汶川地震发震前,增长了约数千帕~数十千帕,即使初期库仑应力为负,在未来某时刻库仑应力仍可能由负转正,并不断增大.通过计算汶川地震震源及其附近区域内相对于多种可能断层倾角的库仑应力,发现库仑应力随断层倾角增大而增加.因此整体来看,紫坪铺水库蓄水对龙门山断裂带起加载作用,有可能触发地震.对紫坪铺库区周围的小震分析也显示,蓄水以来小震明显增多、且不同时段小震发生密集区与水库距离逐渐增大,与孔隙压力扩散趋势一致.以上结果表明紫坪铺水库的蓄水增加了汶川地震的危险性.     

Development of a national VNIR soil-spectral library for soil classification and prediction of organic matter concentrations

Soil visible-near infrared diffuse reflectance spectroscopy(vis-NIR DRS)has become an important area of research in the fields of remote and proximal soil sensing.The technique is considered to be particularly useful for acquiring data for soil digital mapping,precision agriculture and soil survey.In this study,1581 soil samples were collected from 14 provinces in China,including Tibet,Xinjiang,Heilongjiang,and Hainan.The samples represent 16 soil groups of the Genetic Soil Classification of China.After air-drying and sieving,the diffuse reflectance spectra of the samples were measured under laboratory conditions in the range between 350 and 2500 nm using a portable vis-NIR spectrometer.All the soil spectra were smoothed using the Savitzky-Golay method with first derivatives before performing multivariate data analyses.The spectra were compressed using principal components analysis and the fuzzy k-means method was used to calculate the optimal soil spectral classification.The scores of the principal component analyses were classified into five clusters that describe the mineral and organic composition of the soils.The results on the classification of the spectra are comparable to the results of other similar research.Spectroscopic predictions of soil organic matter concentrations used a combination of the soil spectral classification with multivariate calibration using partial least squares regression(PLSR).This combination significantly improved the predictions of soil organic matter(R2=0.899;RPD=3.158)compared with using PLSR alone(R2=0.697;RPD=1.817).     

Atmospheric correction of HJ-1 A/B CCD over land: Land surface reflectance calculation for geographical information product

This paper proposed a method to retrieve the land surface reflectance from the HJ-1A/B CCD data. The aerosol optical depth（AOD）, the most important factor affecting the atmospheric correction of CCD images at all bands, is proposed to retrieve from the CCD imagery by the approach of dense dark vegetation（DDV） method. A look-up table in terms of the transmittances, the path radiances and the atmospheric spherical albedo as functions of the AOD was established for a variety of sun-sensor geometry and aerosol loadings. The atmospheric correction is then achieved with the look-up table and the MODIS surface reflectance output（MOD09） as the priori datasets. Based on the retrieved AOD and the look-up table of atmospheric correction coefficients, the land surface reflectance was retrieved for the HJ-1A/B data according to the atmospheric radiative transfer equation. Some in-situ measurement Data for Yanzhou of Shandong province in East China and MODIS land surface reflectance products MOD09 are used to preliminarily validate the proposed method. The results show that the proposed method can remove effectively the atmospheric contributions, and the overall accuracy of the retrieval land surface reflectance can be improved substantially.     

简单生物圈模式 (SiB2) 中湍流能量通量对近地层两类阻抗系数的敏感性研究

陆—气之间的能量交换是通过近地层湍流热量和水汽通量来实现的.以往的研究表明近地层阻抗对湍流能量通量有着不同程度的影响,但是关于陆面模式中阻抗系数的取值范围却始终没有统一的标准.为了深入了解简单生物圈模式 (SiB₂) 中近地层阻抗系数取值变化对湍流能量通量的影响,我们以那曲站为例,分别采用传统的逐个因子分析法和考虑参数间相互作用的部分因子分析法定量地研究了夏季该观测站近地层湍流能量通量分别对冠层阻抗系数C₁和地表阻抗系数C₂的敏感性响应.结果表明,感热通量对地表阻抗系数C₂更为敏感,而潜热通量则对冠层阻抗系数C₁较为敏感;感热通量随C₁增加而增大,随C₂增加而减小,而潜热通量则随C₁或C₂的增加而减小;不管是感热通量还是潜热通量,它们对阻抗系数的敏感度随阻抗系数的增大而减小,而对阻抗系数的相对敏感度则随阻抗系数的增大而增大.最后,结合那曲站夏季下垫面稀疏短草的分布特点分析了造成感热通量和潜热通量敏感变化各异的原因.