Nd isotopic variations of Phanerozoic paleoceans

The size of various Phanerozoic paleoceans, determined from paleogeographic maps, provides a framework for assigningε_Nd measurements of primary authigenic marine minerals to individual paleoceans and for establishing theirε_Nd versus age curves. We report newε_Nd data (in the range from 80 to 750 Ma) for the largest Phanerozoic and Late Precambrian paleocean, the Pacific/Panthalassa Ocean (PPO), and limited data on some small Paleozoic paleoceans. These new data and all other published data on paleoceanε_Nd values are used to determine the trends through time ofε_Nd in seven Phanerozoic paleoceans. The variations in paleoceanε_Nd curves are due to changes in the rate of crustal additions and in paleoceanic geometry. From these data we can calculate a meanε_Nd curve and a mean Nd model age curve for Phanerozoic seawater. These curves show large variations with age and reflect the variations ofε_Nd in erosion products from the continents as a function of age. The high erosion rates of young, uplifted areas enhanceε_Nd and Nd model age variations of Phanerozoic seawater. This makes it difficult to quantitatively recover crustal addition rates from these curves; however, it makes the episodic nature of continental additions more evident.     

A crystal-chemical basis for Pb retention and fission-track annealing systematics in U-bearing minerals, with implications for geochronology

This study develops an empirical crystal-chemical framework for systematizing the kinetics of Pb loss and fission-track annealing in U-bearing minerals. Ionic porosity, Z (the fraction of a mineral's unit-cell volume not occupied by ions) potentially accounts for kinetic behavior by monitoring mean metal-oxygen bond length/strength. Various tests of a general kinetics-porosity relationship are presented, based upon diverse mineral data including: (1) Pb diffusion parameters; (2) measured closure temperatures (T_C) for fission-track annealing and (3) retentivities of both Pb and fission tracks, from apparent-age data. Every kinetic parameter (including T_C and mineral age for both the U/Pb and fission-track systems) is inversely correlated with Z within the sub-assemblage: zircon (Z ≈ 29%), titanite ( 34%) and apatite ( 38%). Assuming a diffusional closure model, Pb isotopic transport phenomena are described by a T_C-Z scale “calibrated” with field-based T_C data for titanite (≥ 680 ± 20°C) and apatite ( 500°C). Extrapolation of this scale yields T_C estimates for the following minerals: staurolite (T_C ≥ 1060°C, Z ≈ 25%); garnet (≥ 1010°C, 26.5%); zircon (≥900°C); monazite, xenotime, and epidote (≥ 750°C, 32%); and Ca-clinopyroxene (≥ 670 ± 30°C, 34 ± 1%, depending on composition). These empirical results imply that a (U/)Pb/Pb date for staurolite or garnet records the time of mineral growth, not post-growth isotopic closure, as also concluded in recent field studies. Because Z systematizes fission-track annealing, this recrystallization process, like volume-diffusion, must also be rate-limited by the strength of chemical bonds. The extent to which other recrystallization processes are likewise rate-limited is important to U/Pb geochronology because they potentially compete with diffusion as mechanisms for Pb-isotopic resetting in nature.     

Repeated thermal resetting of phengites in the Mulhacen Complex (Betic Zone, southeastern Spain) shown byAr/Ar step heating and single grain laser probe dating

This study reports the results of the first⁴⁰Ar/³⁹Ar combined induction furnace and laser probe dating of phengites from the Mulhacen HP/LT metamorphic complex in the Betic Cordilleras, southern Spain. Laser step heating and spot fusion analyses on different halves of a split single grain were made with a continuous laser probe. Spot fusion analysis resulted in ages of about 30–31 Ma in the core and ages as low as 25–26 Ma in the rim. Laser step heating on the other half of the grain gave a spectrum with apparent ages increasing from about 25 Ma to 29.5 Ma. The age spectrum and the decreasing ages towards the rim of the grain may imply that resetting essentially occurred by volume diffusion of radiogenic⁴⁰Ar due to late stage reheating resulting from extensional tectonics. Ages around 30 Ma in the core of the grain are interpreted as minimum estimates of the cooling age of the main tectono-metamorphic phaseD₂.

Induction furnace step heating on phengite separates from mica schists and one gneiss resulted in two types of age spectra. Type I spectra show monotonously rising apparent ages from14.5 ± 1.9 Ma to20.7 ± 0.2 Ma, and in a second sample from16.9 ± 0.7 to29.7 ± 0.2 Ma. Type II spectra are characterized by plateaus of14.4 ± 0.1 Ma (the gneiss sample),17.3 ± 0.1 Ma and17.6 ± 0.1 Ma. Type II spectra show low temperature apparent ages significantly below the plateau age, implying resetting subsequent to initial cooling. Modelling of the age spectra demonstrated that the plateau ages are possibly the result of strong resetting (75–85% of argon loss) of an older isotope system. Total fusion of a number of phengite single grains from marbles taken close to type II mica schists yielded ages of15.4 ± 1.2 Ma and17.0 ± 0.7 Ma. The observed repeated resetting is coeval with major volcanic activity in basins adjacent to the metamorphic ranges, pointing to a resetting by advective fluid transport related to volcanism.     

Deuterium tracing for the estimation of transpiration from trees Part 3. Measurements of transpiration from Eucalyptus plantation, India

Measurements of transpiration from individual trees of Eucalyptus from plantations at four different sites in Karnataka, Southern India, are presented. These show large (as much as tenfold) differences in the transpiration between premonsoon and postmonsoon periods, a reflection of the effects of soil-moisture stress in the premonsoon periods. For trees with diameters at breast height (DBH) less than 10 cm the transpiration rate of individual trees is proportional to the square of the DBH. For trees which are not experiencing soil-water stress the daily transpiration rate of individual trees, q, is well represented by the relation: q = (6.6 ± 0.3)g (m³ day⁻¹ where g (m²) is the tree basal area. On a unit ground area basis the transpiration rate, expressed as a depth per day is given by the relation: E_t = (0.66 ± 0.03)G (mm day⁻¹ where g(m²ha⁻¹) is the total basal area per hectare. For all the sites studied, although there is evidence for the ‘mining’ of soil water as roots penetrate deeper depths in the soil each year, there is no evidence for direct abstraction from the water table.     

Experimental determination of U and Th partitioning between clinopyroxene and natural and synthetic basaltic liquid

Clinopyroxene-silicate liquid partition coefficients for U and Th have been determined by particle track radiography from 1 atm crystallization experiments at controlledfO₂. Two natural basaltic and one synthetic composition were studied atfO₂ values from the NiNiO oxygen buffer to 1 log unit more oxidizing than FeFeO (IW+ 1). Over the range offO₂ values and compositions studied,D_U^cpx/liq = 0.0034–0.015,D_Th^cpx/liq = 0.008–0.036, andD_ThD_U= 3.4–1.1. With decreasingfO₂,D_ThD_U can decrease by up to a factor of 3 for a given composition, primarily from an increase inD_U^cpx/liq, which we interpret as resulting from an increase in the proportion of tetravalent U in the system with decreasingfO₂. This demonstrates that crystal-liquid UTh fractionation isfO₂ dependent and that U in terrestrial magmas is not entirely tetravalent.D_Th^cpx/liq appears to decrease in the two basalts at the lowestfO₂, possibly as a result of changes in composition withfO₂.

Our data show the sense of UTh fractionation by clinopyroxene-liquid partitioning is consistent with previous experimental determinations, in thatD_Th^cpxD_U^cpx> 1 in all cases. This indicates that, during partial melting, the liquid will have aTh/U ratio less than the clinopyroxene in the source. The observed²³⁸U²³⁰Th disequilibrium in MORB requires that the partial melt should have aTh/U ratio greater than the bulk source, and, therefore, cannot result from clinopyroxene-liquid partitioning. Further, the magnitudes of the measured partition coefficients are too small to generate significant UTh fractionation in either direction. Assuming that clinopyroxene contains the bulk of the U and Th in MORB source, our results indicate that²³⁸U²³⁰Th disequilibrium in MORB may not be caused by partial melting at all.     

The Impact of the Herbicide Diuron on Photosynthesis in Three Species of Tropical Seagrass

The impact and recovery from exposure to the herbicide diuron [DCMU; 3-(3′,4′-dichlorophenyl)-1,1-dimethylurea] was assessed for three tropical seagrasses, maintained in outdoor aquaria over a 10-day period. Photosynthetic stress was detected using chlorophyll a fluorescence, measured with a Diving-PAM (pulse amplitude modulated fluorometer). Exposure to 10 and 100 μg l⁻¹ diuron resulted in a decline in effective quantum yield (ΔF/F_m^′) within 2 h of herbicide exposure in Cymodocea serrulata, Halophila ovalis and Zostera capricorni. Effective quantum yield also declined over the first 24 h of exposure in H. ovalis at even lower diuron concentrations (0.1 and 1.0 μg l⁻¹). Effective quantum yield in H. ovalis and Z. capricorni was significantly depressed at all diuron concentrations (0.1–100 μg l⁻¹) after 5 days exposure, whereas effective quantum yield in C. serrulata was only significantly lower in plants exposed to highest diuron concentrations (10 and 100 μg l⁻¹). Effective quantum yield depression was present 5 days after plants exposed to 10 and 100 μg l⁻¹ diuron were returned to fresh seawater. These results indicate that exposure to herbicide concentrations present in nearshore Queensland sediments present a potential risk to seagrass functioning.     

Late cretaceous to early paleogene paleomagnetic results from Sikhote Alin, far eastern Russia: implications for deformation of East Asia

Welded tuffs in the Bogopol and Sijanov groups were sampled at 27 sites from 12 caldera formations in the Sikhote Alin mountain range around Kavalerovo (44.3°N, 135.0°E) for chronological and paleomagnetic studies. KAr age dates show that the welded tuffs erupted between 66 Ma and 46 Ma. All sites yield reliable paleomagnetic directions, with unblocking temperatures higher than 560°C. The high-temperature component at 12 sites and the medium-temperature component at 3 sites in the Bogopol Group show reversed polarity (D = 193.7°, I = −57.6°,₉₅ = 8.1°). The high-temperature component at 11 sites in the Sijanov Group showed both reversed and normal polarities and its mean direction reveals no detectable deflection from north (D = −2.9°, I = 59.6°,₉₅ = 11.2°). The combined paleomagnetic direction of the two groups yields a paleomagnetic pole of 250.5°E, 84.1°N (A₉₅ = 8.8°), which falls near Cretaceous paleomagnetic poles from Outer Mongolia, Inner Mongolia, the North China Block and the South China Block. The Sikhote Alin area appears not to have been subjected to detectable motion with respect to East Asia since about 50 Ma. This implies that the Sikhote Alin area behaved as an integral part of East Asia during the opening of the Japan Sea at about 15 Ma. However, significant separation between the paleomagnetic poles of East Asia and Europe during the Jurassic-Paleogene implies a major relative movement between these two blocks since the Paleogene.     

Mantle thermal structure and active upwelling during continental breakup in the North Atlantic

Seismic reflection and refraction data acquired on four transects spanning the Southeast Greenland rifted margin and Greenland–Iceland Ridge (GIR) provide new constraints on mantle thermal structure and melting processes during continental breakup in the North Atlantic. Maximum igneous crustal thickness varies along the margin from >30 km in the near-hotspot zone (<500 km from the hotspot track) to 18 km in the distal zone (500–1100 km). Magmatic productivity on summed conjugate margins of the North Atlantic decreases through time from 1800±300 to 600±50 km³/km/Ma in the near-hotspot zone and from 700±200 to 300±50 km³/km/Ma in the distal zone. Comparison of our data with the British/Faeroe margins shows that both symmetric and asymmetric conjugate volcanic rifted margins exist. Joint consideration of crustal thickness and mean crustal seismic velocity suggests that along-margin changes in magmatism are principally controlled by variations in active upwelling rather than mantle temperature. The thermal anomaly (ΔT) at breakup was modest (100–125°C), varied little along the margin, and transient. Data along the GIR indicate that the potential temperature anomaly (125±50°C) and upwelling ratio (4 times passive) of the Iceland hotspot have remained roughly constant since 56 Ma. Our results are consistent with a plume–impact model, in which (1) a plume of radius 300 km and ΔT of 125°C impacted the margin around 61 Ma and delivered warm material to distal portions of the margin; (2) at breakup (56 Ma), the lower half of the plume head continued to feed actively upwelling mantle into the proximal portion of the margin; and (3) by 45 Ma, both the remaining plume head and the distal warm layer were exhausted, with excess magmatism thereafter largely confined to a narrow (<200 km radius) zone immediately above the Iceland plume stem. Alternatively, the warm upper mantle layer that fed excess magmatism in the distal portion of the margin may have been a pre-existing thermal anomaly unrelated to the plume.     

Helium isotope geochemistry of mid-ocean ridge basalts from the South Atlantic

We report new helium isotope results for 49 basalt glass samples from the Mid-Atlantic Ridge between 1°N and 47°S.³He/⁴He in South Atlantic mid-ocean ridge basalts (MORB) varies between 6.5 and 9.0 R_A (R_A is the atmospheric ratio of1.39 × 10⁻⁶), encompassing the range of previously reported values for MORB erupted away from high³He/⁴He hotspots such as Iceland. He, Sr and Pb isotopes show systematic relationships along the ridge axis. The ridge axis is segmented with respect to geochemical variations, and local spike-like anomalies in³He/⁴He, Pb and Sr isotopes, and trace element ratios such as(La/Sm)_N are prevalent at the latitudes of the islands of St. Helena, Tristan da Cunha and Gough to the east of the ridge. The isotope systematics are consistent with injection beneath the ridge of mantle “blobs” enriched in radiogenic He, Pb and Sr, derived from off-axis hotspot sources. The variability in³He/⁴He along the ridge can be used to refine the hotspot source-migrating-ridge sink model.

MORB from the 2–7°S segment are systematically the least radiogenic samples found along the mid-ocean ridge system to date. Here the depleted mantle source is characterized by⁸⁷Sr/⁸⁶Sr of 0.7022, Pb isotopes close to the geochron and with²⁰⁶Pb/²⁰⁴Pb of 17.7, and³He/⁴He of 8.6–8.9 R_A. The “background contamination” of the subridge mantle, by radiogenic helium derived from off-ridge hotspots, displays a maximum between 20 and 24°S. The HePb and HeSr isotope relations along the ridge indicate that the³He/⁴He ratios are lower for the hotspot sources of St. Helena, Tristan da Cunha and Gough than for the MORB source, consistent with direct measurements of³He/⁴He ratios in the island lavas. Details of the HeSrPb isotope systematics between 12 and 22°S are consistent with early, widespread dispersion of the St. Helena plume into the asthenosphere, probably during flattening of the plume head beneath the thick lithosphere prior to continental breakup. The geographical variation in theHe/Pbratio deduced from the isotope systematics suggests only minor degassing of the plume during this stage. Subsequently, it appears that the plume component reaching the mid-Atlantic ridge was partially outgassed of He during off-ridge hotspot volcanism and related melting activity.

Overall, the similar behavior of He and Pb isotopes along the ridge indicates that the respective mantle sources have evolved under conditions which produced related He and Pb isotope variations.     

An integrated paleomagnetic study of Rio Grande de Santiago volcanic succession (trans-Mexican volcanic belt): revisited

We carried out an integrated paleomagnetic, rock-magnetic and paleointensity study of Miocene volcanic succession from the trans-Mexican volcanic belt (TMVB) north of Guadalajara. A total of 37 consecutive basaltic lava flows (326 oriented standard paleomagnetic cores) were collected at Lazo locality. Continuous susceptibility measurements with temperature and hysteresis experiments yield in most cases reasonably reversible curves with Curie points close to that of pseudo-single-domain magnetite. Two geomagnetic reversals were observed in the 300 m thick composite section. Paleosecular variation was lower than the one observed in general during Miocene. It appears that the volcanic units have been emplaced during a relatively short time span of about 1 Ma. The mean paleomagnetic directions obtained from this study do not differ significantly from that expected for the middle Miocene. The mean paleomagnetic direction calculated from all data is I=31.1°, D=354.6°, k=124 and ₉₅=2.1°, N=37. Seventy-two samples with apparently preserved primary magnetic mineralogy and without secondary magnetization, mostly belonging to reverse polarity chron were pre-selected for Thellier paleointensity determination. The flow-mean paleointensity values are ranging from 22.4±3.4 to 53.8±6.0 μT and the corresponding virtual dipole moments (VDMs) are ranging from (5.4±0.8) to (12.0±1.4)×10²² A m². This corresponds to mean value of (7.7±2.2)×10²² A m², which is close to present day geomagnetic field strength. Altogether, our data suggest the existence of relatively high geomagnetic field strength undergoing low fluctuations.     

Inclination shallowing and preferred transitional VGP paths

Transitional VGP paths recorded in sediments cluster into two antipodal preferred longitude bands that tend to lie 90° away from their site longitudes, the latter also being clustered. VGP paths obtained from lava flow sequences, though much fewer, appear not to show these biases, suggesting a rock-magnetic influence on VGP paths recorded in sediments. Inclination shallowing of detrital magnetic remanence, enhanced under low transitional field strengths, is the most likely candidate.

We illustrate the effects of inclination shallowing by applying a simple shallowing model (tan I_R = f tan I_A, where I_A is the inclination of the magnetic remanence and I_A is the inclination of the ambient field) with field variation to hypothetical data sets. Shallowing-induced clustering increases as f decreases and becomes extreme as f approaches 0.1.

We have used the model to ‘de-shallow’ the available set of transitional VGP sediment records for various values of f. The probability that the observations arise from inclination shallowing of a uniform random distribution of paths increases as f decreases. When f drops to 0.13 there is a 50% chance of getting at least as much grouping as observed. To decide if inclination shallowing is a dominant factor in the clustering, we need to know whether such extreme shallowing is widespread in sedimentary records under transitional field conditions. Field and laboratory redeposition data are not yet adequate to resolve this question.     

REE diffusion in calcite

Chemical diffusion of four rare-earth elements (La, Nd, Dy and Yb) has been measured in natural calcite under anhydrous conditions, using rare-earth carbonate powders as the source of diffusants. Experiments were run in sealed silica capsules along with finely ground calcite to ensure stability of the single-crystal samples during diffusion anneals. Rutherford backscattering spectroscopy (RBS) was used to measure diffusion profiles. The following Arrhenius relations were obtained over the temperature range 600–850°C: D_La =2.6×10⁻¹⁴ exp(−147±14 kJ mol⁻¹/RT) m² s⁻¹, D_Nd =2.4×10⁻¹⁴ exp(−150±13 kJ mol⁻¹/RT) m² s⁻¹, D_Dy =2.9×10⁻¹⁴ exp(−145±25 kJ mol⁻¹/RT) m² s⁻¹, D_Yb =3.9×10⁻¹² exp(−186±23 kJ mol⁻¹/RT) m² s⁻¹. In contrast to previous findings for refractory silicates (e.g. zircon), differences in transport rates among the REE are not pronounced over the range of temperature conditions investigated in this study. Diffusion of the REE is significantly slower than diffusion of the divalent cations Sr and Pb and slower than transport of Ca and C at temperatures above 650°C. Fine-scale zoning and isotopic and REE chemical signatures may be retained in calcites under many conditions if diffusion is the dominant process affecting alteration.     

Dating earthquakes with high-precision thorium-230 ages of very young corals

Mass spectrometric techniques have recently yielded significant increases in precision and sensitivity over previous methods for measuring²³⁰Th abundance in corals. To assess the accuracy of²³⁰Th ages, three corals from Vanuatu, whose ages were known from counting annual growth bands, were analyzed. For each sample, the date of growth determined by²³⁰Th analysis (A.D.1969 ± 3, 1932 ± 5,and1806 ± 5; 2σ uncertainties based on analytical error) was indistinguishable from the date determined by counting bands (A.D. 1971–1973, 1935–1937, and 1804–1810), indicating that the²³⁰Th dates are accurate.²³⁰Th dates were also determined for two adjacent emerged heads from Santo Is., Vanuatu, which were thought to have died when they were raised above sea level during coseismic uplift. The dates (A.D.1864 ± 4, 1866 ± 4) were the same, indicating that the heads died at the same time and consistent with the idea that they were killed by coseismic emergence around A.D. 1865. The difference between this date and the date of the only major historically documented earthquake that caused uplift (A.D. 1973,M_s = 7.5), suggests a seismic recurrence interval of108 ± 4y for Santo. Analogous emerged corals from Malekula Is., Vanuatu yielded²³⁰Th dates that were similar to each other (A.D.1729 ± 3, 1718 ± 5) and are inferred to have died during coseismic emergence around A.D. 1729. In conjunction with the date of the only large historically documented earthquake that caused uplift (A.D. 1965,M_s = 7.5), the recurrence interval for Malekula is236 ± 3y. If similar emerged corals can be found, this appraoch may be extended back in time and to other localities because it appears that such features can now be dated both accurately and precisely.     

Glacial–interglacial changes in Subantarctic sea surface temperature and δO-water using foraminiferal Mg

Laboratory culturing experiments with living Globigerina bulloides indicate that Mg/Ca is primarily a function of seawater temperature and suggest that Mg/Ca of fossil specimens is an effective paleotemperature proxy. Using culturing results and a core-top Neogloboquadrina pachyderma calibration, we have estimated glacial–interglacial changes in sea surface temperature (SST) using planktonic Mg/Ca records from core RC11-120 in the Subantarctic Indian Ocean (43°S, 80°E) and core E11-2 in the Subantarctic Pacific Ocean (56°S, 115°W). Our results suggest that glacial SST was about 4°C cooler in the Subantarctic Indian Ocean and 2.5°C cooler in the Subantarctic Pacific. Comparison of SST and planktonic δ¹⁸O records indicates that changes in SST lead changes in δ¹⁸O by on average 1–3 kyr. The glacial–interglacial temperature change indicated by the Subantarctic Mg/Ca records suggests that temperature accounts for 40–60% of the foraminiferal δ¹⁸O change. We have used the Mg/Ca-based SST estimates and δ¹⁸O determinations to generate site-specific seawater δ¹⁸O records, which suggest that seawater δ¹⁸O was on average 1‰ more positive during glacial episodes compared with interglacial episodes.     

Calibrations for benthic foraminiferal Mg/Ca paleothermometry and the carbonate ion hypothesis

Benthic foraminiferal magnesium/calcium ratios were determined on one hundred and forty core-top samples from the Atlantic Ocean, the Norwegian Sea, the Indian Ocean, the Arabian Sea and the Pacific Ocean, mostly at sites with bottom water temperatures below 5 °C. Mg/Ca ratios are consistently lower, by  0.2 mmol/mol, in samples cleaned using oxidative and reductive steps than using oxidative cleaning. Differences between Cibicidoides species have been identified: Mg/Ca of Cibicidoides robertsonianus > Cibicidoides kullenbergi > Cibicidoides wuellerstorfi. Comparison with bottom water temperatures support observations of lowered Mg/Ca of C. wuellerstorfi at temperature below  3 °C compared with values predicted by published calibrations and from other Cibicidoides species. Hydrographic data shows that carbonate ion saturation (Δ[CO₃²⁻]) decreases rapidly below this temperature. An empirical sensitivity of Δ[CO₃²⁻] on Mg/Ca has been established for C. wuellerstorfi of 0.0086 ± 0.0006 mmol/mol/μmol/kg. A novel application using modern temperatures and Last Glacial Maximum temperatures derived via pore fluid modelling supports a carbonate ion saturation state effect on Mg incorporation. This may significantly affect calculated δ¹⁸O_seawater obtained from foraminiferal δ¹⁸O and Mg/Ca temperature.     

Mathematical modeling of PCB bioaccumulation in Perna viridis

In the present work, we built a mathematical model of polychlorinated biphenyl (PCB) bioaccumulation in Perna viridis, namely, a one-compartment model with a time dependent incorporation rate R (μg g⁻¹ lipid per ppb water per day), with positive substrate cooperativity as the underlying physical mechanism. The temporal change of the PCB concentration Q (μg g⁻¹ lipid) in the soft tissues of the mussel depends on the competition of the input rate RW and the output rate kQ, where W is the concentration of PCB in water (ppb water) and k is the elimination rate (per day). From our experimental data, k=0.181±0.017 d⁻¹. The critical concentration in water W_c for positive substrate cooperativity was found to be 2.4 ppb. Below W_c, R is a constant. For a water concentration of 0.5 ppb Aroclor 1254, R=24.0±2.4 μg g⁻¹ lipid ppb⁻¹ d⁻¹. Above W_c, positive substrate cooperativity comes into effect and R becomes a function of time and dependent on the concentration Q in a form R=γQ/(Q+δ). This is the case for a water concentration of 5 ppb Aroclor 1254, where γ=15.1 μg g⁻¹ lipid ppb⁻¹ d¹ and δ≈200 μg g⁻¹ lipid. From this model, the uptake is exponentially increasing when the PCB concentration in the mussel is small compared to 200 μg g⁻¹ lipid, and hyperbolically increasing when the concentration is large compared to 200 μg g⁻¹ lipid, which are consistent with the experimental data. The model is useful for understanding the true processes taking place during the bioaccumulation and for risk assessment with higher confidence. Future experimental data which challenge the present model are anticipated and in fact desirable for improvement and perfection of the model.     

Preliminary stratigraphic and structural architecture of Bhutan: Implications for the along strike architecture of the Himalayan system

High-precision in-situ ion microprobe (SIMS) oxygen isotope analysis of zircons from two diorite intrusions associated with the late Caledonian Lochnagar pluton in Scotland has revealed large differences in the degree of heterogeneity in zircon δ¹⁸O between the diorites. Zircon crystals from the Cul nan Gad diorite (CnG) show a unimodal distribution of oxygen isotope values (δ¹⁸O = 6.0 ± 0.6‰ (2σ)) and no or only minor grain-scale variation. Those from the Allt Darrarie diorite (AD1) show a large range in δ¹⁸O and an apparent bimodal distribution with modes of 6.6 ± 0.4‰ and 7.3 ± 0.4‰. Variations of up to 1.2‰ occur between and within grains; both an increase and decrease in δ¹⁸O with zircon growth has been observed. The δ¹⁸O composition of growing zircon can only change if open-system processes affect the magma composition, i.e. if material of contrasting δ¹⁸O composition is added to the magma. The variability in AD1 is interpreted to represent a cryptic record of magma mixing. A ‘deep crustal hot zone’ is a likely site for generation of the dioritic magmas which developed by mixing of residual melts and crustal partial melts or by melting of mafic lower crustal rocks. The overall small number of zircons with mantle-like δ¹⁸O values (5.3 ± 0.6‰ (2σ)) in the Lochnagar diorites is largely the product of crustal differentiation rather than crustal growth.

The δ¹⁸O of quartz from the CnG and AD1 diorites shows only minor variation (CnG: 10.9 ± 0.5‰ (2σ), AD1: 11.7 ± 0.6‰ (2σ)) within single populations, with no evidence of mixing. Quartz–zircon isotopic disequilibrium is consistent with later crystallisation of quartz from late magmatic fluids, and in case of the AD1 diorite after the inferred magma mixing from a homogenised, higher δ¹⁸O melt.

High-precision SIMS oxygen isotope analysis of zircon provides a new approach to identifying and resolving previously undetected early-stage magma mixing and constraining the compositions and origins of the component magmas. A combination of zircon, quartz and whole-rock data has proven to be a powerful tool in reconstructing the petrogenetic evolution of diorite from early crystallisation to late alteration.     

Sulfate reduction using methane in sediments beneath a bathyal “cold seep” giant clam community off Hatsushima Island, Sagami Bay, Japan

Two sandy sediment cores (Cores D227-120 and D380) were collected from inside a deep-sea giant clam (Calyptogena soyoae) community off Hatsushima Island, western Sagami Bay, central Japan (35°59.9′N, 139°13.6′E; 1160 m deep) and a muddy sediment core (Core D227-202) was obtained from outside the community by the submersibleShinkai 2000. The chloride concentration of the pore waters is constant vertically and sulfate reduction using sedimentary organic matter occurs in Core D227-202 (21 cm long). The chloride concentrations are lower by 7% at the 7.5–9 cm depth in Core D227-120 (9 cm long) and by 3% at the 11–12 cm depth in Core D380 (16 cm long) than those of the overlying bottom waters in the cores from inside of the community. Sulfate concentration decreases remarkably and dissolved inorganic carbon, alkalinity, ammonium-N, and hydrogen sulfide concentrations increase significantly with increasing depth in Core D380.δ³⁴S values of sulfate ions increase from +20.5 to +35.3‰ andδ¹³C values of dissolved inorganic carbon decrease drastically from −7.0 to −45‰ with increasing depth from the top to the bottom of the core, although theδ¹³C values of the organic carbon of the sediments are−23.7 ± 0.9‰ in Core D380. These results indicate that sulfate reduction using methane is active within the sediments just beneath the living clams and that the hydrogen sulfide produced can be used by endosymbiotic sulfur oxidizing bacteria living in the gills ofC. soyoae in the community.     

Coral bleaching on high-latitude marginal reefs at Sodwana Bay,South Africa

Coral bleaching, involving the expulsion of symbiotic zooxanthellae from the host cells, poses a major threat to coral reefs throughout their distributional range. The role of temperature in coral bleaching has been extensively investigated and is widely accepted. A bleaching event was observed on the marginal high-latitude reefs of South Africa located at Sodwana Bay during the summer months of 2000. This was associated with increased sea temperatures with high seasonal peaks in summer and increased radiation in exceptionally clear water. The bleaching was limited to Two-mile Reef and Nine-mile Reef at Sodwana Bay and affected <12% of the total living cover on Two-mile Reef. Montipora spp., Alveopora spongiosa and Acropora spp. were bleached, as well as some Alcyoniidae (Sinularia dura, Lobophytum depressum, L. patulum). A cyclical increase in sea temperature (with a period of 5–6 years) was recorded during 1998–2000 in addition to the regional temperature increase caused by the El Niño Southern Oscillation phenomenon. The mean sea temperature increased at a rate of 0.27 °C year⁻¹ from May 1994 to April 2000. High maximum temperatures were measured (>29 °C). The lowest mean monthly and the mean maximum monthly temperatures at which coral bleaching occurred were 27.5 and 28.8 °C, respectively, while the duration for which high temperatures occurred in 2000 was 67 days at 27.5 °C (4 days at 28.8 °C). Increased water clarity and radiation appeared to be a synergistic cause in the coral bleaching encountered at Sodwana Bay.     

The equation of state of CaSiO3 perovskite to 108 GPa at 300 K

Pressure–volume measurements have been performed for CaSiO₃ perovskite to 108 GPa at 300 K using NaCl and argon pressure-transmitting media, and energy dispersive X-ray diffraction (EDXD) in a diamond-anvil cell (DAC). By determining a parameter that is the product of the elastic anisotropy, S, and the uniaxial stress component, t, for each data point, we define the stress condition of the sample. For different points at the same pressure in a temperature-quenched sample, the St value can differ by as much as a factor of 5, indicating heterogeneity in the stress condition. This may be responsible for the large scatter of earlier P–V measurements in the DAC which in general used a large diameter X-ray beam. Also, the St value provides insight into the elastic anisotropy, S, of CaSiO₃ perovskite and platinum. The sign of S (positive) for CaSiO₃ perovskite agrees with first principles calculations but the magnitude may be inconsistent. A new compression curve at 300 K was obtained for CaSiO₃ perovskite by using those data points which represent the most nearly hydrostatic conditions. It is observed that the data points with high St values yield larger volumes than the points with small St values at a given pressure. By selecting the data points having low St values (St≤0.005), combining with lower pressure large volume press (LVP) measurements and fitting to third order Birch–Murnaghan equation of state (EOS), we find that CaSiO₃ perovskite is more compressible (V₀=45.58±0.05 ų, K_T0=236±4 GPa, and K_T0′=3.9±0.2 GPa) than suggested by previous studies. The density and bulk modulus of CaSiO₃ perovskite at lower mantle pressures and 300 K are 1–3% greater and 5–15% smaller, respectively, than found in previous studies. This study demonstrates that defining the stress state of the sample is crucial to obtain an accurate 300 K compression curve for unquenchable high-pressure phases.