Validation of ocean model syntheses against hydrography using a new web application

Results are presented from a new web application called OceanDIVA - Ocean Data Intercomparison and Visualization Application. This tool reads hydrographic profiles and ocean model output and presents the data on either depth levels or isotherms for viewing in Google Earth, or as probability density functions (PDFs) of regional model-data misfits. As part of the CLIVAR Global Synthesis and Observations Panel, an intercomparison of water mass properties of various ocean syntheses has been undertaken using OceanDIVA. Analysis of model-data misfits reveals significant differences between the water mass properties of the syntheses, such as the ability to capture mode water properties.

Modelling the global coastal ocean

Shelf and coastal seas are regions of exceptionally high biological productivity, high rates of biogeochemical cycling and immense socio-economic importance. They are, however, poorly represented by the present generation of Earth system models, both in terms of resolution and process representation. Hence, these models cannot be used to elucidate the role of the coastal ocean in global biogeochemical cycles and the effects global change (both direct anthropogenic and climatic) are having on them. Here, we present a system for simulating all the coastal regions around the world (the Global Coastal Ocean Modelling System) in a systematic and practical fashion. It is based on automatically generating multiple nested model domains, using the Proudman Oceanographic Laboratory Coastal Ocean Modelling System coupled to the European Regional Seas Ecosystem Model. Preliminary results from the system are presented. These demonstrate the viability of the concept, and we discuss the prospects for using the system to explore key areas of global change in shelf seas, such as their role in the carbon cycle and climate change effects on fisheries.

A study of the sensitivity of ocean overturning circulation and climate to freshwater input in different regions of the North Atlantic

The Atlantic meridional overturning circulation (AMOC) is an important component of the climate system. Models indicate that the AMOC can be perturbed by freshwater forcing in the North Atlantic. Using an ocean-atmosphere general circulation model, we investigate the dependence of such a perturbation of the AMOC, and the consequent climate change, on the region of freshwater forcing. A wide range of changes in AMOC strength is found after 100 years of freshwater forcing. The largest changes in AMOC strength occur when the regions of deepwater formation in the model are forced directly, although reductions in deepwater formation in one area may be compensated by enhanced formation elsewhere. North Atlantic average surface air temperatures correlate linearly with the AMOC decline, but warming may occur in localised regions, notably over Greenland and where deepwater formation is enhanced. This brings into question the representativeness of temperature changes inferred from Greenland ice-core records.

On the climate response of the low-latitude Pacific Ocean to changes in the global freshwater cycle

Under global warming, the predicted intensification of the global freshwater cycle will modify the net freshwater flux at the ocean surface. Since the freshwater flux maintains ocean salinity structures, changes to the density-driven ocean circulation are likely. A modified ocean circulation could further alter the climate, potentially allowing rapid changes, as seen in the past. The relevant feedback mechanisms and timescales are poorly understood in detail, however, especially at low latitudes where the effects of salinity are relatively subtle. In an attempt to resolve some of these outstanding issues, we present an investigation of the climate response of the low-latitude Pacific region to changes in freshwater forcing. Initiated from the present-day thermohaline structure, a control run of a coupled ocean-atmosphere general circulation model is compared with a perturbation run in which the net freshwater flux is prescribed to be zero over the ocean. Such an extreme experiment helps to elucidate the general adjustment mechanisms and their timescales. The atmospheric greenhouse gas concentrations are held constant, and we restrict our attention to the adjustment of the upper 1,000 m of the Pacific Ocean between 40°N and 40°S, over 100 years. In the perturbation run, changes to the surface buoyancy, near-surface vertical mixing and mixed-layer depth are established within 1 year. Subsequently, relative to the control run, the surface of the low-latitude Pacific Ocean in the perturbation run warms by an average of 0.6°C, and the interior cools by up to 1.1°C, after a few decades. This vertical re-arrangement of the ocean heat content is shown to be achieved by a gradual shutdown of the heat flux due to isopycnal (i.e. along surfaces of constant density) mixing, the vertical component of which is downwards at low latitudes. This heat transfer depends crucially upon the existence of density-compensating temperature and salinity gradients on isopycnal surfaces. The timescale of the thermal changes in the perturbation run is therefore set by the timescale for the decay of isopycnal salinity gradients in response to the eliminated freshwater forcing, which we demonstrate to be around 10-20 years. Such isopycnal heat flux changes may play a role in the response of the low-latitude climate to a future accelerated freshwater cycle. Specifically, the mechanism appears to represent a weak negative sea surface temperature feedback, which we speculate might partially shield from view the anthropogenically-forced global warming signal at low latitudes. Furthermore, since the surface freshwater flux is shown to play a role in determining the ocean's thermal structure, it follows that evaporation and/or precipitation biases in general circulation models are likely to cause sea surface temperature biases.

Climate entropy budget of the HadCM3 atmosphere-ocean general circulation model and of FAMOUS, its low resolution version

