Speed of reaction-diffusion fronts in spatially heterogeneous media

The front speed problem for nonuniform reaction rate and diffusion coefficient is studied by using singular perturbation analysis, the geometric approach of Hamilton-Jacobi dynamics, and the local speed approach. Exact and perturbed expressions for the front speed are obtained in the limit of large times. For linear and fractal heterogeneities, the analytic results have been compared with numerical results exhibiting a good agreement. Finally we reach a general expression for the speed of the front in the case of smooth and weak heterogeneities

Seismic attenuation and velocity dispersion in heterogeneous partially saturated porous rocks

Using a numerical approach, we explore wave-induced fluid flow effects in partially saturated porous rocks in which the gas-water saturation patterns are governed by mesoscopic heterogeneities associated with the dry frame properties. The link between the dry frame properties and the gas saturation is defined by the assumption of capillary pressure equilibrium, which in the presence of heterogeneity implies that neighbouring regions can exhibit different levels of saturation. To determine the equivalent attenuation and phase velocity of the synthetic rock samples considered in this study, we apply a numerical upscaling procedure, which permits to take into account mesoscopic heterogeneities associated with the dry frame properties as well as spatially continuous variations of the pore fluid properties. The multiscale nature of the fluid saturation is taken into account by locally computing the physical properties of an effective fluid, which are then used for the larger-scale simulations. We consider two sets of numerical experiments to analyse such effects in heterogeneous partially saturated porous media, where the saturation field is determined by variations in porosity and clay content, respectively. In both cases we also evaluate the seismic responses of corresponding binary, patchy-type saturation patterns. Our results indicate that significant attenuation and modest velocity dispersion effects take place in this kind of media for both binary patchy-type and spatially continuous gas saturation patterns and in particular in the presence of relatively small amounts of gas. The numerical experiments also show that the nature of the gas distribution patterns is a critical parameter controlling the seismic responses of these environments, since attenuation and velocity dispersion effects are much more significant and occur over a broader saturation range for binary patchy-type gas-water distributions. This analysis therefore suggests that the physical mechanisms governing partial saturation should be accounted for when analysing seismic data in a poroelastic framework. In this context, heterogeneities associated with the dry frame properties, which do not play important roles in wave-induced fluid flow processes per se, should be taken into account since they may determine the kind of gas distribution pattern taking place in the porous rock.

Linking life-history traits, ecology, and niche breadth evolution in north american eriogonoids (polygonaceae).

Abstract Macroevolutionary and microevolutionary studies provide complementary explanations of the processes shaping the evolution of niche breadth. Macroevolutionary approaches scrutinize factors such as the temporal and spatial environmental heterogeneities that drive differentiation among species. Microevolutionary studies, in contrast, focus on the processes that affect intraspecific variability. We combine these perspectives by using macroevolutionary models in a comparative study of intraspecific variability. We address potential differences in rates of evolution of niche breadth and position in annual and perennial plants of the Eriogonoideae subfamily of the Polygonaceae. We anticipated higher rates of evolution in annuals than in perennials owing to differences in generation time that are paralleled by rates of molecular evolution. Instead, we found that perennial eriogonoid species present greater environmental tolerance (wider climate niche) than annual species. Niche breadth of perennial species has evolved two to four times faster than in annuals, while niche optimum has diversified more rapidly among annual species than among perennials. Niche breadth and average elevation of species are correlated. Moreover, niche breadth increases more rapidly with mean species elevation in perennials than in annuals. Our results suggest that both environmental gradients and life-history strategy influence rates and patterns of niche breadth evolution.

Méthodes sismiques actives adaptées à l'étude de milieux hétérogènes : application aux glissements de terrain

