Modelling nitrogen fluxes from the land surface to the coastal zone in European systems: the perspective of the INCA project

This paper describes the results and conclusions of the INCA (Integrated Nitrogen Model for European CAtchments) project and sets the findings in the context of the ELOISE (European Land-Ocean Interaction Studies) programme. The INCA project was concerned with the development of a generic model of the major factors and processes controlling nitrogen dynamics in European river systems, thereby providing a tool (a) to aid the scientific understanding of nitrogen transport and retention in catchments and (b) for river-basin management and policy-making. The findings of the study highlight the heterogeneity of the factors and processes controlling nitrogen dynamics in freshwater systems. Nonetheless, the INCA model was able to simulate the in-stream nitrogen concentrations and fluxes observed at annual and seasonal timescales in Arctic, Continental and Maritime-Temperate regimes. This result suggests that the data requirements and structural complexity of the INCA model are appropriate to simulate nitrogen fluxes across a wide range of European freshwater environments. This is a major requirement for the production of coupled fiver-estuary-coastal shelf models for the management of our aquatic environment. With regard to river-basin management, to achieve an efficient reduction in nutrient fluxes from the land to the estuarine and coastal zone, the model simulations suggest that management options must be adaptable to the prevailing environmental and socio-economic factors in individual catchments: 'Blanket approaches' to environmental policy appear too simple. (c) 2004 Elsevier B.V. All rights reserved.

Dynamics of QoI sensitivity in Mycosphaerella fijiensis in Costa Rica during 2000 to 2003

From 1997 onward, the strobilurin fungicide azoxystrobin was widely used in the main banana-production zone in Costa Rica against Mycosphaerella fijiensis var. difformis causing black Sigatoka of banana. By 2000, isolates of M. fijiensis with resistance to the quinolene oxidase inhibitor fungicides were common on some farms in the area. The cause was a single point mutation from glycine to alanine in the fungal target protein, cytochrome b gene. An amplification refractory mutation system Scorpion quantitative polymerase chain reaction assay was developed and used to determine the frequency of G 143A allele in samples of M. fijiensis. Two hierarchical surveys of spatial variability, in 2001 and 2002,found no significant variation in frequency on spatial scales <10 in. This allowed the frequency of G143A alleles on a farm to be estimated efficiently by averaging single samples taken at two fixed locations. The frequency of G 143A allele in bulk samples from I I farms throughout Costa Rica was determined at 2-month intervals. There was no direct relationship between the number of spray applications and the frequency of G143A on individual farms. Instead, the frequency converged toward regional averages, presumably due to the large-scale mixing of ascospores dispersed by wind. Using trap plants in an area remote from the main producing area, immigration of resistant ascospores was detected as far as 6 km away both with and against the prevailing wind.

Ultra-distal tephra deposits from super-eruptions: examples from Toba, Indonesia and Taupo Volcanic Zone, New Zealand

Voluminous rhyolitic eruptions from Toba, Indonesia, and Taupo Volcanic Zone (TVZ), New Zealand, have dispersed volcanic ash over vast areas in the late Quaternary. The ~74 ka Youngest Toba Tuff (YTT) eruption deposited ash over the Bay of Bengal and the Indian subcontinent to the west. The ~340 ka Whakamaru eruption (TVZ) deposited the widespread Rangitawa Tephra, dominantly to the southeast (in addition to occurrences northwest of vent), extending across the landmass of New Zealand, and the South Pacific Ocean and Tasman Sea, with distal terrestrial exposures on the Chatham Islands. These super-eruptions involved ~2500 km^3 and ~1500 km3 of magma (dense-rock equivalent; DRE), respectively. Ultra-distal terrestrial exposures of YTT at two localities in India, Middle Son Valley, Madhya Pradesh, and Jurreru River Valley, Andhra Pradesh, at distances of >2000 km from the source caldera, show a basal ‘primary’ ashfall unit ~4 cm thick, although deposits containing reworked ash are up to ~3 m in total thickness. Exposures of Rangitawa Tephra on the Chatham Islands, >900 km from the source caldera, are ~15-30 cm thick. At more proximal localities (~200 km from source), Rangitawa Tephra is ~55-70 cm thick and characterized by a crystal-rich basal layer and normal grading. Both distal tephra deposits are characterized by very-fine ash (with high PM10 fractions) and are crystal-poor. Glass chemistry, stratigraphy and grain-size data for these distal tephra deposits are presented with comparisons of their correlation, dispersal and preservation. Using field observations, ash transport and deposition were modeled for both eruptions using a semi-analytical model (HAZMAP), with assumptions concerning average wind direction and strength during eruption, column shape and vent size. Model outputs provide new insights into eruption dynamics and better estimates of eruption volumes associ- ated with tephra fallout. Modeling based on observed YTT distal tephra thicknesses indicate a relatively low (<40 km high), very turbulent eruption column, consistent with deposition from a co-ignimbrite cloud extending over a broad region. Similarly, the Whakamaru eruption was modeled as producing a predominantly Plinian column (~45 km high), with dispersal to the southeast by strong prevailing winds. Significant ash fallout of the main dispersal direction, to the northwest of source, cannot be replicated in this modeling. The widespread dispersal of large volumes of fine ash from both eruptions may have had global environmental consequences, acutely affecting areas up to thousands of kilometers from vent.

