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.     

Tephrochronology and tephrostratigraphy of two Pleistocene continental fossiliferous successions from central Italy

The stratigraphical context of two Middle Pleistocene fossiliferous palaeosols from Central Italy (Abruzzo and Tuscany) have been studied. Small mammals and molluscs occur in both palaeosols, which are covered by tephra layers that were analysed using an interdisciplinary approach. Application of fission‐track dating to apatites separated from the Case Picconetto tephra (Pescara, Abruzzo), yielded an age of 0.48 ± 0.04 Ma, indistinguishable from those previously determined for the Campani Quarry (Lower Valdarno, Tuscany) (0.46 ± 0.05 Ma and 0.48 ± 0.05 Ma). Geochemical and petrographic investigations indicate that these tephra originated from different volcanoes, the Alban Hills Volcanic Complex and the Vico Volcano (Latium) respectively. Small mammal and mollusc assemblages indicate different palaeoclimatic and palaeoenvironmental conditions for the Case Picconetto and Campani Quarry palaeosols. Warm and humid conditions can be inferred for the Campani Quarry site, whereas open and cold conditions can be inferred for Case Picconetto. On the basis of faunal data, fission‐track dates and attribution of tephra to specific volcanic eruptions, we suggest a correlation of these faunas with marine oxygen isotope stage 14 (Case Picconetto) and with marine oxygen isotope stage 11 (Campani Quarry), respectively. Copyright © 2003 John Wiley & Sons, Ltd.     

Chemical alteration of tephra in the depositional environment: theoretical stability modelling

The study of the chemical stability of vitreous material in aqueous media is well‐established. There has to date been little consideration of the implications of variations in the chemical durability of tephra in Quaternary tephrochronology. Chemical alteration can take the form of cationic leaching from the matrix, or complete destruction of the silica network, either of which could constrain the ability to chemically identify distal tephra. Here we apply established models of vitreous durability to the published chemical analyses of a large number of Icelandic tephras in order to predict their relative durabilities under equivalent conditions. This suggests that some important tephras have relatively poor chemical stability, and that rhyolitic tephras are, in general, more stable than basaltic. We conclude that tephras should be expected to show predictable differential chemical stability in the post‐depositional environment. Copyright © 2003 John Wiley & Sons, Ltd.     

A high resolution,multiproxy Late‐glacial record of climate change and intrasystem responses in northwest England

A lacustrine carbonate sequence from Hawes Water, Lancashire, UK, has been studied using stable isotopic, lithological, pollen and mineral magnetic analysis. The data reveal four abrupt climatic oscillations in the Late‐glacial Interstadial leading up to the onset of the Loch Lomond Stadial. The data also point to climatic warming relatively early within the stadial, ca. 12 500 GRIP yr, prior to the onset of the Holocene. The oxygen isotope record is taken as a signature of climate forcing against which the response of the lake‐system can be monitored. By adopting this approach it is revealed that the response of the biological system to the rapid climatic oscillations is non‐linear and primarily a function of the antecedent conditions. A significant end‐Devensian isotopic excursion (A) is matched by only minor changes in the cold‐adapted floras and faunas. During the warmer interstadial, the response of the biological ecosystem (events B–D) is clearly influenced by thresholds: major changes in the catchment vegetation associated with relatively minor oscillations in the isotopic signature. The stratigraphical patterns reveal significant lag effects between the onset of climate deterioration and resulting changes in vegetation. Copyright © 2002 John Wiley & Sons, Ltd.     

The influence of a non‐associated flow rule on the calculation of the factor of safety of soil slopes

This paper presents a method that incorporates a non‐associated flow rule into the limit analysis to investigate the influence of the dilatancy angle on the factor of safety for the slope stability analysis. The proposed method retain's the advantage of the upper bound method, which is simple and has no stress involvement in the calculation of the energy dissipation and the factor of safety. Copyright © 2001 John Wiley & Sons, Ltd.     

THE SIMULATED STORAGE CAPACITY OF FLOOD AND WATERLOGGING IN THE TYPICAL AGRICULTURAL REGION IN JIANGHAN PLAIN