Résumé Les glissements de terrain représentent un des principaux risques naturels dans les régions montagneuses. En Suisse, chaque année les glissements de terrains causent des dégâts qui affectent les infrastructures et ont des coûts financiers importants. Une bonne compréhension des mécanismes des glissements peut permettre d'atténuer leur impact. Celle-ci passe notamment par la connaissance de la structure interne du glissement, la détermination de son volume et de son ou ses plans de glissement. Dans un glissement de terrain, la désorganisation et la présence de fractures dans le matériel déplacé engendre un changement des paramètres physiques et en particulier une diminution des vitesses de propagation des ondes sismiques ainsi que de la densité du matériel. Les méthodes sismiques sont de ce fait bien adaptées à l'étude des glissements de terrain. Parmi les méthodes sismiques, l'analyse de la dispersion des ondes de surface est une méthode simple à mettre en oeuvre. Elle présente l'avantage d'estimer les variations des vitesses de cisaillement avec la profondeur sans avoir spécifiquement recours à l'utilisation d'une source d'onde S et de géophones horizontaux. Sa mise en oeuvre en trois étapes implique la mesure de la dispersion des ondes de surface sur des réseaux étendus, la détermination des courbes de dispersion pour finir par l'inversion de ces courbes. Les modèles de vitesse obtenus à partir de cette procédure ne sont valides que lorsque les milieux explorés ne présentent pas de variations latérales. En pratique cette hypothèse est rarement vérifiée, notamment pour un glissement de terrain dans lequel les couches remaniées sont susceptibles de présenter de fortes hétérogénéités latérales. Pour évaluer la possibilité de déterminer des courbes de dispersion à partir de réseaux de faible extension des mesures testes ont été effectuées sur un site (Arnex, VD) équipé d'un forage. Un profil sismique de 190 m de long a été implanté dans une vallée creusée dans du calcaire et remplie par des dépôts glacio-lacustres d'une trentaine de mètres d'épaisseur. Les données acquises le long de ce profil ont confirmé que la présence de variations latérales sous le réseau de géophones affecte l'allure des courbes de dispersion jusqu'à parfois empêcher leur détermination. Pour utiliser l'analyse de la dispersion des ondes de surface sur des sites présentant des variations latérales, notre approche consiste à déterminer les courbes de dispersions pour une série de réseaux de faible extension, à inverser chacune des courbes et à interpoler les différents modèles de vitesse obtenus. Le choix de la position ainsi que de l'extension des différents réseaux de géophones est important. Il tient compte de la localisation des hétérogénéités détectées à partir de l'analyse de sismique réfraction, mais également d'anomalies d'amplitudes observées sur des cartes qui représentent dans le domaine position de tir - position du récepteur, l'amplitude mesurée pour différentes fréquences. La procédure proposée par Lin et Lin (2007) s'est avérée être une méthode efficace permettant de déterminer des courbes de dispersion à partir de réseaux de faible extension. Elle consiste à construire à partir d'un réseau de géophones et de plusieurs positions de tir un enregistrement temps-déports qui tient compte d'une large gamme de distances source-récepteur. Au moment d'assembler les différentes données une correction de phase est appliquée pour tenir compte des hétérogénéités situées entre les différents points de tir. Pour évaluer cette correction nous suggérons de calculer pour deux tir successif la densité spectrale croisée des traces de même offset: Sur le site d'Arnex, 22 courbes de dispersions ont été déterminées pour de réseaux de géophones de 10 m d'extension. Nous avons également profité du forage pour acquérir un profil de sismique verticale en ondes S. Le modèle de vitesse S déduit de l'interprétation du profil de sismique verticale est utilisé comme information à priori lors l'inversion des différentes courbes de dispersion. Finalement, le modèle en deux dimension qui a été établi grâce à l'analyse de la dispersion des ondes de surface met en évidence une structure tabulaire à trois couches dont les limites coïncident bien avec les limites lithologiques observées dans le forage. Dans celui-ci des argiles limoneuses associées à une vitesse de propagation des ondes S de l'ordre de 175 m/s surmontent vers 9 m de profondeur des dépôts de moraine argilo-sableuse caractérisés par des vitesses de propagation des ondes S de l'ordre de 300 m/s jusqu'à 14 m de profondeur et supérieur ou égal à 400 m/s entre 14 et 20 m de profondeur. Le glissement de la Grande Combe (Ballaigues, VD) se produit à l'intérieur du remplissage quaternaire d'une combe creusée dans des calcaires Portlandien. Comme dans le cas du site d'Arnex les dépôts quaternaires correspondent à des dépôts glacio-lacustres. Dans la partie supérieure la surface de glissement a été localisée à une vingtaine de mètres de profondeur au niveau de l'interface qui sépare des dépôts de moraine jurassienne et des dépôts glacio-lacustres. Au pied du glissement 14 courbes de dispersions ont été déterminées sur des réseaux de 10 m d'extension le long d'un profil de 144 m. Les courbes obtenues sont discontinues et définies pour un domaine de fréquence de 7 à 35 Hz. Grâce à l'utilisation de distances source-récepteur entre 8 et 72 m, 2 à 4 modes de propagation ont été identifiés pour chacune des courbes. Lors de l'inversion des courbes de dispersion la prise en compte des différents modes de propagation a permis d'étendre la profondeur d'investigation jusqu'à une vingtaine de mètres de profondeur. Le modèle en deux dimensions permet de distinguer 4 couches (Vs1 < 175 m/s, 175 m/s < Vs2 < 225 m/s, 225 m/s < Vs3 < 400 m/s et Vs4 >.400 m/s) qui présentent des variations d'épaisseur. Des profils de sismiques réflexion en ondes S acquis avec une source construite dans le cadre de ce travail, complètent et corroborent le modèle établi à partir de l'analyse de la dispersion des ondes de surface. Un réflecteur localisé entre 5 et 10 m de profondeur et associé à une vitesse de sommation de 180 m/s souligne notamment la géométrie de l'interface qui sépare la deuxième de la troisième couche du modèle établi à partir de l'analyse de la dispersion des ondes de surface. Abstract Landslides are one of the main natural hazards in mountainous regions. In Switzerland, landslides cause damages every year that impact infrastructures and have important financial costs. In depth understanding of sliding mechanisms may help limiting their impact. In particular, this can be achieved through a better knowledge of the internal structure of the landslide, the determination of its volume and its sliding surface or surfaces In a landslide, the disorganization and the presence of fractures in the displaced material generate a change of the physical parameters and in particular a decrease of the seismic velocities and of the material density. Therefoe, seismic methods are well adapted to the study of landslides. Among seismic methods, surface-wave dispersion analysis is a easy to implement. Through it, shearwave velocity variations with depth can be estimated without having to resort to an S-wave source and to horizontal geophones. Its 3-step implementation implies measurement of surface-wave dispersion with long arrays, determination of the dispersion curves and finally inversion of these curves. Velocity models obtained through this approach are only valid when the investigated medium does not include lateral variations. In practice, this assumption is seldom correct, in particular for landslides in which reshaped layers likely include strong lateral heterogeneities. To assess the possibility of determining dispersion curves from short array lengths we carried out tests measurements on a site (Arnex, VD) that includes a borehole. A 190 m long seismic profile was acquired in a valley carved into limestone and filled with 30 m of glacio-lacustrine sediments. The data acquired along this profile confirmed that the presence of lateral variations under the geophone array influences the dispersion-curve shape so much that it sometimes preventes the dispersion curves determination. Our approach to use the analysis of surface-wave dispersion on sites that include lateral variations consists in obtaining dispersion curves for a series of short length arrays; inverting each so obtained curve and interpolating the different obtained velocity model. The choice of the location as well as the geophone array length is important. It takes into account the location of the heterogeneities that are revealed by the seismic refraction interpretation of the data but also, the location of signal amplitude anomalies observed on maps that represent, for a given frequency, the measured amplitude in the shot position - receiver position domain. The procedure proposed by Lin and Lin (2007) turned out to be an efficient one to determine dispersion curves using short extension arrays. It consists in building a time-offset from an array of geophones with a wide offset range by gathering seismograms acquired with different source-to-receiver offsets. When assembling the different data, a phase correction is applied in order to reduce static phase error induced by lateral variation. To evaluate this correction, we suggest to calculate, for two successive shots, the cross power spectral density of common offset traces. On the Arnex site, 22 curves were determined with 10m in length geophone-arrays. We also took advantage of the borehole to acquire a S-wave vertical seismic profile. The S-wave velocity depth model derived from the vertical seismic profile interpretation is used as prior information in the inversion of the dispersion-curves. Finally a 2D velocity model was established from the analysis of the different dispersion curves. It reveals a 3-layer structure in good agreement with the observed lithologies in the borehole. In it a clay layer with a shear-wave of 175 m/s shear-wave velocity overlies a clayey-sandy till layer at 9 m depth that is characterized down to 14 m by a 300 m/s S-wave velocity; these deposits have a S-wave velocity of 400 m/s between depths of 14 to 20 m. The La Grand Combe landslide (Ballaigues, VD) occurs inside the Quaternary filling of a valley carved into Portlandien limestone. As at the Arnex site, the Quaternary deposits correspond to glaciolacustrine sediments. In the upper part of the landslide, the sliding surface is located at a depth of about 20 m that coincides with the discontinuity between Jurassian till and glacio-lacustrine deposits. At the toe of the landslide, we defined 14 dispersion curves along a 144 m long profile using 10 m long geophone arrays. The obtained curves are discontinuous and defined within a frequency range of 7 to 35 Hz. The use of a wide range of offsets (from 8 to 72 m) enabled us to determine 2 to 4 mode of propagation for each dispersion curve. Taking these higher modes into consideration for dispersion curve inversion allowed us to reach an investigation depth of about 20 m. A four layer 2D model was derived (Vs1< 175 m/s, 175 m/s <Vs2< 225 m/s, 225 m/s < Vs3 < 400 m/s, Vs4> 400 m/s) with variable layer thicknesses. S-wave seismic reflection profiles acquired with a source built as part of this work complete and the velocity model revealed by surface-wave analysis. In particular, reflector at a depth of 5 to 10 m associated with a 180 m/s stacking velocity image the geometry of the discontinuity between the second and third layer of the model derived from the surface-wave dispersion analysis.

