Cokriging nonstationary data

Universal cokriging is used to obtain predictions when dealing with multivariate random functions. An important type of nonstationarity is defined in terms of multivariate random functions with increments which are stationary of orderk. The covariance between increments of different variables is modeled by means of the pseudo-cross-covariance function. Criteria are formulated to which the parameters of pseudo-cross-covariance functions must comply so as to ensure positive-definiteness. Cokriging equations and the induced cokriging equations are given. The study is illustrated by an example from soil science.     

On least squares estimation of generalized covariance functions

Least squares estimation (LSE) is theoretically related to quadratic unbiased estimation of variance components. It is argued that these methods of estimation of variance components essentially generalize LSE though they are not formally equivalent.     

Palaeotemperature reconstructions of the European permafrost zone during marine oxygen isotope Stage 3 compared with climate model results

A palaeotemperature reconstruction based on periglacial phenomena in Europe north of approximately 51 °N, is compared with high‐resolution regional climate model simulations of the marine oxygen isotope Stage 3 (Stage 3) palaeoclimate. The experiments represent Stage 3 warm (interstadial), Stage 3 cold (stadial) and Last Glacial Maximum climatic conditions. The palaeotemperature reconstruction deviates considerably for the Stage 3 cold climate experiments, with mismatches up to 11 °C for the mean annual air temperature and up to 15 °C for the winter temperature. However, in this reconstruction various factors linking climate and permafrost have not been taken into account. In particular a relatively thin snow cover and high climatic variability of the glacial climate could have influenced temperature limits for ice‐wedge growth. Based on modelling the 0 °C mean annual ground temperature proves to be an appropriate upper temperature limit. Using this limit, mismatches with the Stage 3 cold climate experiments have been reduced but still remain. We therefore assume that the Stage 3 ice wedges were generated during short (decadal time‐scale) intervals of extreme cold climate, below the mean temperatures indicated by the Stage 3 cold climate model simulations. Copyright © 2003 John Wiley & Sons, Ltd.     

The Preboreal climate reversal and a subsequent solar‐forced climate shift

Accurate chronologies are essential for linking palaeoclimate archives. Carbon‐14 wiggle‐match dating was used to produce an accurate chronology for part of an early Holocene peat sequence from the Borchert (The Netherlands). Following the Younger Dryas–Preboreal transition, two climatic shifts could be inferred. Around 11 400 cal. yr BP the expansion of birch (Betula) forest was interrupted by a dry continental phase with dominantly open grassland vegetation, coeval with the PBO (Preboreal Oscillation), as observed in the GRIP ice core. At 11 250 cal. yr BP a sudden shift to a humid climate occurred. This second change appears to be contemporaneous with: (i) a sharp increase of atmospheric ¹⁴C; (ii) a temporary decline of atmospheric CO₂; and (iii) an increase in the GRIP ¹⁰Be flux. The close correspondence with excursions of cosmogenic nuclides points to a decline in solar activity, which may have forced the changes in climate and vegetation at around 11 250 cal. yr BP. Copyright © 2004 John Wiley & Sons, Ltd.     

An optimized method for automated analysis of algal pigments by HPLC

A recent development in algal pigment analysis by high-performance liquid chromatography (HPLC) is the application of automation. An optimization of a complete sampling and analysis protocol applied specifically in automation has not yet been performed. In this paper we show that automation can only be successful if the various methodological aspects of the sampling and analysis protocol are considered in coherence. We introduce an optimized protocol that involves freeze-drying of the sample, subsequent extraction in 90% acetone and the application of water-packing during analysis. The method was evaluated on both natural plankton populations and a broad spectrum of microalgal cultures: Thalassiosira weisflogii (Bacillariophyceae), Emiliania huxleyi (Prymnesiophyceae), Phaeocystis globosa and Phaeocystis antarctica (Prymnesiophyceae) and Pyramimonas sp. (Prasinophyceae). Whereas pigment extracts were unstable in methanol, with recorded chlorophyll a losses from 10% to 60% per day, pigment degradation rates in acetone were generally less than 1% over 18 h storage in the autosampler (4 °C). In addition, it was found that the extraction efficiency of acetone significantly increased upon freeze-drying prior to extraction. Increases as high as 50–60% were measured in P. antarctica. The application of water-packing of the sample during injection resulted in improved peak shape and peak separation, without diluting the pigment concentrations. Automation is especially beneficial for application in the field, when mixed algal assemblages and low biomass put a high demand on the sensitivity as well as reproducibility of the method.     

Transient biouptake flux and accumulation by microorganisms: The case of two types of sites with Langmuir adsorption

The uptake of a chemical species by an aquatic microorganism is modelled considering two kinds of sites where Langmuirian adsorption is followed by first order internalisation kinetics. Simpler models, such as only one internalisation route (while most of the adsorption takes place on non-internalising sites) or a linear isotherm for adsorption on one or both sites, become limiting cases of this double-Langmuirian model. The model considers the sites located on the spherical (or semi-spherical) surface of the organism, and takes diffusion from the medium into explicit account. The numerical solution for the internalisation flux shows a maximum. We provide an estimate for the time needed to reach a certain proximity to steady state. The transient solution confirms that the analytical expressions for the steady-state flux are usually valid and that the accumulated amounts reflect the impact of the short-time uptake. The Instantaneous Steady-State Approximation (ISSA), where an intercept of the linear regression of accumulated amount as a function of time is interpreted as an adsorbed amount, can be critically assessed with the transient numerical code for two cases: (i) when the total burden of metal on the cell is the input data and (ii) when an extraction procedure provides further information on the adsorbed and internalised amount.     

Diet shifts of Caribbean grunts (Haemulidae) and snappers (Lutjanidae) and the relation with nursery-to-coral reef migrations

