Benthic foraminifera of late Quaternary sediments in the northern North Atlantic

Four long sediment cores from locations in the Framstrait, the Norwegian-Greenland Seas and the northern North Atlantic were analysed in a high resolution sampling mode (1 - 2 cm density) for their benthic foraminiferal content. In particular the impact of the intense climatic changes at glacial/interglacial transitions (terminations I and II) on the benthic community have been of special interest. The faunal data were investigated by means of multivariate analysis and represented in their chronological occurence. The most prominent species of benthic foraminifera in the Norwegian-Greenland Seas are Oridorsalis umbonatus, Cibicidoides wuellerstorfi, the group of Cassidulina, Pyrgo rotalaria, Globocassidulina subglobosa and fragmented tubes of arenaceous species. The climatic signal of termination I as well as termination II is recorded in the fossil foraminiferal tests as divided transition from glacial to interglacial. The elder INDAR maximum (individuals accumulation rate = individuals/sq cm * 1.000 y; Norwegian-Greenland Seas: average 3.000 - 6.000 individuals/sq cm * 1.000 y; northern North Atlantic: average 150 individuals/sq cm * 1.000 y) is followed by a period of decreased values. The second, younger maximum reaches comparable values as the elder maximum. The interglacial INDAR are in average 700 individuals/sq cm * 1.000 y in the Norwegian-Greenland Seas and 200 individuals/sq cm * 1.000 y in average in the northern North Atlantic. The occurence of the elder INDAR maximum shows a distinct chronological transgressivity between the northern North Atlantic (12.400 ybp.) and the Framstrait (8.900 ybp.). The time shift from south to north amounts 3.500 yrs., the average expanding velocity 0,78 km per year. Within the Norwegian-Greenland Seas the average expanding velocity amounts 0,48 km per year. This chronological transgressivity is interpreted as impact of the progressive expanding of the North Atlantic and the Norwegian Current during the deglaciation. The dynamic of the faunal development is defined as increasing INDAR per time. The elder INDAR maximum shows in both glacial/interglacial transitions an exponential increase from south to north. Termination II is characterized by a general higher dynamic as termination I. By means of the high resolution sampling density the impact of regional isotopic recognized melt-water events is recognized by an increase of endobenthic and t-ubiquitous species in the Norwegian-Greenland Seas sediments. During termination I the relative minimum between both INDAR maxima occur chronological with an decrease of calculated sea surface temperatures. This is interpreted as indication of the close pelagic - benthic coupling. The climatic signal in the northern North Atlantic recorded in the fossil benthic foraminiferal community shows a lower amplitude as in the Norwegian-Greenland Seas. The occurence of the epibenthic Cibicidoides wuellersforfi allows to evaluate the variability of the bottom water mass. In general at all core locations increasing lateral bottom currents are recognized with the occurence of the second younger INDAR maximum. In comparison with various paleo-climatological data sets fossil benthic foraminifers show a distinct koherence with changes of the atmospheric temperatures, the SSTs and the postglacial sea level increase. The benthic foraminiferal fauna is bound indirectly on and indicative for regional climatic changes, but principal dependent upon global climatic changes.

Searching for convergence in phylogenetic Markov chain Monte Carlo

Markov chain Monte Carlo (MCMC) is a methodology that is gaining widespread use in the phylogenetics community and is central to phylogenetic software packages such as MrBayes. An important issue for users of MCMC methods is how to select appropriate values for adjustable parameters such as the length of the Markov chain or chains, the sampling density, the proposal mechanism, and, if Metropolis-coupled MCMC is being used, the number of heated chains and their temperatures. Although some parameter settings have been examined in detail in the literature, others are frequently chosen with more regard to computational time or personal experience with other data sets. Such choices may lead to inadequate sampling of tree space or an inefficient use of computational resources. We performed a detailed study of convergence and mixing for 70 randomly selected, putatively orthologous protein sets with different sizes and taxonomic compositions. Replicated runs from multiple random starting points permit a more rigorous assessment of convergence, and we developed two novel statistics, delta and epsilon, for this purpose. Although likelihood values invariably stabilized quickly, adequate sampling of the posterior distribution of tree topologies took considerably longer. Our results suggest that multimodality is common for data sets with 30 or more taxa and that this results in slow convergence and mixing. However, we also found that the pragmatic approach of combining data from several short, replicated runs into a metachain to estimate bipartition posterior probabilities provided good approximations, and that such estimates were no worse in approximating a reference posterior distribution than those obtained using a single long run of the same length as the metachain. Precision appears to be best when heated Markov chains have low temperatures, whereas chains with high temperatures appear to sample trees with high posterior probabilities only rarely. [Bayesian phylogenetic inference; heating parameter; Markov chain Monte Carlo; replicated chains.]

