Rapid Lateglacial tree population dynamics and ecosystem changes in the eastern Baltic region

A growing body of evidence implies that the concept of 'treeless tundra' in eastern and northern Europe fails to explain the rapidity of Lateglacial and postglacial tree population dynamics of the region, yet the knowledge of the geographic locations and shifting of tree populations is fragmentary. Pollen, stomata and plant macrofossil stratigraphies from Lake Kurjanovas in the poorly studied eastern Baltic region provide improved knowledge of ranges of north‐eastern European trees during the Lateglacial and subsequent plant population responses to the abrupt climatic changes of the Lateglacial/Holocene transition. The results prove the Lateglacial presence of tree populations (Betula, Pinus and Picea) in the eastern Baltic region. Particularly relevant is the stomatal and plant macrofossil evidence showing the local presence of reproductive Picea populations during the Younger Dryas stadial at 12 900–11 700 cal. a BP, occurring along with Dryas octopetala and arctic herbs, indicating semi‐open vegetation. The spread of Pinus–Betula forest at ca. 14 400 cal. a BP, the rise of Picea at ca. 12 800 cal. a BP and the re‐establishment of Pinus–Betula forest at ca. 11 700 cal. a BP within a span of centuries further suggest strikingly rapid, climate‐driven ecosystem changes rather than gradual plant succession on a newly deglaciated land. Copyright © 2009 John Wiley & Sons, Ltd.     

The second order characteristic condition for plastic flow

In this paper the second order characteristic (discontinuous bifurcation) condition is derived for the granular flow (fully plastic) equations. This second order bifurcation equation is shown to be formally identical to the first order localization requirement during steady elastoplastic deformation provided the elastic compliance tensor is substituted for the product of the plastic multiplier with the flow Hessian. For isotropic yield and flow functions the invariant form of the characteristic condition is given in detail, as well as an alternative expression in adapted co‐ordinates. The characteristic condition can be regarded as defining a hardening function which is maximized to identify the critical angles. When the method is applied to 3D Coulomb flow, Mohr's 3D fracture plane conditions are obtained uniquely. Copyright © 2003 John Wiley & Sons, Ltd.     

River sediments,great floods and centennial‐scale Holocene climate change

A new analysis of all 346 published ¹⁴C dated Holocene alluvial units in Britain offers a unique insight into the regional impacts of global change and shows how surprisingly sensitive British rivers have been to relatively modest but repeated changes in climate. Fourteen major but probably brief periods of flooding are identified bracketed within the periods 400–1070, 1940–3940, 7520–8100 and at ca. 10 420 cal. yr BP. There is a strong correspondence between climatic deteriorations inferred from mire wet shifts and major periods of flooding, especially at ca. 8000 cal. yr BP and since ca. 4000 cal. yr BP. The unusually long and complete British record also demonstrates that alterations in land cover have resulted in a step change in river basin sensitivity to variations in climate. This has very important implications for assessing and mitigating the impact of increasing severe flooding. In small and medium‐sized river basins land use is likely to play a key role in either moderating or amplifying the climatic signal. Copyright © 2003 John Wiley & Sons, Ltd.     

塔中4号油藏渗透率分布的非均质性及渗流方向研究

研究油藏内部油水渗流方向时,主要考虑渗透率影响,引入两个水平方向渗透率比、水平与垂向渗透率比参数来分析。前者反映样品所在储层平面上两个水平方向渗流能力的差异,后者反映样品所在储层水平与垂向渗流能力的差异,并经小柱样测试和全岩心测试两种方法确定两个参数。通过对塔中4号油藏CⅠ油组9口井2306个样品两种渗透率比值分布直方图分析,研究了CⅠ油组渗透率的非均质性及其渗流方向:CI油组储层在水平方向的渗流能力有很大差异,具有明显的方向性。储层垂向渗流能力明显小于水平方向的渗流能力。     

A new stereo‐analytical method for determination of removal blocks in discontinuous rock masses

