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.     

Numerical modelling of the impact crater depth-diameter dependence in an acoustically fluidized target

2D numerical modelling of impact cratering has been utilized to quantify an important depth-diameter relationship for different crater morphologies, simple and complex. It is generally accepted that the final crater shape is the result of a gravity-driven collapse of the transient crater, which is formed immediately after the impact. Numerical models allow a quantification of the formation of simple craters, which are bowl-shaped depressions with a lens of rock debris inside, and complex craters, which are characterized by a structural uplift. The computation of the cratering process starts with the first contact of the impactor and the planetary surface and ends with the morphology of the final crater. Using different rheological models for the sub-crater rocks, we quantify the influence on crater mechanics. To explain the formation of complex craters in accordance to the threshold diameter between simple and complex craters, we utilize the Acoustic Fluidization model. We carried out a series of simulations over a broad parameter range with the goal to fit the observed depth/diameter relationships as well as the observed threshold diameters on the Moon, Earth and Venus.     

Exposure dating and reinterpretation of coarse debris accumulations (‘rock glaciers’) in the Cairngorm Mountains,Scotland

Relict rock glaciers have considerable potential for contributing to palaeoclimatic reconstruction, but this potential is often undermined by lack of dating control and problems of interpretation. Here we reinvestigate and date four proposed ‘rock glaciers’ in the Cairngorm Mountains and show that the morphology of only one of these appears consistent with that of a true rock glacier produced by creep of underlying ice or ice‐rich sediment. All four features comprise rockslide or rock avalanche runout debris, and the possibility that all four represent unmodified runout accumulations cannot be discounted. Surface exposure dating of the four debris accumulations using cosmogenic ¹⁰Be produced uncertainty‐weighted mean ages of 15.4 ± 0.8 ka, 16.2 ± 1.0 ka, 12.1 ± 0.6 ka and 12.7 ± 0.8 ka. All four ages imply emplacement under cold stadial conditions, two prior to the Windermere Interstade of ca. 14.5–12.9 cal. ka BP and two during the Loch Lomond Stade of ca. 12.9–11.5 cal. ka BP. The above ages indicate that paraglacial rock‐slope failure on granite rockwalls occurred within a few millennia after deglaciation. The mean exposure ages obtained for runout debris at two sites – Strath Nethy (16.2 ± 1.0 ka) and Lairig Ghru (15.4 ± 0.8 ka) – are consistent with basal radiocarbon ages from Loch Etteridge, 22 km to the southwest (mean = 15.6 ± 0.3 cal. ka BP) and imply widespread deglaciation of the Cairngorms and adjacent valleys before 15 ka and possibly 16 ka. Copyright © 2008 John Wiley & Sons, Ltd.     

Biostratigraphic and aminostratigraphic constraints on the age of the Middle Pleistocene glacial succession in north Norfolk,UK

