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.     

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.     

Restricted Three-Body Problem: An Approximation of its General Solution – Part One – the Manifold of Symmetric Periodic Solutions

This paper gives the results of a programme attempting to exploit ‘la seule bréche’ (Poincaré, 1892, p. 82) of non-integrable systems, namely to develop an approximate general solution for the three out of its four component-solutions of the planar restricted three-body problem. This is accomplished by computing a large number of families of ‘solutions précieuses’ (periodic solutions) covering densely the space of initial conditions of this problem. More specifically, we calculated numerically and only for μ = 0.4, all families of symmetric periodic solutions (1st component of the general solution) existing in the domain D:(x ₀ ∊ [−2,2],C ∊ [−2,5]) of the (x ₀, C) space and consisting of symmetric solutions re-entering after 1 up to 50 revolutions (see graph in Fig. 4). Then we tested the parts of the domain D that is void of such families and established that they belong to the category of escape motions (2nd component of the general solution). The approximation of the 3rd component (asymmetric solutions) we shall present in a future publication. The 4th component of the general solution of the problem, namely the one consisting of the bounded non-periodic solutions, is considered as approximated by those of the 1st or the 2nd component on account of the `Last Geometric Theorem of Poincaré' (Birkhoff, 1913). The results obtained provoked interest to repeat the same work inside the larger closed domain D:(x ₀ ∊ [−6,2], C ∊ [−5,5]) and the results are presented in Fig. 15. A test run of the programme developed led to reproduction of the results presented by Hénon (1965) with better accuracy and many additional families not included in the sited paper. Pointer directions construed from the main body of results led to the definition of useful concepts of the basic family of order n, n = 1, 2,… and the completeness criterion of the solution inside a compact sub-domain of the (x ₀, C) space. The same results inspired the ‘partition theorem’, which conjectures the possibility of partitioning an initial conditions domain D into a finite set of sub-domains D i that fulfill the completeness criterion and allow complete approximation of the general solution of this problem by computing a relatively small number of family curves. The numerical results of this project include a large number of families that were computed in detail covering their natural termination, the morphology, and stability of their member solutions. Zooming into sub-domains of D permitted clear presentation of the families of symmetric solutions contained in them. Such zooming was made for various values of the parameter N, which defines the re-entrance revolutions number, which was selected to be from 50 to 500. The areas generating escape solutions have being investigated. In Appendix A we present families of symmetric solutions terminating at asymptotic solutions, and in Appendix B the morphology of large period symmetric solutions though examples of orbits that re-enter after from 8 to 500 revolutions. The paper concludes that approximations of the general solution of the planar restricted problem is possible and presents such approximations, only for some sub-domains that fulfill the completeness criterion, on the basis of sufficiently large number of families.     

Small-Amplitude Disturbances in a Radiating and Scattering Grey MediumIII. Gravitational Effects on the Solutions of Given Real Wave Number k

We study the fundamental modes of radiation hydrodynamic linear waves that arise from one-dimensional small-amplitude initial fluctuations with wave number k in a radiating and scattering grey medium by taking into account the gravitational effects. The equation of radiative acoustics is derived from three hydrodynamic equations, Poisson’s equation, and two moment equations of radiation, by assuming a spherical symmetry for the matter and radiation and by using the Eddington approximation. We solve the dispersion relation as a quintic function of angular frequency ω, the wave number k being a real parameter. Numerical results reveal that wave patterns of five solutions are distinguished into three types: the radiation-dominated, type 1, and type 2 matter-dominated cases. In the case of no gravitaional effects (Kaneko et al., 2005), the following wave modes appear: radiation wave, conservative radiation wave, entropy wave, Newtonian-cooling wave, opacity-damped and cooling-damped waves, constant-volume and constant-pressure diffusions, adiabatic sound wave, cooling-damped and drag-force-damped isothermal sound waves, isentropic radiation-acoustic wave, and gap mode. Meanwhile, the gravitaional effects being taken into account, the growing gravo-diffusion mode newly arises from the constant-pressure diffusion at the point that k agrees with Jeans’ wave number specified by the isothermal sound speed. This mode changes to the growing radiation-acoustic gravity mode near the point that k becomes Jeans’ wave number specified by the isentropic radiation-acoustic speed. In step with a transition between them, the isentropic radiation-acoustic wave splits into the damping radiation-acoustic gravity mode and constant-volume diffusion. The constant-volume diffusion emerges twice if the gravitational effects are taken into account. Since analytic solutions are derived for all wave modes, we discuss their physical significance. The critical conditions are given which distinguish between radiation-dominated and type 1 matter-dominated cases, and between type 1 and type 2 matter-dominated cases. Waves in a self-gravitating scattering grey medium are also analyzed, which provides us some hints for the effects of energy and momentum exchange between matter and radiation.     

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.     

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.