Three dimensional simulation of fluid flow in X‐ray CT images of porous media

A numerical scheme is developed in order to simulate fluid flow in three dimensional (3‐D) microstructures. The governing equations for steady incompressible flow are solved using the semi‐implicit method for pressure‐linked equations (SIMPLE) finite difference scheme within a non‐staggered grid system that represents the 3‐D microstructure. This system allows solving the governing equations using only one computational cell. The numerical scheme is verified through simulating fluid flow in idealized 3‐D microstructures with known closed form solutions for permeability. The numerical factors affecting the solution in terms of convergence and accuracy are also discussed. These factors include the resolution of the analysed microstructure and the truncation criterion. Fluid flow in 2‐D X‐ray computed tomography (CT) images of real porous media microstructure is also simulated using this numerical model. These real microstructures include field cores of asphalt mixes, laboratory linear kneading compactor (LKC) specimens, and laboratory Superpave gyratory compactor (SGC) specimens. The numerical results for the permeability of the real microstructures are compared with the results from closed form solutions. Copyright © 2004 John Wiley & Sons, Ltd.     

Reliability assessment of granular filters in embankment dams

Empirical criteria have been used successfully to design filters of most embankment large dam projects throughout the world. However, these empirical rules are only applicable to a particular range of soils tested in laboratory and do not take into account the variability of the base material and filter particle sizes. In addition, it is widely accepted that the safety of fill dams is mainly dependent on the reliability of their filter performance. The work herein presented consists in a new general method for assessing the probability of fulfilling any empirical filter design criteria accounting for base and filter heterogeneity by means of first‐order reliability methods (FORM), so that reliability indexes and probabilities of fulfilling any particular criteria are obtained. This method will allow engineers to estimate the safety of existing filters in terms of probability of fulfilling their design criteria and might also be used as a decision tool on sampling needs and material size tolerances during construction. In addition, sensitivity analysis makes possible to analyse how reliabilities are influenced by different sources of input data. Finally, in case of a portfolio risk assessment, this method will allow engineers to compare the safety of several existing dams in order to prioritize safety investments and it is expected to be a very useful tool to evaluate probabilities of failure due to internal erosion. Copyright © 2006 John Wiley & Sons, Ltd.     

Recurrent Red-tides in the Southampton Water Estuary Caused by the Phototrophic CiliateMesodinium rubrum

In the Southampton Water estuary (southern England, U.K.), red-tides caused by the planktonic, phototrophic ciliateMesodinium rubrum(=Myrionecta rubra) occur during most summers and sometimes in autumn. These events were investigated in detail between 1985 and 1987 and were characterized by levels of chlorophylla(chl a) of over 100 μg l⁻¹, cell numbers ofM. rubrumof over 1×10³ ml⁻¹, oxygen saturations of around 150%, and depleted numbers of macrozooplankton. Initiation of red-water did not appear to be triggered by irradiance or nutrients, but coincided with an increase in temperature and water column stability. This enhanced stability was promoted by increased surface to bottom gradients of both temperature and salinity, and by reduced mixing during neap tides. Development of red-water was accompanied by removal of most of the dissolved NH⁺₄from the water column, whereas some NO⁻₃persisted, presumably maintained by freshwater input. NO⁻₃and NH⁺₄gradually returned to pre-bloom concentrations as the red-water declined in late summer. Maximal biomass ofM. rubrumappeared to be limited by irradiance, and self-shading probably imposed an upper limit of around 300 mg chl a m⁻²within the water column. At the observed levels of chl a, irradiance values within the population maximum between 1 and 3 m depth were only just of the order (≈15 μmol photons m⁻² s⁻¹) required to balance estimated respiratory demands. Oxygen concentration became undersaturated during the late bloom phase, with minimal values of 20–30% saturation recorded in deeper waters; however, despite this and reduced numbers of macrozooplankton, direct deleterious effects on other organisms were not observed.     

Depositional and post‐depositional history of warm stage deposits at Knocknacran,Co. Monaghan,Ireland: implications for preservation of Irish last interglacial deposits

