THE INFRARED COLORS OF THE SUN

Solar infrared colors provide powerful constraints on the stellar effective temperature scale, but they must be measured with both accuracy and precision in order to do so. We fulfill this requirement by using line-depth ratios to derive in a model-independent way the infrared colors of the Sun, and we use the latter to test the zero point of the Casagrande et al. effective temperature scale, confirming its accuracy. Solar colors in the widely used Two Micron All Sky Survey (2MASS) JHK(s) and WISE W1-4 systems are provided: (V - J)(circle dot) = 1.198, (V - H)(circle dot) = 1.484, (V - K-s)(circle dot) = 1.560, (J - H)(circle dot) = 0.286, (J - K-s)(circle dot) = 0.362, (H - K-s)(circle dot) = 0.076, (V - W1)(circle dot) = 1.608, (V - W2)(circle dot) = 1.563, (V - W3)(circle dot) = 1.552, and (V - W4)(circle dot) = 1.604. A cross-check of the effective temperatures derived implementing 2MASS or WISE magnitudes in the infrared flux method confirms that the absolute calibration of the two systems agrees within the errors, possibly suggesting a 1% offset between the two, thus validating extant near-and mid-infrared absolute calibrations. While 2MASS magnitudes are usually well suited to derive T-eff, we find that a number of bright, solar-like stars exhibit anomalous WISE colors. In most cases, this effect is spurious and can be attributed to lower-quality measurements, although for a couple of objects (3%+/- 2% of the total sample) it might be real, and may hint at the presence of warm/hot debris disks.

U-Pb zircon in situ dating with LA-MC-ICP-MS using a mixed detector configuration

The LA-MC-ICP-MS method applied to U-Pb in situ dating is still rapidly evolving due to improvements in both lasers and ICP-MS. To test the validity and reproducibility of the method, 5 different zircon samples, including the standard Temora-2, ranging in age between 2.2 Ga and 246 Ma, were dated using both LA-MC-ICP-MS and SHRIMP. The selected zircons were dated by SHRIMP and, after gentle polishing, the laser spot was driven to the same site or on the same zircon phase with a 213 nm laser microprobe coupled to a multi-collector mixed system. The data were collected with a routine spot size of 25 μm and, in some cases, of 15 and 40 μm. A careful cross-calibration using a diluted U-Th-Pb solution to calculate the Faraday reading to counting rate conversion factors and the highly suitable GJ-1 standard zircon for external calibrations were of paramount importance for obtaining reliable results. All age results were concordant within the experimental errors. The assigned age errors using the LA-MC-ICP-MS technique were, in most cases, higher than those obtained by SHRIMP, but if we are not faced with a high resolution stratigraphy, the laser technique has certain advantages.

Updated status and performance for the cosmic origins spectrograph onboard the Hubble Space Telescope

The Cosmic Origins Spectrograph (COS) was installed on the Hubble Space Telescope in May 2009. Although COS was initially designed to perform high-sensitivity medium- and low-resolution spectroscopy of astronomical objects in the 1150-3200 Å wavelength range, new wavelength settings have recently become available that allow medium-resolution spectroscopy down to 900 Å, at effective areas comparable to those of FUSE. Here we provide an update on the implementation of the new short wavelength settings G130M/1222, 1096, and 1055. We discuss changes to the Far-Ultraviolet (FUV) and Near-Ultraviolet (NUV) dark rates, FUV pulse height filtering, new and improved flux calibrations for FUV Lifetime Positions 1 and 2, changes in sensitivity for both the NUV and FUV channels, and give a general overview of the calibration projects undertaken in Cycles 19 and 20.

Troposphere Calibration Techniques for Deep Space Probe Tracking

