Application of low temperature thermochronology to understanding cooling and exhumation in the Central Swiss Alps

THESIS ABSTRACT : Low-temperature thermochronology relies on application of radioisotopic systems whose closure temperatures are below temperatures at which the dated phases are formed. In that sense, the results are interpreted as "cooling ages" in contrast to "formation ages". Owing to the low closure-temperatures, it is possible to reconstruct exhumation and cooling paths of rocks during their residence at shallow levels of the crust, i.e. within first ~10 km of depth. Processes occurring at these shallow depths such as final exhumation, faulting and relief formation are fundamental for evolution of the mountain belts. This thesis aims at reconstructing the tectono-thermal history of the Aar massif in the Central Swiss Alps by means of zircon (U-Th)/He, apatite (U-Th)/He and apatite fission track thermochronology. The strategy involved acquisition of a large number of samples from a wide range of elevations in the deeply incised Lötschen valley and a nearby NEAT tunnel. This unique location allowed to precisely constrain timing, amount and mechanisms of exhumation of the main orographic feature of the Central Alps, evaluate the role of topography on the thermochronological record and test the impact of hydrothermal activity. Samples were collected from altitudes ranging between 650 and 3930 m and were grouped into five vertical profiles on the surface and one horizontal in the tunnel. Where possible, all three radiometric systems were applied to each sample. Zircon (U-Th)/He ages range from 5.1 to 9.4 Ma and are generally positively correlated with altitude. Age-elevation plots reveal a distinct break in slope, which translates into exhumation rate increasing from ~0.4 to ~3 km/Ma at 6 Ma. This acceleration is independently confirmed by increased cooling rates on the order of 100°C/Ma constrained on the basis of age differences between the zircon (U-Th)/He and the remaining systems. Apatite fission track data also plot on a steep age-elevation curve indicating rapid exhumation until the end of the Miocene. The 6 Ma event is interpreted as reflecting tectonically driven uplift of the Aar massif. The late Miocene timing implies that the increase of precipitation in the Pliocene did not trigger rapid exhumation in the Aar massif. The Messinian salinity crisis in the Mediterranean could not directly intensify erosion of the Aar but associated erosional output from the entire Alps may have tapered the orogenic wedge and caused reactivation of thrusting in the Aar massif. The high exhumation rates in the Messinian were followed by a decrease to ~1.3 km/Ma as evidenced by ~8 km of exhumation during last 6 Ma. The slowing of exhumation is also apparent from apatite (U-Th)1He age-elevation data in the northern part of the Lötschen valley where they plot on a ~0.5km/Ma line and range from 2.4 to 6.4 Ma However, from the apatite (U-Th)/He and fission track data from the NEAT tunnel, there is an indication of a perturbation of the record. The apatite ages are youngest under the axis of the valley, in contrast to an expected pattern where they would be youngest in the deepest sections of the tunnel due to heat advection into ridges. The valley however, developed in relatively soft schists while the ridges are built of solid granitoids. In line with hydrological observations from the tunnel, we suggest that the relatively permeable rocks under the valley floor, served as conduits of geothermal fluids that caused reheating leading to partial Helium loss and fission track annealing in apatites. In consequence, apatite ages from the lowermost samples are too young and the calculated exhumation rates may underestimate true values. This study demonstrated that high-density sampling is indispensable to provide meaningful thermochronological data in the Alpine setting. The multi-system approach allows verifying plausibility of the data and highlighting sources of perturbation. RÉSUMÉ DE THÈSE : La thermochronologie de basse température dépend de l'utilisation de systèmes radiométriques dont la température de fermeture est nettement inférieure à la température de cristallisation du minéral. Les résultats obtenus sont par conséquent interprétés comme des âges de refroidissement qui diffèrent des âges de formation obtenus par le biais d'autres systèmes de datation. Grâce aux températures de refroidissement basses, il est aisé de reconstruire les chemins de refroidissement et d'exhumation des roches lors de leur résidence dans la croute superficielle (jusqu'à 10 km). Les processus qui entrent en jeu à ces faibles profondeurs tels que l'exhumation finale, la fracturation et le faillage ainsi que la formation du relief sont fondamentaux dans l'évolution des chaînes de montagne. Ces dernières années, il est devenu clair que l'enregistrement thermochronologique dans les orogènes peut être influencé par le relief et réinitialisé par l'advection de la chaleur liée à la circulation de fluides géothermaux après le refroidissement initial. L'objectif de cette thèse est de reconstruire l'histoire tectono-thermique du massif de l'Aar dans les Alpes suisses Centrales à l'aide de trois thermochronomètres; (U-Th)/He sur zircon, (U-Th)/He sur apatite et les traces de fission sur apatite. Afin d'atteindre cet objectif, nous avons récolté un grand nombre d'échantillons provenant de différentes altitudes dans la vallée fortement incisée de Lötschental ainsi que du tunnel de NEAT. Cette stratégie d'échantillonnage nous a permis de contraindre de manière précise la chronologie, les quantités et les mécanismes d'exhumation de cette zone des Alpes Centrales, d'évaluer le rôle de la topographie sur l'enregistrement thermochronologique et de tester l'impact de l'hydrothermalisme sur les géochronomètres. Les échantillons ont été prélevés à des altitudes comprises entre 650 et 3930m selon 5 profils verticaux en surface et un dans le tunnel. Quand cela à été possible, les trois systèmes radiométriques ont été appliqués aux échantillons. Les âges (U-Th)\He obtenus sur zircons sont compris entre 5.l et 9.4 Ma et sont corrélés de manière positive avec l'altitude. Les graphiques représentant l'âge et l'élévation montrent une nette rupture de la pente qui traduisent un accroissement de la vitesse d'exhumation de 0.4 à 3 km\Ma il y a 6 Ma. Cette accélération de l'exhumation est confirmée par les vitesses de refroidissement de l'ordre de 100°C\Ma obtenus à partir des différents âges sur zircons et à partir des autres systèmes géochronologiques. Les données obtenues par traces de fission sur apatite nous indiquent également une exhumation rapide jusqu'à la fin du Miocène. Nous interprétons cet évènement à 6 Ma comme étant lié à l'uplift tectonique du massif de l'Aar. Le fait que cet évènement soit tardi-miocène implique qu'une augmentation des précipitations au Pliocène n'a pas engendré cette exhumation rapide du massif de l'Aar. La crise Messinienne de la mer méditerranée n'a pas pu avoir une incidence directe sur l'érosion du massif de l'Aar mais l'érosion associée à ce phénomène à pu réduire le coin orogénique alpin et causer la réactivation des chevauchements du massif de l'Aar. L'exhumation rapide Miocène a été suivie pas une diminution des taux d'exhumation lors des derniers 6 Ma (jusqu'à 1.3 km\Ma). Cependant, les âges (U-Th)\He sur apatite ainsi que les traces de fission sur apatite des échantillons du tunnel enregistrent une perturbation de l'enregistrement décrit ci-dessus. Les âges obtenus sur les apatites sont sensiblement plus jeunes sous l'axe de la vallée en comparaison du profil d'âges attendus. En effet, on attendrait des âges plus jeunes sous les parties les plus profondes du tunnel à cause de l'advection de la chaleur dans les flancs de la vallée. La vallée est creusée dans des schistes alors que les flancs de celle-ci sont constitués de granitoïdes plus durs. En accord avec les observations hydrologiques du tunnel, nous suggérons que la perméabilité élevée des roches sous l'axe de la vallée à permi l'infiltration de fluides géothermaux qui a généré un réchauffement des roches. Ce réchauffement aurait donc induit une perte d'Hélium et un recuit des traces de fission dans les apatites. Ceci résulterait en un rajeunissement des âges apatite et en une sous-estimation des vitesses d'exhumation sous l'axe de la vallée. Cette étude à servi à démontrer la nécessité d'un échantillonnage fin et précis afin d'apporter des données thermochronologiques de qualité dans le contexte alpin. Cette approche multi-système nous a permi de contrôler la pertinence des données acquises ainsi que d'identifier les sources possibles d'erreurs lors d'études thermochronologiques. RÉSUMÉ LARGE PUBLIC Lors d'une orogenèse, les roches subissent un cycle comprenant une subduction, de la déformation, du métamorphisme et, finalement, un retour à la surface (ou exhumation). L'exhumation résulte de la déformation au sein de la zone de collision, menant à un raccourcissement et un apaissessement de l'édifice rocheux, qui se traduit par une remontée des roches, création d'une topographie et érosion. Puisque l'érosion agit comme un racloir sur la partie supérieure de l'édifice, des tentatives de corrélation entre les épisodes d'exhumation rapide et les périodes d'érosion intensive, dues aux changements climatiques, ont été effectuées. La connaissance de la chronologie et du lieu précis est d'une importance capitale pour une quelconque reconstruction de l'évolution d'une chaîne de montagne. Ces critères sont donnés par un retraçage des changements de la température de la roche en fonction du temps, nous donnant le taux de refroidissement. L'instant auquel les roches ont refroidit, passant une certaine température, est contraint par l'application de techniques de datation par radiométrie. Ces méthodes reposent sur la désintégration des isotopes radiogéniques, tels que l'uranium et le potassium, tous deux abondants dans les roches de la croûte terrestre. Les produits de cette désintégration ne sont pas retenus dans les minéraux hôtes jusqu'au moment du refroidissement de la roche sous une température appelée 'de fermeture' , spécifique à chaque système de datation. Par exemple, la désintégration radioactive des atomes d'uranium et de thorium produit des atomes d'hélium qui s'échappent d'un cristal de zircon à des températures supérieures à 200°C. En mesurant la teneur en uranium-parent, l'hélium accumulé et en connaissant le taux de désintégration, il est possible de calculer à quel moment la roche échantillonnée est passée sous la température de 200°C. Si le gradient géothermal est connu, les températures de fermeture peuvent être converties en profondeurs actuelles (p. ex. 200°C ≈ 7km), et le taux de refroidissement en taux d'exhumation. De plus, en datant par système radiométrique des échantillons espacés verticalement, il est possible de contraindre directement le taux d'exhumation de la section échantillonnée en observant les différences d'âges entre des échantillons voisins. Dans les Alpes suisses, le massif de l'Aar forme une structure orographique majeure. Avec des altitudes supérieures à 4000m et un relief spectaculaire de plus de 2000m, le massif domine la partie centrale de la chaîne de montagne. Les roches aujourd'hui exposées à la surface ont été enfouies à plus de 10 km de profond il y a 20 Ma, mais la topographie actuelle du massif de l'Aar semble surtout s'être développée par un soulèvement actif depuis quelques millions d'années, c'est-à-dire depuis le Néogène supérieur. Cette période comprend un changement climatique soudain ayant touché l'Europe il y a environ 5 Ma et qui a occasionné de fortes précipitations, entraînant certainement une augmentation de l'érosion et accélérant l'exhumation des Alpes. Dans cette étude, nous avons employé le système de datation (U-TH)/He sur zircon, dont la température de fermeture de 200°C est suffisamment basse pour caractériser l'exhumation du Néogène sup. /Pliocène. Les échantillons proviennent du Lötschental et du tunnel ferroviaire le plus profond du monde (NEAT) situé dans la partie ouest du massif de l'Aar. Considérés dans l'ensemble, ces échantillons se répartissent sur un dénivelé de 3000m et des âges de 5.1 à 9.4 Ma. Les échantillons d'altitude supérieure (et donc plus vieux) documentent un taux d'exhumation de 0.4 km/Ma jusqu'à il y a 6 Ma, alors que les échantillons situés les plus bas ont des âges similaires allant de 6 à 5.4 Ma, donnant un taux jusqu'à 3km /Ma. Ces données montrent une accélération dramatique de l'exhumation du massif de l'Aar il y a 6 Ma. L'exhumation miocène sup. du massif prédate donc le changement climatique Pliocène. Cependant, lors de la crise de salinité d'il y a 6-5.3 Ma (Messinien), le niveau de la mer Méditerranée est descendu de 3km. Un tel abaissement de la surface d'érosion peut avoir accéléré l'exhumation des Alpes, mais le bassin sud alpin était trop loin du massif de l'Aar pour influencer son érosion. Nous arrivons à la conclusion que la datation (U-Th)/He permet de contraindre précisément la chronologie et l'exhumation du massif de l'Aar. Concernant la dualité tectonique-érosion, nous suggérons que, dans le cas du massif de l'Aar, la tectonique prédomine.