Organic‐rich deposits, uncovered during overburden removal from mantled gypsum karst at Knocknacran opencast gypsum mine, Co. Monaghan, are the best candidate to date for a last interglacial record in Ireland. The two till and organic‐rich deposits (preserved at different quarry elevations) were emplaced on to a Tertiary dolerite surface during high‐energy flood events and subsequently folded and faulted by movement towards sinkholes in underlying gypsum. Uranium–thorium disequilibrium dating suggests that the organic‐rich deposits in the upper section were hydrologically isolated at ca. 41 ka and those in the lower section at ca. 86 ka. Interpretation of the pollen content, although tentative because of the depositional and post‐depositional history of the material, suggests that the organic material originated in a warm stage possibly warmer than the post‐Eemian interstadials. The unusual setting of preservation may indicate that in situ, last interglacial deposits have generally been removed by erosion in Ireland. Copyright © 2004 John Wiley & Sons, Ltd.     

Simulating methane emission from a Chinese rice field as influenced by fertilizer and water level

Natural and agricultural wetlands are considered to be the major sources of global atmospheric methane (CH₄). A one‐dimensional model was developed to simulate methane emission and used to examine the influence of various physical processes on the rate of methane emission. Three processes involved in the methane emission are implemented in the model: production, reoxidation and transport. Three transport pathways were considered: diffusion across water–air or soil–air interfaces, ebullition and diffusion through plants. These pathways are influenced by soil properties, plant growth, water‐table conditions, temperature and external inputs (e.g. fertilizer). The model was used to examine the seasonal variation of the methane emission at a rice field in Hunan, China, which was observed during a field experiment for consecutive (early and late) rice seasons in 1992. The observed seasonal variations of methane emission, and role of plants in transporting methane to the atmosphere, are captured by the model simulation. Further model applications were conducted to simulate effects of fertilizer and water‐level condition on the methane emission. The results indicate that unfermented organic fertilizer produces a higher methane emission rate than mineral fertilizer. The simulations with treatments of a deep‐water covering and constant moisture reduced the methane emission. The rice field study provides a framework for further development of the model towards simulations based on spatially distributed variables (e.g. water table, soil temperature and vegetation) at a regional scale. Copyright © 2003 John Wiley & Sons, Ltd.     

Estimation of snow water equivalent using microwave radiometry over Arctic first‐year sea ice

The magnitude and spatial distribution of snow on sea ice are both integral components of the ocean–sea‐ice–atmosphere system. Although there exists a number of algorithms to estimate the snow water equivalent (SWE) on terrestrial surfaces, to date there is no precise method to estimate SWE on sea ice. Physical snow properties and in situ microwave radiometry at 19, 37 and 85 GHz, V and H polarization were collected for a 10‐day period over 20 first‐year sea ice sites. We present and compare the in situ physical, electrical and microwave emission properties of snow over smooth Arctic first‐year sea ice for 19 of the 20 sites sampled. Physical processes creating the observed vertical patterns in the physical and electrical properties are discussed. An algorithm is then developed from the relationship between the SWE and the brightness temperature measured at 37 GHz (55°) H polarization and the air temperature. The multiple regression between these variables is able to account for over 90% of the variability in the measured SWE. This algorithm is validated with a small in situ data set collected during the 1999 field experiment. We then compare our data against the NASA snow thickness algorithm, designed as part of the NASA Earth Enterprise Program. The results indicated a lack of agreement between the NASA algorithm and the algorithm developed here. This lack of agreement is attributed to differences in scale between the Special Sensor Microwave/Imager and surface radiometers and to differences in the Antarctic versus Arctic snow physical and electrical properties. Copyright © 2003 John Wiley & Sons, Ltd.     

One hundred years of mountain hydrology in Switzerland by the WSL