Ground-based Earth troposphere calibration systems play an important role in planetary exploration, especially to carry out radio science experiments aimed at the estimation of planetary gravity fields. In these experiments, the main observable is the spacecraft (S/C) range rate, measured from the Doppler shift of an electromagnetic wave transmitted from ground, received by the spacecraft and coherently retransmitted back to ground. If the solar corona and interplanetary plasma noise is already removed from Doppler data, the Earth troposphere remains one of the main error sources in tracking observables. Current Earth media calibration systems at NASA’s Deep Space Network (DSN) stations are based upon a combination of weather data and multidirectional, dual frequency GPS measurements acquired at each station complex. In order to support Cassini’s cruise radio science experiments, a new generation of media calibration systems were developed, driven by the need to achieve the goal of an end-to-end Allan deviation of the radio link in the order of 3×〖10〗^(-15) at 1000 s integration time. The future ESA’s Bepi Colombo mission to Mercury carries scientific instrumentation for radio science experiments (a Ka-band transponder and a three-axis accelerometer) which, in combination with the S/C telecommunication system (a X/X/Ka transponder) will provide the most advanced tracking system ever flown on an interplanetary probe. Current error budget for MORE (Mercury Orbiter Radioscience Experiment) allows the residual uncalibrated troposphere to contribute with a value of 8×〖10〗^(-15) to the two-way Allan deviation at 1000 s integration time. The current standard ESA/ESTRACK calibration system is based on a combination of surface meteorological measurements and mathematical algorithms, capable to reconstruct the Earth troposphere path delay, leaving an uncalibrated component of about 1-2% of the total delay. In order to satisfy the stringent MORE requirements, the short time-scale variations of the Earth troposphere water vapor content must be calibrated at ESA deep space antennas (DSA) with more precise and stable instruments (microwave radiometers). In parallel to this high performance instruments, ESA ground stations should be upgraded to media calibration systems at least capable to calibrate both troposphere path delay components (dry and wet) at sub-centimetre level, in order to reduce S/C navigation uncertainties. The natural choice is to provide a continuous troposphere calibration by processing GNSS data acquired at each complex by dual frequency receivers already installed for station location purposes. The work presented here outlines the troposphere calibration technique to support both Deep Space probe navigation and radio science experiments. After an introduction to deep space tracking techniques, observables and error sources, in Chapter 2 the troposphere path delay is widely investigated, reporting the estimation techniques and the state of the art of the ESA and NASA troposphere calibrations. Chapter 3 deals with an analysis of the status and the performances of the NASA Advanced Media Calibration (AMC) system referred to the Cassini data analysis. Chapter 4 describes the current release of a developed GNSS software (S/W) to estimate the troposphere calibration to be used for ESA S/C navigation purposes. During the development phase of the S/W a test campaign has been undertaken in order to evaluate the S/W performances. A description of the campaign and the main results are reported in Chapter 5. Chapter 6 presents a preliminary analysis of microwave radiometers to be used to support radio science experiments. The analysis has been carried out considering radiometric measurements of the ESA/ESTEC instruments installed in Cabauw (NL) and compared with the requirements of MORE. Finally, Chapter 7 summarizes the results obtained and defines some key technical aspects to be evaluated and taken into account for the development phase of future instrumentation.

Calcolo della posizione di reti di stazioni permanenti

