Spatial variability in growth and reproduction of the Pacific oyster Crassostrea gigas (Thunberg, 1793) along the west European coast

The Pacific oyster Crassostrea gigas was introduced in Europe for commercial purposes in the mid 1960s. It was initially thought that low winter temperatures would restrain this species' reproduction and settlement; however, its present distribution in areas where no introduction has taken place suggests that natural invasion and expansion has occurred. Along the European coast, wild populations of Pacific oysters are already found from northern Germany to southern Portugal. Whether C. gigas will continue to further expand through northern waters will depend on its physiological performance. In this study, the performance of wild oyster populations has been studied in terms of growth and reproduction at three stations: La Rochelle (France; 46°N), Yerseke (Oosterschelde estuary, The Netherlands, 51°N), and Texel (Wadden Sea estuary, The Netherlands, 53°N). The French population had the lowest somatic-shell mass ratio and an increase in maximum shell length, somatic and gonadal mass was observed from France to the Netherlands. In addition, mean oocyte diameter decreased significantly from south to north. The combination of increasing gonadal mass and decreasing oocyte volume suggests an increasing reproductive output in terms of egg numbers from France to The Netherlands. Differences in temperature between locations will at least be partly responsible for the observed patterns; however, other environmental factors (such as food availability, predation pressure, sediment type and/or seston concentration) cannot be excluded. Since smaller eggs (oocytes) are thought to have a longer development time, the environmental conditions along the Dutch coast may result in increased larval dispersal and possibly in further population expansion.     

Growth and reproduction of the bivalve Spisula subtruncata (da Costa) in Dutch coastal waters

The bivalve Spisula subtruncata is usually abundant in shallow coastal waters along the Dutch coast. However, its biomass has been decreasing since 1995. In order to assess whether reproductive failure may be the cause of the observed decline over the last decades, the energy investment in reproduction of a population of S. subtruncata from central Dutch coastal waters was studied. The population studied consisted of individuals of up to four years old. Shell length reached maximum values of around 32 mm and individual total body, somatic and gonadal ash-free dry mass reached maximum values of about 278 mg AFDM, 252 mg AFDM and 76 mg AFDM, respectively. A clear seasonal cycle in somatic and gonadal mass was observed. Somatic and gonadal mass indices increased in early spring and reached maximum values during summer, followed by a decrease to minimum values at the beginning of the following year. Spawning was in June–July and settlement of spat seems to have occurred in July–August. Mean oocyte diameter was 57.43 ± 0.03 μm, corresponding to a volume of 98972 μm³. These results suggested that reproductive failure was not the cause of the current population decline. Most likely, unsuccessful settlement of spat and/or severe predation during the first months of life were responsible for the observed patterns.     

The effect of submarine groundwater discharge on nutrient and salinity regimes in a coastal lagoon off Perth, Western Australia

Submarine groundwater discharge (SGD) into a coastal lagoon off Perth, Western Australia, contains nitrate and silicate in concentrations two orders of magnitude higher than those of the receiving waters. This discharge delivers enough nitrate to replace that dissolved in the lagoon water mass about every eight days and enough silicate to replace the lagoon silicate in about 48 days. The delivery rate of nitrate nitrogen by SGD is equal to about 48% of that required for observed growth rates of lagoon macrophytes. Surface salinity is lower close to the shore as a result of SGD. During calm conditions a salinity front was observed in the lagoon, with a nearshore pool of nutrient-enriched water floating above the more saline ocean water.     

Th sorption and export models in the water column: A review