１ＩＮＴＲＯＤＵＣＴＩＯＮＪｉａｎｇｈａｎＰｌａｉｎｉｓｏｎｅｏｆｔｈｅｉｍｐｏｒｔａｎｔｂａｓｅｓｏｆＣｈｉｎａ′ｓａｇｒｉｃｕｌｔｕｒａｌｐｒｏｄｕｃｔｉｖｉｔｙｌｏｃａｔｅｄｉｎｔｈｅｍｉｄｄｌｅＣｈａｎｇｊｉａｎｇＲｉｖｅｒｖａｌｌｅｙｏｆｃｅｎｔｒａｌＣｈｉｎａ．ＴｈｅＪｉａｎｇｈａｎｐｌａｉｎａｐｐｒｏｘｉｍａｔｅｌｙ６．６×１０４ｋｍ２ｉｎａｒｅａａｎｄｏｖｅｒ３０ｍｉｌｌｉｏｎｉｎｐｏｐｕｌａｔｉｏｎｉｓｔｈｅｒｅｎｏｗｎｅｄ“ｃｏｕｎｔｒｙｏｆｆｉｓｈａｎｄｒｉｃｅ＂ｉｎ…     

陕西秦巴山区农业自然资源开发利用

农业自然资源据其在农业利用过程中的地位可分为两类:一类是生成农业自然资源的环境,另一类为可资直接利用的自然生成物(即狭义的农业自然资源)。陕西秦巴山区自然环境的主要特点是:复合山体高峻庞大,地跨暖温带和北亚热带,垂直分异明显。区内有:森林、草场、野生生物、林特产与农作物等多种农业自然资源。据此,提出了本区农业自然资源的开发利用途径。     

Structures in the phase space of a four dimensional symplectic map

We numerically investigate the projections of non periodic orbits in a 4-dimensional (4-D) symplectic map composed of two coupled 2-dimensional (2-D) maps. We describe in detail the structures that are produced in different planes of projection and we find how the morphology of the 4-D orbits is influenced by the features of the 2-D maps as the coupling parameter increases. We give an empirical law that describes this influence.     

AREAL DIFFERENTIATION AND CONSTRUCTION OF OASIS AGRICULTURAL ECOSYSTEM IN THE SHIYANG RIVER BASIN

On the basis of the study on areal differentiation of the natural environment of oasis agriculture ecosystems in the Shiyang River Basin, this paper comparatively analyzes the natural productivities, water economic benefits, production efficiency, ecological stabilities and developmental conditions of the Wuwei Oasis agricultural ecosystem in the middle reaches of the river basin and the Minqin Oasis agricultural ecosystem in the lower reaches. Under a same management level and investment of . material and energy, primary productiveness and economic benefits of the former are higher than those of the latter. Construction directions of Wuwei and Minqin oases should be different in order to alleviate the water- use contradiction between the middle and lower reaches. The construction objective of Wuwei Oasis should be efficient irrigated farming production system and Minqin Oasis should become a mixed forestry-pastoral-farming ecosystem taking ecological protection as its major function.     

Sediment transport in analogue flume models compared with real‐world sedimentary systems: a new look at scaling evolution of sedimentary systems in a flume

Evolution of sedimentary systems at large temporal and spatial scales cannot be scaled down to laboratory dimensions by conventional hydraulic Froude scaling. Therefore, many researchers question the validity of experiments aiming to simulate this evolution. Yet, it has been shown that laboratory experiments yield stratigraphic responses to allocyclic forcing that are remarkably similar to those in real‐world prototypes, hinting at scale independency with strong dependence on boundary conditions but weak dependence on the actual sediment transport dynamics. This paper addresses the dilemma by contrasting sediment transport rules that apply in the laboratory with those that apply in real‐world geological systems. It is demonstrated that the generation of two‐dimensional stratigraphy in a flume can be simulated numerically by the non‐linear diffusion equation. Sediment transport theory is used to demonstrate that only suspension‐dominated meandering rivers should be simulated with linear diffusion. With increasing grain‐size (coarse sand to gravel) and shallowness of river systems, the prediction of long‐term transport must be simulated by non‐linear, slope‐dependent diffusion to allow for increasing transport rates and thus change in stratigraphic style. To point out these differences in stratigraphic style, three stages in infill of accommodation have been defined here: (i) a start‐up stage, when the system is prograding to base level (e.g. the shelf edge) with no sediment flux beyond the base‐level point; (ii) a fill‐up stage, when the system is further aggrading while progressively more sediment is bypassing base level with the progression of the infill; and (iii) a keep‐up stage, when more than 90% of the input is bypassing the base level and less than 10% is used for filling the accommodation. By plotting the rate of change in flux for various degrees of non‐linearity (varying the exponent in the diffusion equation) it was found that the error between model and real‐world prototype is largest for the suspension‐dominated prototypes, although never more than 30% and only at the beginning of the fill‐up stage. The error reduces to only 10% for the non‐linear sandy‐gravelly and gravelly systems. These results are very encouraging and open up ways to calibrate numerical models of sedimentary system evolution by such experiments.