Seismic attenuation in heterogeneous partially saturated rocks

We study wave-induced fluid flow effects in porous rocks partially saturated with gas and water, where the saturation patterns are governed by mesoscopic heterogeneities associated with the dry frame properties. The link between the dry frame properties and the gas saturation is defined by the assumption of capillary pressure equilibrium, which in the presence of heterogeneity implies that neighboring regions can exhibit different levels of saturation. In order to determine the equivalent attenuation and phase velocity of the synthetic rock samples considered in this study, we apply a numerical upscaling procedure, which permits to take into account mesoscopic heterogeneities associated with the dry frame properties as well as spatially continuous variations of the pore fluid properties. We consider numerical experiments to analyze such effects in heterogeneous partially saturated porous media, where the saturation field is determined by realistic variations in porosity. Our results indicate that the spatially continuous nature of gas saturation inherent to this study is a critical parameter controlling the seismic response of these environments, which in turn suggests that the physical mechanisms governing partial saturation should be accounted for when analyzing seismic data in a poro-elastic context.

Has the Euro-Mediterranean partnership affected Mediterranean business cycles?

We date turning points of the reference cycle for 19 Mediterranean countries andanalyze their structure and interdependences. Fluctuations are volatile and not highlycorrelated across countries; recessions are deep but asynchronous making average outputlosses in the area limited. Heterogeneities across countries and regions are substantial.Mediterranean cycles are time varying but their evolution is not linked withthe Euro-Mediterranean partnership process. The concordance of cyclical fluctuationsis poorly related to trade and financial linkages and to their evolution over time.