Over the past few decades, the radioisotope pair of ²³⁸U / ²³⁴Th has been widely and increasingly used to describe particle dynamics and particle export fluxes in a variety of aquatic systems. The present paper is one of five review articles dedicated to ²³⁴Th. It is focused on the models associated with ²³⁴Th whereas the companion papers (same issue) are focused on present and future methodologies and techniques (Rutgers van der Loeff et al.), C / ²³⁴Th ratios (Buesseler et al.), ²³⁴Th speciation (Santschi et al.) and present and future applications of ²³⁴Th [Waples, J.T., Benitez-Nelson, C.R., Savoye, N., Rutgers van der Loeff, M., Baskaran, M., Gustafsson, Ö., this issue. An Introduction to the application and future use of ²³⁴Th in aquatic systems. Marine Chemistry, FATE special issue]. In this paper, we review current ²³⁴Th scavenging models and discuss the relative importance of the non-steady state and physical terms associated with the most commonly used model to estimate ²³⁴Th flux. Based on this discussion we recommend that for future work the use of models should be accompanied by a discussion of the effect that model and data uncertainty have on the model results. We also suggest that future field work incorporate repeat occupations of sample sites on time scales of 1–4 weeks in order to evaluate steady state versus non-steady state estimates of ²³⁴Th export, especially during high flux events (> ca. 800 dpm m^− 2 d^− 1). Finally, knowledge of the physical oceanography of the study area is essential, particularly in ocean margins and in areas of established upwelling (e.g., Equatorial Pacific). These suggestions will greatly enhance the application of ²³⁴Th as a tracer of particle dynamics and flux in more complicated regimes.     

Solar rotation velocity as determined from sunspot drawings of J. Hevelius in the 17th century

Using two original copies of Hevelius' Selenographia and reducing spot positions with two different methods, we found that the solar angular rotation velocity at the beginning of the Maunder minimum was about the same as today. The gradient of the differential rotation was slightly steeper than given in modern reductions, but not significantly different. These findings are in contradiction to those published by Eddy et al. (1976).     

Regional record of a global oceanic anoxic event: OAE1a on the Apulia Platform margin, Gargano Promontory, southern Italy

The Gargano Promontory of southern Italy, located on the eastern margin of the Apulia Platform, represents a peculiar Tethyan area where the transition between carbonate platform and adjacent basins is exposed on land. The Aptian stratigraphic record, represented in shallow-water, slope and deep-water deposits, provides a good opportunity to investigate the regional response to the worldwide documented climatic, biotic and palaeoceanographic changes related to Oceanic Anoxic Event 1a (OAE1a). A synthesis of data previously published (Ischitella and Coppitella sections), together with original data (Val Carbonara and Coppa della Guardia sections), from four stratigraphic sections from different depositional settings (proximal to distal) is provided, using an integrated, high-resolution micropalaeontological (planktonic foraminifera and calcareous nannofossils) and, for one section, geochemical (stable carbon and oxygen isotopes) approach. Organic matter preservation is confined to the more distal areas and consists of two thin intervals of black shales in the Aptian portion of the Marne a Fucoidi Formation. Biostratigraphic data assign the older black shale (5 cm thick) to the Selli Level equivalent (OAE1a, Lower Aptian); this carbon-rich interval is immediately followed by another black shale (7–10 cm thick) of early Late Aptian age.OAE1a is generally interpreted as a high-productivity event during a warming interval, followed by a cooling trend. In the Gargano Promontory, although the oxygen isotope curve indicates the above-mentioned climatic evolution, the micropalaeontological data do not support high fertility of the surface water, whereas micropalaeontological and geochemical data for the younger black shale do record high productivity (radiolarian increase in overall abundance, low nannofossil and foraminiferal species richness, increase in abudance of nannofossil fertility indices) associated with a cooling trend. The carbon and oxygen isotope record is in line with evidence from the curves documented elsewhere, whereas, among the biotic events, only the “nannoconid crisis” preceding OAE1a is revealed to be globally correlable.Environmental models for the two episodes of organic matter preservation are proposed, taking into account both global and local controlling factors.     

Outlier detection algorithms and their performance in GOCE gravity field processing

The satellite missions CHAMP, GRACE, and GOCE mark the beginning of a new era in gravity field determination and modeling. They provide unique models of the global stationary gravity field and its variation in time. Due to inevitable measurement errors, sophisticated pre-processing steps have to be applied before further use of the satellite measurements. In the framework of the GOCE mission, this includes outlier detection, absolute calibration and validation of the SGG (satellite gravity gradiometry) measurements, and removal of temporal effects. In general, outliers are defined as observations that appear to be inconsistent with the remainder of the data set. One goal is to evaluate the effect of additive, innovative and bulk outliers on the estimates of the spherical harmonic coefficients. It can be shown that even a small number of undetected outliers (<0.2 of all data points) can have an adverse effect on the coefficient estimates. Consequently, concepts for the identification and removal of outliers have to be developed. Novel outlier detection algorithms are derived and statistical methods are presented that may be used for this purpose. The methods aim at high outlier identification rates as well as small failure rates. A combined algorithm, based on wavelets and a statistical method, shows best performance with an identification rate of about 99%. To further reduce the influence of undetected outliers, an outlier detection algorithm is implemented inside the gravity field solver (the Quick-Look Gravity Field Analysis tool was used). This results in spherical harmonic coefficient estimates that are of similar quality to those obtained without outliers in the input data.     

