科尔沁沙地土壤水分动态分析及其概率密度函数模拟

基于Laio土壤水分动态随机模型(Laio模型),利用2006-2010年5~9月土壤水分连续监测数据及日降水资料,分析科尔沁沙地固定沙丘和沙质草地生长季根系层土壤水分动态及其与降水格局的关系,研究点尺度土壤水分概率密度函数,并对Laio模型涉及的13个参数进行了敏感性分析。结果表明:① 研究区年降水的季节分配极不均匀,主要集中在4~10月的生长季,占全年降水量的93%;0~5 mm降水事件占全年降水事件的73%,但其降水量只占全年降水量的25%;降水间隔期以0~10 d为主,占全年无降水期的38%,其频数最高,占全年间隔期频数的87%。② 固定沙丘和沙质草地根系层厚度分别为0~100 cm和0~70 cm,沙质草地根系层土壤水分显著高于固定沙丘;两类沙地7月份的土壤水分都显著高于生长季其他月份。③ 两类沙地生长季根系层土壤水分均服从正态分布;通过Laio模型得到了两类沙地生长季根系层土壤水分概率密度函数p(s),其峰值及峰值出现的位置和峰的阔度均与观测结果很接近,说明Laio模型能对科尔沁沙地土壤水分概率密度函数进行较好的模拟。④ Laio模型涉及的13个参数中,对p(s)最为敏感的参数是降水频率λ、平均降水量α、最大蒸散量E_max、水分胁迫点s*和凋萎系数s_w,主要影响p(s)曲线的峰值。     

三维地质建模软件发展现状及问题探讨

分析了三维地质建模软件可以明显地提高工作效率与工程质量,建立勘探区的三维地质模型,检查原始数据的正确性,通过各种“虚拟切割”手段检查各种地质要素间的相互关系。探讨了三维地质建模软件开发中所存在的几个问题,包括软件开发团队应适时转为公司机制,审慎选择开源工具,三维地质建模的三维模型数据标准,在勘探规范中增加三维地质建模要求,二、三维联动编辑,选准产品的用户群的与盈利模式等。     

Incremental growth of normal faults: Insights from a laser-equipped analog experiment

We conducted a laser-equipped analog experiment aimed at quasi-continuously monitoring the growth of a dense population of normal faults in homogeneous conditions. To further understand the way geological faults progressively gain in slip and length as they accumulate more strain, we measured with great precision the incremental slip and length changes that the analog faults sustain as they grow. These measurements show that the analog faults share common features with the natural ones. In particular, during their growth, the faults develop and maintain cumulative slip profiles that are generally triangular and asymmetric. The growth takes place through two distinct phases: an initial, short period of rapid lateral lengthening, followed by a longer phase of slip accumulation with little or no lengthening. The incremental slip is found to be highly variable in both space (along the faults) and time, resulting in variable slip rates. In particular, ‘short- and long-term’ slip rates are markedly different. We also find that slip measurements at local points on fault traces do not contain clear information on the slip increment repeat mode. Finally, while the fault growth process is highly heterogeneous when considered at the scale of a few slip events, it appears homogeneous and self-similar at longer time scales which integrate many slip increments. This is likely to be the result of a feedback between stress heterogeneities and slip development. The long-term scale homogeneity also implies that the long-term faulting process is primarily insensitive to the short-term heterogeneities that are rapidly smoothed or redistributed. We propose a new conceptual scenario of fault growth that integrates the above observations and we suggest that faults grow in a bimodal way as a result of a self-driven and self-sustaining process.     

Hydrogeochemistry and strontium isotopes in the Arno River Basin (Tuscany, Italy): Constraints on natural controls by statistical modeling

In this paper the chemistry of major and trace elements and, for the first time, strontium isotopic ratios measured in running waters from the Arno River Basin (Tuscany, central-northern Italy) and thermal springs discharging in the same hydrographic system are presented and discussed. Classical graphical methods (e.g. mixing diagrams) have here been improved to identify, in a correct statistical sample space, extreme chemical compositions attributable to the action of geochemical processes and/or inherited from specific lithologies (namely contributions or components 1, 2 and 3) to be used in inverse modeling procedures, due to the absence of clear end-members. A linear least squares problem, with non-negativity constraints and distances, as required for compositional data (convex linear mixing) was solved by considering the contribution of the most discriminant variables given by the ⁸⁷Sr/⁸⁶Sr isotopic ratios and the concentration of Ca²⁺, Mg²⁺, Sr and Rb. Following such approach, it can be assumed (p < 0.05) that component 1, characterized by a Ca²⁺(Mg²⁺)– radiogenic-rich (⁸⁷Sr/⁸⁶Sr = 0.71274; Rb/Sr = 0.039), represents a dominant feature at basin scale, with a weight ranging from 69% to 100%. Much lower percentages are related to component 2, represented by a Ca²⁺(Mg²⁺)– facies with intermediate ⁸⁷Sr/⁸⁶Sr (0.70874) and low Rb/Sr (2.8 × 10⁻⁴) ratios and component 3, identified by Ca²⁺– facies with less radiogenic ⁸⁷Sr/⁸⁶Sr (0.70827) and low Mg²⁺/Ca²⁺ (0.011) ratios. These components are mainly dominated by the dissolution of evaporitic rocks and/or mixing with thermal waters in the southern part of the Arno River Basin and by dissolution of the carbonatic fraction, kinetically favored with respect to that of the silicatic minerals, in the upper reaches of the main course and its tributaries, respectively.     