Considerable debate surrounds the age of the Middle Pleistocene glacial succession in East Anglia following some recent stratigraphical reinterpretations. Resolution of the stratigraphy here is important since it not only concerns the glacial history of the region but also has a bearing on our understanding of the earliest human occupation of north‐western Europe. The orthodox consensus that all the tills were emplaced during the Anglian (Marine Isotope Stage (MIS) 12) has recently been challenged by a view assigning each major till to a different glacial stage, before, during and after MIS 12. Between Trimingham and Sidestrand on the north Norfolk coast, datable organic sediments occur immediately below and above the glacial succession. The oldest glacial deposit (Happisburgh Till) directly overlies the ‘Sidestrand Unio‐bed’, here defined as the Sidestrand Hall Member of the Cromer Forest‐bed Formation. Dating of these sediments therefore has a bearing on the maximum age of the glacial sequence. This paper reviews the palaeobotany and describes the faunal assemblages recovered from the Sidestrand Unio‐bed, which accumulated in a fluvial environment in a fully temperate climate with regional deciduous woodland. There are indications from the ostracods for weakly brackish conditions. Significant differences are apparent between the Sidestrand assemblages and those from West Runton, the type site of the Cromerian Stage. These differences do not result from contrasting facies or taphonomy but reflect warmer palaeotemperatures at Sidestrand and a much younger age. This conclusion is suggested by the higher proportion of thermophiles at Sidestrand and the occurrence of a water vole with unrooted molars (Arvicola) rather than its ancestor Mimomys savini with rooted molars. Amino acid racemisation data also indicate that Sidestrand is significantly younger than West Runton. These data further highlight the stratigraphical complexity of the ‘Cromerian Complex’ and support the conventional view that the Happisburgh Till was emplaced during the Anglian rather than the recently advanced view that it dates from MIS 16. Moreover, new evidence from the Trimingham lake bed (Sidestrand Cliff Formation) above the youngest glacial outwash sediments (Briton's Lane Formation) indicates that they also accumulated during a Middle Pleistocene interglacial – probably MIS 11. All of this evidence is consistent with a short chronology placing the glacial deposits within MIS 12, rather than invoking multiple episodes of glaciation envisaged in the ‘new glacial stratigraphy’ during MIS 16, 12, 10 and 6. Copyright © 2009 John Wiley & Sons, Ltd.     

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.     

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.     

Ice‐cover variability on shallow lakes at high latitudes: model simulations and observations

A one‐dimensional thermodynamic model for simulating lake‐ice phenology is presented and evaluated. The model can be driven with observed daily or hourly atmospheric forcing of air temperature, relative humidity, wind speed, cloud amount and snowfall. In addition to computing the energy balance components, key model output includes the temperature profile at an arbitrary number of levels within the ice/snow (or the water temperature if there is no ice) and ice thickness (clear ice and snow‐ice) on a daily basis, as well as freeze‐up and break‐up dates. The lake‐ice model is used to simulate ice‐growth processes on shallow lakes in arctic, sub‐arctic, and high‐boreal forest environments. Model output is compared with field and remote sensing observations gathered over several ice seasons. Simulated ice thickness, including snow‐ice formation, compares favourably with field measurements. Ice‐on and ice‐off dates are also well simulated when compared with field and satellite observations, with a mean absolute difference of 2 days. Model simulations and observations illustrate the key role that snow cover plays on the seasonal evolution of ice thickness and the timing of spring break‐up. It is also shown that lake morphometry, depth in particular, is a determinant of ice‐off dates for shallow lakes at high latitudes. Copyright © 2003 John Wiley & Sons, Ltd.     

A micromechanical finite element model for linear and damage‐coupled viscoelastic behaviour of asphalt mixture

This study presents a finite element (FE) micromechanical modelling approach for the simulation of linear and damage‐coupled viscoelastic behaviour of asphalt mixture. Asphalt mixture is a composite material of graded aggregates bound with mastic (asphalt and fine aggregates). The microstructural model of asphalt mixture incorporates an equivalent lattice network structure whereby intergranular load transfer is simulated through an effective asphalt mastic zone. The finite element model integrates the ABAQUS user material subroutine with continuum elements for the effective asphalt mastic and rigid body elements for each aggregate. A unified approach is proposed using Schapery non‐linear viscoelastic model for the rate‐independent and rate‐dependent damage behaviour. A finite element incremental algorithm with a recursive relationship for three‐dimensional (3D) linear and damage‐coupled viscoelastic behaviour is developed. This algorithm is used in a 3D user‐defined material model for the asphalt mastic to predict global linear and damage‐coupled viscoelastic behaviour of asphalt mixture. For linear viscoelastic study, the creep stiffnesses of mastic and asphalt mixture at different temperatures are measured in laboratory. A regression‐fitting method is employed to calibrate generalized Maxwell models with Prony series and generate master stiffness curves for mastic and asphalt mixture. A computational model is developed with image analysis of sectioned surface of a test specimen. The viscoelastic prediction of mixture creep stiffness with the calibrated mastic material parameters is compared with mixture master stiffness curve over a reduced time period. In regard to damage‐coupled viscoelastic behaviour, cyclic loading responses of linear and rate‐independent damage‐coupled viscoelastic materials are compared. Effects of particular microstructure parameters on the rate‐independent damage‐coupled viscoelastic behaviour are also investigated with finite element simulations of asphalt numerical samples. Further study describes loading rate effects on the asphalt viscoelastic properties and rate‐dependent damage behaviour. Copyright © 2006 John Wiley & Sons, Ltd.