Testing for convergence clubs in income per-capita: A predictive density approach

The paper proposes a technique to jointly test for groupings of unknown size in the cross sectional dimension of a panel and estimates the parameters of each group, and applies it to identifying convergence clubs in income per-capita. The approach uses the predictive density of the data, conditional on the parameters of the model. The steady state distribution of European regional data clusters around four poles of attraction with different economic features. The distribution of incomeper-capita of OECD countries has two poles of attraction and each grouphas clearly identifiable economic characteristics.

Quantitative comparison of simulations of seismic wave propagation in heterogeneous poro-elastic media involving fluid-solid interfaces and in equivalent visco-elastic solids

There is increasing evidence to suggest that the presence of mesoscopic heterogeneities constitutes the predominant attenuation mechanism at seismic frequencies. As a consequence, centimeter-scale perturbations of the subsurface physical properties should be taken into account for seismic modeling whenever detailed and accurate responses of the target structures are desired. This is, however, computationally prohibitive since extremely small grid spacings would be necessary. A convenient way to circumvent this problem is to use an upscaling procedure to replace the heterogeneous porous media by equivalent visco-elastic solids. In this work, we solve Biot's equations of motion to perform numerical simulations of seismic wave propagation through porous media containing mesoscopic heterogeneities. We then use an upscaling procedure to replace the heterogeneous poro-elastic regions by homogeneous equivalent visco-elastic solids and repeat the simulations using visco-elastic equations of motion. We find that, despite the equivalent attenuation behavior of the heterogeneous poro-elastic medium and the equivalent visco-elastic solid, the seismograms may differ due to diverging boundary conditions at fluid-solid interfaces, where there exist additional options for the poro-elastic case. In particular, we observe that the seismograms agree for closed-pore boundary conditions, but differ significantly for open-pore boundary conditions. This is an interesting result, which has potentially important implications for wave-equation-based algorithms in exploration geophysics involving fluid-solid interfaces, such as, for example, wave field decomposition.