An empirical model of event scale cohesive bank profile evolution

The highly stochastic nature of riverbank erosion has driven the need for spatially explicit empirical models. Detailed bank profile surveys along a meander bend of the Brandywine Creek in Pennsylvania, USA, before and after 28 high flow events over a 2·5 year period are used to develop an empirical model of cohesive bank profile erosion. Two hundred and thirty‐six bank erosion observations are classified as hydraulic erosion or subaerial erosion. Threshold conditions required to initiate bank erosion cannot be defined based on field measurements. Using the near‐bank velocity and the number of freeze–thaw cycles as predictors, regression equations are derived for hydraulic erosion that specify the length, thickness, and location on the bank face of eroded blocks. An empirical discriminant function defines the critical geometry of overhang failures, and the volumes removed by overhang failures are computed using another regression equation. All the regression equations are significant, but have low correlation coefficients, suggesting that cohesive bank erosion has a strong stochastic component. Individual events typically remove small masses of soil (average volume 0·084 m³/m) a few centimeters thick (median = 0·057 m) and a few decimeters in length (median = 0·50 m) from the lower third of the bank. Hydraulic erosion is responsible for 87% of all erosion. When applied to three survey sites not used in its development, the profile model predicts the total volume of erosion with errors of 23%, 5% and 1%. Twenty‐four percent of computed erosion volumes for single events are within 50% of observed volumes at these three sites. Extending the approach to decadal timescales and to entire bends will require three‐dimensional observations of bank failure, and spatially and temporally explicit methods to account for the influence of individual large trees on bank failures and near‐bank hydraulic processes. Copyright © 2009 John Wiley & Sons, Ltd.     

The effects of interdune vegetation changes on eolian dune field evolution: a numerical‐modeling case study at Jockey's Ridge,North Carolina,USA

Changes in vegetation cover within dune fields can play a major role in how dune fields evolve. To better understand the linkage between dune field evolution and interdune vegetation changes, we modified Werner's (Geology, 23, 1995: 1107–1110) dune field evolution model to account for the stabilizing effects of vegetation. Model results indicate that changes in the density of interdune vegetation strongly influence subsequent trends in the height and area of eolian dunes. We applied the model to interpreting the recent evolution of Jockey's Ridge, North Carolina, where repeat LiDAR surveys and historical aerial photographs and maps provide an unusually detailed record of recent dune field evolution. In the absence of interdune vegetation, the model predicts that dunes at Jockey's Ridge evolve towards taller, more closely‐spaced, barchanoid dunes, with smaller dunes generally migrating faster than larger dunes. Conversely, the establishment of interdune vegetation causes dunes to evolve towards shorter, more widely‐spaced, parabolic forms. These results provide a basis for understanding the increase in dune height at Jockey's Ridge during the early part of the twentieth century, when interdune vegetation was sparse, followed by the decrease in dune height and establishment of parabolic forms from 1953‐present when interdune vegetation density increased. These results provide a conceptual model that may be applicable at other sites with increasing interdune vegetation cover, and they illustrate the power of using numerical modeling to model decadal variations in eolian dune field evolution. We also describe model results designed to test the relative efficacy of alternative strategies for mitigating dune migration and deflation. Installing sand‐trapping fences and/or promoting vegetation growth on the stoss sides of dunes are found to be the most effective strategies for limiting dune advance, but these strategies must be weighed against the desire of many park visitors to maintain the natural state of the dunes. Copyright © 2009 John Wiley & Sons, Ltd.     

Enhanced application of root‐reinforcement algorithms for bank‐stability modeling

