Evaluating the century C model using long-term fertilizer trials in the Indo-Gangetic Plains, India

The GEFSOC Project developed a system for estimating soil carbon (C) stocks and changes at the national and sub-national scale. As part of the development of the system, the Century ecosystem model was evaluated for its ability to simulate soil organic C (SOC) changes in environmental conditions in the Indo-Gangetic Plains, India (IGP). Two long-term fertilizer trials (LTFT), with all necessary parameters needed to run Century, were used for this purpose: a jute (Corchorus capsularis L.), rice (Oryza sativa L.) and wheat (Triticum aestivum L.) trial at Barrackpore, West Bengal, and a rice-wheat trial at Ludhiana, Punjab. The trials represent two contrasting climates of the IGP, viz. semi-arid, dry with mean annual rainfall (MAR) of < 800 mm and humid with > 1600 turn. Both trials involved several different treatments with different organic and inorganic fertilizer inputs. In general, the model tended to overestimate treatment effects by approximately 15%. At the semi-arid site, modelled data simulated actual data reasonably well for all treatments, with the control and chemical N + farm yard manure showing the best agreement (RMSE = 7). At the humid site, Century performed less well. This could have been due to a range of factors including site history. During the study, Century was calibrated to simulate crop yields for the two sites considered using data from across the Indian IGP. However, further adjustments may improve model performance at these sites and others in the IGP. The availability of more longterm experimental data sets (especially those involving flooded lowland rice and triple cropping systems from the IGP) for testing and validation is critical to the application of the model's predictive capabilities for this area of the Indian sub-continent. (C) 2007 Elsevier B.V. All rights reserved.

Terrestrial implications for the maritime geoarchaeological resource: A view from the Lower Palaeolithic

Stone tools and faunal remains have been recovered from the English Channel and the North Sea through trawling, dredging for aggregates, channel clearance, and coring. These finds highlight the potential for a maritime Lower Palaeolithic archaeological resource. It is proposed here that any Lower Palaeolithic artefacts, faunal remains, and sediments deposited in the maritime zone during dry, low-stand phases were once (and may still be) contextually similar to their counterparts in the terrestrial Lower Palaeolithic records of north-western Europe. Given these similarities, can interpretive models and analytical frameworks developed for terrestrial archaeology be profitably applied to an assessment of the potential value of any maritime resource? The terrestrial geoarchaeological resource for the Lower Palaeolithic is dominated by artefacts and ecofacts that have been fluvially reworked. The spatio-temporal resolution of these data varies from entire river valleys and marine isotope stages to river channel gravel bar surfaces and decadal timescales, thus supporting a variety of questions and approaches. However, the structure of the terrestrial resource also highlights two fundamental limitations in current maritime knowledge that can restrict the application of terrestrial approaches to any potential maritime resource: (i) how have the repetitive transgressions and regressions of the Middle and Late Pleistocene modified the terrace landforms and sediments associated with the river systems of the English Channel and southern North Sea basins?; and (ii) do the surviving submerged terrace landforms and fluvial sedimentary deposits support robust geochronological models, as is the case with the classical terrestrial terrace sequences? This paper highlights potential approaches to these questions, and concludes that the fluvial palaeogeography, Pleistocene fossils, and potential Lower Palaeolithic artefacts of the maritime geoarchaeological resource can be profitably investigated in future as derived, low-resolution data sets, facilitating questions of colonisation, occupation, demography, and material culture.

Marine04 marine radiocarbon age calibration, 0-26 cal kyr BP