Con il trascorrere del tempo, le reti di stazioni permanenti GNSS (Global Navigation Satellite System) divengono sempre più un valido supporto alle tecniche di rilevamento satellitare. Esse sono al tempo stesso un’efficace materializzazione del sistema di riferimento e un utile ausilio ad applicazioni di rilevamento topografico e di monitoraggio per il controllo di deformazioni. Alle ormai classiche applicazioni statiche in post-processamento, si affiancano le misure in tempo reale sempre più utilizzate e richieste dall’utenza professionale. In tutti i casi risulta molto importante la determinazione di coordinate precise per le stazioni permanenti, al punto che si è deciso di effettuarla tramite differenti ambienti di calcolo. Sono stati confrontati il Bernese, il Gamit (che condividono l’approccio differenziato) e il Gipsy (che utilizza l’approccio indifferenziato). L’uso di tre software ha reso indispensabile l’individuazione di una strategia di calcolo comune in grado di garantire che, i dati ancillari e i parametri fisici adottati, non costituiscano fonte di diversificazione tra le soluzioni ottenute. L’analisi di reti di dimensioni nazionali oppure di reti locali per lunghi intervalli di tempo, comporta il processamento di migliaia se non decine di migliaia di file; a ciò si aggiunge che, talora a causa di banali errori, oppure al fine di elaborare test scientifici, spesso risulta necessario reiterare le elaborazioni. Molte risorse sono quindi state investite nella messa a punto di procedure automatiche finalizzate, da un lato alla preparazione degli archivi e dall’altro all’analisi dei risultati e al loro confronto qualora si sia in possesso di più soluzioni. Dette procedure sono state sviluppate elaborando i dataset più significativi messi a disposizione del DISTART (Dipartimento di Ingegneria delle Strutture, dei Trasporti, delle Acque, del Rilevamento del Territorio - Università di Bologna). E’ stato così possibile, al tempo stesso, calcolare la posizione delle stazioni permanenti di alcune importanti reti locali e nazionali e confrontare taluni fra i più importanti codici scientifici che assolvono a tale funzione. Per quanto attiene il confronto fra i diversi software si è verificato che: • le soluzioni ottenute dal Bernese e da Gamit (i due software differenziati) sono sempre in perfetto accordo; • le soluzioni Gipsy (che utilizza il metodo indifferenziato) risultano, quasi sempre, leggermente più disperse rispetto a quelle degli altri software e mostrano talvolta delle apprezzabili differenze numeriche rispetto alle altre soluzioni, soprattutto per quanto attiene la coordinata Est; le differenze sono però contenute in pochi millimetri e le rette che descrivono i trend sono comunque praticamente parallele a quelle degli altri due codici; • il citato bias in Est tra Gipsy e le soluzioni differenziate, è più evidente in presenza di determinate combinazioni Antenna/Radome e sembra essere legato all’uso delle calibrazioni assolute da parte dei diversi software. E’ necessario altresì considerare che Gipsy è sensibilmente più veloce dei codici differenziati e soprattutto che, con la procedura indifferenziata, il file di ciascuna stazione di ciascun giorno, viene elaborato indipendentemente dagli altri, con evidente maggior elasticità di gestione: se si individua un errore strumentale su di una singola stazione o se si decide di aggiungere o togliere una stazione dalla rete, non risulta necessario il ricalcolo dell’intera rete. Insieme alle altre reti è stato possibile analizzare la Rete Dinamica Nazionale (RDN), non solo i 28 giorni che hanno dato luogo alla sua prima definizione, bensì anche ulteriori quattro intervalli temporali di 28 giorni, intercalati di sei mesi e che coprono quindi un intervallo temporale complessivo pari a due anni. Si è così potuto verificare che la RDN può essere utilizzata per l’inserimento in ITRF05 (International Terrestrial Reference Frame) di una qualsiasi rete regionale italiana nonostante l’intervallo temporale ancora limitato. Da un lato sono state stimate le velocità ITRF (puramente indicative e non ufficiali) delle stazioni RDN e, dall’altro, è stata effettuata una prova di inquadramento di una rete regionale in ITRF, tramite RDN, e si è verificato che non si hanno differenze apprezzabili rispetto all’inquadramento in ITRF, tramite un congruo numero di stazioni IGS/EUREF (International GNSS Service / European REference Frame, SubCommission for Europe dello International Association of Geodesy).

Towards a Reconstruction of Thermal Properties of Light Nuclei from Fusion - Evaporation reactions

This thesis work has been developed in the framework of a new experimental campaign, proposed by the NUCL-EX Collaboration (INFN III Group), in order to progress in the understanding of the statistical properties of light nuclei, at excitation energies above particle emission threshold, by measuring exclusive data from fusion-evaporation reactions. The determination of the nuclear level density in the A~20 region, the understanding of the statistical behavior of light nuclei with excitation energies ~3 A.MeV, and the measurement of observables linked to the presence of cluster structures of nuclear excited levels are the main physics goals of this work. On the theory side, the contribution to this project given by this work lies in the development of a dedicated Monte-Carlo Hauser-Feshbach code for the evaporation of the compound nucleus. The experimental part of this thesis has consisted in the participation to the measurement 12C+12C at 95 MeV beam energy, at Laboratori Nazionali di Legnaro - INFN, using the GARFIELD+Ring Counter(RCo) set-up, from the beam-time request to the data taking, data reduction, detector calibrations and data analysis. Different results of the data analysis are presented in this thesis, together with a theoretical study of the system, performed with the new statistical decay code. As a result of this work, constraints on the nuclear level density at high excitation energy for light systems ranging from C up to Mg are given. Moreover, pre-equilibrium effects, tentatively interpreted as alpha-clustering effects, are put in evidence, both in the entrance channel of the reaction and in the dissipative dynamics on the path towards thermalisation.

