推求土壤水分运动参数的简单入渗法 I.理论分析

预报非饱和土壤水分运动必须首先获得土壤水分运动参数。土壤水分运动参数包括土壤水分特征曲线和导水率。本文使用积分方法求解了一维水平非饱和土壤水分运动问题，根据其解建立了推求非饱和土壤水分运动参数的简单入渗法，用以推求ｖａｎ Ｇｅｎｕｃｈｔｅｎ特征曲线模型中的参数α和ｎ。α和ｎ是根据湿润区的特征长度、吸渗率和土壤的饱和导水率（ｋｓ）来确定的，而非饱和导水率可由α、ｎ和Ｋｓ确定。这一新的简单入渗法是基于Ｒｉｃｈａｒｄｓ方程和土壤导水特征的闭合型方程。简单入渗法提供了利用瞬态水流方法来确定土壤水分特征曲线而替代通常的平衡法。简单入渗法是一个全新的、简捷的确定土壤导水特性的方法。

根系吸水机理研究进展

植物根系除支撑和固定其地上部这一力学功能外的一个重要功能就是从土壤中吸收水分和养分以满足植物地上部生长所需.表征植物根吸水能力的一个重要的水力学参数是水力导度(用单位时间单位面积的水流速率来表示),可在细胞(细胞水力导度)、单根和整株根系水平上来表达,其中单根导度可分为径向导度和轴向导度,仅径向导度反映了单根吸收水分的能力,而轴向导度则反映了植物根系输导水分的能力,但在整株根系水平上则以通过整个根系的水流通量与根木质部和根表土壤间的水势差之比来表示,既包括径向导度也包括轴向导度.近年来,对植物根系吸水或根水力导度研究已取得了许多重要进展,这对阐明根系吸水机理和地上地下部关系起到了重要作用.下面做一简要综述.1　根系吸收水分的主要部位根系吸水的部位一般认为在距根尖10～100ｍｍ的区域内,这种看法是基于解剖上的证据.从根表面到根中心,依次为根表皮、下表皮、皮层、内皮层、中柱,其中根表皮是有最高吸收活性的根区,但一般仅可存活几天,而内皮层将皮层和中柱分开形成了根内侧的一个界面,老根一般有周皮或栓质化的内皮层,有很强的不透水性(即阻力很大).但Ｓａｎｄｅｒｓｏｎ[1]对大麦的研究发现,虽然老根区内重度栓化的内皮层已...

压力室测定根系导水率方法探讨

用压力室连续测定了玉米根系升压和降压过程的导水率。结果表明 ,降压过程测得的根系导水率显著大于用升压过程的 ,并且前者的相关系数大于后者。这种差异是由于这两个过程中质外体途径细胞壁空间充水量不同造成的。开始升压时 ,由于细胞壁空间含水量低 ,质外体途径阻力大 ,导致非结构阻力 ;随着压力的升高 ,细胞壁空间含水量增大 ,质外体途径导度增大 ,减小甚至可以消除非结构阻力。降压法可以使根系快速复水 ,消除传统方法因长时间复水所致根结构的改变。建议用降压法测定根系导水率

细沟股流剥蚀率与载沙量以及沟长的耦合关系

土壤剥蚀是指由侵蚀动力引起的土壤颗粒从土壤母质移动的过程。细沟剥蚀土粒随着细沟股流中含沙量的增加而减少 ,已有的一些侵蚀模型 (如 WEPP)均提到了这一点。用黄土高原一种典型的粉壤土 ,在 5种坡度、3种流量下进行了细沟侵蚀模拟试验。对试验结果进行了回归 ,分析了黄土高原斜坡及陡坡地、细沟股流剥蚀率随含沙量以及沟长变化的函数关系。这对细沟侵蚀动力过程的研究深入 ,以及对侵蚀过程的预测预报提供了有力的参考依据

陡坡细沟含沙水流剥蚀率的试验研究及其计算方法