New radiocarbon calibration curves, IntCal04 and Marine04, have been constructed and internationally ratified to replace the terrestrial and marine components of IntCal98. The new calibration data sets extend an additional 2000 yr, from 0-26 cal kyr BP (Before Present, 0 cal. BP = AD 1950), and provide much higher resolution, greater precision, and more detailed structure than IntCal98. For the Marine04 curve, dendrochronologically-dated tree-ring samples, converted with a box diffusion model to marine mixed-layer ages, cover the period from 0-10.5 call kyr BR Beyond 10.5 cal kyr BP, high-resolution marine data become available from foraminifera in varved sediments and U/Th-dated corals. The marine records are corrected with site-specific C-14 reservoir age information to provide a single global marine mixed-layer calibration from 10.5-26.0 cal kyr BR A substantial enhancement relative to IntCal98 is the introduction of a random walk model, which takes into account the uncertainty in both the calendar age and the C-14 age to calculate the underlying calibration curve (Buck and Blackwell, this issue). The marine data sets and calibration curve for marine samples from the surface mixed layer (Marine04) are discussed here. The tree-ring data sets, sources of uncertainty, and regional offsets are presented in detail in a companion paper by Reimer et al. (this issue).

Sensitivity analysis of a catchment-scale nitrogen model

There are now considerable expectations that semi-distributed models are useful tools for supporting catchment water quality management. However, insufficient attention has been given to evaluating the uncertainties inherent to this type of model, especially those associated with the spatial disaggregation of the catchment. The Integrated Nitrogen in Catchments model (INCA) is subjected to an extensive regionalised sensitivity analysis in application to the River Kennet, part of the groundwater-dominated upper Thames catchment, UK The main results are: (1) model output was generally insensitive to land-phase parameters, very sensitive to groundwater parameters, including initial conditions, and significantly sensitive to in-river parameters; (2) INCA was able to produce good fits simultaneously to the available flow, nitrate and ammonium in-river data sets; (3) representing parameters as heterogeneous over the catchment (206 calibrated parameters) rather than homogeneous (24 calibrated parameters) produced a significant improvement in fit to nitrate but no significant improvement to flow and caused a deterioration in ammonium performance; (4) the analysis indicated that calibrating the flow-related parameters first, then calibrating the remaining parameters (as opposed to calibrating all parameters together) was not a sensible strategy in this case; (5) even the parameters to which the model output was most sensitive suffered from high uncertainty due to spatial inconsistencies in the estimated optimum values, parameter equifinality and the sampling error associated with the calibration method; (6) soil and groundwater nutrient and flow data are needed to reduce. uncertainty in initial conditions, residence times and nitrogen transformation parameters, and long-term historic data are needed so that key responses to changes in land-use management can be assimilated. The results indicate the general, difficulty of reconciling the questions which catchment nutrient models are expected to answer with typically limited data sets and limited knowledge about suitable model structures. The results demonstrate the importance of analysing semi-distributed model uncertainties prior to model application, and illustrate the value and limitations of using Monte Carlo-based methods for doing so. (c) 2005 Elsevier B.V. All rights reserved.

National and sub-national assessments of soil organic carbon stocks and changes: The GEFSOC modelling system

Soil organic carbon (SOC) plays a vital role in ecosystem function, determining soil fertility, water holding capacity and susceptibility to land degradation. In addition, SOC is related to atmospheric CO, levels with soils having the potential for C release or sequestration, depending on land use, land management and climate. The United Nations Convention on Climate Change and its Kyoto Protocol, and other United Nations Conventions to Combat Desertification and on Biodiversity all recognize the importance of SOC and point to the need for quantification of SOC stocks and changes. An understanding of SOC stocks and changes at the national and regional scale is necessary to further our understanding of the global C cycle, to assess the responses of terrestrial ecosystems to climate change and to aid policy makers in making land use/management decisions. Several studies have considered SOC stocks at the plot scale, but these are site specific and of limited value in making inferences about larger areas. Some studies have used empirical methods to estimate SOC stocks and changes at the regional scale, but such studies are limited in their ability to project future changes, and most have been carried out using temperate data sets. The computational method outlined by the Intergovernmental Panel on Climate Change (IPCC) has been used to estimate SOC stock changes at the regional scale in several studies, including a recent study considering five contrasting eco regions. This 'one step' approach fails to account for the dynamic manner in which SOC changes are likely to occur following changes in land use and land management. A dynamic modelling approach allows estimates to be made in a manner that accounts for the underlying processes leading to SOC change. Ecosystem models, designed for site scale applications can be linked to spatial databases, giving spatially explicit results that allow geographic areas of change in SOC stocks to be identified. Some studies have used variations on this approach to estimate SOC stock changes at the sub-national and national scale for areas of the USA and Europe and at the watershed scale for areas of Mexico and Cuba. However, a need remained for a national and regional scale, spatially explicit system that is generically applicable and can be applied to as wide a range of soil types, climates and land uses as possible. The Global Environment Facility Soil Organic Carbon (GEFSOC) Modelling System was developed in response to this need. The GEFSOC system allows estimates of SOC stocks and changes to be made for diverse conditions, providing essential information for countries wishing to take part in an emerging C market, and bringing us closer to an understanding of the future role of soils in the global C cycle. (C) 2007 Elsevier B.V. All rights reserved.