In 1903 the Swiss Federal Research Institute WSL started its first forest hydrology measurements with the aim to deliver a sound scientific basis for the implementation of new forest legislation introduced in Switzerland in 1876. This legislation was triggered by several large floods that occurred in Switzerland, for which a major cause was widely seen as the poor condition of forests at that time. Consequently, hydrologic research at WSL first focused on the influence of forests on floods. In the second half of the 20th century, other hydrological issues such as water quality, snow hydrology and sediment transport complemented the hydrologic research at WSL. Some recent results of this work are presented in three papers joining this introductory paper to mark the 100th anniversary of hydrologic research at WSL. Copyright © 2006 John Wiley & Sons, Ltd.     

Sensitivity analysis of a deterministic water temperature model to forest canopy and soil temperature in Catamaran Brook (New Brunswick,Canada)

A coupled deterministic hydrological and water temperature model, CEQUEAU, was modified to include soil temperature and crown closure in its calculation of local advective terms in the heat budget. The modified model was than tested to verify its sensitivity to these modifications. An analysis of the heat budget of a small forested catchment in eastern Canada revealed that the advective term related to interflow plays a significant role in the daily water heat budget, providing on average 28% of the local advective budget (which also includes advective heat terms from surface runoff and groundwater) and nearly 14% of the total heat budget (which includes all radiative terms at the water surface, convection and evaporation, as well as the local advective terms). Relative sensitivity indices (RSIs) were used to verify the impact of the newly introduced parameters and variables. Among them, parameters related to the forest cover (crown closure and leaf area index) have a maximum RSI of ?0·6; i.e. a 100% increase in value produces a 60% decrease in the local advective term. Parameters with the greatest influence are the volume of water contributing to interflow and the amplitude of the net radiative flux at the soil surface, which, if doubled, would double the contribution of the local interflow advective term to the heat budget. Copyright © 2003 John Wiley & Sons, Ltd.     

Modelling the effect of low soil temperatures on transpiration by Scots pine

For ecosystem modelling of the Boreal forest it is important to include processes associated with low soil temperature during spring‐early summer, as these affect the tree water uptake. The COUP model, a physically based SVAT model, was tested with 2 years of soil and snow physical measurements and sap flow measurements in a 70‐year‐old Scots pine stand in the boreal zone of northern Sweden. During the first year the extent and duration of soil frost was manipulated in the field. The model was successful in reproducing the timing of the soil warming after the snowmelt and frost thaw. A delayed soil warming, into the growing season, severely reduced the transpiration. We demonstrated the potential for considerable overestimation of transpiration by the model if the reduction of the trees' capacity to transpire due to low soil temperatures is not taken into account. We also demonstrated that the accumulated effect of aboveground conditions could be included when simulating the relationship between soil temperature and tree water uptake. This improved the estimated transpiration for the control plot and when soil warming was delayed into the growing season. The study illustrates the need of including antecedent conditions on root growth in the model in order to catch these effects on transpiration. The COUP model is a promising tool for predicting transpiration in high‐latitude stands. Copyright © 2006 John Wiley & Sons, Ltd.     

全球海洋与大气界面的熵通量估算

给出并讨论了海洋系统中由于热量交换、物质交换以及外力做功引起的熵通量的表达式;同时,根据新近调整的全球海 气界面的气候平均资料,估计了全球海洋与大气界面气候年平均的熵通量(熵流),这一熵通量主要取决于海 气界面的热量和物质交换的空间不均匀性。计算结果显示,由热量交换分量引起的熵通量对海洋系统的总熵通量起到主要的贡献作用,其值约为-555.6mW·(m2·K)-1;由海表风应力做功引起的熵通量相对较小,约为-0.09mW·(m2·K)-1;由物质交换引起的熵通量最小,仅为-0.02mW·(m2·K)-1。总的来说,海洋系统从外界获得的总熵通量为-555.7mW·(m2·K)-1,这也就意味着在气候平均定态下,海洋系统内部的熵产生在量值上等于系统的熵交换,即为555.7mW·(m2·K)-1。海洋系统的负熵流与其内部的各种不可逆过程引起的熵产生取得平衡,确保了全球海洋系统处于非平衡热力学定态,并维持着海洋系统中各种尺度的时 空有序现象的消长过程。