细沟剥蚀土粒随着细沟股流中含沙量的增加而减少 ,这一概念已在一些侵蚀模型 (如 WEPP)中得到应用。用黄土高原一种典型的粉壤土 ,在 5种坡度 (5°,10°,15°,2 0°,2 5°) ,3种流量 (2 ,4,8L/m in)条件下进行了细沟侵蚀模拟试验 ,试验沟长 0 .5～ 8m。通过 40 5次试验 ,确定了不同坡度、入流量条件下 ,侵蚀产沙量与细沟长度的定量函数关系。在假定细沟径流和土壤侵蚀沿细沟的行为相同条件下 ,提出了一种计算含沙水剥蚀率的方法 ,并进一步表达了细沟剥蚀率随含沙量以及沟长变化的函数关系。实验结果在 15°,2 0°,2 5°时表现出很好的显著性

Observations of Kuroshio intrusion into the South China Sea

To discuss the intrusion of the Kuroshio into the SCS, we examined the mixing between the North Pacific and South China Sea (SCS) waters based on in-situ CTD data collected in August and September 2008 and the moored ADCP data taken from mid September 2008 to early July 2009. The CTD survey included four meridional sections from 119A degrees E to 122A degrees E around the Luzon Strait, during which pressure, temperature, and salinity were measured. The CTD data show that the isopycnal surface tilted from the SCS to the North Pacific; and it was steeper in the lower layers than in the upper ones. Meanwhile, we found strong vertical mixing taken place in the areas near 121A degrees E. The Kuroshio in high temperature and salinity intruded westward through Luzon Strait. The frequency of buoyancy was one order of magnitude greater than that of the common ones in the ocean, suggesting stronger stratification in the northeastern SCS. On the other hand, the long-term ADCP data show that before late October 2008, the direction of water flow in the SCS was eastward, and from November 2008 to late February 2009, it turned northwestward in the layers shallower than 150 m, while remained unchanged in deep layers from 200 to 450 m. From March to June 2009, the direction shifted with increasing depth from northward to southward, akin to the Ekman spiral. EOF analysis of the current time series revealed dominant empirical modes: the first mode corresponded to the mean current and showed that the Kuroshio intrusion occurred in the upper layers only from late December to early March. The temporal coefficient of the first and the second mode indicated clearly a dominant signal in a quasi-seasonal cycle.

地下水流动数值模拟的并行计算方法研究