Characterisation and field deployment of a novel quantitative Time-of-Flight Aerosol Mass Spectrometer (ToF-AMS)

Das Time-of-Flight Aerosol Mass Spectrometer (ToF-AMS) der Firma Aerodyne ist eine Weiterentwicklung des Aerodyne Aerosolmassenspektrometers (Q-AMS). Dieses ist gut charakterisiert und kommt weltweit zum Einsatz. Beide Instrumente nutzen eine aerodynamische Linse, aerodynamische Partikelgrößenbestimmung, thermische Verdampfung und Elektronenstoß-Ionisation. Im Gegensatz zum Q-AMS, wo ein Quadrupolmassenspektrometer zur Analyse der Ionen verwendet wird, kommt beim ToF-AMS ein Flugzeit-Massenspektrometer zum Einsatz. In der vorliegenden Arbeit wird anhand von Laborexperimenten und Feldmesskampagnen gezeigt, dass das ToF-AMS zur quantitativen Messung der chemischen Zusammensetzung von Aerosolpartikeln mit hoher Zeit- und Größenauflösung geeignet ist. Zusätzlich wird ein vollständiges Schema zur ToF-AMS Datenanalyse vorgestellt, dass entwickelt wurde, um quantitative und sinnvolle Ergebnisse aus den aufgenommenen Rohdaten, sowohl von Messkampagnen als auch von Laborexperimenten, zu erhalten. Dieses Schema basiert auf den Charakterisierungsexperimenten, die im Rahmen dieser Arbeit durchgeführt wurden. Es beinhaltet Korrekturen, die angebracht werden müssen, und Kalibrationen, die durchgeführt werden müssen, um zuverlässige Ergebnisse aus den Rohdaten zu extrahieren. Beträchtliche Arbeit wurde außerdem in die Entwicklung eines zuverlässigen und benutzerfreundlichen Datenanalyseprogramms investiert. Dieses Programm kann zur automatischen und systematischen ToF-AMS Datenanalyse und –korrektur genutzt werden.

Entwicklung, Aufbau und Charakterisierung eines neuartigen Ionenfallen-Massenspektrometers für Aerosolpartikel (AIMS)