Temperature-dependent elasticity of microtubules

Central to the biological function of microtubules is their ability to modify their length which occurs by addition and removal of subunits at the ends of the polymer, both in vivo and in vitro. This dynamic behavior is strongly influenced by temperature. Here, we show that the lateral interaction between tubulin subunits forming microtubule is strongly temperature dependent. Microtubules deposited on prefabricated substrates were deformed in an atomic force microscope during imaging, in two different experimental geometries. Microtubules were modeled as anisotropic, with the Young's modulus corresponding to the resistance of protofilaments to stretching and the shear modulus describing the weak interaction between the protofilaments. Measurements involving radial compression of microtubules deposited on flat mica confirm that microtubule elasticity depends on the temperature. Bending measurements performed on microtubules deposited on lithographically fabricated substrates show that this temperature dependence is due to changing shear modulus, implying that the lateral interaction between the protofilaments is strongly determined by the temperature. These measurements are in good agreement with previously reported measurements of the disassembly rate of microtubules, demonstrating that the mechanical and dynamic properties of microtubules are closely related.

Comparative winter thermoregulation and body temperature in three sympatric Apodemus species (A-alpicola, A-flavicollis, and A-sylvaticus)

When living in sympatry with Apodemus sylvaticus and A. flavicollis, A. alpicola dominates numerically at higher altitudes. A more efficient winter thermal isolation or a higher winter thermogenic capacity procuring a physiological advantage could explain at least part of this domination. We therefore measured body temperature (Tb), oxygen consumption (VO2), wet minimal thermal conductance (C) and non shivering thermogenesis (NST) at different ambient temperatures (Ta) on winter acclimated mice of the three species, and this for the first time in A. alpicola. NST was high and C low in the three species. No significant difference could be noticed either in Tb between 5 and -10 degrees C, in VO2 measurements at a Ta of -10 degrees C or in C. The NST measurements represent, respectively, 135.2% for A. sylvaticus, 142.8% for A. flavicollis and 140.5% for A. alpicola of the expected values, the values for A. sylvaticus being significantly lower than for the other two species. The basal metabolic rates (BMR) represent 169.4% for A. sylvaticus, 161.6% for A. flavicollis and 138.3% for A. alpicola of the expected values. Having removed the effect of body weight, the BMR value was significantly lower in A. alpicola than in A. flavicollis, but no difference could be noticed between A. sylvaticus and the other two species. In conclusion, the three species of mice have very similar acclimated thermoregulatory characteristics, well adapted to cold ambient conditions. One discriminating and advantageous factor could be the lower basal metabolic rate measured in A. alpicola compared to the other two species.