The paper provides a new stereo‐analytical method, which is a combination of the stereographic method and analytical methods, to separate finite removable blocks from the infinite and tapered blocks in discontinuous rock masses. The methodology has applicability to both convex and concave blocks. Application of the methodology is illustrated through examples. Addition of this method to the existing block theory procedures available in the literature improves the capability of block theory in solving practical problems in rock engineering. Copyright © 2003 John Wiley & Sons, Ltd.     

Infiltration at the Tai rain forest (Ivory Coast): Measurements and modelling

Infiltration experiments have been performed at three sites along a well-known catena under virgin tropical rain forest using a portable sprinkling infiltrometer. Experimentally determined infiltration curves are presented. Infiltration curves are also simulated on the basis of the Mein-Larson equation. The parameters for this model have been obtained from the infiltration curves (saturated conductivity) and simple soil moisture determinations (fillable porosity). The agreement between experimentally determined and modelled infiltration is reasonable, provided (a) saturated conductivity as derived from the experimental data is corrected, (b) a storage parameter, also derived from the experimental data, is added to the Mein-Larson model, and (c) the decline in soil porosity with depth is either small or occurs abruptly at shallow depth. Comparison of observed infiltration rates with rainfall intensity shows that Horton Overland Flow has to occur naturally at least on the middle and lower section of the catena. Despite the fact that most parameters can be estimated in principle from basic soil data, it remains advisable to obtain sprinkling infiltrometer field measurements, because of soil variability due to dynamic surface conditions, macroporosity, air entrapment, and irregularity of the wetting front.     

Hysteresis‐based analysis of overland metal transport

Introducing a concept of equivalent mass depth of flow, this study describes the phenomenon of non‐point source pollutant (metal) transport for pavement (or overland) flow in analogy with wave propagation in wide open channels. Hysteretic and normal mass rating curves are developed for runoff rate and mass of 12 dissolved and particulate‐bound metal elements (pollutants) using the rainfall‐runoff and water quality data of the 15 × 20 m² instrumented pavement in Cincinnati, USA. Normal mass rating curves developed for easy computation of pollutant load are found to be of a form similar to Manning's, which is valid for open channel flows. Based on the hysteresis analysis, wave types for dissolution and mixing of particulate‐bound metals are identified. The analysis finds that the second‐order partial‐differential equation normally used for metal transport does not have the efficacy to describe fully the strong non‐linear phenomena such as is described for various metal elements by dynamic waves. In addition, the proportionality concept of the popular SCS‐CN concept is extended for determining the potential maximum metal mass M_p of all the 12 elements transported by a rain storm and related to the antecedent dry period (ADP). For the primary metal zinc element, M_p is found to increase with the ADP and vice versa. 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.     

Fast iterative solution of large undrained soil‐structure interaction problems

In view of rapid developments in iterative solvers, it is timely to re‐examine the merits of using mixed formulation for incompressible problems. This paper presents extensive numerical studies to compare the accuracy of undrained solutions resulting from the standard displacement formulation with a penalty term and the two‐field mixed formulation. The standard displacement and two‐field mixed formulations are solved using both direct and iterative approaches to assess if it is cost‐effective to achieve more accurate solutions. Numerical studies of a simple footing problem show that the mixed formulation is able to solve the incompressible problem ‘exactly’, does not create pressure and stress instabilities, and obviate the need for an ad hoc penalty number. In addition, for large‐scale problems where it is not possible to perform direct solutions entirely within available random access memory, it turns out that the larger system of equations from mixed formulation also can be solved much more efficiently than the smaller system of equations arising from standard formulation by using the symmetric quasi‐minimal residual (SQMR) method with the generalized Jacobi (GJ) preconditioner. Iterative solution by SQMR with GJ preconditioning also is more elegant, faster, and more accurate than the popular Uzawa method. Copyright © 2003 John Wiley & Sons, Ltd.     

Analytical solutions of one‐dimensional infiltration before and after ponding

This paper is a continuation of previous research, which obtained a convenient solution for arbitrary surface fluxes before ponding. By means of Fourier Transformation this has been extended to derive analytical solutions of a linearized Richards' equation for arbitrary input fluxes after surface saturation. Copyright © 2003 John Wiley & Sons, Ltd.