Quantifying rates of landscape evolution and tectonic processes by thermochronology and numerical modeling of crustal heat transport using PECUBE

PECUBE is a three-dimensional thermal-kinematic code capable of solving the heat production-diffusion-advection equation under a temporally varying surface boundary condition. It was initially developed to assess the effects of time-varying surface topography (relief) on low-temperature thermochronological datasets. Thermochronometric ages are predicted by tracking the time-temperature histories of rock-particles ending up at the surface and by combining these with various age-prediction models. In the decade since its inception, the PECUBE code has been under continuous development as its use became wider and addressed different tectonic-geomorphic problems. This paper describes several major recent improvements in the code, including its integration with an inverse-modeling package based on the Neighborhood Algorithm, the incorporation of fault-controlled kinematics, several different ways to address topographic and drainage change through time, the ability to predict subsurface (tunnel or borehole) data, prediction of detrital thermochronology data and a method to compare these with observations, and the coupling with landscape-evolution (or surface-process) models. Each new development is described together with one or several applications, so that the reader and potential user can clearly assess and make use of the capabilities of PECUBE. We end with describing some developments that are currently underway or should take place in the foreseeable future. (C) 2012 Elsevier B.V. All rights reserved.

Triassic stratigraphic evolution of the Arabian-Greater India embayment of the southern Tethys margin

An exceptional, tectonically remarkably unaffected, nearly 200 m-thick continuous section of hemipelagic and turbiditic sediments, covering most of the Triassic is described from the Batain Complex of north-eastern Oman. According to conodont and radiolarian data the sequence spans the late Scythian to the early Norian, a time period of nearly 30 M. Coupled with a high resolution stratigraphy, the lithostratigraphy, sedimentology, as well as sequence and isotope stratigraphy of the section are documented. For the Triassic of the Batain Plain we propose the new name Sal Formation, which replaces the formerly used Matbat Formation, and subdivide it into three new members. The Sal Formation was deposited on the proximal continental margin of northeastern Arabia and records various depositional environments. The lower member is interpreted as the distal part of a homoclinal ramp which evolves to a distally steepened ramp during time of deposition of the middle member. The upper member displays a toe of slope position which is indicated by an increase of proximal turbidites. These sediments form part of a segment of the Neo-Tethyan embayment between Arabia and India. The stratigraphic analysis indicates highly varying sedimentation rates from a minimum of 2 m/M gamma around the Anisian/Ladinian boundary up to 15 m/M gamma during the Lower and Upper Triassic. Sequence-stratigraphically, the Sal section is subdivided into six third order cycles which are biochronologically well integrated into the global Triassic cycle chart. The mixed siliciclastic-calcareous upper member of the Sal Formation typically shows highstand related carbonate shedding. It is, therefore, an important test case for sequence-stratigraphic controlled carbonate export to mixed basin fills. The well developed sequence stratigraphic cycles are mirrored in the isotope patterns. Additionally, the carbon and oxygen isotope data from the Sal Formation record the same chemostratigraphic marker at the Spathian/Anisian boundary known from other Tethyan sections.

Carbon isotope excursions and microfacies changes in marine Permian-Triassic boundary sections in Hungary

Several Permian-Triassic boundary sections occur in various structural units within Hungary. These sections represent different facies zones of the western Palaeotethys margin. The Gardony core in the NE part of the Transdanubian Range typically represents the inner ramp, while the Balvany section in the Bukk Mountains of northern Hungary represents an outer ramp setting. The two sections have different patterns for their delta(13)C values. The Balvany section shows a continuous change towards more negative delta(13)C values starting at the first biotic decline, followed by a sharp, quasi-symmetric negative peak at the second decline. The appearance of the delta(13)C peak has no relationship to the lithology and occurs within a shale with low overall carbonate content, indicating that the peak is not related to diagenesis or other secondary influences. Instead, the shift and the peak reflect primary processes related to changes in environmental conditions. The continuous shift in delta(13)C values is most probably related to a decrease in bioproductivity, whereas the sharp peak can be attributed to an addition of C strongly depleted in (13)C to the ocean-atmosphere system. The most plausible model is a massive release of methane-hydrate. The quasi-symmetric pattern suggests a rapid warming-cooling cycle or physical unroofing of sediments through slope-failure and releasing methane-hydrate. The Gidony-1 core shows a continuous negative delta(13)C shift starting below the P-T boundary. However, the detailed analyses revealed a sharp delta(13)C peak in the boundary interval, just below the major biotic decline, although its magnitude doesn't reach that observed in the Balvany section. Based on careful textural examination and high-resolution stable isotope microanalyses we suggest that the suppression of the delta(13)C peak that is common in the oolitic boundary sections is due to combined effects of condensed sedimentation, sediment reworking and erosion, as well as perhaps diagenesis. (c) 2005 Elsevier B.V All rights reserved.