Linking temperature estimates and microstructures in deformed polymineralic mantle rocks

To constrain deformation temperatures of mantle shear zones, we studied a strike-slip shear zone (Hilti massif, Semail ophiolite, Oman) and focused on the interaction between microstructural mechanisms and chemical equilibration processes. Quantitative microfabric analysis on harzburgites with different deformation intensity (porphyroclastic tectonite, mylonite, and ultramylonite) was combined with orthopyroxene geothermometry. The average grain size of all phases decreases with decreasing shear zone thickness. Dynamic recrystallization of porphyroclasts in combination with dissolution-precipitation and nucleation result in small-sized, chemically equilibrated pyroxenes. The composition of orthopyroxene was used to calculate deformation temperatures. In the case of the porphyroclastic tectonites, the chemical composition of orthopyroxene has been reset by diffusion yielding temperature estimates of 880-900 degrees C. The mylonites were deformed by dislocation creep of olivine and show a broad range of calculated temperatures, which result from a combination of grain size reduction and inheritance of equilibrium compositions from earlier high-temperature events and diffusion. In mylonites, diffusion profiles combined with geothermometry and grain size analysis indicate a mylonitic deformation temperature of 800-900 degrees C possibly followed by diffusion. In ultramylonites, the smallest grains (<30 mu m) reveal equilibration at temperatures of similar to 700 degrees C during the last stages of ductile deformation, which was dominated by diffusion creep of olivine. Our results provide a crucial link between temperature and evolution of microstructures from dislocation creep to diffusion creep in mantle shear zones.

