SURFACE ROUGHNESS EVOLUTION OF SOILS CONTAINING ROCK FRAGMENTS

Soil surface roughness is a dynamic property which determines, to a large extent, erosion and infiltration rates. Although soils containing rock fragments are widespread in the Mediterranean region, the effect of the latter on surface roughness evolution is yet poorly understood. Therefore, laboratory experiments were conducted in order to investigate the effect of rock fragment content, rock fragment size and initial moisture content of the fine earth on the evolution of interrill surface roughness during simulated rainfall. Surface elevations of simulated plough layers along transects of 50 cm length were measured before and after simulated rainfall (totalling 192.5 mm, I = 70 mm h⁻¹) with a laser microreliefmeter. The results were used to investigate whether systematic variations in interrill surface roughness along stony hillslopes in southeastern Spain could be attributed to rock fragment cover and rock fragment size. Soil surface elevations were measured along the contour lines (50 cm long transects) with a contact microreliefmeter. Roughness was expressed by two parameters related to the height and frequency of roughness elements, respectively: standard deviation of de-trended surface elevations (random roughness: RR), and correlation length (L) derived from exponential fits of the autocorrelation functions. The frequently used assumption that surface roughness (RR) of cultivated topsoils decreases exponentially with cumulative rain is not valid for soil surfaces covered by rock fragments. The RR of soils containing small rock fragments (1.7–2.7 cm) increased with cumulative rainfall after an initial decrease during the first 17.5 mm of rainfall. For soils containing large rock fragments (7.7 cm), RR increased with rainfall above a threshold rock fragment content by mass of 52 per cent. For a given rainfall application, RR increased non-linearly with rock fragment content. The correlation length for soils containing small rock fragments decreases with rock fragment content and is significantly lower than for soils with large rock fragments. Soils covered with small rock fragments (large RR and small L) are thus well protected against raindrop impact by a water film in the depressions between the rock fragments. On abandoned agricultural fields along hillslopes in southeastern Spain, rock fragments cover increases non-linearly with slope owing to selective erosion of finer particles on steep slopes. The increase of surface cover by large rock fragments (>25 mm) is even more pronounced. The simultaneous increase of rock fragment cover and rock fragment size with slope explains the non-linear increase of RR with slope. These relationships differ for soils covered by platy misaschists and those covered with cubic andesites. The variations in correlation length along the hillslopes are not clear, probably owing to a simultaneous increase in rock fragment cover and rock fragment size. These findings may provide a better prediction of soil surface roughness of interrill areas covered by rock fragments using slope angle and lithology.     

ESTIMATING INFILTRATION PARAMETERS FROM BASIC SOIL PROPERTIES

Infiltration data were collected on two rectangular grids with 25 sampling points each. Both experimental grids were located in tropical rain forest (Guyana), the first in an Arenosol area and the second in a Ferralsol field. Four different infiltration models were evaluated based on their performance in describing the infiltration data. The model parameters were estimated using non-linear optimization techniques. The infiltration behaviour in the Ferralsol was equally well described by the equations of Philip, Green–Ampt, Kostiakov and Horton. For the Arenosol, the equations of Philip, Green–Ampt and Horton were significantly better than the Kostiakov model. Basic soil properties such as textural composition (percentage sand, silt and clay), organic carbon content, dry bulk density, porosity, initial soil water content and root content were also determined for each sampling point of the two grids. The fitted infiltration parameters were then estimated based on other soil properties using multiple regression. Prior to the regression analysis, all predictor variables were transformed to normality. The regression analysis was performed using two information levels. The first information level contained only three texture fractions for the Ferralsol (sand, silt and clay) and four fractions for the Arenosol (coarse, medium and fine sand, and silt and clay). At the first information level the regression models explained up to 60% of the variability of some of the infiltration parameters for the Ferralsol field plot. At the second information level the complete textural analysis was used (nine fractions for the Ferralsol and six for the Arenosol). At the second information level a principal components analysis (PCA) was performed prior to the regression analysis to overcome the problem of multicollinearity among the predictor variables. Regression analysis was then carried out using the orthogonally transformed soil properties as the independent variables. Results for the Ferralsol data show that the parameters of the Green–Ampt and Kostiakov model were estimated relatively accurately (maximum R² = 0.76). For the Arenosol, use of the second information level together with PCA produced regression models with an R² value ranging from 0.38 to 0.68. For the Ferralsol, most of the variance was explained by the root content and organic matter content. In the Arenosol plot, the fractions medium and fine sand explained most of the observed variance.     