The entropy budget is calculated of the coupled atmosphere–ocean general circulation model HadCM3. Estimates of the different entropy sources and sinks of the climate system are obtained directly from the diabatic heating terms, and an approximate estimate of the planetary entropy production is also provided. The rate of material entropy production of the climate system is found to be ∼50 mW m−2 K−1, a value intermediate in the range 30–70 mW m−2 K−1 previously reported from different models. The largest part of this is due to sensible and latent heat transport (∼38 mW m−2 K−1). Another 13 mW m−2 K−1 is due to dissipation of kinetic energy in the atmosphere by friction and Reynolds stresses. Numerical entropy production in the atmosphere dynamical core is found to be about 0.7 mW m−2 K−1. The material entropy production within the ocean due to turbulent mixing is ∼1 mW m−2 K−1, a very small contribution to the material entropy production of the climate system. The rate of change of entropy of the model climate system is about 1 mW m−2 K−1 or less, which is comparable with the typical size of the fluctuations of the entropy sources due to interannual variability, and a more accurate closure of the budget than achieved by previous analyses. Results are similar for FAMOUS, which has a lower spatial resolution but similar formulation to HadCM3, while more substantial differences are found with respect to other models, suggesting that the formulation of the model has an important influence on the climate entropy budget. Since this is the first diagnosis of the entropy budget in a climate model of the type and complexity used for projection of twenty-first century climate change, it would be valuable if similar analyses were carried out for other such models.

Modelling the movement and survival of the root-feeding clover weevil, Sitona lepidus, in the root-zone of white clover

White clover (Trifolium repens) is an important pasture legume but is often difficult to sustain in a mixed sward because, among other things, of the damage to roots caused by the soil-dwelling larval stages of S. lepidus. Locating the root nodules on the white clover roots is crucial for the survival of the newly hatched larvae. This paper presents a numerical model to simulate the movement of newly hatched S. lepidus larvae towards the root nodules, guided by a chemical signal released by the nodules. The model is based on the diffusion-chemotaxis equation. Experimental observations showed that the average speed of the larvae remained approximately constant, so the diffusion-chernotaxis model was modified so that the larvae respond only to the gradient direction of the chemical signal but not its magnitude. An individual-based lattice Boltzmann method was used to simulate the movement of individual larvae, and the parameters required for the model were estimated from the measurement of larval movement towards nodules in soil scanned using X-ray microtomography. The model was used to investigate the effects of nodule density, the rate of release of chemical signal, the sensitivity of the larvae to the signal, and the random foraging of the larvae on the movement and subsequent survival of the larvae. The simulations showed that the most significant factors for larval survival were nodule density and the sensitivity of the larvae to the signal. The dependence of larval survival rate on nodule density was well fitted by the Michealis-Menten kinetics. (c) 2005 Elsevier B.V All rights reserved.

Carbon (C-13) and nitrogen (N-15) translocation in a maize-Striga hermonthica association

The translocation of C and N in a maize-Striga hermonthica association was investigated at three rates of nitrogen application in a glasshouse experiment. The objectives were to measure the transfer of C and N from maize to S. hermonthica and to determine whether the amount of N in the growing medium affected the proportions of C and N transferred. Young plants of maize were labelled in a (CO2)-C-13 atmosphere and leaf tips were immersed in ((NH4)-N-15)(2)SO4 Solution. The Striga x N interaction was not significant for any of the responses measured. Total dry matter for infected maize was significantly smaller than for uninfected maize from 43 to 99 days after planting, but N application increased total dry matter at all sampling times. Infected maize plants partitioned 39-45 % of their total dry matter to the roots compared with 28-31 % for Uninfected maize. Dry matter of S. hermonthica was not affected by the rate of N applied. S. hermonthica derived 100 % of its carbon from maize before emergence, decreasing to 22-59 % thereafter; the corresponding values for nitrogen were up to 59 % pre-emergence and Lip to 100 % after emergence. The relative proportions of nitrogen depleted from the host (up to 10 %) were greater than those of carbon (maximum 1.2 %) at all times of sampling after emergence of the parasite. The results show that the parasite was more dependent on the host for nitrogen than for carbon.

Dynamics of potassium leaching on a hillslope grassland soil

There have been only a few studies of potassium (K) losses from grassland systems, and little is known about their dynamics, especially in relation to nitrogen (N) management. A study was performed during the autumn and winter of 1999 and 2000 to understand the effects of N and drainage on the dynamics of K leaching on a hillslope grassland soil in southwestern England. Two N application rates were studied (0 and 280 kg N ha(-1) yr(-1)), both with and without the drainage. Treatments receiving N also received farmyard manure (FM). Higher total K losses and K concentrations in the leachates were found in the N + FM treatments (150 and 185% higher than in 0 N treatments), which were related to K additions in the FM. Drainage reduced K losses by 35% because of an increase in dry matter production and a reduction in overland and preferential flow. The pattern of change in K concentration in the leachates was associated with preferential flow at the beginning of the drainage season and with matrix flow later in winter, and was best described by a double exponential curve. Rainfall intensity and the autumn application of FM were the main determinants of K losses by leaching. The study provided new insights into the relationships between soil hydrology, rainfall, and K leaching and its implications for grassland systems.