Evaluating thermal treeline indicators based on air and soil temperature using an air-to-soil temperature transfer model

In recent research, both soil (root-zone) and air temperature have been used as predictors for the treeline position worldwide. In this study, we intended to (a) test the proposed temperature limitation at the treeline, and (b) investigate effects of season length for both heat sum and mean temperature variables in the Swiss Alps. As soil temperature data are available for a limited number of sites only, we developed an air-to-soil transfer model (ASTRAMO). The air-to-soil transfer model predicts daily mean root-zone temperatures (10cm below the surface) at the treeline exclusively from daily mean air temperatures. The model using calibrated air and root-zone temperature measurements at nine treeline sites in the Swiss Alps incorporates time lags to account for the damping effect between air and soil temperatures as well as the temporal autocorrelations typical for such chronological data sets. Based on the measured and modeled root-zone temperatures we analyzed. the suitability of the thermal treeline indicators seasonal mean and degree-days to describe the Alpine treeline position. The root-zone indicators were then compared to the respective indicators based on measured air temperatures, with all indicators calculated for two different indicator period lengths. For both temperature types (root-zone and air) and both indicator periods, seasonal mean temperature was the indicator with the lowest variation across all treeline sites. The resulting indicator values were 7.0 degrees C +/- 0.4 SD (short indicator period), respectively 7.1 degrees C +/- 0.5 SD (long indicator period) for root-zone temperature, and 8.0 degrees C +/- 0.6 SD (short indicator period), respectively 8.8 degrees C +/- 0.8 SD (long indicator period) for air temperature. Generally, a higher variation was found for all air based treeline indicators when compared to the root-zone temperature indicators. Despite this, we showed that treeline indicators calculated from both air and root-zone temperatures can be used to describe the Alpine treeline position.