Im Rahmen der Arbeit wurde ein neuartiges Aerosol-Ionenfallen-Massenspektrometer (AIMS) aufgebaut und umfassend charakterisiert. Mit dem AIMS kann die chemische Zusammensetzung der verdampfbaren Komponente (bei etwa 600 °C) von Aerosolpartikeln quantitativ und on-line bestimmt werden. Die Durchmesser der Teilchen, die analysiert werden können, liegen zwischen etwa 30 und 500 nm. Der experimentelle Aufbau greift auf ein bereits gut charakterisiertes Einlasssystem des Aerodyne Aerosol-Massenspektrometers (AMS) zurück, das einen Partikeleinlass, bestehend aus einer kritischen Düse und einer aerodynamischen Linse, einen Verdampfer für die Aerosolteilchen und eine Elektronenstoß-Ionenquelle enthält. Das kommerzielle AMS verwendet entweder ein lineares Quadrupol-Massenfilter (Q-AMS) oder ein Flugzeit-Massenspektrometer (ToF-AMS). Im AIMS hingegen wird eine dreidimensionale Ionenfalle als Massenanalysator eingesetzt. Dadurch eröffnen sich unter anderem Möglichkeiten zur Durchführung von MSn-Studien und Ionen/Molekül-Reaktionsstudien. Das Massenspektrometer und wichtige Teile der Steuerungselektronik wurden am Max-Planck-Institut für Chemie in Mainz entworfen und hergestellt. Das AIMS wird von einem PC und einer Software, die in der Programmiersprache LabVIEW verfasst ist, gesteuert. Aufgrund seiner Kompaktheit ist das Instrument auch für den Feldeinsatz geeignet. Mit der Software Simion 7.0 wurden umfangreiche Simulationsstudien durchgeführt. Diese Studien beinhalten Simulationen zur Ermittlung der optimalen Spannungseinstellungen für den Ionentransfer von der Ionenquelle in die Ionenfalle und eine Abschätzung der Sammeleffizienz der Ionenfalle, die gut mit einem gemessenen Wert übereinstimmt. Charakterisierungsstudien zeigen einige instrumentelle Merkmale des AIMS auf. Es wurde beispielsweise ein Massenauflösungsvermögen von 807 für m/z 121 gefunden, wenn eine Analyserate von 1780 amu/s verwendet wird. Wird die Analyserate verringert, dann lässt sich das Massenauflösungsvermögen noch erheblich steigern. Bei m/z 43 kann dann ein Wert von > 1500 erzielt werden, wodurch sich Ionenfragmente wie C2H3O+ (m/z 43.0184) und C3H7+ (m/z 43.0548) voneinander trennen lassen. Der Massenbereich des AIMS lässt sich durch resonante Anregung erweitern; dies wurde bis zu einer Masse von 1000 amu getestet. Kalibrationsmessungen mit laborgenerierten Partikeln zeigen eine hervorragende Linearität zwischen gemessenen Signalstärken und erzeugten Aerosol-Massenkonzentrationen. Diese Studien belegen im Zusammenhang mit den gefundenen Nachweisgrenzen von Nitrat (0.16 μg/m³) und Sulfat (0.65 μg/m³) aus Aerosolpartikeln, dass das AIMS für quantitative Messungen von atmosphärischem Aerosol geeignet ist. Ein Vergleich zwischen dem AIMS und dem Q-AMS für Nitrat in städtischem Aerosol zeigt eine gute Übereinstimmung der gefundenen Messwerte. Für laborgenerierte Polystyren-Latexpartikel wurde eine MS/MS-Studie unter der Anwendung von collision induced dissociation (CID) durchgeführt. Das Verhältnis von Fragmentionen zu Analytionen wurde zu einem Wert von > 60% bestimmt. In der Zukunft können ähnliche MS/MS-Studien auch für atmosphärische Aerosolpartikel angewandt werden, wodurch sich neue Perspektiven für die Speziation von Aerosolbestandteilen eröffnen. Dann sollen vor allem Kondensationsprozesse, das heißt die Bildung von sekundärem Aerosol, detailliert untersucht werden.

Development of a model-based transient calibration process for diesel engine electronic control module tables – Part 1: data requirements, processing, and analysis

This is the first part of a study investigating a model-based transient calibration process for diesel engines. The motivation is to populate hundreds of parameters (which can be calibrated) in a methodical and optimum manner by using model-based optimization in conjunction with the manual process so that, relative to the manual process used by itself, a significant improvement in transient emissions and fuel consumption and a sizable reduction in calibration time and test cell requirements is achieved. Empirical transient modelling and optimization has been addressed in the second part of this work, while the required data for model training and generalization are the focus of the current work. Transient and steady-state data from a turbocharged multicylinder diesel engine have been examined from a model training perspective. A single-cylinder engine with external air-handling has been used to expand the steady-state data to encompass transient parameter space. Based on comparative model performance and differences in the non-parametric space, primarily driven by a high engine difference between exhaust and intake manifold pressures (ΔP) during transients, it has been recommended that transient emission models should be trained with transient training data. It has been shown that electronic control module (ECM) estimates of transient charge flow and the exhaust gas recirculation (EGR) fraction cannot be accurate at the high engine ΔP frequently encountered during transient operation, and that such estimates do not account for cylinder-to-cylinder variation. The effects of high engine ΔP must therefore be incorporated empirically by using transient data generated from a spectrum of transient calibrations. Specific recommendations on how to choose such calibrations, how many data to acquire, and how to specify transient segments for data acquisition have been made. Methods to process transient data to account for transport delays and sensor lags have been developed. The processed data have then been visualized using statistical means to understand transient emission formation. Two modes of transient opacity formation have been observed and described. The first mode is driven by high engine ΔP and low fresh air flowrates, while the second mode is driven by high engine ΔP and high EGR flowrates. The EGR fraction is inaccurately estimated at both modes, while EGR distribution has been shown to be present but unaccounted for by the ECM. The two modes and associated phenomena are essential to understanding why transient emission models are calibration dependent and furthermore how to choose training data that will result in good model generalization.