Stress field control of eruption dynamics at a rift volcano: Nyamuragira, D.R.Congo

Using the record of 30 flank eruptions over the last 110 years at Nyamuragira, we have tested the relationship between the eruption dynamics and the local stress field. There are two groups of eruptions based on their duration (< 80days >) that are also clustered in space and time. We find that the eruptions fed by dykes parallel to the East African Rift Valley have longer durations (and larger volumes) than those eruptions fed by dykes with other orientations. This is compatible with a model for compressible magma transported through an elastic-walled dyke in a differential stress field from an over-pressured reservoir (Woods et al., 2006). The observed pattern of eruptive fissures is consistent with a local stress field modified by a northwest-trending, right lateral slip fault that is part of the northern transfer zone of the Kivu Basin rift segment. We have also re-tested with new data the stochastic eruption models for Nyamuragira of Burt et al. (1994). The time-predictable, pressure-threshold model remains the best fit and is consistent with the typically observed declining rate of sulphur dioxide emission during the first few days of eruption with lava emission from a depressurising, closed, crustal reservoir. The 2.4-fold increase in long-term eruption rate that occurred after 1977 is confirmed in the new analysis. Since that change, the record has been dominated by short-duration eruptions fed by dykes perpendicular to the Rift. We suggest that the intrusion of a major dyke during the 1977 volcano-tectonic event at neighbouring Nyiragongo volcano inhibited subsequent dyke formation on the southern flanks of Nyamuragira and this may also have resulted in more dykes reaching the surface elsewhere. Thus that sudden change in output was a result of a changed stress field that forced more of the deep magma supply to the surface. Another volcano-tectonic event in 2002 may also have changed the magma output rate at Nyamuragira.

Adaptation-induced collective dynamics of a single-cell protozoan

We investigate the behavior of a single-cell protozoan in a narrow tubular ring. This environment forces them to swim under a one-dimensional periodic boundary condition. Above a critical density, single-cell protozoa aggregate spontaneously without external stimulation. The high-density zone of swimming cells exhibits a characteristic collective dynamics including translation and boundary fluctuation. We analyzed the velocity distribution and turn rate of swimming cells and found that the regulation of the turing rate leads to a stable aggregation and that acceleration of velocity triggers instability of aggregation. These two opposing effects may help to explain the spontaneous dynamics of collective behavior. We also propose a stochastic model for the mechanism underlying the collective behavior of swimming cells.

The ocean’s role in setting the mean position of the Inter-Tropical Convergence Zone

Through study of observations and coupled climate simulations, it is argued that the mean position of the Inter-Tropical Convergence Zone (ITCZ) north of the equator is a consequence of a northwards heat transport across the equator by ocean circulation. Observations suggest that the hemispheric net radiative forcing of climate at the top of the atmosphere is almost perfectly symmetric about the equator, and so the total (atmosphere plus ocean) heat transport across the equator is small (order 0.2 PW northwards). Due to the Atlantic ocean’s meridional overturning circulation, however, the ocean carries significantly more heat northwards across the equator (order 0.4 PW) than does the coupled system. There are two primary consequences. First, atmospheric heat transport is southwards across the equator to compensate (0.2 PW southwards), resulting in the ITCZ being displaced north of the equator. Second, the atmosphere, and indeed the ocean, is slightly warmer (by perhaps 2 °C) in the northern hemisphere than in the southern hemisphere. This leads to the northern hemisphere emitting slightly more outgoing longwave radiation than the southern hemisphere by virtue of its relative warmth, supporting the small northward heat transport by the coupled system across the equator. To conclude, the coupled nature of the problem is illustrated through study of atmosphere–ocean–ice simulations in the idealized setting of an aquaplanet, resolving the key processes at work.

Late Quaternary hydrological dynamics in the Middle Kalahari: forcing and feedbacks