Riparian vegetation is known to exert a number of mechanical and hydrologic controls on bank stability. In particular, plant roots provide mechanical reinforcement to a soil matrix due to the different responses of soils and roots to stress. Root reinforcement is largely a function of the strength of the roots crossing potential shear planes, and the number and diameter of such roots. However, previous bank stability models have been constrained by limited field data pertaining to the spatial and temporal variability of root networks within stream banks. In this paper, a method is developed to use root‐architecture data to derive parameters required for modeling temporal and spatial changes in root reinforcement. Changes in root numbers over time were assumed to follow a sigmoidal curve, which commonly represents the growth rates of organisms. Regressions for numbers of roots crossing potential shear planes over time showed small variations between species during the juvenile growth phase, but extrapolation led to large variations in root numbers by the time the senescent phase of the sigmoidal growth curve had been reached. In light of potential variability in the field data, the mean number of roots crossing a potential shear plane at each year of tree growth was also calculated using data from all species and an additional sigmoidal regression was run. After 30 years the mean number of roots predicted to cross a 1 m shear plane was 484, compared with species‐specific curves whose values ranged from 240 roots for black willow trees to 890 roots for western cottonwood trees. In addition, the effect of spatial variations in rooting density with depth on stream‐bank stability was modeled using the bank stability and toe erosion model (BSTEM). Three root distributions, all approximating the same average root reinforcement (5 kPa) over the top 1 m of the bank profile, were modeled, but with differing vertical distributions (concentrated near surface, non‐linear decline with depth, uniform over top meter). It was found that stream‐bank F_S varied the most when the proportion of the failure plane length to the depth of the rooting zone was greatest. Copyright © 2008 John Wiley & Sons, Ltd.     

基于三维模拟的海洋CSEM资料处理

海洋可控源电磁法已经成为海洋油气勘探一个重要工具,但是其资料处理和解释还处于定性和一维模拟阶段.在积分方程三维模拟的基础上对Troll油田实测数据进行了处理,采用人机交互三维模拟寻找背景模型和异常体初始模型,最后对异常体电阻率采用准线性近似快速反演,取得了定量的结果.同时,说明对于二维测线和二维模型依然可以用三维来模拟,其结果优于二维反演.在电子计算机技术快速发展的今天,可以预计三维反演将成为资料处理解释的主流.     

An easily implemented agro-hydrological procedure with dynamic root simulation for water transfer in the crop–soil system: Validation and application

Models for water transfer in the crop–soil system are key components of agro-hydrological models for irrigation, fertilizer and pesticide practices. Many of the hydrological models for water transfer in the crop–soil system are either too approximate due to oversimplified algorithms or employ complex numerical schemes. In this paper we developed a simple and sufficiently accurate algorithm which can be easily adopted in agro-hydrological models for the simulation of water dynamics. We used a dual crop coefficient approach proposed by the FAO for estimating potential evaporation and transpiration, and a dynamic model for calculating relative root length distribution on a daily basis. In a small time step of 0.001 d, we implemented algorithms separately for actual evaporation, root water uptake and soil water content redistribution by decoupling these processes. The Richards equation describing soil water movement was solved using an integration strategy over the soil layers instead of complex numerical schemes. This drastically simplified the procedures of modeling soil water and led to much shorter computer codes. The validity of the proposed model was tested against data from field experiments on two contrasting soils cropped with wheat. Good agreement was achieved between measurement and simulation of soil water content in various depths collected at intervals during crop growth. This indicates that the model is satisfactory in simulating water transfer in the crop–soil system, and therefore can reliably be adopted in agro-hydrological models. Finally we demonstrated how the developed model could be used to study the effect of changes in the environment such as lowering the groundwater table caused by the construction of a motorway on crop transpiration.     

Modeling the dynamic process of tsunami earthquake by liquid-solid coupling model

Tsunami induced by earthquake is an interaction problem between liquid and solid.Shallow-water wave equation is often used to modeling the tsunami,and the boundary or initial condition of the problem is determined by the displacement or velocity field from the earthquake under sea floor,usually no interaction between them is consid-ered in pure liquid model.In this study,the potential flow theory and the finite element method with the interaction between liquid and solid are employed to model the dynamic processes of the earthquake and tsunami.For model-ing the earthquake,firstly the initial stress field to generate the earthquake is set up,and then the occurrence of the earthquake is simulated by suddenly reducing the elastic material parameters inside the earthquake fault.It is dif-ferent from seismic dislocation theory in which the relative slip on the fault is specified in advance.The modeling results reveal that P,SP and the surface wave can be found at the sea surface besides the tsunami wave.The surface wave arrives at the distance of 600 km from the epicenter earlier than the tsunami 48 minutes,and its maximum amplitude is 0.55 m,which is 2 times as large as that of the sea floor.Tsunami warning information can be taken from the surface wave on the sea surface,which is much earlier than that obtained from the seismograph stations on land.The tsunami speed on the open sea with 3 km depth is 175.8 m/s,which is a little greater than that pre-dicted by long wave theory,(gh)1/2=171.5 m,and its wavelength and amplitude in average are 32 km and 2 m,respectively.After the tsunami propagates to the continental shelf,its speed and wavelength is reduced,but its amplitude become greater,especially,it can elevate up to 10 m and run 55 m forward in vertical and horizontal directions at sea shore,respectively.The maximum vertical accelerations at the epicenter on the sea surface and on the earthquake fault are 5.9 m/s2 and 16.5 m/s2,respectively,the later is 2.8 times the former,and therefore,sea water is a good shock