A methodology to estimate renewal time scales in estuaries: the Tagus Estuary case

High-resolution hydrodynamic models are a common tool to simulate water dynamics in estuaries. Results from these models are, however, difficult to interpret without the aid of additional parameters to integrate the information. In this paper a methodology to understand the transport patterns in the Tagus Estuary is proposed. It is based on the computation of two renewal time scales: residence time and integrated water fraction. This last parameter is used to build a dependency matrix that gives the integrated influence of each region of the estuary at a selected point. The parameters are computed using a Lagrangian transport model coupled to the hydrodynamic model. Results show that Tagus Estuary has two different types of regions: the central part of the estuary, with low renewal efficiency, and three regions with higher renewal efficiency. Renewal mechanisms are, however, different for each region as shown by the dependency matrix. Comparison of renewal time scales with results from a water-quality model revealed that residence time is not a limiting parameter for primary production in the Tagus Estuary.Responsible Editor: Hans Burchard     

Mineral-scale Sr-isotope constraints on magma evolution and chamber dynamics in the Rum layered intrusion,Scotland

Sr isotopic zoning within single plagioclase crystals from rocks from Unit 9 of the Rum layered intrusion is used to infer events during crystal growth in a magma undergoing contamination. The ⁸⁷Sr/⁸⁶Sr diversity among minerals and between cores and rims of plagioclase crystals increases as the boundary between unit 9 and the overlying Unit 10 peridotite is approached. Models of near-solidus interaction of the cumulate with a fluid or melt, or large scale textural re-equilibration, cannot easily account for the systematic differences in ⁸⁷Sr/⁸⁶Sr between small crystals and the rims of larger crystals.We propose a simple interpretation in which crystal growth is concentrated along the cool margins of the reservoir. Crystals are subsequently advected to a site of accumulation at the base of the reservoir, probably by episodic plume-like dense downwellings allowing mixing of isotopically zoned and unzoned crystals.If the core-rim isotope variations are inherited from primary magmatic growth, then the small distances over which they are now preserved (1–2 mm) place constraints on the minimum cooling rate of the intrusion. Although the length scale of diffusive equilibration is influenced by a number of poorly-constrained variables (starting temperature, feldspar composition, temperature-time path) cooling was clearly very rapid with cooling to effective closure (~1,000 °C) within a few thousand years.Editorial responsibility: I. Parsons     

Regional-scale hydrothermal alteration in the Central Blake River Group,western Abitibi subprovince,Canada: implications for VMS prospectivity

The Late Archean Blake River Group is a thick succession of predominantly mafic volcanic rocks within the southern zone of the Abitibi greenstone belt. It contains a number of silicic volcanic centers of different size, including the large Noranda volcanic complex, which is host to 17 past-producing volcanogenic massive sulfide deposits. The Noranda complex consists of a 7- to 9-km-thick succession of bimodal mafic and felsic volcanic rocks erupted during five major cycles of volcanism. Massive sulfide formation coincided with a period of intense magmatic activity (cycle III) and the formation of the Noranda cauldron. Hydrothermal alteration in these rocks is interpreted to reflect large-scale hydrothermal fluid flow associated with rapid crustal extension and rifting of the volcanic complex. The alteration includes abundant albite, chlorite, epidote and quartz (silicification), which exhibit broad stratigraphic and structural control and correlate with previously mapped whole-rock oxygen isotope zonation. The Mine Sequence volcanic rocks are characterized by abundant iron-rich chlorite (Fe/Fe+Mg >0.5), hydrothermal amphibole (ferroactinolite) and coarse-grained epidote of clinozoisite composition (<10 wt% Fe ₂O ₃). Volcanic rocks of the pre-cauldron sequences, which contain only subeconomic stringer mineralization, are characterized by less abundant chlorite and mainly fine-grained epidote (>10 wt% Fe ₂O ₃) lacking the clinozoisite solid solution. Alteration in the Mine Sequence volcanic rocks persists along strike well beyond the limits of the main ore deposits (as far as several tens of kilometers) and can be readily distinguished from greenschist facies metamorphic assemblages at a regional scale. The lack of similar alteration in the pre-cauldron sequences is consistent with limited ¹⁸O-depletion and suggests that the early history of the volcanic complex did not support large-scale, high-temperature fluid flow in these rocks. Comparisons with a much smaller, barren volcanic complex in nearby Ben Nevis township reveal important differences in the alteration mineralogy between volcanoes of different size, with implications for area selection during regional-scale mineral exploration. The Ben Nevis Complex consists of a 3- to 4-km-thick succession of mafic, intermediate and felsic volcanic rocks centered on a small subvolcanic intrusion. Alteration of the volcanic rocks comprises mainly low-temperature assemblages of prehnite, pumpellyite, magnesium-rich chlorite (Fe/Fe+Mg <0.5), iron-rich epidote (>10 wt% Fe ₂O ₃) and calcite. Actinolite ± magnetite alteration occurs proximal to the intrusive core of the complex, but the limited extent of this alteration indicates only local high-temperature fluid circulation adjacent to the intrusion. A distal zone of carbonate alteration is located 4–6 km from the center of the volcano. Although iron-bearing carbonates are present locally within this zone, the absence of siderite argues against a high-temperature origin for this alteration. These observations do not offer positive encouragement for the existence of a fossil geothermal system of sufficient size or intensity to have produced a large massive sulfide deposit.     

A review of mortality maps included in the annual reports on vital statistics for Ireland, 1864-1998