Age of cichlids: new dates for ancient lake fish radiations

Timing divergence events allow us to infer the conditions under which biodiversity has evolved and gain important insights into the mechanisms driving evolution. Cichlid fishes are a model system for studying speciation and adaptive radiation, yet, we have lacked reliable timescales for their evolution. Phylogenetic reconstructions are consistent with cichlid origins prior to Gondwanan landmass fragmentation 121-165 MYA, considerably earlier than the first known fossil cichlids (Eocene). We examined the timing of cichlid evolution using a relaxed molecular clock calibrated with geological estimates for the ages of 1) Gondwanan fragmentation and 2) cichlid fossils. Timescales of cichlid evolution derived from fossil-dated phylogenies of other bony fishes most closely matched those suggested by Gondwanan breakup calibrations, suggesting the Eocene origins and marine dispersal implied by the cichlid fossil record may be due to its incompleteness. Using Gondwanan calibrations, we found accumulation of genetic diversity within the radiating lineages of the African Lakes Malawi, Victoria and Barombi Mbo, and Palaeolake Makgadikgadi began around or after the time of lake basin formation. These calibrations also suggest Lake Tanganyika was colonized independently by the major radiating cichlid tribes that then began to accumulate genetic diversity thereafter. These results contrast with the widely accepted theory that diversification into major lineages took place within the Tanganyika basin. Together, this evidence suggests that ancient lake habitats have played a key role in generating and maintaining diversity within radiating lineages and also that lakes may have captured preexisting cichlid diversity from multiple sources from which adaptive radiations have evolved.

DEVELOPMENT OF A PREDICTIVE COMBUSTION MODEL OF A SPARK IGNITED ENGINE WITH GASOLINE DIRECT INJECTION, VARIABLE VALVE TIMING, DURATION AND LIFT TECHNOLOGIES

There is a need by engine manufactures for computationally efficient and accurate predictive combustion modeling tools for integration in engine simulation software for the assessment of combustion system hardware designs and early development of engine calibrations. This thesis discusses the process for the development and validation of a combustion modeling tool for Gasoline Direct Injected Spark Ignited Engine with variable valve timing, lift and duration valvetrain hardware from experimental data. Data was correlated and regressed from accepted methods for calculating the turbulent flow and flame propagation characteristics for an internal combustion engine. A non-linear regression modeling method was utilized to develop a combustion model to determine the fuel mass burn rate at multiple points during the combustion process. The computational fluid dynamic software Converge ©, was used to simulate and correlate the 3-D combustion system, port and piston geometry to the turbulent flow development within the cylinder to properly predict the experimental data turbulent flow parameters through the intake, compression and expansion processes. The engine simulation software GT-Power © is then used to determine the 1-D flow characteristics of the engine hardware being tested to correlate the regressed combustion modeling tool to experimental data to determine accuracy. The results of the combustion modeling tool show accurate trends capturing the combustion sensitivities to turbulent flow, thermodynamic and internal residual effects with changes in intake and exhaust valve timing, lift and duration.

IMPACT OF ENGINE CALIBRATION ON PM OXIDATION IN A CATALYZED PARTICULATE FILTER OVER A TRANSIENT CYCLE: A MODELLING STUDY

The combustion strategy in a diesel engine has an impact on the emissions, fuel consumption and the exhaust temperatures. The PM mass retained in the CPF is a function of NO2 and PM concentrations in addition to the exhaust temperatures and the flow rates. Thus the engine combustion strategy affects exhaust characteristics which has an impact on the CPF operation and PM mass retained and oxidized. In this report, a process has been developed to simulate the relationship between engine calibration, performance and HC and PM oxidation in the DOC and CPF respectively. Fuel Rail Pressure (FRP) and Start of Injection (SOI) sweeps were carried out at five steady state engine operating conditions. This data, along with data from a previously carried out surrogate HD-FTP cycle [1], was used to create a transfer function model which estimates the engine out emissions, flow rates, temperatures for varied FRP and SOI over a transient cycle. Four different calibrations (test cases) were considered in this study, which were simulated through the transfer function model and the DOC model [1, 2]. The DOC outputs were then input into a model which simulates the NO2 assisted and thermal PM oxidation inside a CPF. Finally, results were analyzed as to how engine calibration impacts the engine fuel consumption, HC oxidation in the DOC and the PM oxidation in the CPF. Also, active regeneration for various test cases was simulated and a comparative analysis of the fuel penalties involved was carried out.