Analysis of a boundary protein delimiting chromatin domains in yeast

SUMMARY The expression state of a eukaryotic gene depends in part on its location in the chromosome. This position effect results from the organization of eukaryotic genomes into discrete functional domains, defined by local differences in chromatin structure. The expression of genes within each domain appears to be defined and maintained by the concerted action of regulatory elements such as promoters, enhancers, silencers and locus control regions. Individual domains may be bordered by boundary elements that separate regions of permissive and silent chromatin. When located next to chromosomal elements such as telomeres, genes can be subjected to epigenetic silencing. In yeast, this is mediated by the propagation of the SIR proteins from telomeres towards more centromeric regions. Particular transcription factors can protect downstream genes from silencing when tethered between the gene and the telomere, and they may thus act as chromatin domain boundaries. Here we have studied one of these transcription factors, CTF-1, that binds directly histone H3. A deletion mutagenesis localized the barrier activity to CTF-1 histone-binding domain. A saturating point mutagenesis of this domain identified several amino-acid substitutions that similarly inhibited the boundary and histone-binding activities. Chromatin immunoprecipitation experiments indicated that the barrier protein efficiently prevents the spreading of SIR proteins, and that it separates domains of hypoacetylated and hyperacetylated histones. Together, these results suggest a mechanism by which proteins such as CTF-1 may interact directly with histone H3 to prevent the propagation of a silent chromatin structure, thereby defining boundaries of permissive and silent chromatin domains. RESUME L'expression des gènes eucaryotes dépend en partie de leur localisation sur les chromosomes. Cet effet de position résulte de l'organisation des génomes eucaryotes en domaines fonctionnels, définis par des changements locaux au niveau de la structure de la chromatine. Dans chacun de ces domaines, l'expression des gènes est définie et maintenue par l'action concertée de différents éléments régulateurs tels que les promoteurs, les amplificateurs, les silenceurs et les locus control régions. Ces domaines peuvent être entourés par des éléments barrière, séparant les régions de chromatine répressive des régions permissive pour l'expression des gènes. Lorsqu'ils se situent à proximité d'éléments chromosomiques comme les telomères, les gènes peuvent être réprimés de manière épigénétique. Chez la levure, cette répression est établie par la propagation des protéines SIR depuis les télomères vers les régions centromériques. Certains facteurs de transcription peuvent empêcher la répression d'un gène, lorsqu'ils sont placés entre ce gène et le télomère. Nous avons étudié un de ces facteurs, CTF-1, qui a la particularité de lier directement l'histone H3. La délétion de certaines parties de CTF-1 a permis de déterminer que la région responsable de l'activité barrière correspond au domaine d'interaction avec H3. Plusieurs mutations points effectuées dans ce domaine inhibent à la fois l'activité barrière et la capacité de lier H3. Des expériences d'immuno-précipitation de la chromatine indiquent que la protéine barrière CTF-1 prévient efficacement la propagation des protéines SIR et sépare des domaines contenant des histones hypo-acétylées de ceux constitués d'histones hyper-acétylées. Ces résultats suggèrent que CTF-1 interagit directement avec l'histone H3 pour empêcher la propagation de la chromatine répressive, délimitant ainsi des domaines de chromatine permissive et des domaines de chromatine silencieuse.