Numerical modeling of groundwater is very important for understanding groundwater flow and solving hydrogeological problem. Today, groundwater studies require massive model cells and high calculation accuracy, which are beyond single-CPU computer’s capabilities. With the development of high performance parallel computing technologies, application of parallel computing method on numerical modeling of groundwater flow becomes necessary and important. Using parallel computing can improve the ability to resolve various hydro-geological and environmental problems. In this study, parallel computing method on two main types of modern parallel computer architecture, shared memory parallel systems and distributed shared memory parallel systems, are discussed. OpenMP and MPI (PETSc) are both used to parallelize the most widely used groundwater simulator, MODFLOW. Two parallel solvers, P-PCG and P-MODFLOW, were developed for MODFLOW. The parallelized MODFLOW was used to simulate regional groundwater flow in Beishan, Gansu Province, which is a potential high-level radioactive waste geological disposal area in China. 1. The OpenMP programming paradigm was used to parallelize the PCG (preconditioned conjugate-gradient method) solver, which is one of the main solver for MODFLOW. The parallel PCG solver, P-PCG, is verified using an 8-processor computer. Both the impact of compilers and different model domain sizes were considered in the numerical experiments. The largest test model has 1000 columns, 1000 rows and 1000 layers. Based on the timing results, execution times using the P-PCG solver are typically about 1.40 to 5.31 times faster than those using the serial one. In addition, the simulation results are the exact same as the original PCG solver, because the majority of serial codes were not changed. It is worth noting that this parallelizing approach reduces cost in terms of software maintenance because only a single source PCG solver code needs to be maintained in the MODFLOW source tree. 2. P-MODFLOW, a domain decomposition–based model implemented in a parallel computing environment is developed, which allows efficient simulation of a regional-scale groundwater flow. The basic approach partitions a large model domain into any number of sub-domains. Parallel processors are used to solve the model equations within each sub-domain. The use of domain decomposition method to achieve the MODFLOW program distributed shared memory parallel computing system will process the application of MODFLOW be extended to the fleet of the most popular systems, so that a large-scale simulation could take full advantage of hundreds or even thousands parallel processors. P-MODFLOW has a good parallel performance, with the maximum speedup of 18.32 (14 processors). Super linear speedups have been achieved in the parallel tests, indicating the efficiency and scalability of the code. Parallel program design, load balancing and full use of the PETSc were considered to achieve a highly efficient parallel program. 3. The characterization of regional ground water flow system is very important for high-level radioactive waste geological disposal. The Beishan area, located in northwestern Gansu Province, China, is selected as a potential site for disposal repository. The area includes about 80000 km2 and has complicated hydrogeological conditions, which greatly increase the computational effort of regional ground water flow models. In order to reduce computing time, parallel computing scheme was applied to regional ground water flow modeling. Models with over 10 million cells were used to simulate how the faults and different recharge conditions impact regional ground water flow pattern. The results of this study provide regional ground water flow information for the site characterization of the potential high-level radioactive waste disposal.

Recursos hidricos no Brasil e no mundo.

2001

Generalized semi-geostrophic theory on a sphere

Douglas, Robert; Cullen, M.J.P.; Roulston, I.; Sewell, M.J., (2005) 'Generalized semi-geostrophic theory on a sphere', Journal of Fluid Mechanics 531 pp.123-157 RAE2008

Clergy, crusified within the body of Christ - a challenge to pastoral care

M.A. Thesis / University of Pretoria / Department of Practical Theology / Advised by Prof M Masango

A Marine Environmental Recorder

An inexpensive Marine Environmental Recorder is described. The instrument system is small, lightweight and of low-power consumption. It is flexible, simple to operate and economical. It can be used remotely in a moored, buoyed or towed instrument system, recording measurements continuously for up to 24 h, or intermittently for 1 min every hour, for a period of up to 60 d. It has been used extensively in the Continuous Plankton Recorder and the Undulating Oceanographic Recorder to measure temperature, depth and occasionally chlorophyll and radiant energy; as a temperature recorder, it has a resolution of 0.1 Co, an uncertainty of measurement of ±0.1 Co and a stability of calibration within ±0.1 Co over a period of several months. With optional additional sensors for pitch, roll, vibration, acceleration and water-flow, the instrument system has been used to measure the performance of underwater towed vehicles and plankton samplers. The Marine Environmental Recorder is being incorporated into an instrument system in a data buoy, for automatically monitoring the marine environment in estuaries around the British Isles.

Assessing the sensitivity of blue mussel beds to pressures associated with human activities.