Functional diversity decreases with temperature in high elevation ant fauna

1. Severe environmental conditions filter community species compositions, forming clines of functional diversity along environmental gradients. Here, the changes in functional diversity in ant assemblages with severe environmental conditions in the Swiss Alps were investigated. 2. Eight sites were sampled along an elevation gradient (1800-2550 m). The variation in functional diversity was analysed along an elevation gradient considering four traits: social structure (monogynous vs. polygynous), worker size, pupal development, and nest structure. 3. Ant species richness and functional diversity decreased with decreasing temperature. Species found in colder habitats tended to live in subterranean nests rather than in mounds and exhibit a polymorphism in queen number, either within or across populations. The phylogenetic diversity did not decrease at colder temperature: Formicinae and Myrmicinae occupied the full range of elevations investigated. 4. An insulation experiment indicated that mounds are more thermally insulated against the cold compared with soil. The absence of a mound-building ant from high elevations probably results from a reduction in the amount of vegetal materials provided by coniferous trees. 5. More severe abiotic conditions at higher elevations act as a filter on ant assemblages, directly through physiological tolerances to the abiotic conditions and indirectly as the vegetation necessary for nest building shifts with elevation.</list-item

High elevation Plantago lanceolata plants are less resistant to herbivory than their low elevation conspecifics: is it just temperature?