Siberian river run-off and Late Quaternary glaciation in the southern Kara Sea, Arctic Ocean: preliminary results

The extent of the Barents-Kara Ice Sheet during the eastern Last Glacial Maximum (LGM) is not yet fully known. A detailed echo-sounding survey performed during the Boris Petrov Expedition 2001 permitted the detailed mapping of part of it. Based on the profiling results, a southern connection between the LGM Barents-Kara Ice Sheet and a local ice sheet on Taymyr Peninsula appears to be unlikely. Based on sediment core data and profiling results, most of the terrigenous river-derived material accumulated in the estuaries during late Holocene times, whereas during early Holocene times of lowered sea level major amounts were transported further offshore and accumulated on the shelf. During the post-glacial sea level rise, the main depocentre migrated southward, reaching its present position no earlier than about 6 cal. Ky BP (or 5.2 Kya). Future studies of accelerator mass spectrometry (AMS) ¹⁴C-dated sediment cores will allow a detailed reconstruction of the variability of fluvial sediment discharge and the history of glaciation in the Kara Sea during late Quaternary times.     

Melt/biotite B/B isotopic fractionation and the boron local environment in the structure of volcanic glasses

This paper is focused on the role of boron coordination in determining the ¹¹B/¹⁰B isotopic fractionation between melt/glass and biotite at magmatic temperatures. For this purpose, three evolved volcanic rocks from Roccastrada, Mt. Amiata, and Mt. Cimini belonging to the Neogene-Quaternary magmatism of central Italy were studied. In these samples, the measured boron biotite-glass partition coefficient ranges between 0.004 and 0.011, indicating that boron behaves as an incompatible element during biotite crystallization. The ¹¹B magic-angle spinning nuclear magnetic resonance (NMR) spectra reveal the presence of trigonal BO_3/2 units, tetrahedral BO_4/2⁻ sites, and three-coordinated BO_2/2O⁻ species containing one nonbridging oxygen. The relative contributions of these different boron sites were estimated by spectral deconvolution, and it was observed that the fraction of trigonally coordinated boron decreases with increasing K₂O concentration in the glass. The ¹¹B/¹⁰B isotopic fractionation between biotite and melt/glass was observed to be large even at magmatic temperatures and was found to be 1.0066 (Roccastrada sample), 1.00535 (Mt. Amiata sample), and 1.00279 (Mt. Cimini sample). Fractionation is mostly related to the relative amount of trigonal and tetrahedral boron sites in the glass network rather than to other processes, including the speciation of hydrous species in the glass structure. The measured α values are significantly higher than the calculated ones obtained using the reduced partition function ratios (RPFRs) for B(OH)₃ and B(OH)₄⁻ as reported by Kakihana et al. (1977) and the abundance of trigonal and tetrahedral boron obtained by ¹¹B NMR spectra. Furthermore, a nonlinear relationship is observed between the percentage of BO₄ in the glass structure and the measured 1000lnα, suggesting that the approximation of monomeric B(OH)₃ and B(OH)₄⁻ species contributions through ideal mixing in calculating the RPFRs in polyanions (Oi et al., 1989) probably does not apply to silicate glasses.The large B isotopic fractionation measured between glass and biotite and its dependence on the boron coordination in the glass are a limitation to the use of δ¹¹B in the mineral to characterize magmas. Nonetheless, the high incompatible behavior of boron in the most common magmatic minerals rules out that fractional crystallization significantly modified the B isotopic composition of the melt.     

The forward problem of electromagnetic induction: accurate finite-difference approximations for two-dimensional discrete boundaries with arbitrary geometry

We describe finite-difference approximations to the equations of 2-D electromagnetic induction that permit discrete boundaries to have arbitrary geometrical relationships to the nodes. This allows finite-difference modelling with the flexibility normally ascribed to finite-element modelling. Accuracy is demonstrated by comparison with finite-element computations. We also show that related approximations lead to substantially improved accuracy in regions of steep, but not discontinuous, conductivity gradient.