Potassium budgets in grassland systems as affected by nitrogen and drainage

The main inputs, outputs and transfers of potassium (K) in soils and swards under typical south west England conditions were determined during 1999/00 and 2000/01 to establish soil and field gate K budgets under different fertilizer nitrogen (N) (0 and 280 kg ha(-1) yr(-1)) and drainage (undrained and drained) treatments. Plots receiving fertilizer N also received farmyard manure (FYM). Potassium soil budgets ranged, on average for the two years, from -5 (+N, drained) to +9 (no N and undrained) kg K ha(-1) yr(-1) and field gate budgets from +23 (+N, drained) to +89 (+N, undrained). The main inputs and outputs to the soil K budgets were fertilizer application (65%) and plant uptake (93%). Animals had a minor effect on K export but a major impact on K recycling. Nitrogen fertilizer application and drainage increased K uptake by the grass and, with it, the efficiency of K used. It also depleted easily available soil K, which could be associated with smaller K losses by leaching.

Factors affecting potassium leaching in different soils

A set of lysimeter based experiments was carried out during 2000/01 to evaluate the impact of soil type and grassland management on potassium (K) leaching. The effects of (1) four soil textures (sand, loam, loam over chalk and clay), (2) grazing and cutting (with farmyard manure application), and (3) K applied as inorganic fertilizer, dairy slurry or a mixture of both sources were tested. Total K losses in the clay soil were more than twice those in the sand soil (13 and 6 kg K ha(-1), respectively) because of the development of preferential flow in the clay soil. They were also greater in the cut treatment than in the grazed treatment (82 and 51 kg K ha(-1), respectively; P less than or equal to0.01), associated with a 63% increase of K concentration in the leachates from the former (6.7 +/- 0.28 and 4.1 +/- 0.22 mg K L-1 for cut and grazed, respectively; P less than or equal to0.01) because of the K input from the farmyard manure. The source of fertilizer did not affect total K losses or the average K concentration in the leachates (P >0.05), but it changed the pattern of these over time.

Determination of the impact of continuous defoliation of Lolium perenne and Trifolium repens on bacterial and fungal community structure in rhizosphere soil

To determine the effects of defoliation on microbial community structure, rhizosphere soil samples were taken pre-, and post-defoliation from the root tip and mature root regions of Trifolium repens L. and Lolium perenne L. Microbial DNA isolated from samples was used to generate polymerase chain reaction-denaturing gradient gel electrophoresis molecular profiles of bacterial and fungal communities. Bacterial plate counts were also obtained. Neither plant species nor defoliation affected the bacterial and fungal community structures in both the root tip and mature root regions, but there were significant differences in the bacterial and fungal community profiles between the two root regions for each plant. Prior to defoliation, there was no difference between plants for bacterial plate counts of soils from the root tip regions; however, counts were greater in the mature root region of L. perenne than T. repens. Bacterial plate counts for T. repens were higher in the root tip than the mature root region. After defoliation, there was no effect of plant type, position along the root or defoliation status on bacterial plate counts, although there were significant increases in bacterial plate counts with time. The results indicate that a general effect existed during maturation in the root regions of each plant, which had a greater impact on microbial community structure than either plant type or the effect of defoliation. In addition there were no generic consequences with regard to microbial populations in the rhizosphere as a response to plant defoliation.

Exudation of alcohol and aldehyde sugars from roots of defoliated Lolium perenne L. grown under sterile conditions

Root exudates were collected over a 27 day period from defoliated and non-defoliated Lolium perenne L. plants grown under sterile conditions in microlysimeters. Eleven individual sugars, including both aldehyde and alcohol sugars, were identified and quantified with a gas chromatograph-mass spectrometer (GC-MS). There was no change in the number of sugars present between 7 and 27 days, but the exudation of alcohol sugars decreased rapidly at about day 12. Xylose and glucose were present in the largest amounts. Defoliation initially increased the total amount of sugars in the exudates, but continuous defoliation reduced total sugar exudation by 16% and induced changes in the exudation patterns of individual sugars. Defoliation enhanced exudation of erythritol, threitol, and xylitol, reduced exudation of glucose and arabitol, but had little effect on the amounts of other sugars exuded. The more complex 6 C, 5 OH aldehyde sugars, especially glucose, showed changes earlier and to a greater extent (17 days), than the 5 C, 4 OH (xylose and ribose) and 6 C 4 OH (fucose) aldehyde groups. These findings confirm the general finding that repeated defoliation reduces the quantity of total sugars exuded, but the pattern of release of individual sugars is complex and variable.