Labor market heterogeneity and the aggregate matching function

The matching function -a key building block in models of labor market frictions- impliesthat the job finding rate depends only on labor market tightness. We estimate such amatching function and find that the relation, although remarkably stable over 1967-2007,broke down spectacularly after 2007. We argue that labor market heterogeneities are notfully captured by the standard matching function, but that a generalized matching functionthat explicitly takes into account worker heterogeneity and market segmentation is fullyconsistent with the behavior of the job finding rate. The standard matching function canbreak down when, as in the Great Recession, the average characteristics of the unemployedchange too much, or when dispersion in labor market conditions -the extent to which somelabor markets fare worse than others- increases too much.

Numerical simulation of ambient seismic wavefield modification caused by pore-fluid effects in an oil reservoir

We have modeled numerically the seismic response of a poroelastic inclusion with properties applicable to an oil reservoir that interacts with an ambient wavefield. The model includes wave-induced fluid flow caused by pressure differences between mesoscopic-scale (i.e., in the order of centimeters to meters) heterogeneities. We used a viscoelastic approximation on the macroscopic scale to implement the attenuation and dispersion resulting from this mesoscopic-scale theory in numerical simulations of wave propagation on the kilometer scale. This upscaling method includes finite-element modeling of wave-induced fluid flow to determine effective seismic properties of the poroelastic media, such as attenuation of P- and S-waves. The fitted, equivalent, viscoelastic behavior is implemented in finite-difference wave propagation simulations. With this two-stage process, we model numerically the quasi-poroelastic wave-propagation on the kilometer scale and study the impact of fluid properties and fluid saturation on the modeled seismic amplitudes. In particular, we addressed the question of whether poroelastic effects within an oil reservoir may be a plausible explanation for low-frequency ambient wavefield modifications observed at oil fields in recent years. Our results indicate that ambient wavefield modification is expected to occur for oil reservoirs exhibiting high attenuation. Whether or not such modifications can be detected in surface recordings, however, will depend on acquisition design and noise mitigation processing as well as site-specific conditions, such as the geologic complexity of the subsurface, the nature of the ambient wavefield, and the amount of surface noise.