Moving on. Where should the university map collection be going?

This paper argues for the importance of retaining a map library presence on UK university campuses at a time when many are under threat of closure, and access to geospatial data is increasingly moving to web-based services. It is suggested that the need for local expertise is undiminished and map curators need to redefine themselves as geoinformation specialists, preserving their paper map collections, but also meeting some of the challenges of GIS, and contributing to national developments in the construction of distributed geolibraries and the provision of metadata, especially with regard to local data sets.

Dating the volcanic eruption at Thera

The eruption of the volcano at Thera (Santorini) in the Aegean Sea undoubtedly had a profound influence on the civilizations of the surrounding region. The date of the eruption has been a subject of much controversy because it must be linked into the established and intricate archaeological phasings of both the prehistoric Aegean and the wider east Mediterranean. Radiocarbon dating of material from the volcanic destruction layer itself can provide some evidence for the date of the eruption, but because of the shape of the calibration curve for the relevant period, the value of such dates relies on there being no biases in the data sets. However, by dating the material from phases earlier and later than the eruption, some of the problems of the calibration data set can be circumvented and the chronology for the region can be resolved with more certainty. In this paper, we draw together the evidence we have accumulated so far, including new data on the destruction layer itself and for the preceding cultural horizon at Thera, and from associated layers at Miletos in western Turkey. Using Bayesian models to synthesize the data and to identify outliers, we conclude from the most reliable C-14 evidence (and using the INTCAL98 calibration data set) that the eruption of Thera occurred between 1663 and 1599 BC.

IntCal04 terrestrial radiocarbon age calibration, 0-26 cal kyr BP

A new calibration curve for the conversion of radiocarbon ages to calibrated (cal) ages has been constructed and internationally ratified to replace ImCal98, which extended from 0-24 cal kyr BP (Before Present, 0 cal BP = AD 1950). The new calibration data set for terrestrial samples extends from 0-26 cal kyr BP, but with much higher resolution beyond 11.4 cal kyr BP than ImCal98. Dendrochronologically-dated tree-ring samples cover the period from 0-12.4 cal kyr BP. Beyond the end of the tree rings, data from marine records (corals and foraminifera) are converted to the atmospheric equivalent with a site-specific marine reservoir correction to provide terrestrial calibration from 12.4-26.0 cal kyr BP. A substantial enhancement relative to ImCal98 is the introduction of a coherent statistical approach based on a random walk model, which takes into account the uncertainty in both the calendar age and the C-14 age to calculate the underlying calibration curve (Buck and Blackwell, this issue). The tree-ring data sets, sources of uncertainty, and regional offsets are discussed here. The marine data sets and calibration curve for marine samples from the surface mixed layer (Marine 04) are discussed in brief, but details are presented in Hughen et al. (this issue a). We do not make a recommendation for calibration beyond 26 cal kyr BP at this time; however, potential calibration data sets are compared in another paper (van der Plicht et al., this issue).