Structure et évolution thermique de la croûte supérieure des chaînes de subduction : approches thermologique et modélisation numérique, dans les North Cascades (U.S.A) et de la zone de faille Motagua (Guatemala)

Cette thèse cible l'étude de la structure thermique de la croûte supérieure (<10km) dans les arcs magmatiques continentaux, et son influence sur l'enregistrement thermochronologique de leur exhumation et de leur évolution topographique. Nous portons notre regard sur deux chaînes de montagne appartenant aux Cordillères Américaines : Les Cascades Nord (USA) et la zone de faille Motagua (Guatemala). L'approche utilisée est axée sur la thermochronologie (U-Th-Sm)/He sur apatite et zircon, couplée avec la modélisation numérique de la structure thermique de la croûte. Nous mettons en évidence la variabilité à la fois spatiale et temporelle du gradient géothermique, et attirons l'attention du lecteur sur l'importance de prendre en compte la multitude des processus géologiques perturbant la structure thermique dans les chaînes de type cordillère, c'est à dire formées lors de la subduction océanique sous un continent.Une nouvelle approche est ainsi développée pour étudier et contraindre la perturbation thermique autour des chambres magmatiques. Deux profiles âge-elevation (U-Th-Sm)/He sur apatite et zircon, ont été collectées 7 km au sud du batholithe de Chilliwack, Cascades Nord. Les résultats montrent une variabilité spatiale et temporelle du gradient géothermique lors de l'emplacement magmatique qui peut être contrainte et séparé de l'exhumation. Durant l'emplacement de l'intrusion, la perturbation thermique y atteint un état d'équilibre (-80-100 °C/km) qui est fonction du flux de magma et de ia distance à la source du magma, puis rejoint 40 °C/km à la fin du processus d'emplacement magmatique.Quelques nouvelles données (U-Th)/He, replacées dans une compilation des données existantes dans les Cascades Nord, indiquent une vitesse d'exhumation constante (-100 m/Ma) dans le temps et l'espace entre 35 Ma et 2 Ma, associée à un soulèvement uniforme de la chaîne contrôlé par l'emplacement de magma dans la croûte durant toute l'activité de l'arc. Par contre, après ~2 Ma, le versant humide de la chaîne est affecté par une accélération des taux d'exhumation, jusqu'à 3 km de croûte y sont érodés. Les glaciations ont un triple effet sur l'érosion de cette chaîne: (1) augmentation des vitesses d'érosion, d'exhumation et de soulèvement la où les précipitations sont suffisantes, (2) limitation de l'altitude contrôlé par la position de Γ Ε LA, (3) élargissement du versant humide et contraction du versant aride de la chaîne.Les modifications des réseaux de drainage sont des processus de surface souvent sous-estimés au profil d'événements climatiques ou tectoniques. Nous proposons une nouvelle approche couplant une analyse géomorphologique, des données thermochronologiques de basse température ((U-Th-Sm)/He sur apatite et zircon), et l'utilisation de modélisation numérique thermo-cinématique pour les mettre en évidence et les dater; nous testons cette approche sur la gorge de la Skagit river dans les North Cascades.De nouvelles données (U-Th)/He sur zircons, complétant les données existantes, montrent que le déplacement horizontal le long de la faille transformante continentale Motagua, la limite des plaques Caraïbe/Amérique du Nord, a juxtaposé un bloc froid, le bloc Maya (s.s.), contre un bloque chaud, le bloc Chortis (s.s.) originellement en position d'arc. En plus de donner des gammes d'âges thermochronologiques très différents des deux côtés de la faille, le déplacement horizontal rapide (~2 cm/a) a produit un fort échange thermique latéral, résultant en un réchauffement du côté froid et un refroidissement du côté chaud de la zone de faille de Motagua.Enfin des données (U-Th-Sm)/He sur apatite témoignent d'un refroidissement Oligocène enregistré uniquement dans la croûte supérieure de la bordure nord de la zone de faille Motagua. Nous tenterons ultérieurement de reproduire ce découplage vertical de la structure thermique par la modélisation de la formation d'un bassin transtensif et de circulation de fluides le long de la faille de Motagua. - This thesis focuses on the influence of the dynamic thermal structure of the upper crust (<10km) on the thermochronologic record of the exhumational and topographic history of magmatic continental arcs. Two mountain belts from the American Cordillera are studied: the North Cascades (USA) and the Motagua fault zone (Guatemala). I use a combined approach coupling apatite and zircon (U-Th-Sm}/He thermochronology and thermo- kinematic numerical modelling. This study highlights the temporal and spatial variability of the geothermal gradient and the importance to take into account the different geological processes that perturb the thermal structure of Cordilleran-type mountain belts (i.e. mountain belts related to oceanic subduction underneath a continent}.We integrate apatite and zircon (U-Th)/He data with numerical thermo-kinematic models to study the relative effects of magmatic and surface processes on the thermal evolution of the crust and cooling patterns in the Cenozoic North Cascades arc (Washington State, USA). Two age-elevation profiles that are located 7 km south of the well-studied Chiliiwack intrusions shows that spatial and temporal variability in geothermal gradients linked to magma emplacement can be contrained and separated from exhumation processes. During Chiliiwack batholith emplacement at -35-20 Ma, the geothermal gradient of the country rocks increased to a very high steady-state value (80-100°C/km), which is likely a function of magma flux and the distance from the magma source area. Including temporally varying geothermal gradients in the analysis allows quantifying the thermal perturbation around magmatic intrusions and retrieving a relatively simple denudation history from the data.The synthesis of new and previously published (U-Th)/He data reveals that denudation of the Northern Cascades is spatially and temporally constant at -100 m/Ma between ~32 and ~2 Ma, which likely reflects uplift due to magmatic crustal thickening since the initiation of the Cenozoic stage of the continental magmatic arc. In contrast, the humid flank of the North Cascades is affected by a ten-fold acceleration in exhumation rate at ~2 Ma, which we interpret as forced by the initiation of glaciations; around 3 km of crust have been eroded since that time. Glaciations have three distinct effects on the dynamics of this mountain range: (1) they increase erosion, exhumation and uplift rates where precipitation rates are sufficient to drive efficient glacial erosion; (2) they efficiently limit the elevation of the range; (3) they lead to widening of the humid flank and contraction of the arid flank of the belt.Drainage reorganizations constitute an important agent of landscape evolution that is often underestimated to the benefit of tectonic or climatic events. We propose a new method that integrates geomorphology, low-temperature thermochronometry (apatite and zircon {U-Th-Sm)/He), and 3D numerical thermal-kinematic modelling to detect and date drainage instability producing recent gorge incision, and apply this approach to the Skagit River Gorge, North Cascades.Two zircon (U-Th)/He age-elevation profiles sampled on both sides of the Motagua Fault Zone (MFZ), the boundary between the North American and the Caribbean plates, combined with published thermochronological data show that strike-slip displacement has juxtaposed the cold Maya block (s.s.) against the hot, arc derived, Chortis block (s.s ), producing different age patterns on both sides of the fault and short-wavelength lateral thermal exchange, resulting in recent heating of the cool side and cooling of the hot side of the MFZ.Finally, an apatite (U-Th-Sm)/He age-elevation profile records rapid cooling at -35 Ma localized only in the upper crust along the northern side of the Motagua fault zone. We will try to reproduce these data by modeling the thermal perturbation resulting from the formation of a transtensional basin and of fluid flow activity along a crustal- scale strike-slip fault.