Diffusion-weighted spectroscopy: a novel approach to determine macromolecule resonances in short-echo time 1H-MRS.

Quantification of short-echo time proton magnetic resonance spectroscopy results in >18 metabolite concentrations (neurochemical profile). Their quantification accuracy depends on the assessment of the contribution of macromolecule (MM) resonances, previously experimentally achieved by exploiting the several fold difference in T(1). To minimize effects of heterogeneities in metabolites T(1), the aim of the study was to assess MM signal contributions by combining inversion recovery (IR) and diffusion-weighted proton spectroscopy at high-magnetic field (14.1 T) and short echo time (= 8 msec) in the rat brain. IR combined with diffusion weighting experiments (with δ/Δ = 1.5/200 msec and b-value = 11.8 msec/μm(2)) showed that the metabolite nulled spectrum (inversion time = 740 msec) was affected by residuals attributed to creatine, inositol, taurine, choline, N-acetylaspartate as well as glutamine and glutamate. While the metabolite residuals were significantly attenuated by 50%, the MM signals were almost not affected (< 8%). The combination of metabolite-nulled IR spectra with diffusion weighting allows a specific characterization of MM resonances with minimal metabolite signal contributions and is expected to lead to a more precise quantification of the neurochemical profile.

Do seismic waves and fluid flow sense the same permeability?

We study the discrepancy between the effective flow permeability and the effective seismic permeability, that is, the effective permeability controlling seismic attenuation due to wave-induced fluid flow, in 2D rock samples having mesoscopic heterogeneities and in the presence of strong permeability fluctuations. In order to do so, we employ a numerical oscillatory compressibility test to determine attenuation and velocity dispersion due to wave-induced fluid flow in these kinds of media and compare the responses with those obtained by replacing the heterogeneous permeability field by constant values, including the average permeability as well as the effective flow permeability of the sample. The latter is estimated in a separate upscaling procedure by solving the steady-state flow equation in the rock sample under study. Numerical experiments let us verify that attenuation levels are less significant and the attenuation peak gets broader in the presence of such strong permeability fluctuations. Moreover, we observe that for very low frequencies the effective seismic permeability is similar to the effective flow permeability, while for very high frequencies it approaches the arithmetic average of the permeability field.

Seismic attenuation due to wave-induced fluid flow at microscopic and macroscopic scales

Wave-induced fluid flow at microscopic and mesoscopic scales arguably constitutes the major cause of intrinsic seismic attenuation throughout the exploration seismic and sonic frequency ranges. The quantitative analysis of these phenomena is, however, complicated by the fact that the governing physical processes may be dependent. The reason for this is that the presence of microscopic heterogeneities, such as micro-cracks or broken grain contacts, causes the stiffness of the so-called modified dry frame to be complex-valued and frequency-dependent, which in turn may affect the viscoelastic behaviour in response to fluid flow at mesoscopic scales. In this work, we propose a simple but effective procedure to estimate the seismic attenuation and velocity dispersion behaviour associated with wave-induced fluid flow due to both microscopic and mesoscopic heterogeneities and discuss the results obtained for a range of pertinent scenarios.

Actively maintained lipid nanodomains in biomembranes

Lipid rafts, defined as domains rich in cholesterol and sphingolipids, are involved in many important plasma membrane functions. Recent studies suggest that the way cells handle membrane cholesterol is fundamental in the formation of such lateral heterogeneities. We propose to model the plasma membrane as a nonequilibrium phase-separating system where cholesterol is dynamically incorporated and released. The model shows how cellular regulation of membrane cholesterol may determine the nanoscale lipid organization when the lipid mixture is close to a phase separation boundary, providing a plausible mechanism for raft formation in vivo.