The behaviour of the rare earth elements during estuarine mixing-revisited

Complete rare earth element (except Eu) and Y concentrations from the estuarine mixing zone (salinity =0.2 to 33) of Elimbah Creek, Queensland, Australia, were measured by quadrupole ICP-MS without preconcentration. High sampling density in the low salinity regime along with high quality data allow accurate tracing of the development of the typical marine rare earth element anomalies as well as Y/Ho fractionation. Over the entire estuary, the rare earth elements are strongly removed relative to a freshwater endmember (60-80% removal). This large overall removal occurs despite a strong remineralisation peak (190% for La, 130% for Y relative to the freshwater endmember) in the mid-salinity zone. Removal and remineralisation are accompanied by fractionation of the original (freshwater) rare earth element pattern, resulting in light rare earth element depletion. Estuarine fractionation generates a large positive La anomaly and a superchondritic Y/Ho ratio. Conversely, we observe no evidence to support the generation of the negative Ce anomaly in the estuary. With the exception of Ce, the typical marine rare earth element features can thus be attributed to estuarine mixing processes. The persistence of these features in hydrogenous sediments for at least 3.71 Ga highlights the importance of estuarine processes for marine chemistry on geological timescales. (c) 2005 Elsevier B.V. All rights reserved.

Contrast sensitivity for complex and random gratings

This thesis studied the effect of (i) the number of grating components and (ii) parameter randomisation on root-mean-square (r.m.s.) contrast sensitivity and spatial integration. The effectiveness of spatial integration without external spatial noise depended on the number of equally spaced orientation components in the sum of gratings. The critical area marking the saturation of spatial integration was found to decrease when the number of components increased from 1 to 5-6 but increased again at 8-16 components. The critical area behaved similarly as a function of the number of grating components when stimuli consisted of 3, 6 or 16 components with different orientations and/or phases embedded in spatial noise. Spatial integration seemed to depend on the global Fourier structure of the stimulus. Spatial integration was similar for sums of two vertical cosine or sine gratings with various Michelson contrasts in noise. The critical area for a grating sum was found to be a sum of logarithmic critical areas for the component gratings weighted by their relative Michelson contrasts. The human visual system was modelled as a simple image processor where the visual stimuli is first low-pass filtered by the optical modulation transfer function of the human eye and secondly high-pass filtered, up to the spatial cut-off frequency determined by the lowest neural sampling density, by the neural modulation transfer function of the visual pathways. The internal noise is then added before signal interpretation occurs in the brain. The detection is mediated by a local spatially windowed matched filter. The model was extended to include complex stimuli and its applicability to the data was found to be successful. The shape of spatial integration function was similar for non-randomised and randomised simple and complex gratings. However, orientation and/or phase randomised reduced r.m.s contrast sensitivity by a factor of 2. The effect of parameter randomisation on spatial integration was modelled under the assumption that human observers change the observer strategy from cross-correlation (i.e., a matched filter) to auto-correlation detection when uncertainty is introduced to the task. The model described the data accurately.

Denoising Manifolds for Dimension

2000 Mathematics Subject Classification: 68T01, 62H30, 32C09.

Undersampled Critical Branching Processes on Small-World and Random Networks Fail to Reproduce the Statistics of Spike Avalanches