Influences of space, soil, nematodes and plants on microbial community composition of chalk grassland soils

Microbial communities respond to a variety of environmental factors related to resources (e.g. plant and soil organic matter), habitat (e.g. soil characteristics) and predation (e.g. nematodes, protozoa and viruses). However, the relative contribution of these factors on microbial community composition is poorly understood. Here, we sampled soils from 30 chalk grassland fields located in three different chalk hill ridges of Southern England, using a spatially explicit sampling scheme. We assessed microbial communities via phospholipid fatty acid (PLFA) analyses and PCR-denaturing gradient gel electrophoresis (DGGE) and measured soil characteristics, as well as nematode and plant community composition. The relative influences of space, soil, vegetation and nematodes on soil microorganisms were contrasted using variation partitioning and path analysis. Results indicate that soil characteristics and plant community composition, representing habitat and resources, shape soil microbial community composition, whereas the influence of nematodes, a potential predation factor, appears to be relatively small. Spatial variation in microbial community structure was detected at broad (between fields) and fine (within fields) scales, suggesting that microbial communities exhibit biogeographic patterns at different scales. Although our analysis included several relevant explanatory data sets, a large part of the variation in microbial communities remained unexplained (up to 92% in some analyses). However, in several analyses, significant parts of the variation in microbial community structure could be explained. The results of this study contribute to our understanding of the relative importance of different environmental and spatial factors in driving the composition of soil-borne microbial communities.

Extraction of tidal channel networks from aerial photographs alone and combined with laser altimetry

Tidal channel networks play an important role in the intertidal zone, exerting substantial control over the hydrodynamics and sediment transport of the region and hence over the evolution of the salt marshes and tidal flats. The study of the morphodynamics of tidal channels is currently an active area of research, and a number of theories have been proposed which require for their validation measurement of channels over extensive areas. Remotely sensed data provide a suitable means for such channel mapping. The paper describes a technique that may be adapted to extract tidal channels from either aerial photographs or LiDAR data separately, or from both types of data used together in a fusion approach. Application of the technique to channel extraction from LiDAR data has been described previously. However, aerial photographs of intertidal zones are much more commonly available than LiDAR data, and most LiDAR flights now involve acquisition of multispectral images to complement the LiDAR data. In view of this, the paper investigates the use of multispectral data for semiautomatic identification of tidal channels, firstly from only aerial photographs or linescanner data, and secondly from fused linescanner and LiDAR data sets. A multi-level, knowledge-based approach is employed. The algorithm based on aerial photography can achieve a useful channel extraction, though may fail to detect some of the smaller channels, partly because the spectral response of parts of the non-channel areas may be similar to that of the channels. The algorithm for channel extraction from fused LiDAR and spectral data gives an increased accuracy, though only slightly higher than that obtained using LiDAR data alone. The results illustrate the difficulty of developing a fully automated method, and justify the semi-automatic approach adopted.

Use of a portable topographic mapping millimetre wave radar at an active lava flow

A ground-based millimetre wave radar, AVTIS (All-weather Volcano Topography Imaging Sensor), has been developed for topographic monitoring. The instrument is portable and capable of measurements over ranges up to similar to 7 km through cloud and at night. In April and May 2005, AVTIS was deployed at Arenal Volcano, Costa Rica, in order to determine topographic changes associated with the advance of a lava flow. This is the first reported application of mm-wave radar technology to the measurement of lava flux rates. Three topographic data sets of the flow were acquired from observation distances of similar to 3 km over an eight day period, during which the flow front was detected to have advanced similar to 200 m. Topographic differences between the data sets indicated a flow thickness of similar to 10 m, and a dense rock equivalent lava flux of similar to 0.20 +/- 0.08 m(3) s(-1).

Assessing the uncertainty associated with intermittent rainfall fields