The Kalahari region has become a major source of Quaternary palaeoenvironmental data derived primarily from the analysis of geomorphological proxies of environmental change. One suite of data, from palaeolacustrine landforms, has recently provided a new record of major hydrological changes in the last 150 ka [Burrough, S. L., Thomas, D. S. G., Bailey, R. M., 2009. Mega-Lake in the Kalahari: A Late Pleistocene record of the Palaeolake Makgadikgadi system. Quaternary Science Reviews, in press.]. Here we present an improved analysis of the drivers and feedbacks of lake level change, utilising information from three main sources: data from the lake system itself, from analyses of other late Quaternary records within the region and from climate modelling. Simulations using the Hadley Centre coupled climate model, HadCM3, suggest that once triggered, the lake body was large enough to potentially affect both local and regional climates. Surface waters and their interactions with the climate are therefore an important component of environmental dynamics during the late Quaternary. Through its capacity to couple Middle Kalahari environments to distant forcing mechanisms and to itself force environmental change, we demonstrate that the existence or absence of megalake Makgadikgadi adds a new level of complexity to the interpretations of environmental proxy records in southern Africa's summer rainfall zone.

Granular flow in the marginal ice zone

The region of sea ice near the edge of the sea ice pack is known as the marginal ice zone (MIZ), and its dynamics are complicated by ocean wave interaction with the ice cover, strong gradients in the atmosphere and ocean and variations in sea ice rheology. This paper focuses on the role of sea ice rheology in determining the dynamics of the MIZ. Here, sea ice is treated as a granular material with a composite rheology describing collisional ice floe interaction and plastic interaction. The collisional component of sea ice rheology depends upon the granular temperature, a measure of the kinetic energy of flow fluctuations. A simplified model of the MIZ is introduced consisting of the along and across momentum balance of the sea ice and the balance equation of fluctuation kinetic energy. The steady solution of these equations is found to leading order using elementary methods. This reveals a concentrated region of rapid ice flow parallel to the ice edge, which is in accordance with field observations, and previously called the ice jet. Previous explanations of the ice jet relied upon the existence of ocean currents beneath the ice cover. We show that an ice jet results as a natural consequence of the granular nature of sea ice.

Geophysical characterization of the complex dynamics of groundwater and seawater exchange in a highly stressed aquifer system linked to a coastal lagoon (SE Spain)

Anthropogenic pressure influences the two-way interactions between shallow aquifers and coastal lagoons. Aquifer overexploitation may lead to seawater intrusion, and aquifer recharge from rainfall plus irrigation may, in turn, increase the groundwater discharge into the lagoon. We analyse the evolution, since the 1950s up to the present, of the interactions between the Campo de Cartagena Quaternary aquifer and the Mar Menor coastal lagoon (SE Spain). This is a very heterogeneous and anisotropic detrital aquifer, where aquifer–lagoon interface has a very irregular geometry. Using electrical resistivity tomography, we clearly identified the freshwater–saltwater transition zone and detected areas affected by seawater intrusion. Severity of the intrusion was spatially variable and significantly related to the density of irrigation wells in 1950s–1960s, suggesting the role of groundwater overexploitation. We distinguish two different mechanisms by which water from the sea invades the land: (a) horizontal advance of the interface due to a wide exploitation area and (b) vertical rise (upconing) caused by local intensive pumping. In general, shallow parts of the geophysical profiles show higher electrical resistivity associated with freshwater mainly coming from irrigation return flows, with water resources mostly from deep confined aquifers and imported from Tagus river, 400 km north. This indicates a likely reversal of the former seawater intrusion process.

Interacting controls on innate sources of CO2 efflux from a calcareous arid zone soil under experimental acidification and wetting

More than half of global soil carbon is stored as carbonates, primarily in arid and semi-arid zones. Climate change models predict more frequent and severe rainfall events in some parts of the globe, many of which are dominated by calcareous soils. Such events trigger substantial increases in soil CO2 efflux. We hypothesised that the primary source of CO2 emissions from calcareous, arid zone soil during a single wetting event is abiotic and that soil acidification and wetting have a positive, potentially interacting, effect. We manipulated soil pH, soil moisture, and controlled soil respiration by gamma irradiating half of an 11 day incubation experiment. All manipulated experimental treatments had a rapid and enormous effect on CO2 emission. Respiration contributed ca. 5% of total CO2 efflux; the major source (carbonate buffering) varied depending on the extent of acidification and wetting. Maximum CO2 efflux occurred when pH was lowest and at intermediate matric potential. CO2 efflux was lowest at native pH when soil was air dry. Our data suggest that there may be an underestimate of soil-atmosphere carbon fluxes in arid ecosystems with calcareous soils. There is also a clear potential that these soils may become net carbon sources depending on changes in rainfall patterns, rainfall acidity, and future land management. Our findings have major implications for carbon cycling in arid zone soil and further study of carbon dynamics in these terrestrial systems at a landscape level will be required if we are to improve global climate and carbon cycling models.