The Joint Nature Conservation Committee (JNCC) commissioned this project to generate an improved understanding of the sensitivities of blue mussel (Mytilus edulis) beds, found in UK waters, to pressures associated with human activities in the marine environment. The work will provide an evidence base that will facilitate and support management advice for Marine Protected Areas, development of UK marine monitoring and assessment, and conservation advice to offshore marine industries. Blue mussel beds are identified as a Habitat of Principle Importance (HPI) under the Natural Environment and Rural Communities (NERC) Act 2006, as a Priority Marine Feature (PMF) under the Marine (Scotland) Act 2010, and included on the OSPAR (Annex V) list of threatened and declining species and habitats. The purpose of this project was to produce sensitivity assessments for the blue mussel biotopes included within the HPI, PMF and OSPAR habitat definitions, and clearly document the supporting evidence behind the assessments and any differences between them. A total of 20 pressures falling in five categories - biological, hydrological, physical damage, physical loss, and pollution and other chemical changes - were assessed in this report. The review examined seven blue mussel bed biotopes found on littoral sediment and sublittoral rock and sediment. The assessments were based on the sensitivity of M. edulis rather than associated species, as M. edulis was considered the most important characteristic species in blue mussel beds. To develop each sensitivity assessment, the resistance and resilience of the key elements are assessed against the pressure benchmark using the available evidence gathered in this review. The benchmarks were designed to provide a ‘standard’ level of pressure against which to assess sensitivity. Blue mussel beds were highly sensitive to a few human activities: • introduction or spread of non-indigenous species (NIS); • habitat structure changes - removal of substratum (extraction); and • physical loss (to land or freshwater habitat). Physical loss of habitat and removal of substratum are particularly damaging pressures, while the sensitivity of blue mussel beds to non-indigenous species depended on the species assessed. Crepidula fornicata and Crassostrea gigas both had the potential to outcompete and replace mussel beds, so resulted in a high sensitivity assessment. Mytilus spp. populations are considered to have a strong ability to recover from environmental disturbance. A good annual recruitment may allow a bed to recovery rapidly, though this cannot always be expected due to the sporadic nature of M. edulis recruitment. Therefore, blue mussel beds were considered to have a 'Medium' resilience (recovery within 2-10 years). As a result, even where the removal or loss of proportion of a mussel bed was expected due to a pressure, a sensitivity of 'Medium' was reported. Hence, most of the sensitivities reported were 'Medium'. It was noted, however, that the recovery rates of blue mussel beds were reported to be anywhere between two years to several decades. In addition, M. edulis is considered very tolerant of a range of physical and chemical conditions. As a result, blue mussel beds were considered to be 'Not sensitive' to changes in temperature, salinity, de-oxygenation, nutrient and organic enrichment, and substratum type, at the benchmark level of pressure. The report found that no distinct differences in overall sensitivity exist between the HPI, PMF and OSPAR definitions. Individual biotopes do however have different sensitivities to pressures, and the OSPAR definition only includes blue mussel beds on sediment. These differences were determined by the position of the habitat on the shore and the sediment type. For example, the infralittoral rock biotope (A3.361) was unlikely to be exposed to pressures that affect sediments. However in the case of increased water flow, mixed sediment biotopes were considered more stable and ‘Not sensitive’ (at the benchmark level) while the remaining biotopes were likely to be affected.

Using a clearly documented, evidence-based approach to create sensitivity assessments allows the assessment basis and any subsequent decision making or management plans to be readily communicated, transparent and justifiable. The assessments can be replicated and updated where new evidence becomes available ensuring the longevity of the sensitivity assessment tool. For every pressure where sensitivity was previously assessed as a range of scores in MB0102, the assessments made by the evidence review have supported one of the MB0102 assessments. The evidence review has reduced the uncertainty around assessments previously undertaken in the MB0102 project (Tillin et al., 2010) by assigning a single sensitivity score to the pressures as opposed to a range. Finally, as blue mussel bed habitats also contribute to ecosystem function and the delivery of ecosystem services, understanding the sensitivity of these biotopes may also support assessment and management in regard to these. Whatever objective measures are applied to data to assess sensitivity, the final sensitivity assessment is indicative. The evidence, the benchmarks, the confidence in the assessments and the limitations of the process, require a sense-check by experienced marine ecologists before the outcome is used in management decisions.

Anatomy and facies association of a drumlin in Co. Down, Northern Ireland, from seismic and electrical resistivity surveys