[1] In many practical situations where spatial rainfall estimates are needed, rainfall occurs as a spatially intermittent phenomenon. An efficient geostatistical method for rainfall estimation in the case of intermittency has previously been published and comprises the estimation of two independent components: a binary random function for modeling the intermittency and a continuous random function that models the rainfall inside the rainy areas. The final rainfall estimates are obtained as the product of the estimates of these two random functions. However the published approach does not contain a method for estimation of uncertainties. The contribution of this paper is the presentation of the indicator maximum likelihood estimator from which the local conditional distribution of the rainfall value at any location may be derived using an ensemble approach. From the conditional distribution, representations of uncertainty such as the estimation variance and confidence intervals can be obtained. An approximation to the variance can be calculated more simply by assuming rainfall intensity is independent of location within the rainy area. The methodology has been validated using simulated and real rainfall data sets. The results of these case studies show good agreement between predicted uncertainties and measured errors obtained from the validation data.

Evaluation of the sea ice simulation in a new coupled atmosphere-ocean climate model (HadGEM1)

[ 1] A rapid increase in the variety, quality, and quantity of observations in polar regions is leading to a significant improvement in the understanding of sea ice dynamic and thermodynamic processes and their representation in global climate models. We assess the simulation of sea ice in the new Hadley Centre Global Environmental Model (HadGEM1) against the latest available observations. The HadGEM1 sea ice component uses elastic-viscous-plastic dynamics, multiple ice thickness categories, and zero-layer thermodynamics. The model evaluation is focused on the mean state of the key variables of ice concentration, thickness, velocity, and albedo. The model shows good agreement with observational data sets. The variability of the ice forced by the North Atlantic Oscillation is also found to agree with observations.

Indian Ocean-monsoon coupled interactions and impending monsoon droughts

Monsoon droughts over the Indian subcontinent emanate from failures in the seasonal (June-September) monsoon rains. While prolonged dry-spells ("monsoon-breaks'') pervade on sub-seasonal/intra-seasonal time-scales, the underlying causes for these long-lasting anomalies remain elusive. Based on analyses of a suite of observed data sets, we report an ocean-atmosphere dynamical coupling on intra-seasonal time-scales, in the tropical Indian Ocean, which is pivotal in forcing extended monsoon-breaks and causing droughts over the subcontinent. This coupling involves a feedback between the monsoonal flow and thermocline depth in the Equatorial Eastern Indian Ocean (EEIO), in which an anomaly of the summer monsoon circulation induces downwelling and maintains a higher-than-normal heat-content. The near-equatorial anomalies induce strong and sustained suppression of monsoon rainfall over the subcontinent. It is concluded that the intra-seasonal evolution of the ocean-monsoon coupled system is a vital key to unlocking the dynamics of monsoon droughts.

Associations between stratospheric variability and tropospheric blocking

There is widely believed to be a link between stratospheric flow variability and stationary, persistent “blocking” weather systems, but the precise nature of this link has proved elusive. Using data from the ERA-40 Reanalysis and an atmospheric general circulation model (GCM) with a well-resolved stratosphere (HadGAM), it is shown that there are in fact several different highly significant associations, with blocking in different regions being related to different patterns of stratospheric variability. This is true in both hemispheres and in both data sets. The associations in HadGAM are shown to be very similar to those in ERA-40, although the model has a tendency to underestimate both European blocking and the wave number 2 stratospheric variability to which this is related. Although the focus is on stratospheric variability in general, several of the blocking links are seen to occur in association with the major stratospheric sudden warmings. In general, the direction of influence appears to be upward, as blocking anomalies are shown to modify the planetary stationary waves, leading to an upward propagation of wave activity into the stratosphere. However, significant correlations are also apparent with the zonal mean flow in the stratosphere leading the occurrence of blocking at high latitudes. Finally, the underestimation of blocking is an enduring problem in GCMs, and an example has recently been given in which improving the resolution of the stratosphere improved the representation of blocking. Here, however, another example is given, in which increasing the stratospheric resolution unfortunately does not lead to an improvement in blocking.