Traits that mediate species interactions are evolutionarily shaped by biotic and abiotic drivers, yet we know relatively little about the relative importance of these factors. Herbivore pressure, along with resource availability and third-party' mutualists, are hypothesized to play a major role in the evolution of plant defence traits. Here, we used the model system Plantago lanceolata, which grows along steep elevation gradients in the Swiss Alps, to investigate the effect of elevation, herbivore pressure, mycorrhizal inoculation and temperature on plant resistance. Over a 1200 m elevation gradient, the levels of herbivory and iridoid glycosides (IGs) declined with increasing elevation. By planting seedlings at three different elevations, we further showed that both low-elevation growing conditions and mycorrhizal inoculation resulted in increased plant resistance to herbivores. Finally, using a temperature-controlled experiment comparing high- and low-elevation ecotypes, we showed that high-elevation ecotypes are less resistant to herbivory, and that lower temperatures impair IGs deployment after herbivore attack. We thus propose that both lower herbivore pressure, and colder temperatures relax the defense syndrome of high elevation plants.

A pseudo-spectral method for the simulation of poro-elastic seismic wave propagation in 2D polar coordinates using domain decomposition

We present a novel numerical approach for the comprehensive, flexible, and accurate simulation of poro-elastic wave propagation in 2D polar coordinates. An important application of this method and its extensions will be the modeling of complex seismic wave phenomena in fluid-filled boreholes, which represents a major, and as of yet largely unresolved, computational problem in exploration geophysics. In view of this, we consider a numerical mesh, which can be arbitrarily heterogeneous, consisting of two or more concentric rings representing the fluid in the center and the surrounding porous medium. The spatial discretization is based on a Chebyshev expansion in the radial direction and a Fourier expansion in the azimuthal direction and a Runge-Kutta integration scheme for the time evolution. A domain decomposition method is used to match the fluid-solid boundary conditions based on the method of characteristics. This multi-domain approach allows for significant reductions of the number of grid points in the azimuthal direction for the inner grid domain and thus for corresponding increases of the time step and enhancements of computational efficiency. The viability and accuracy of the proposed method has been rigorously tested and verified through comparisons with analytical solutions as well as with the results obtained with a corresponding, previously published, and independently bench-marked solution for 2D Cartesian coordinates. Finally, the proposed numerical solution also satisfies the reciprocity theorem, which indicates that the inherent singularity associated with the origin of the polar coordinate system is adequately handled.

Flowering date of taxonomic families predicts phenological sensitivity to temperature: Implications for forecasting the effects of climate change on unstudied taxa.

PREMISE OF THE STUDY: Numerous long-term studies in seasonal habitats have tracked interannual variation in first flowering date (FFD) in relation to climate, documenting the effect of warming on the FFD of many species. Despite these efforts, long-term phenological observations are still lacking for many species. If we could forecast responses based on taxonomic affinity, however, then we could leverage existing data to predict the climate-related phenological shifts of many taxa not yet studied. METHODS: We examined phenological time series of 1226 species occurrences (1031 unique species in 119 families) across seven sites in North America and England to determine whether family membership (or family mean FFD) predicts the sensitivity of FFD to standardized interannual changes in temperature and precipitation during seasonal periods before flowering and whether families differ significantly in the direction of their phenological shifts. KEY RESULTS: Patterns observed among species within and across sites are mirrored among family means across sites; early-flowering families advance their FFD in response to warming more than late-flowering families. By contrast, we found no consistent relationships among taxa between mean FFD and sensitivity to precipitation as measured here. CONCLUSIONS: Family membership can be used to identify taxa of high and low sensitivity to temperature within the seasonal, temperate zone plant communities analyzed here. The high sensitivity of early-flowering families (and the absence of early-flowering families not sensitive to temperature) may reflect plasticity in flowering time, which may be adaptive in environments where early-season conditions are highly variable among years.

Chemical labelling of oyster shells used for time-calibrated high-resolution Mg/Ca ratios: A tool for estimation of past seasonal temperature variations