LISEM: A SINGLE-EVENT,PHYSICALLY BASED HYDROLOGICAL AND SOIL EROSION MODEL FOR DRAINAGE BASINS. II: SENSITIVITY ANALYSIS,VALIDATION AND APPLICATION

A new hydrological and soil erosion model has been developed and tested: LISEM, the Limburg soil erosion model. The model uses physically based equations to describe interception, infiltration and soil water transport, storage in surface depressions, splash and flow detachment, transport capacity and overland and channel flow. From the validation results it is clear that, although the model has several advantages over other models, the results of LISEM 1.0 are far from perfect. Based on the sensitivity analysis and field observations, the main reasons for these differences seems to be the spatial and temporal variability of the soil hydraulic conductivity and the initial pressure head at the basin scale. Another reason for the differences between measured and simulated results is our lack or understanding of the theory of hydrological and soil erosion processes.     

COMPUTING FIELD STATICS WITH THE HELP OF SEISMIC TOMOGRAPHY*

Field static corrections in general need be applied to all onshore seismic reflection data to eliminate the disturbing effects a weathering layer or near-surface low velocity zone has on the continuity of deep seismic reflections. The traveltimes of waves refracted at the bottom of the low velocity zone (or intermediate refracting interfaces) can often be observed as first breaks on shot records and used to develop a laterally inhomogeneous velocity model for this layer, from which the field static corrections can then be obtained. A simple method is described for computing accurate field statics from first breaks. It is based on a linearization principal for traveltimes and leads to the algorithms that are widely and successfully applied within the framework of seismic tomography. We refine an initial model for the low velocity layer (estimated by a standard traveltime inversion technique) by minimizing the errors between the observed first arrivals on field records and those computed by ray theory through an initial model of the low velocity layer. Thus, one can include more lateral velocity variations within the low velocity layers, which are important to obtain good field static corrections. Traditional first break traveltime inversion methods cannot, in general, provide such refined velocity values. The technique is successfully applied to seismic data from the Amazon Basin. It is based on a simple model for the low velocity layer that consists of an undulating earth surface and one planar horizontal refractor overlain by a laterally changing velocity field.     

旋转流场中的格子波耳兹曼模型

在大尺度的旋转流场中,由于存在哥氏力,使流体的流动出现了一系列复杂的动力学现象.在原格子Boltzmann模型研究的基础上,引入了哥氏力效应,发展了一个旋转流场中的格子Boltzmann模型.从该模型出发可导出地球流体力学方程,用这一模型对理想边界条件下的北半球大气环流进行了数值计算.数值结果很好地再现了大尺度地转流的流动特征.从理论和数值实验上验证了该模型的适用性.     

井间电磁成像系统应用研究

井间电磁（EM）成像系统的技术目标是实现井间电学特性的直接测量,并提供井间电阻率的二维及三维图像.从1997年开始,胜利油田与美国EMI公司合作,开展了EM技术的应用与研究,主要是通过大型、系统的现场试验,在大井间距的条件下,进一步验证仪器的性能,进行油藏研究的适用性和可行性分析;以系统的现场实测数据为依据,加快成像处理方法和软件的开发;开展穿透金属套管井的EM试验,探讨金属套管条件下的电阻率成像方法。两年来,利用EMI提供的XBH2000型测量系统,对胜利油田典型的低电阻率砂泥岩剖面,成功地进行了3对井、10个井次的试验。试验证明了仪器性能的可靠,取得了测量重复性好、系统完整的井间EM数据.并反演得到电阻率成像,在分析井间油气分布和油层开采动态方面,取得了良好的地质效果.     

INTERRELATIONSHIP BETWEEN WATER QUALITY AND GROUNDWATER FLOW DYNAMICS IN A SMALL WETLAND SYSTEM ALONG A SANDY HILL RIDGE

