Spatial isolation and genetic differentiation in naturally fragmented plant populations of the Swiss Alps

Aims The effect Of anthropogenic landscape fragmentation on the genetic diversity and adaptive potential of plant populations is a major issue in conservation biology. However, little is known about the partitioning of genetic diversity in alpine species, which occur in naturally fragmented habitats. Here, we, investigate molecular patterns of three alpine plants (Epilobium fleischeri, Geum reptans and Campanula thyrsoides) across Switzerland and ask whether Spatial isolation has led to high levels of populations differentiation, increasing over distance, and a decrease of within-population variability. We further hypothesize that file contrasting potential for long-distance dispersal (LDD) of Seed in these Species will considerably influence and explain diversity partitioning. Methods For each study species, we Sampled 20-23 individuals from each of 20-32 populations across entire Switzerland. We applied Random Amplified Polymorphic Dimorphism markers to assess genetic diversity within (Nei's expected heterozygosity, H-e; percentage of polymorphic hands, P-P) and among (analysis of molecular variance, Phi(st)) populations and correlated population size and altitude with within-populalion diversity. Spatial patterns of genetic relatedness were investigated using Mantel tests and standardized major axis regression as well as unweighted pair group method with arithmetic mean cluster analyses and Monmonier's algorithm. To avoid known biases, We standardized the numbers of populations, individuals and markers using multiple random reductions. We modelled LDD with a high alpine wind data set using the terminal velocity and height of seed release as key parameters. Additionally, we assessed a number of important life-history traits and factors that potentially influence genetic diversity partitioning (e.g. breeding system, longevity and population size). Important findings For all three species, We found a significant isolation-by-distance relationship but only a moderately high differentiation among populations (Phi(st): 22.7, 48 and 16.8%, for E. fleischeri, G. reptans and C. thyrsoides, respectively). Within-population diversity (H-c: 0.19-0.21, P-p: 62-75%) was not reduced in comparison to known results from lowland species and even small populations with < 50 reproductive individuals contained high levels of genetic diversity. We further found no indication that a high long-distance seed dispersal potential enhances genetic connectivity among populations. Gene flow seems to have a strong stochastic component causing large dissimilarity between population pairs irrespective of the spatial distance. Our results suggest that other life-history traits, especially the breeding System, may play an important role in genetic diversity partitioning. We conclude that spatial isolation in the alpine environment has a strong influence on population relatedness but that a number of factors can considerably influence the strength of this relationship.

Cellular and network mechanisms of rhythm generation in spinal cord circuits

The generation of rhythmic electrical activity is a prominent feature of spinal cord circuits that is used for locomotion and also for circuit refinement during development. The mechanisms involved in rhythm generation in spinal cord networks are not fully understood. It is for example not known whether spinal cord rhythms are driven by pacemaker neurons and if yes, which neurons are involved in this function. We studied the mechanisms involved in rhythm generation in slice cultures from fetal rats that were grown on multielectrode arrays (MEAs). We combined multisite extracellular recordings from the MEA electrodes with intracellular patch clamp recordings from single neurons. We found that spatially restricted oscillations of activity appeared in most of the cultures spontaneously. Such activity was based on intrinsic activity in a percentage of the neurons that could activate the spinal networks through recurrent excitation. The local oscillator networks critically involved NMDA, AMPA and GABA / glycine receptors at subsequent phases of the oscillation cycle. Intrinsic spiking in individual neurons (in the absence of functional synaptic coupling) was based on persistent sodium currents. Intrinsic firing as well as persistent sodium currents were increased by 5-HT through 5-HT2 receptors. Comparing neuronal activity to muscle activity in co-cultures of spinal cord slices with muscle fibers we found that a percentage of the intrinsically spiking neurons were motoneurons. These motoneurons were electrically coupled among each other and they could drive the spinal networks through cholinergic recurrent excitation. These findings open the possibility that during development rhythmic activity in motoneurons is not only involved in circuit refinement downstream at the neuromuscular endplates but also upstream at the level of spinal cord circuits.

Light curing through glass ceramics with a second- and a third-generation LED curing unit: effect of curing mode on the degree of conversion of dual-curing resin cements

Objectives The aim of this study was to measure the degree of conversion (DC) of five dual-curing resin cements after different curing modes with a second- and a third-generation light-emitting diode (LED) curing unit. Additionally, irradiance of both light curing units was measured at increasing distances and through discs of two glass ceramics for computer-aided design/manufacturing (CAD/CAM). Materials and methods Irradiance and spectra of the Elipar FreeLight 2 (Standard Mode (SM)) and of the VALO light curing unit (High Power Mode (HPM) and Xtra Power Mode (XPM)) were measured with a MARC radiometer. Irradiance was measured at increasing distances (control) and through discs (1.5 to 6 mm thickness) of IPS Empress CAD and IPS e.max CAD. DC of Panavia F2.0, RelyX Unicem 2 Automix, SpeedCEM, BisCem, and BeautiCem SA was measured with an attenuated total reflectance–Fourier transform infrared spectrometer when self-cured (negative control) or light cured in SM for 40 s, HPM for 32 s, or XPM for 18 s. Light curing was performed directly (positive control) or through discs of either 1.5- or 3-mm thickness of IPS Empress CAD or IPS e.max CAD. DC was analysed with Kruskal–Wallis tests followed by pairwise Wilcoxon rank sum tests (α = 0.05). Results Maximum irradiances were 1,545 mW/cm2 (SM), 2,179 mW/cm2 (HPM), and 4,156 mW/cm2 (XPM), and all irradiances decreased by >80 % through discs of 1.5 mm, ≥95 % through 3 mm, and up to >99 % through 6 mm. Generally, self-curing resulted in the lowest DC. For some cements, direct light curing did not result in higher DC compared to when light cured through ceramic discs. For other cements, light curing through ceramic discs of 3 mm generally reduced DC. Conclusions Light curing was favourable for dual-curing cements. Some cements were more susceptible to variations in curing mode than others. Clinical relevance When light curing a given cement, the higher irradiances of the third-generation LED curing unit resulted in similar DC compared to the second-generation one, though at shorter light curing times.

Spatial variability of recent sedimentation in Lake Ohrid (Albania/Macedonia)

Abstract. Lake Ohrid is likely of Pliocene age and thus commonly referred to as the oldest existing lake in Europe. In this study spatial variability of recent sediment composition is assessed using >50 basin wide distributed surface sediment samples. Analysis of biogeochemical bulk parameters, selected metals, pigment concentrations as well as grain size distributions revealed a significant spatial heterogeneity in surface sediment composition. It implies that sedimentation in Lake Ohrid is controlled by an interaction of multiple natural and anthropogenic factors and processes. Major factors controlling surface sediment composition are related to differences in geological catchment characteristics, anthropogenic land use, and a counterclockwise rotating surface water current. In some instances processes controlling sediment composition also seem to impact distribution patterns of biodiversity, which suggests a common interaction of processes responsible for both patterns.

Deciphering the driving forces of erosion rates on millennial to million-year timescales in glacially impacted landscapes: An example from the Western Alps

In many regions, tectonic uplift is the main driver of erosion over million-year (Myr) timescales, but climate changes can markedly affect the link between tectonics and erosion, causing transient variations in erosion rates. Here we study the driving forces of millennial to Myr-scale erosion rates in the French Western Alps, as estimated from in situ produced cosmogenic 10Be and a newly developed approach integrating detrital and bedrock apatite fission-track thermochronology. Millennial erosion rates from 10Be analyses vary between ~0.27 and ~1.33 m/kyr, similar to rates measured in adjacent areas of the Alps. Significant positive correlations of millennial erosion rates with geomorphic measures, in particular with the LGM ice thickness, reveal a strong transient morphological and erosional perturbation caused by repeated Quaternary glaciations. The perturbation appears independent of Myr-scale uplift and erosion gradients, with the effect that millennial erosion rates exceed Myr-scale erosion rates only in the internal Alps where the latter are low (<0.4 km/Myr). These areas, moreover, exhibit channels that clearly plot above a general linear positive relation between Myr-scale erosion rates and normalized steepness index. Glacial erosion acts irrespective of rock uplift and thus not only leads to an overall increase in erosion rates but also regulates landscape morphology and erosion rates in regions with considerable spatial gradients in Myr-scale tectonic uplift. Our study demonstrates that climate change, e.g., through occurrence of major glaciations, can markedly perturb landscape morphology and related millennial erosion rate patterns, even in regions where Myr-scale erosion rates are dominantly controlled by tectonics.

Verb generation in children with spina bifida.

We investigated verb generation in children with spina bifida meningomyelocele (SBM; n = 55) and in typically developing controls (n = 32). Participants completed 6 blocks (40 trials each) of a task requiring them to produce a semantically related verb in response to a target noun and an additional 40 trials on which they were simply required to read target nouns aloud. After controlling for reading response time, groups did not differ significantly in verb generation response time or learning. Children with SBM produced more non-verb errors than controls and tended to repeat their mistakes over blocks. Verb generation performance was associated with brain volume measures in participants with SBM. Congenital cerebellar dysmorphology is associated with impaired performance in verb generation accuracy, although not with increased response times to produce verbs

Geochemical synthesis for the Effingen Member in boreholes at Oftringen, Gösgen and Küttigen

The Effingen Member is a low-permeability rock unit of Oxfordian age (ca. 160 Ma) that occurs across northern Switzerland. It comprises sandy calcareous marls and (argillaceous) limestones. This report describes the hydrogeochemistry, mineralogy and supporting physical properties of the Effingen Member in three boreholes in the Jura-Südfuss area: Oftringen, Gösgen and Küttigen, where it is 220–240 m thick. The top of the Effingen Member is at 420, 66 and 32 m depths at the three sites. Core materials are available from Oftringen and Gösgen, whereas information from Küttigen is limited to cuttings, in-situ hydrogeological testing and geophysical logging. Hydrogeological boundaries of the Effingen Member vary between locations. Ground-water flows were identified during drilling at the top (Geissberg Member), but not at the base, of the Effingen Member at Oftringen, at the base (Hauptrogenstein Formation) of the Effingen Member at Gösgen, and in a limestone layer (Gerstenhübel unit) within the Effingen Member at Küttigen. The marls and limestones of the Effingen Member have carbonate contents of 46–91 wt.-% and clay-mineral contents of 5–37 wt.-%. Pyrite contents are up to 1.6 wt.-%, but no sulphate minerals were detected by routine analyses. Clay minerals are predominantly mixed-layer illite-smectite, illite and kaolinite, with sporadic traces of chlorite and smectite. Veins filled with calcite ± celestite occur through the Effingen Member at Oftringen but not at Gösgen or Küttigen. They formed at 50–70 ºC from externally derived fluids, probably of Miocene age. Water contents are 0.7–4.2 wt.-%, corresponding to a water-loss porosity range of 1.9–10.8 vol.-%. Specific surface areas, measured by the BET method, are 2–30 m2/g, correlating with clay-mineral contents. Water activity has been measured and yielded surprisingly low values down to 0.8. These cannot be explained by pore-water salinity alone and include other effects, such as changes in the fabric due to stress release or partial saturation. Observed variations in measurements are not fully understood. Cation exchange capacity (CEC) and exchangeable cation populations have been studied by the Ni-en method. CEC, derived from the consumption of the index cation Ni, is 9–99 meq/kgrock at a solid:liquid ratio of 1, correlating with the clay-mineral content. Cation concentrations in Ni-en extract solutions are in the order Na+≥Ca2+>Mg2+>K+>Sr2+. However, the analytical results from the Ni-en extractions have additional contributions from cations originating from pore water and from mineral dissolution reactions that occurred during extraction, and it was not possible to reliably quantify these contributions. Therefore, in-situ cation populations and selectivity coefficients could not be derived. A suite of methods have been used for characterising the chemical compositions of pore waters in the Effingen Member. Advective displacement was used on one sample from each Oftringen and Gösgen and is the only method that produces results that approach complete hydrochemical compositions. Aqueous extraction was used on core samples from these two boreholes and gives data only for Cl- and, in some cases, Br-. Out-diffusion was used on core samples from Oftringen and similarly gives data for Cl- and Br- only. For both aqueous extraction and out-diffusion, reaction of the experimental water with rock affected concentrations of cations, SO42 and alkalinity in experimental solutions. Another method, centrifugation, failed to extract pore water. Stable isotope ratios (δ18O and δ2H) of pore waters in core samples from Oftringen were analysed by the diffusive exchange method and helium contents of pore water in Oftringen samples were extracted for mass spectrometric analysis by quantitative outgassing of preserved core samples. Several lines of evidence indicate that drillcore samples might not have been fully saturated when opened and subsampled in the laboratory. These include comparisons of water-loss porosities with physical porosities, water-activity measurements, and high contents of dissolved gas as inferred from ground-water samples. There is no clear proof of partial saturation and it is unclear whether this might represent in-situ conditions or is due to exsolution of gas due to the pressure release since drilling. Partial saturation would have no impact on the recalculation of pore-water compositions from aqueous extraction experiments using water-loss porosity data. The largest uncertainty in the pore-water Cl- concentrations recalculated from aqueous extraction and out-diffusion experiments is the magnitude of the anion-accessible fraction of water-loss porosity. General experience of clay-mineral rich formations suggests that the anion-accessible porosity fraction is very often about 0.5 and generally in a range of 0.3 to 0.6 and tends to be inversely correlated with clay-mineral contents. Comparisons of the Cl- concentration in pore water obtained by advective displacement with that recalculated from aqueous extraction of an adjacent core sample suggests a fraction of 0.27 for an Oftringen sample, whereas the same procedure for a Gösgen sample suggests a value of 0.64. The former value for anion-accessible porosity fraction is presumed to be unrepresentative given the local mineralogical heterogeneity at that depth. Through-diffusion experiments with HTO and 36Cl- suggest that the anion-accessible porosity fraction in the Effingen Member at Oftringen and Gösgen is around 0.5. This value is proposed as a typical average for rocks of the Effingen Member, bearing in mind that it varies on a local scale in response to the heterogeneity of lithology and pore-space architecture. The substantial uncertainties associated with the approaches to estimating anion-accessible porosity propagate into the calculated values of in-situ pore-water Cl- concentrations. On the basis of aqueous extraction experiments, and using an anion-accessible porosity fraction of 0.5, Cl- concentrations in the Effingen Member at Oftringen reach a maximum of about 14 g/L in the centre. Cl- decreases upwards and downwards from that, forming a curved depth profile. Cl- contents in the Effingen Member at Gösgen increase with depth from about 3.5 g/L to about 14 g/L at the base of the cored profile (which corresponds to the centre of the formation). Out-diffusion experiments were carried out on four samples from Oftringen, distributed through the Effingen Member. Recalculated Cl- concentrations are similar to those from aqueous extraction for 3 out of the 4 samples, and somewhat lower for one sample. Concentrations of other components, i.e. Na+, K+, Ca2+, Mg2+, Sr2+, SO42- and HCO3- cannot be obtained from the aqueous extraction and out-diffusion experimental data because of mineral dissolution and cation exchange reactions during the experiments. Pore-water pH also is not constrained by those extraction experiments. The only experimental approach to obtain complete pore-water compositions for samples from Oftringen and Gösgen is advective displacement of pore water. The sample from Oftringen used for this experiment is from 445 m depth in the upper part of the Effingen Member and gave eluate with 16.5 g/L Cl- whereas aqueous extraction from a nearby sample indicated about 9 g/L Cl-. The sample from Gösgen used for advective displacement is from 123 m depth in the centre of the Effingen Member sequence and gave eluate with about 9 g/L Cl- whereas aqueous extraction gave 11.5 g/L Cl-. In both cases the pore waters have Na-(Ca)-Cl compositions and SO42- concentrations of about 1.1 g/L. The Gösgen sample has a Br/Cl ratio similar to that of sea water, whereas this ratio is lower for the Oftringen sample. Taking account of uncertainties in the applied experimental approaches, it is reasonable to place an upper limit of ca. 20 g/L on Cl- concentration for pore water in the Effingen Member in this area. There are major discrepancies between pore-water SO42- concentrations inferred from aqueous extraction or out-diffusion experiments and those obtained from advective displacement in both the Oftringen and Gösgen cases. A general conclusion is that all or at least part of the discrepancies are attributable to perturbation of the sulphur system and enhancement of SO42- by sulphate mineral dissolution and possibly minor pyrite oxidation during aqueous extraction and out-diffusion. Therefore, data for SO42- calculated from those pore-water sampling methods are considered not to be representative of in-situ conditions. A reference pore-water composition was defined for the Effingen Member in the Jura Südfuss area. It represents the probable upper limits of Cl- contents and corresponding anion and cation concentrations that are reasonably constrained by experimental data. Except for Cl- and possibly Na+ concentrations, this composition is poorly constrained especially with respect to SO42- and Ca2+ concentrations, and pH and alkalinity. Stable isotope compositions, δ18O and δ2H, of pore waters in the Effingen Member at Oftringen plot to the right of the meteoric water line, suggesting that 18O has been enriched by water-rock exchange, which indicates that the pore waters have a long residence time. A long residence time of pore water is supported by the level of dissolved 4He that has accumulated in pore water of the Effingen Member at Oftringen. This is comparable with, or slightly higher than, the amounts of 4He in the Opalinus Clay at Benken. Ground waters were sampled from flowing zones intersected by boreholes at the three locations. The general interpretation is that pore waters and ground-water solutes may have similar origins in Mesozoic and Cenozoic brackish-marine formations waters, but ground-water solutes have been diluted rather more than pore waters by ingress of Tertiary and Quaternary meteoric waters. The available hydrochemical data for pore waters from the Effingen Member at these three locations in the Jura-Südfuss area suggest that the geochemical system evolved slowly over geological periods of time, in which diffusion was an important mechanism of solute transport. The irregularity of Cl- and δ18O profiles and spatial variability of advective ground-water flows in the Malm-Dogger system suggests that palaeohydrogeological and hydrochemical responses to changing tectonic and surface environmental conditions were complex.

Transport of 234U in the Opalinus Clay on centimetre to decimetre scales

The Opalinus Clay formation in North Switzerland is a potential host rock for a deep underground radioactive waste repository. The distribution of U-238, U-234 and Th-230 was studied in rock samples of the Opalinus Clay from an exploratory borehole at Benken (Canton of Zurich) using MC-ICP-MS. The aim of U-234 was to assess the in situ, long-term migration behaviour in this rock. Very low hydraulic conductivities of the Opalinus Clay, reducing potential of the pore water and its chemical equilibrium with the host rock are expected to render both U-238 and Th-230 immobile. If U is heterogeneously distributed in the Opalinus Clay, gradients in the supply of U-234 from the rock matrix to the pore water by the decay of U-238 will be established. Diffusive redistribution separates U-234 from its immobile parent U-238 resulting in bulk rock U-234/U-238 activity disequilibria. These may provide a means of estimating the mobility of U-234 in the rock if the diffusion rate of U-234 is significant compared to its decay rate. Sampling was carried out on two scales. Drilling of cm-spaced samples from the drill-core was done to study mobility over short distances and elucidate possible small-scale lithological control. Homogenized 25-cm-long portions of a 2-m-long drill-core section were prepared to provide information on transport over a longer distance. Variations in U and/or Th content on the cm-scale between clays and carbonate-sandy layers are revealed by beta-scanning, which shows that the (dominant) clay is richer in both elements. Samples were digested using aqua regia followed by total HF dissolution, yielding two fractions. in all studied samples U was found to be concentrated in the HF digestion fraction. It has a high U/Th ratio and a study by SEM-EDS points to sub-mu m up to several mu m in size zircon grains as the main U-rich phase. This fraction consistently has U-234/U-238 activity ratios below unity. The minute zircon grains constitute the major reservoir of U in the rock and act as constant rate suppliers of U-234 into the rock matrix and the pore water. The aqua regia leach fraction was found to be enriched in Th, and complementary to the HF fraction, having U-234/U-238 activity ratios above unity. It is believed that these U activity ratios reflect the surplus of having U-234 delivered from the zircon grains. Some cm-spaced samples show bulk rock U-234/U-238 activity ratios that are markedly out of equilibrium. In most of them a striking negative correlation between the total U content and the bulk rock U-234/U-238 activity ratios is observed. This is interpreted to indicate net U-234 transfer from regions of higher supply of U-234 towards those of lower supply which is, in most cases, equivalent to transfer from clayey towards carbonate/sandy portions of the rock. In contrast, the 25 cm averaged samples all have uniform bulk rock U-234/U-238 activity ratios in equilibrium, indicating U immobility in the last 1-1.5 Ma on this spatial scale. It is concluded that the small-scale lithological variations which govern U spatial distribution in the Opalinus Clay are the major factor determining U-234 in situ supply rates, regulating its diffusive fluxes and controlling the observed bulk rock U-234/U-238 activity ratios. A simple box-model is presented to simulate the measured bulk rock U-234/U-238 activity ratios and to give an additional insight into the studied system. (C) 2008 Elsevier Ltd. All rights reserved.

Precipitation over the past four centuries in the Dieshan Mountains as inferred from tree rings: An introduction to an HHT-based method

To improve our understanding of the Asian monsoon system, we developed a hydroclimate reconstruction in a marginal monsoon shoulder region for the period prior to the industrial era. Here, we present the first moisture sensitive tree-ring chronology, spanning 501 years for the Dieshan Mountain area, a boundary region of the Asian summer monsoon in the northeastern Tibetan Plateau. This reconstruction was derived from 101 cores of 68 old-growth Chinese pine (Pinus tabulaeformis) trees. We introduce a Hilbert–Huang Transform (HHT) based standardization method to develop the tree-ring chronology, which has the advantages of excluding non-climatic disturbances in individual tree-ring series. Based on the reliable portion of the chronology, we reconstructed the annual (prior July to current June) precipitation history since 1637 for the Dieshan Mountain area and were able to explain 41.3% of the variance. The extremely dry years in this reconstruction were also found in historical documents and are also associated with El Niño episodes. Dry periods were reconstructed for 1718–1725, 1766–1770 and 1920–1933, whereas 1782–1788 and 1979–1985 were wet periods. The spatial signatures of these events were supported by data from other marginal regions of the Asian summer monsoon. Over the past four centuries, out-of-phase relationships between hydroclimate variations in the Dieshan Mountain area and far western Mongolia were observed during the 1718–1725 and 1766–1770 dry periods and the 1979–1985 wet period.

Filling the Eastern European gap in millennium-long temperature reconstructions

Tree ring–based temperature reconstructions form the scientific backbone of the current global change debate. Although some European records extend into medieval times, high-resolution, long-term, regional-scale paleoclimatic evidence is missing for the eastern part of the continent. Here we compile 545 samples of living trees and historical timbers from the greater Tatra region to reconstruct interannual to centennial-long variations in Eastern European May–June temperature back to 1040 AD. Recent anthropogenic warming exceeds the range of past natural climate variability. Increased plague outbreaks and political conflicts, as well as decreased settlement activities, coincided with temperature depressions. The Black Death in the mid-14th century, the Thirty Years War in the early 17th century, and the French Invasion of Russia in the early 19th century all occurred during the coldest episodes of the last millennium. A comparison with summer temperature reconstructions from Scandinavia, the Alps, and the Pyrenees emphasizes the seasonal and spatial specificity of our results, questioning those large-scale reconstructions that simply average individual sites.

Spatial and Temporal Melt Variability at Helheim Glacier, East Greenland, and Its Effect on Ice Dynamics

Understanding the behavior of large outlet glaciers draining the Greenland Ice Sheet is critical for assessing the impact of climate change on sea level rise. The flow of marine-terminating outlet glaciers is partly governed by calving-related processes taking place at the terminus but is also influenced by the drainage of surface runoff to the bed through moulins, cracks, and other pathways. To investigate the extent of the latter effect, we develop a distributed surface-energy-balance model for Helheim Glacier, East Greenland, to calculate surface melt and thereby estimate runoff. The model is driven by data from an automatic weather station operated on the glacier during the summers of 2007 and 2008, and calibrated with independent measurements of ablation. Modeled melt varies over the deployment period by as much as 68% relative to the mean, with melt rates approximately 77% higher on the lower reaches of the glacier trunk than on the upper glacier. We compare melt variations during the summer season to estimates of surface velocity derived from global positioning system surveys. Near the front of the glacier, there is a significant correlation (on >95% levels) between variations in runoff (estimated from surface melt) and variations in velocity, with a 1 day delay in velocity relative to melt. Although the velocity changes are small compared to accelerations previously observed following some calving events, our findings suggest that the flow speed of Helheim Glacier is sensitive to changes in runoff. The response is most significant in the heavily crevassed, fast-moving region near the calving front. The delay in the peak of the cross-correlation function implies a transit time of 12-36 h for surface runoff to reach the bed.

Model-Based Estimates of Calanus Finmarchicus Abundance in the Gulf of Maine

Ocean observing systems and satellites routinely collect a wealth of information on physical conditions in the ocean. With few exceptions, such as chlorophyll concentrations, information on biological properties is harder to measure autonomously. Here, we present a system to produce estimates of the distribution and abundance of the copepod Calanus finmarchicus in the Gulf of Maine. Our system uses satellite-based measurements of sea surface temperature and chlorophyll concentration to determine the developmental and reproductive rates of C. finmarchicus. The rate information then drives a population dynamics model of C. finmarchicus that is embedded in a 2-dimensional circulation field. The first generation of this system produces realistic information on interannual variability in C. finmarchicus distribution and abundance during the winter and spring. The model can also be used to identify key drivers of interannual variability in C. finmarchicus. Experiments with the model suggest that changes in initial conditions are overwhelmed by variability in growth rates after approximately 50 d. Temperature has the largest effect on growth rate. Elevated chlorophyll during the late winter can lead to increased C. finmarchicus abundance during the spring, but the effect of variations in chlorophyll concentrations is secondary to the other inputs. Our system could be used to provide real-time estimates or even forecasts of C. finmarchicus distribution. These estimates could then be used to support management of copepod predators such as herring and right whales.

Regulation of Phytoplankton Carbon to Chlorophyll Ratio By Light, Nutrients and Temperature in the Equatorial Pacific Ocean: A Basin-Scale Model

The complex effects of light, nutrients and temperature lead to a variable carbon to chlorophyll (C:Chl) ratio in phytoplankton cells. Using field data collected in the Equatorial Pacific, we derived a new dynamic model with a non-steady C:Chl ratio as a function of irradiance, nitrate, iron, and temperature. The dynamic model is implemented into a basin-scale ocean circulation-biogeochemistry model and tested in the Equatorial Pacific Ocean. The model reproduces well the general features of phytoplankton dynamics in this region. For instance, the simulated deep chlorophyll maximum (DCM) is much deeper in the western warm pool (similar to 100 m) than in the Eastern Equatorial Pacific (similar to 50 m). The model also shows the ability to reproduce chlorophyll, including not only the zonal, meridional and vertical variations, but also the interannual variability. This modeling study demonstrates that combination of nitrate and iron regulates the spatial and temporal variations in the phytoplankton C:Chl ratio in the Equatorial Pacific. Sensitivity simulations suggest that nitrate is mainly responsible for the high C:Chl ratio in the western warm pool while iron is responsible for the frontal features in the C:Chl ratio between the warm pool and the upwelling region. In addition, iron plays a dominant role in regulating the spatial and temporal variations of the C:Chl ratio in the Central and Eastern Equatorial Pacific. While temperature has a relatively small effect on the C:Chl ratio, light is primarily responsible for the vertical decrease of phytoplankton C:Chl ratio in the euphotic zone.