Seismic refraction and electrical resistivity geophysical techniques were used to reconstruct the internal architecture of a drumlin in Co. Down, Northern Ireland. Geophysical results were both validated and complemented by borehole drilling, ground water flow modelling, and geologic mapping. The geophysical anatomy of the drumlin consists of five successive layers with depth including; topsoil, partially saturated and saturated glacial tills, and weathered and more competent greywacke bedrock. There are numerous, often extensive inclusions of clay, sand, gravel, cobbles, and boulders within the topsoil and the till units. Together geophysical and geotechnical findings imply that the drumlin is part of the subglacial lodgement, melt-out, debris flow, sheet flow facies described by previous authors, and formed by re-sedimentation and streamlining of pre-existing sediments during deglaciation of the Late Devensian ice sheet. Seismic refraction imaging is particularly well suited to delineating layering within the drumlin, and is able to reconstruct depths to interfaces to within ± 0.5 m accuracy. Refraction imaging ascertained that the weathered bedrock layer is continuous and of substantial thickness, so that it acts as a basal aquifer which underdrains the bulk of the drumlin. Electrical resistivity imaging was found to be capable of delineating relative spatial changes in the moisture content of the till units, as well as mapping sedimentary inclusions within the till. The moisture content appeared to be elevated near the margins of the drumlin, which may infer a weakening of the drumlin slopes. Our findings advocate the use of seismic refraction and electrical resistivity methods in future sedimentological and geotechnical studies of internal drumlin architecture and drumlin formation, owing particularly to the superior, 3- D spatial coverage of these methods.

Comparative tracing experiments to investigate epikarst structural and compositional heterogeneity

Comparative tracer testing may be used to evaluate the vulnerability of groundwater to specific contaminants by comparing reactive tracer response to that of a simultaneously injected non-reactive “conservative” substance. Conversely, knowledge of tracer reaction with specific materials permits information about subsurface heterogeneity to be inferred. A series of tests completed in the vadose zone overlying a limestone aquifer employed a cocktail of particles along with reactive and non-reactive solute tracers to investigate transport rates between the ground surface and monitoring points approximately 10 m below ground. Short pulse tests revealed both solutes and particulate contaminants could travel at rates of over 10 m/h. Comparison of particle (microorganisms) and non-reactive solute tracer breakthrough revealed that particle tracers experience pore exclusion resulting in higher peak relative concentrations which arrive earlier than those of the solute. Prolonged tracer injection during subsequent experiments confirmed the response observed and illustrated that over 40 % of flow paths between injection and monitoring points were inaccessible to particles, but could allow solutes to pass through them. Similarly, the difference in response between various reactive tracers demonstrated tracers reached monitoring points via multiple flow paths and suggests geochemical heterogeneity plays an important role in influencing tracer behaviour. The results of this investigation highlight the complexity of water flow through the epikarst and the vulnerability of groundwater in karst aquifers to contamination when soil cover is thin to absent.

Measurement of permeability using a bench-top centrifuge

The commonly used British Standard constant head triaxial permeability test for testing of fine-grained soils is relatively time consuming. A reduction in the required time for soil permeability testing would provide potential cost savings to the construction industry, particularly in the construction quality assurance of landfill clay liners. The purpose of this paper is to evaluate an alternative approach of measuring permeability of fine-grained soils benefiting from accelerated time scaling for seepage flow when testing specimens in elevated gravity conditions provided by a centrifuge. As part of the investigation, an apparatus was designed and produced to measure water flow through soil samples under conditions of elevated gravitational acceleration using a small desktop laboratory centrifuge. A membrane was used to hydrostatically confine the test sample. A miniature data acquisition system was designed and incorporated in the apparatus to monitor and record changes in head and flow throughout the tests. Under enhanced gravity in the centrifuge, the flow through the sample was under ‘variable head' conditions as opposed to ‘constant head' conditions as in the classic constant head permeability tests conducted at 1 g . A mathematical model was developed for analysis of Darcy's coefficient of permeability under conditions of elevated gravitational acceleration and verified using the results obtained. The test data compare well with the results on analogous samples obtained using the classical British Standard constant head permeability tests.