Stratigraphie, morphodynamique, paléoenvironnements des terrains sédimentaires meubles à forte déclivité du domaine périglaciaire alpin

Dans le contexte d'un climat de plus en plus chaud, une étude « géosystémique » de la répartition du pergélisol dans l'ensemble d'un versant périglaciaire alpin, de la paroi rocheuse jusqu'au glacier rocheux, s'avère primordiale. S'insérant dans cette problématique, ce travail de thèse vise comme objectif général l'étude des versants d'éboulis situés à l'intérieur de la ceinture du pergélisol discontinu selon deux volets de recherche différents : une étude de la stratigraphie et de la répartition du pergélisol dans les éboulis de haute altitude et des processus qui lui sont associés ; une reconstitution de l'histoire paléoenvironnementale du domaine périglaciaire alpin pendant le Tardiglaciaire et l'Holocène. La stratigraphie et la répartition spatiale du pergélisol a été étudiée dans cinq éboulis des Alpes Valaisannes (Suisse), dont trois ont fait l'objet de forages profonds, grâce à la prospection géophysique de détail effectuée à l'aide de méthodes thermiques, de résistivité, sismiques et nucléaires. Les mesures effectuées ont permis de mettre en évidence que, dans les cinq éboulis étudiés, la répartition du pergélisol est discontinue et aucun des versants n'est intégralement occupé par du pergélisol. En particulier, il a été possible de prouver de manière directe que, dans un éboulis, le pergélisol est présent dans les parties inférieures du versant et absent dans les parties supérieures. Trois facteurs de contrôle principaux de la répartition du pergélisol déterminée au sein des éboulis étudiés ont été individualisés, pouvant agir seuls ou de manière combinée : la ventilation ascendante, l'augmentation de la granulométrie en direction de l'aval et la redistribution de la neige par le vent et les avalanches. Parmi ceux-ci, la relation ventilation-granulométrie semble être le facteur de contrôle principal permettant d'expliquer la présence de pergélisol dans les parties inférieures d'un éboulis et son absence dans les parties supérieures. Enfin, l'analyse de la structure des éboulis périglaciaires de haute altitude a permis de montrer que la stratigraphie du pergélisol peut être un élément important pour l'interprétation de la signification paléoclimatique de ce type de formes. Pour le deuxième volet de la recherche, grâce aux datations relatives effectuées à l'aide de l'utilisation conjointe de la méthode paléogéographique et du marteau de Schmidt, il a été possible de définir la chrono-stratigraphie du retrait glaciaire et du développement des glaciers rocheux et des versants d'éboulis des quatre régions des Alpes suisses étudiées (régions du Mont Gelé - Mont Fort, des Fontanesses et de Chamosentse, dans les Alpes Valaisannes, et Massif de la Cima di Gana Bianca, dans les Alpes Tessinoises). La compilation de toutes les datations effectuées a permis de montrer que la plupart des glaciers rocheux actifs étudiés se seraient développés soit juste avant et/ou pendant l'Optimum Climatique Holocène de 9.5-6.3 ka cal BP, soit au plus tard juste après cet évènement climatique majeur du dernier interglaciaire. Parmi les glaciers rocheux fossiles datés, la plupart aurait commencé à se former dans la deuxième moitié du Tardiglaciaire et se serait inactivé dans la première partie de l'Optimum Climatique Holocène. Pour les éboulis étudiés, les datations effectuées ont permis d'observer que leur surface date de la période entre le Boréal et l'Atlantique récent, indiquant que les taux d'éboulisation après la fin de l'Optimum Climatique Holocène ont dû être faibles, et que l'intervalle entre l'âge maximal et l'âge minimal est dans la plupart des cas relativement court (4-6 millénaires), indiquant que les taux d'éboulisation durant la période de formation des éboulis ont dû être importants. Grâce au calcul des taux d'érosion des parois rocheuses sur la base du volume de matériaux rocheux pour quatre des éboulis étudiés, il a été possible mettre en évidence l'existence d'une « éboulisation parapériglaciaire » liée à la dégradation du pergélisol dans les parois rocheuses, fonctionnant principalement durant les périodes de réchauffement climatique rapide comme cela a été le cas au début du Bølling, du Préboréal à la fin de l'Atlantique récent et, peut-être, à partir des années 1980. - In the context of a warmer climate, a « geosystemical » study of the permafrost distribution in a whole alpine periglacial hillslope, from the rockwall to the rockglacier, is of great importance. With respect to this problem, the general objective of this PhD thesis is the global study of talus slopes located within the alpine periglacial belt following two different research axes: the analysis of the internal structure and of the permafrost distribution of high altitude talus slopes and of the related processes; the reconstruction of the palaeoenvironmental history of the alpine periglacial belt during the Lateglacial and the Holocene. The stratigraphy and the permafrost distribution were studied in five talus slopes of the Valais Alps (Switzerland) with the analysis of borehole data (on three of the five talus slopes) and other methods of permafrost prospecting: Electrical Resistivity Tomography (ERT), Refraction Seismic Tomography (RST) and nuclear well logging. The collected data shows that, in all of the studied talus slopes, permafrost distribution is discontinuous and that neither of the hillslopes is integrally characterised by permafrost. In particular, this data proves by direct investigations that, in talus slopes, permafrost is present in the lower parts of the hillslope, whereas it is absent in the upper parts. Permafrost distribution in alpine talus slopes is depending of the combination of almost three controlling factors, whose respective importance is variable: the chimney effect, the increase of grain size downslope and the redistribution of snow by avalanches. Depending on the size of the talus and on topographical and geomorphological heterogeneities, various cases are possible: one dominant controlling factor or the combination of various factors. Nevertheless, it would be an error to consider each controlling factor independently, without considering their relationships. Between these controlling factors, the relationship chimney effect/grain size seems to be the most important factor controlling the presence of permafrost in the lowest part of periglacial talus slopes, and its absence in the upper parts. Finally, the analysis of the talus structure shows that the permafrost stratigraphy may be an important element of interpretation of the palaeoclimatic significance of an alpine talus slope. The second research axe focused on the establishment of a chronology of the Lateglacial glacier retreat and the dating of rockglaciers and talus slopes development in four studied regions of the Swiss Alps (Mont Gelé - Mont Fort, Fontanesses and Chamosentse regions, in the Valais Alps, and the Cima di Gana Bianca Massif, in the Ticino Alps). The compilation of the dates acquired through the combination of the palaeogeographical method and of the Schmidt hammer indicates that most of the investigated active rockglaciers started to evolve during the early phases of the Holocene or, at the latest, after the early-to-mid Holocene Climatic Optimum (ending around 6.3 ka cal BP). For the dated relict rockglaciers, most of them started to evolve in the second half of the Lateglacial, and probably became inactive at the beginning of the Holocene Climatic Optimum. For the investigated talus slopes, the relative dating carried out allowed to show that their surface date from the period included between the Boreal and the end of the Atlantic, pointing out that the rockwall retreat after the end of the Holocene Climatic Optimum was weak, and that the interval between maximal and minimal ages is in most cases relatively short (4-6 millennia). Therefore, the rockwall retreat during the development period of the talus slopes must has been considerable. Thanks to the calculation of rockwall erosion rates based on the volume of talus accumulations for four of the investigated hillslopes, it was possible to find evidences of the existence of "paraperiglacial rockfall phases" related to the permafrost degradation in rockwalls. These phases coincide with rapid climate warming periods, as at the beginning of the Bølling, during the Preboreal or, maybe, since 1980.