The spatial size distribution of grunts and snappers have previously indicated the separation of juveniles in nursery habitats from the adults on the coral reef. This implies life cycle migrations from nursery habitats (such as seagrass beds and mangroves) to the coral reef. If diet shifts are related to such migrations, then the diets of these fish must change before or around the fish size at which such migrations take place. A wide size range of juveniles of two grunt species (Haemulon sciurus and Haemulon flavolineatum) and of two snapper species (Lutjanus apodus and Ocyurus chrysurus) were caught in seagrass beds and mangroves, and their gut contents identified and quantified. Regression analysis between fish size and dietary importance of small crustaceans showed a negative relationship in all four species. Positive relations were found for H. sciurus, L. apodus and O. chrysurus between fish length and the dietary importance of decapods, and for L. apodusand O. chrysurus between fish length and prey fish importance. Critical changes in the fish diets with fish size were examined by application of a Canonical Correspondence Analysis (CCA). The CCA yielded three clusters of size-classes of fishes with similar diets, and application of a Mantel test showed that each of these clusters had significantly different diets, and that each cluster diet was significantly specialised. The size at which a fish species ‘switched’ from one cluster to another was compared with size-at-maturity data and with the typical size at which these species migrate from the nursery habitats to the coral reef. H. sciurus and H. flavolineatum may be prompted to migrate from the nursery habitats to coral reef habitats because of dietary changes, or because of the development of the gonads. For L. apodus and O. chrysurus, a dietary changeover forms a more likely explanation for nursery-to-reef migrations than does sexual maturation because these species reach maturity at sizes much larger than the maximum size of individuals found in nursery habitats. Although other factors may theoretically initiate or promote the migration patterns, the results of this study indicate that ontogenetic dietary changes may crucially influence the nursery-to-coral reef migrations of these reef fish species.     

Genetic Relationships Among One Non-Endemic and Two Endemic Mediterranean Triplefin Blennies (Pisces, Blennioidei)

Abstract. Three triplefin blennies occur sympatrically in the Mediterranean Sea; Tripterygion tripteronotus and T. melanurus are endemic, whereas T. delaisi is also found in the Eastern Atlantic. Although very similar in morphology, ecology and behaviour, some striking differences exist among reproductive strategies. Several authors proposed hypotheses on the evolution of these species. In order to enhance insight into the genetic structure of this genus, an enzyme electrophoretic survey was performed. Samples of 40 to 45 specimens of each species were collected near Calvi (Corsica, France). Twenty-four enzyme loci were analysed, of which 17 appeared polymorphic. T. delaisi showed a high level of average observed heterozygosity (H_o) of 0.140, as opposed to very low levels of 0.009 and 0.021 in T. tripteronotus and T. melanurus , respectively. Nei's genetic distance was 0.21 between T. delaisi and T. tripteronotus , 0.35 between T. delaisi and T. melanurus and 0.45 between T. tripteronotus and T. melanurus . We hypothesise that divergence of these species started before the Pleistocene and that the endemic species survived the Pleistocene glaciations in refugia within the Mediterranean. The possibility that T. tripteronotus and T. delaisi diverged through sympatric speciation is discussed.     

The dimensions of sand ripples in full-scale oscillatory flows

New large-scale experiments have been carried out in two oscillatory flow tunnels to study ripple regime sand suspension and net sand transport processes in full-scale oscillatory flows. The paper focuses on ripple dimensions and the new data are combined with existing data to make a large dataset of ripple heights and lengths for flows with field-scale amplitudes and periods. A feature of the new experiments is a focus on the effect of flow irregularity. The combined dataset is analysed to examine the range of hydraulic conditions under which oscillatory flow ripples occur, to examine the effects of flow irregularity and ripple three-dimensionality on ripple dimensions and to test and improve existing methods for predicting ripple dimensions.The following are the main conclusions. (1) The highest velocities in a flow time-series play an important role in determining the type of bedform occurring in oscillatory flow. Bedform regime is well characterised by mobility number based on maximum velocity in the case of regular flow and based on the mean of the highest one tenth peak velocities in the case of irregular flow. (2) For field-scale flows, sand size is the primary factor determining whether equilibrium ripples will be 2D or 3D. 2D ripples occur when the sand D₅₀ ≥ 0.30 mm and 3D ripples occur when D₅₀ ≤ 0.22 mm (except when the flow orbital diameter is low). (3) Ripple type (2D or 3D) is the same for regular and irregular flows and ripple dimensions produced by equivalent regular and irregular flows follow a similar functional dependence on mobility number, with mobility number based on maximum velocity in the case of regular flow and based on the mean of the highest one tenth velocities in the case of irregular flow. For much of the ripple regime, ripple dimensions have weak dependency on mobility number and ripple dimensions are similar for regular and irregular flows with the same flow orbital amplitude. However, differences in ripples produced by equivalent regular and irregular flows become significant at the high mobility end of the ripple regime. (4) Ripple dimensions predicted using the Wiberg and Harris formulae are in poor agreement with measured ripple dimensions from the large-scale experiments. Predictions based on the Mogridge et al. and the Nielsen formulae show better overall agreement with the data but also show systematic differences in cases of 3D ripples and ripples generated by irregular flows. (5) Based on the combined large-scale data, modifications to the Nielsen ripple dimension equations are proposed for the heights and lengths of 2D ripples. The same equations apply to regular and irregular flows, but with mobility number appropriately defined. 3D ripples are generally smaller than 2D ripples and estimates of 3D ripple height and length may be obtained by applying multipliers of 0.55 and 0.73 respectively to the 2D formulae. The proposed modified Nielsen formulae provide an improved fit to the large-scale data, accounting for flow irregularity and ripple three-dimensionality.