Detailed modelling of the hydrological setting of fen meadows appears to be possible provided that detailed information on geomorphology, hydrochemistry and piezometric heads is available for a number of years. In the Laegieskamp, a small wetland reserve located in the central part of The Netherlands, a piezometric monitoring network was sampled for water quality analysis and piezometric heads between 1986 and 1992. Average yearly discharge and recharge periods were used for FLOWNET calculations. First, the models were used to determine, with the help of information on water quality, the hydrological systems in the study area. Secondly, they were used to define the present and past hydrological setting of a fen meadow in the reserve. The hydrological systems and water quality in the study area have changed considerably over the past 65 years. At present the fen meadow is mainly fed by precipitation. The mineral-rich conditions favouring the fen meadow vegetation are thought to be maintained thanks to a clayey peat layer and an oscillating shallow water body that prevents rapid leaching of minerals. The sulphate content in the fen exhibits a pattern of temporal variation, which is related to the severity of the annual drought. Our study showed that groundwater flow is mainly lateral, instead of the assumed vertical infiltration of groundwater in previous regional studies. This led us to the conclusion that conservation and restoration perspectives are much better than previously expected. The polluted middle, deep groundwater is not a major threat to this fen at the moment. © 1997 John Wiley & Sons, Ltd.     

The Merzbacher & Eggler (1984) Geohygrometer: a Cautionary Note on its Suitability for High-K Suites

The use of the popular Merzbacher & Eggler (1984, Geology12, 587–590) experimental geohygrometer for calcalkalineandesites and dacites is critically evaluated and two pitfallsare found. First, calculation of the correct projection parametersis problematic because two endmember calculation schemes arefound in the literature; Baker & Eggler (1983, Journal ofVolcanology and Geothermal Research 18, 387–404) and Walkeret al. (1979, Contributions to Mineralogy and Petrology 70,111–125). Although related, these two schemes have crucialdifferences that can result in very different projection parametersfor the plagioclase component. This is most crucial for high-Kcompositions; H₂O contents estimated using the Walker et al.(1979) scheme can be as much as 100% higher than those estimatedusing the Baker & Eggler (1983) projection. Incorrect projectionparameter calculation has led to overestimation of water contentsin high-K andesites and dacites from the Central Andes. Second,for medium-K and high-K andesitic–dacitic compositionswater contents derived using the Merzbacher & Eggler (1984)geohygrometer deviate considerably from water contents estimatedusing other methods. Experimental data from the literature,and our studies of water contents inferred from melt inclusionsand plagioclase–melt equilibrium for dacites from theAltiplano–Puna Volcanic Complex of the Central Andes indicatethat the Merzbacher & Eggler (1984) geohygrometer shouldnot be applied to compositions with K₂O >1·9 wt %,as originally calibrated. KEY WORDS: magmatic volatiles; pre-eruptive water; experimental geohygrometer; high-K suites; melt inclusion     

The Role of the Solar Cycle in the Relationship Between the North Atlantic Oscillation and Northern Hemisphere Surface Temperatures

The North Atlantic Oscillation （NAO） is one of the leading modes of climate variability in the Northern Hemisphere. It has been shown that it clearly relates to changes in meteorological variables, such as surface temperature, at hemispherical scales. However, recent studies have revealed that the NAO spatial pattern also depends upon solar forcing. Therefore, its effects on meteorological variables must vary depending upon this factor. Moreover, it could be that the Sun affects climate through variability patterns, a hypothesis that is the focus of this study. We find that the relationship between the NAO/AO and hemispheric temperature varies depending upon solar activity. The results show a positive significant correlation only when solar activity is high. Also, the results support the idea that solar activity influences tropospheric climate fluctuations in the Northern Hemisphere via the fluctuations of the stratospheric polar vortex .     

CEMPS: A prototype spatial decision support system to aid in planning emergency evacuations

CEMPS is a prototype spatial decision support system which links the topographical support and analysis provided by a geographic information system, ARC/INFO, with the ability to simulate the dynamics of an evacuation process. CEMPS has been designed to enable emergency planners to experiment with different emergency evacuation plans in order to devise a plan which meets their requirements. ARC/INFO is used to parameterize a dynamic simulation with topographical information and to display its results. The prototype runs on a Sun SPARCStation cluster but could be modified to run on other hardware and software.