Since 1864, some 135 reports have been produced examining marriages, birth and deaths for Ireland. The reports produced during the period prior to World War I are notable for the maps and charts that are included and the level of detail provided. Recent reports have failed to include mortality maps, while those produced since 1950 generally show a declining level of detail.     

An empirical garnet (YAG) – xenotime thermometer

 A pronounced negative correlation between the yttrium concentration in garnet ([Y]_Grt) and temperature has been observed in xenotime (YPO₄)-bearing metapelites from central New England, USA. The [Y]_Grt decreases roughly two orders of magnitude (∼5500 to less than 100 ppm Y) over a 150 °C interval. A regression of ln([Y]_Grt) against estimated reciprocal temperature yields the following relationship:

Evaporation over a Heterogeneous Land Surface: EVA_GRIPS and the LITFASS-2003 Experiment—An Overview

The Evaporation at Grid/Pixel Scale (EVA_GRIPS) project was realised in order to determine the area-averaged evaporation over a heterogeneous land surface at the scale of a grid box of a regional numerical weather prediction or climate model, and at the scale of a pixel of a satellite image. EVA_GRIPS combined surface-based and airborne measurements, satellite data analysis, and numerical modelling activities. A mesoscale field experiment, LITFASS-2003, was carried out in the heterogeneous landscape around the Meteorological Observatory Lindenberg (MOL) of the German Meteorological Service in May and June, 2003. The experiment was embedded in the comprehensive, operational measurement program of the MOL. Experimental determination of surface fluxes on a variety of spatial scales was achieved by employing micrometeorological flux stations, scintillometers, a combination of ground-based remote sensing instruments, and the Helipod, a turbulence probe carried by a helicopter. Surface energy fluxes were also derived from satellite data. Modelling work included the use of different Soil–Vegetation–Atmosphere Transfer schemes, a large-eddy simulation model and three mesoscale atmospheric models. The paper gives an overview on the background of EVA_GRIPS, and on the measurements and meteorological conditions during LITFASS-2003. A few general results are discussed.     

Investigating the amount of zeroing in modern sediments of River Danube, Austria

Sediments attributed to flooding events of River Danube concerning the bleaching of the optical stimulated luminescence (OSL) signal were investigated. It is demonstrated that the OSL signal in both quartz and feldspar is not completely but differentially bleached in the sediment grains. Partial bleaching of the samples is clearly indicated by the scatter of equivalent dose determined for several individual single aliquots. It is also shown that residual ages in feldspars are significantly higher than those calculated for quartz. It is furthermore demonstrated that analysing measurement-time dependent equivalent dose estimates is not a suitable method to identify partial bleaching in the investigated sediment grains. However, the transport and deposition process of the investigated samples was probably disturbed by an artificial input of sediment, and this case study may thus not be representative of undisturbed high flood events in the past.     

Carbon dioxide in igneous petrogenesis: II. Fluid dynamics of mantle metasomatism

Petrographic, fluid inclusion, geochemical and isotopic evidence from xenoliths in alkali basalts suggests that low-viscosity fluids rich in O-H-C, dissolved silicates and especially the incompatible elements may ascend, decompress and precipitate crystalline phases and/or induce partial fusion in the upper mantle. Such mantle metasomatic fluids (MMF) may be important in generating isotopic heterogeneity and in transporting and focusing mantle heat. In order to model the movement of MMF, the ordinary differential equations governing the variation ofP, T, ascent velocity and fluid density of a compressible, viscous, single-phase (H₂O or CO₂) non-reacting fluid ascending through a vertical crack of constant width have been solved. A large number of numerical simulations were carried out in which the significant factors affecting flow behavior (thermodynamic and transport fluid properties, roughness and width of cracks, geothermal gradient, initial conditions, etc.) were systematically varied. The calculations show that: (1) MMF tends to move at uniform rates following a short period of rapid initial acceleration, (2) MMF ascends nearly isothermally, (3) MMF acts as an efficient heat transfer agent; numerical experiments show that transport of heat into regions undergoing metasomatism can lead to partial fusion. The heat transported by movement of MMF averaged over the age of the Earth is sufficient to generate about 0.1 km³ of basaltic magma per year, which is approximately equal to the production rate of alkaline magma. If an intense period of mantle degassing occured early in the history of the Earth, the transport of heat and mass (K, U, Rb, LREE) by migrating fluids might have been important.     

Second order elasticity theory: An improved formulation of the Grüneisen parameter at high pressure

The second order theory of elasticity, in which terms to second order in strain are retained in calculating atomic bond length changes and elastic moduli, is extended to describe thermal vibration of a face-centred cubic crystal. Coupled with equations relating the pressure dependences of elastic constants, this yields a new formulation of the thermal Grüneisen parameter γ in terms of pressure P, incompressibility K and rigidity, μ

$\text{λ}\text{1}\text{2}\text{d}\text{K}\text{d}\text{P}\text{?}\text{1}\text{2}\text{+}\text{1}\text{9}\text{P}\text{K}\text{?}\text{1}\text{3}\text{?}\text{1}\text{9}\text{P}\text{K}\text{f}\text{1?}\text{2}\text{3}\text{P}\text{K}\text{?}\text{2}\text{3}\text{P}\text{K}\text{f}$

where f = 24 (3 K ? 2 P)/(3 K + 115 μ + 90 P). The factor f arises from bond interactions and the case f = 1, representing independent bonds (no interactions), yields the free-volume γ- Since we have shown earlier that the second order elasticity theory provides a convincing explanation of the elasticity of the inner core, we believe that the new formula is appropriate for the inner core. It is, however, inadequate to describe the lower mantle γ, in which atomic bond angle rigidity, not considered here, may be appreciable.