Stable isotope and trace element stratigraphy across the Permian-Triassic transition: A redefinition of the boundary in the Velebit Mountain, Croatia

Stable isotopes of carbonates (delta(13)C(carb), delta(18)O(carb)), organic matter (delta(13)C(org), delta(15)N(org)) and major, trace and rare earth element (REE) compositions of marine carbonate rocks of Late Permian to Early Triassic age were used to establish the position of the Permian-Triassic boundary (PTB) at two continuous sections in the Velebit Mountain, Croatia. The chosen sections - Rizvanusa and Brezimenjaca - are composed of two lithostratigraphic units, the Upper Permian Transitional Dolomite and the overlying Sandy Dolomite. The contact between these units, characterized by the erosional features and sudden occurrence of ooids and siliciclastic grains, was previously considered as the chronostratigraphic PTB. The Sandy Dolomite is characterized by high content of non-carbonate material (up to similar to 30 wt.% insoluble residue), originated from erosion of the uplifted hinterland. A relatively rich assemblage of Permian fossils (including Geinitzina, Globivalvulina, Hemigordius, bioclasts of gastropods, ostracods and brachiopods) was found for the first time in Sandy Dolomite, 5 m above the lithologic boundary in the Rizvanusa section. A rather abrupt negative delta(13)C(carb) excursion in both sections appears in rocks showing no recognizable facies change within the Sandy Dolomite, -2 parts per thousand at Rizvanusa and -1.2 parts per thousand at Brezimenjaca, 11 m and 0.2 m above the lithologic contact, respectively. This level within the lower part of the Sandy Dolomite is proposed as the chemostratigraphic PTB. In the Rizvanusa section, the delta(13)C(org) values decline gradually from similar to-25 parts per thousand in the Upper Permian to similar to-29 parts per thousand in the Lower Triassic. The first negative delta(13)C(org) excursion occurs above the lithologic contact, within the uppermost Permian deposits, and appears to be related to the input of terrigenous material. The release of isotopically light microbial soil-biomass into the shallow-marine water may explain this sudden decrease of delta(13)C(org) values below the PTB. This would support the hypothesis that in the western Tethyan realm the land extinction, triggering a sudden drop of woody vegetation and related land erosion, preceded the marine extinction. The relatively low delta(15)N(org) values at the Permian-Triassic (P-Tr) transition level, close to approximate to 0 parts per thousand, and a secondary negative delta(13)C(org) excursion of -0.5 parts per thousand point to significant terrestrial input and primary contribution of cyanobacteria. The profiles of the concentrations of redox-sensitive elements (Ce, Mn, Fe, V), biogenic or biogenic-scavenged elements (P, Ba, Zn, V), Ce/Ce* values, and normalized trace elements, including Ba/Al, Ba/Fe, Ti/Al, Al/(Al + Fe + Mn) and Mn/Ti show clear excursions at the Transitional Dolomite-Sandy Dolomite lithologic boundary and the chemostratigraphic P-Tr boundary. The stratigraphic variations indicate a major regression phase marking the lithologic boundary, transgressive phases in the latest Permian and a gradual change into shallow/stagnant anoxic marine environment towards the P-Tr boundary level and during the earliest Triassic. (C) 2010 Elsevier B.V. All rights reserved.