Reproducibility of cerebral phenylalanine levels in patients with phenylketonuria determined by 1H-MR spectroscopy

The reproducibility of metabolite content determined by MR spectroscopy (MRS) is usually at best a few percent for the prominent singlets. When studying low-concentration metabolites, like phenylalanine (Phe), where tissue content can be <100 micromol/kg, better reproducibility is paramount-particularly in view of using MRS results for potential individual treatment advice. An optimized, targeted spectroscopy method was established at 1.5T and reproducibility was established in 21 patients with phenylketonuria (PKU) where three spectra were recorded in each of three independent sessions, two of which were in immediate succession to minimize physiologic variation. Intersession variation was found to be only 7 micromol/kg Phe for back-to-back repetition of sessions, in close agreement with the variation of 16 micromol/kg observed for single spectra within a session. Analysis of variance proved the individuality of the blood/brain Phe ratio-though this ratio seems to be influenced by physiologic factors that are not stable in time. The excellent reproducibility was achieved through optimization of various factors, including signal-to-noise ratio, repositioning, and prescan calibrations, but also by enforcing as much prior information as possible (e.g., lineshape and phase from reference scans, constant prior-knowledge-locked baseline). While the application of maximum general prior knowledge is a general method to reduce fluctuations, one should remember that it may introduce systematic errors.

Genetic variation and demographic history of the Haplochromis laparogramma group of Lake Victoria-An analysis based on SINEs and mitochondrial DNA

More than 500 endemic haplochromine cichlid species inhabit Lake Victoria. This striking species diversity is a classical example of recent explosive adaptive radiation thought to have happened within the last similar to 15,000 years. In this study, we examined the population structure and historical demography of 3 pelagic haplochromine cichlid species that resemble in morphology and have similar niche, Haplochromis (Yssichromis) laparogramma, Haplochromis (Y.) pyrrhocephalus, and Haplochromis (Y.) sp. "glaucocephalus". We investigated the sequences of the mitochondrial DNA control region and the insertion patterns of short interspersed elements (SINEs) of 759 individuals. We show that sympatric forms are genetically differentiated in 4 of 6 cases, but we also found apparent weakening of the genetic differentiation in areas with turbid water. We estimated the timings of population expansion and species divergence to coincide with the refilling of the lake at the Pleistocene/Holocene boundary. We also found that estimates can be altered significantly by the choice of the shape of the molecular clock. If we employ the nonlinear clock model of evolutionary rates in which the rates are higher towards the recent, the population expansion was dated at around the event of desiccation of the lake ca. 17,000 YBP. Thus, we succeeded in clarifying the species and population structure of closely related Lake Victoria cichlids and in showing the importance of applying appropriate clock calibrations in elucidating recent evolutionary events. (C) 2009 Elsevier B.V. All rights reserved.

The Accumulation of Reproductive Incompatibilities in African Cichlid Fish

The rate at which different components of reproductive isolation accumulate with divergence time between species has only been studied in a limited, but growing, number of species. We measured premating isolation and hybrid inviability at four different ontogenetic stages from zygotes to adults in interspecific hybrids of 26 pairs of African cichlid species, spanning the entire East African haplochromine radiation. We then used multiple relaxed molecular clock calibrations to translate genetic distances into absolute ages to compare evolutionary rates of different components of reproductive isolation. We find that premating isolation accumulates fast initially but then changes little with increasing genetic distance between species. In contrast, postmating isolation between closely related species is negligible but then accumulates rapidly, resulting in complete hybrid inviability after 4.4/8.5/18.4 million years (my). Thus, the rate at which complete intrinsic incompatibilities arise in this system is orders of magnitude lower than rates of speciation within individual lake radiations. Together these results suggest divergent ecological adaptations may prevent populations from interbreeding and help maintain cichlid species diversity, which may be vulnerable to environmental degradation. By quantifying the capacity to produce viable hybrids between allopatric, distantly related lineages our results also provide an upper divergence time limit for the "hybrid swarm origin" model of adaptive radiation.