The geochemical compositions of biogenic carbonates are increasingly used for palaeoenvironmental reconstructions. The skeletal delta O-18 temperature relationship is dependent on water salinity, so many recent studies have focused on the Mg/Ca and Sr/Ca ratios because those ratios in water do not change significantly on short time scales. Thus, those elemental ratios are considered to be good palaeotemperature proxies in many biominerals, although their use remains ambiguous in bivalve shells. Here, we present the high-resolution Mg/Ca ratios of two modern species of juvenile and adult oyster shells, Crassostrea gigas and Ostrea edulis. These specimens were grown in controlled conditions for over one year in two different locations. In situ monthly Mn-marking of the shells has been used for day calibration. The daily Mg/Ca.ratios in the shell have been measured with an electron microprobe. The high frequency Mg/Ca variation of all specimens displays good synchronism with lunar cycles, suggesting that tides strongly influence the incorporation of Mg/Ca into the shells. Highly significant correlation coefficients (0.70<R<0.83, p<0.0001) between the Mg/Ca ratios and the seawater temperature are obtained only for juvenile C. gigas samples, while metabolic control of Mg/Ca incorporation and lower shell growth rates preclude the use of the Mg/Ca ratio in adult shells as a palaeothermometer. Data from three juvenile C. gigas shells from the two study sites are selected to establish a relationship: T = 3.77Mg/Ca + 1.88, where T is in degrees C and Mg/Ca in mmol/mol. (c) 2012 Elsevier B.V. All rights reserved.

Temperature-induced sex reversal is not responsible for sex-ratio distortions in grayling Thymallus thymallus or brown trout Salmo trutta.

On the basis of the experiments carried out over various years, it was concluded that (1) grayling Thymallus thymallus and brown trout Salmo trutta are resistant to temperature-induced sex reversal at ecologically relevant temperatures, (2) environmental sex reversal is unlikely to cause the persistent sex ratio distortion observed in at least one of the study populations and (3) sex-specific tolerance of temperature-related stress may be the cause of distorted sex ratios in populations of T. thymallus or S. trutta.

Role of the circadian clock gene Per2 in adaptation to cold temperature.

Adaptive thermogenesis allows mammals to resist to cold. For instance, in brown adipose tissue (BAT) the facultative uncoupling of the proton gradient from ATP synthesis in mitochondria is used to generate systemic heat. However, this system necessitates an increase of the Uncoupling protein 1 (Ucp1) and its activation by free fatty acids. Here we show that mice without functional Period2 (Per2) were cold sensitive because their adaptive thermogenesis system was less efficient. Upon cold-exposure, Heat shock factor 1 (HSF1) induced Per2 in the BAT. Subsequently, PER2 as a co-activator of PPARα increased expression of Ucp1. PER2 also increased Fatty acid binding protein 3 (Fabp3), a protein important to transport free fatty acids from the plasma to mitochondria to activate UCP1. Hence, in BAT PER2 is important for the coordination of the molecular response of mice exposed to cold by synchronizing UCP1 expression and its activation.

Broad-scale adaptive genetic variation in alpine plants is driven by temperature and precipitation.

Identifying adaptive genetic variation is a challenging task, in particular in non-model species for which genomic information is still limited or absent. Here, we studied distribution patterns of amplified fragment length polymorphisms (AFLPs) in response to environmental variation, in 13 alpine plant species consistently sampled across the entire European Alps. Multiple linear regressions were performed between AFLP allele frequencies per site as dependent variables and two categories of independent variables, namely Moran's eigenvector map MEM variables (to account for spatial and unaccounted environmental variation, and historical demographic processes) and environmental variables. These associations allowed the identification of 153 loci of ecological relevance. Univariate regressions between allele frequency and each environmental factor further showed that loci of ecological relevance were mainly correlated with MEM variables. We found that precipitation and temperature were the best environmental predictors, whereas topographic factors were rarely involved in environmental associations. Climatic factors, subject to rapid variation as a result of the current global warming, are known to strongly influence the fate of alpine plants. Our study shows, for the first time for a large number of species, that the same environmental variables are drivers of plant adaptation at the scale of a whole biome, here the European Alps.