The stationary distribution of a continuously varying strategy in a class-structured population under mutation-selection-drift balance.

Many traits and/or strategies expressed by organisms are quantitative phenotypes. Because populations are of finite size and genomes are subject to mutations, these continuously varying phenotypes are under the joint pressure of mutation, natural selection and random genetic drift. This article derives the stationary distribution for such a phenotype under a mutation-selection-drift balance in a class-structured population allowing for demographically varying class sizes and/or changing environmental conditions. The salient feature of the stationary distribution is that it can be entirely characterized in terms of the average size of the gene pool and Hamilton's inclusive fitness effect. The exploration of the phenotypic space varies exponentially with the cumulative inclusive fitness effect over state space, which determines an adaptive landscape. The peaks of the landscapes are those phenotypes that are candidate evolutionary stable strategies and can be determined by standard phenotypic selection gradient methods (e.g. evolutionary game theory, kin selection theory, adaptive dynamics). The curvature of the stationary distribution provides a measure of the stability by convergence of candidate evolutionary stable strategies, and it is evaluated explicitly for two biological scenarios: first, a coordination game, which illustrates that, for a multipeaked adaptive landscape, stochastically stable strategies can be singled out by letting the size of the gene pool grow large; second, a sex-allocation game for diploids and haplo-diploids, which suggests that the equilibrium sex ratio follows a Beta distribution with parameters depending on the features of the genetic system.