The power-law size distributions obtained experimentally for neuronal avalanches are an important evidence of criticality in the brain. This evidence is supported by the fact that a critical branching process exhibits the same exponent t~3=2. Models at criticality have been employed to mimic avalanche propagation and explain the statistics observed experimentally. However, a crucial aspect of neuronal recordings has been almost completely neglected in the models: undersampling. While in a typical multielectrode array hundreds of neurons are recorded, in the same area of neuronal tissue tens of thousands of neurons can be found. Here we investigate the consequences of undersampling in models with three different topologies (two-dimensional, small-world and random network) and three different dynamical regimes (subcritical, critical and supercritical). We found that undersampling modifies avalanche size distributions, extinguishing the power laws observed in critical systems. Distributions from subcritical systems are also modified, but the shape of the undersampled distributions is more similar to that of a fully sampled system. Undersampled supercritical systems can recover the general characteristics of the fully sampled version, provided that enough neurons are measured. Undersampling in two-dimensional and small-world networks leads to similar effects, while the random network is insensitive to sampling density due to the lack of a well-defined neighborhood. We conjecture that neuronal avalanches recorded from local field potentials avoid undersampling effects due to the nature of this signal, but the same does not hold for spike avalanches. We conclude that undersampled branching-process-like models in these topologies fail to reproduce the statistics of spike avalanches.

Undersampled Critical Branching Processes on Small-World and Random Networks Fail to Reproduce the Statistics of Spike Avalanches

The power-law size distributions obtained experimentally for neuronal avalanches are an important evidence of criticality in the brain. This evidence is supported by the fact that a critical branching process exhibits the same exponent t~3=2. Models at criticality have been employed to mimic avalanche propagation and explain the statistics observed experimentally. However, a crucial aspect of neuronal recordings has been almost completely neglected in the models: undersampling. While in a typical multielectrode array hundreds of neurons are recorded, in the same area of neuronal tissue tens of thousands of neurons can be found. Here we investigate the consequences of undersampling in models with three different topologies (two-dimensional, small-world and random network) and three different dynamical regimes (subcritical, critical and supercritical). We found that undersampling modifies avalanche size distributions, extinguishing the power laws observed in critical systems. Distributions from subcritical systems are also modified, but the shape of the undersampled distributions is more similar to that of a fully sampled system. Undersampled supercritical systems can recover the general characteristics of the fully sampled version, provided that enough neurons are measured. Undersampling in two-dimensional and small-world networks leads to similar effects, while the random network is insensitive to sampling density due to the lack of a well-defined neighborhood. We conjecture that neuronal avalanches recorded from local field potentials avoid undersampling effects due to the nature of this signal, but the same does not hold for spike avalanches. We conclude that undersampled branching-process-like models in these topologies fail to reproduce the statistics of spike avalanches.

Biologia do caranguejo-uçá, Ucides cordatus (Crustacea: Decapoda: Brachyura), na Baía da Babitonga, Santa Catarina, Brasil

O objetivo deste estudo é caracterizar pela primeira vez alguns aspectos da reprodução do caranguejo-uçá em manguezais da Baía da Babitonga (Santa Catarina). Além disso, a densidade e o tamanho do estoque deste recurso pesqueiro foram também estimados. Os exemplares foram coletados mensalmente, de maio de 2002 a abril de 2003, em duas áreas distintas: Iperoba e Palmital; um total de 2265 espécimes (1623 machos e 642 fêmeas) foi analisado. Os machos com gônadas maturas foram registrados durante todo o ano, enquanto as fêmeas com gônadas maturas ocorreram em apenas cinco meses. As fêmeas ovígeras foram registradas apenas em dezembro e janeiro. O etograma do fenômeno de migração reprodutiva (andada) esteve em concordância com a maior atividade de caranguejos associada às luas cheias e novas, com maior intensidade em dezembro e janeiro, relacionados ao verão austral. A densidade total no Manguezal de Iperoba foi de 2,05 ± 0,97 ind./m², não diferindo significativamente daquela registrada para o Manguezal do Palmital (2,06 ± 1,08 ind./m²) (p < 0,05). A média global para a estimativa de densidade na Baia da Babitonga foi de 2,05 ± 1,00 ind./m², correspondendo a 1,42 ± 0,89 ind./m² com base nas galerias abertas e 0,64 ± 0,63 ind./m² para as galerias fechadas.

Low-density koala (Phascolarctos cinereus) populations in the mulgalands of south-west Queensland. I. Faecal pellet sampling protocol

Koala (Phascolarctos cinereus) populations in eastern Australia are threatened by land clearing for agricultural and urban development. At the same time, conservation efforts are hindered by a dearth of information about inland populations. Faecal deposits offer a source of information that is readily available and easily collected non-invasively. We detail a faecal pellet sampling protocol that was developed for use in a large rangeland biogeographic region. The method samples trees in belt transects, uses a thorough search at the tree base to quickly identify trees with koala pellets under them, then estimates the abundance of faecal pellets under those trees using 1-m(2) quadrats. There was a strong linear relationship between these estimates and a complete enumeration of pellet abundance under the same trees. We evaluated the accuracy of our method in detecting trees where pellets were present by means of a misclassification index that was weighed more heavily for missed trees that had high numbers of pellets under them. This showed acceptable accuracy in all landforms except riverine, where some trees with large numbers of pellets were missed. Here, accuracy in detecting pellet presence was improved by sampling with quadrats, rather than basal searches. Finally, we developed a method to reliably age pellets and demonstrate how this protocol could be used with the faecal-standing-crop method to derive a regional estimate of absolute koala abundance.

Uncertainties in annual riverine phosphorus load estimation: Impact of load estimation methodology, sampling frequency, baseflow index and catchment population density

Models developed to identify the rates and origins of nutrient export from land to stream require an accurate assessment of the nutrient load present in the water body in order to calibrate model parameters and structure. These data are rarely available at a representative scale and in an appropriate chemical form except in research catchments. Observational errors associated with nutrient load estimates based on these data lead to a high degree of uncertainty in modelling and nutrient budgeting studies. Here, daily paired instantaneous P and flow data for 17 UK research catchments covering a total of 39 water years (WY) have been used to explore the nature and extent of the observational error associated with nutrient flux estimates based on partial fractions and infrequent sampling. The daily records were artificially decimated to create 7 stratified sampling records, 7 weekly records, and 30 monthly records from each WY and catchment. These were used to evaluate the impact of sampling frequency on load estimate uncertainty. The analysis underlines the high uncertainty of load estimates based on monthly data and individual P fractions rather than total P. Catchments with a high baseflow index and/or low population density were found to return a lower RMSE on load estimates when sampled infrequently than those with a tow baseflow index and high population density. Catchment size was not shown to be important, though a limitation of this study is that daily records may fail to capture the full range of P export behaviour in smaller catchments with flashy hydrographs, leading to an underestimate of uncertainty in Load estimates for such catchments. Further analysis of sub-daily records is needed to investigate this fully. Here, recommendations are given on load estimation methodologies for different catchment types sampled at different frequencies, and the ways in which this analysis can be used to identify observational error and uncertainty for model calibration and nutrient budgeting studies. (c) 2006 Elsevier B.V. All rights reserved.

Sampling series for determination of probability density of CDMA multiple access interference

Consideration is given to a standard CDMA system and determination of the density function of the interference with and without Gaussian noise using sampling theory concepts. The formula derived provides fast and accurate results and is a simple, useful alternative to other methods

A sampling device to quantify offspring release of sessile marine invertebrates

Quantifying the rate of propagule release is of most importance to estimate reproductive output of natural populations, but simple methods to obtain such data are seldom reported. We designed and tested an inexpensive apparatus capable of reliably measure the release of gametes, eggs or larvae of sessile marine invertebrates in vertical walls. A population of the acom barnacle Chthamalus bisinuatus was sampled with this trap over 68d to obtain a time series of naupliar release. An apparent semilunar trend is shown, indicating the effectiveness of this sampling method.