No description available.

We consider the system of equations arising from mantle dynamics introduced by McKenzie (J. Petrology, 1985). In this multi-phase model, the fluid melt velocity obeys Darcy's law while the deformable "solid" matrix is governed by a highly viscous Stokes equation. The system is then coupled through mass conservation and compaction relations. Together these equations form a coupled Darcy-Stokes system on a continuous single-domain mixture of fluid and matrix. The porosity φ, representing the relative volume of fluid melt to the bulk volume, is assumed to be much smaller than one. When coupled with solute transport and thermal evolution in a time-dependent problem, the model transitions dynamically from a non-porous single phase solid to a two-phase porous medium. Such mixture models have an advantage for numerical approximation since the free boundary between the one and two-phase regions need not be determined explicitly. The equations of mantle dynamics apply to a wide range of applications in deep earth physics such as mid-ocean ridges, subduction zones, and hot-spot volcanism, as well as to glacier dynamics and other two-phase flows in porous media. Mid-ocean ridges form when viscous corner flow of the solid mantle focuses fluid toward a central ridge. Melt is believed to migrate upward until it reaches the lithospheric "tent" where it then moves toward the ridge in a high porosity band. Simulation of this physical phenomenon required confidence in numerical methods to handle highly heterogeneous porosity as well as the single-phase to two-phase transition. In this work we present a standard mixed finite element method for the equations of mantle dynamics and investigate its limitations for vanishing porosity. While stable and optimally convergent for porosity bounded away from zero, the stability estimates we obtain suggest, and numerical results show, the method becomes unstable as porosity approaches zero. Moreover, the fluid pressure is no longer a physical variable when the fluid phase disappears and thus is not a good variable for numerical methods. Inspired by the stability estimates of the standard method, we develop a novel stable mixed method with uniqueness and existence of solutions by studying a linear degenerate elliptic sub-problem akin to the Darcy part of the full model: [mathematical equation], where a and b satisfy a(0)=b(0)=0 and are otherwise positive, and the porosity φ ≥ 0 may be zero on a set of positive measure. Using scaled variables and mild assumptions on the regularity of φ, we develop a practical mass-conservative method based on lowest order Raviart-Thomas finite elements. Finally, we adapt the numerical method for the sub-problem to the full system of equations. We show optimal convergence for sufficiently smooth solutions for a compacting column and mid-ocean ridge-like corner flow examples, and investigate accuracy and stability for less regular problems / text

A tear in the subducting Nazca slab is located between the end of the Pampean flat slab and normally subducting oceanic lithosphere. Tomographic studies suggest mantle material flows through this opening. The best way to probe this hypothesis is through observations of seismic anisotropy, such as shear wave splitting. We examine patterns of shear wave splitting using data from two seismic deployments in Argentina that lay updip of the slab tear. We observe a simple pattern of plate-motion-parallel fast splitting directions, indicative of plate-motion-parallel mantle flow, beneath the majority of the stations. Our observed splitting contrasts previous observations to the north and south of the flat slab region. Since plate-motion-parallel splitting occurs only coincidentally with the slab tear, we propose mantle material flows through the opening resulting in Nazca plate-motion-parallel flow in both the subslab mantle and mantle wedge.

The Western U.S. has experienced widespread extension during the past 10’s of millions of years, largely within the Basin and Range and Rio Grande Rift provinces. Tomography results from previous studies revealed narrow fast seismic velocity anomalies in the mantle on either side of the Rio Grande Rift as well as at the western edge of the Colorado Plateau. The fast mantle anomalies have been interpreted as down-welling that is part of small scale mantle convection at the edge of extending provinces. It was also found that crust was thicker than average ab¬¬ove the possible mantle down-welling, indicating that mantle dynamics may influence crustal flow. We present results from P/S conversion receiver functions using SIEDCAR (Seismic Investigation of Edge Driven Convection Associated with the Rio Grande Rift) data to determine crustal and lithospheric structure beneath the east flank of the Rio Grande Rift. Crustal and lithosphere thickness are estimated using P-to-S and S-to-P receiver functions respectively. Receiver function migration methods were applied to produce images of the crust and lithosphere. The results show variable crustal thickness through the region with an average thickness of 45 km. The crust achieves its maximum thickness of 60km at 105W longitude, between 33.5N and 32.2N latitude. This observation confirms previous receiver function results from Wilson et al, 2005. Body wave tomography (Rocket, 2011; Schmandt and Humphreys, 2010) using similar data to what we used for the receiver function analysis, shows mantle downwelling closely associated with the thickened crust. We believe that the thickened crust might be due to lower crustal flow associated with mantle downwelling or mantle delamination at the edge of the Rio Grande Rift. In this model the sinking mantle pulls the crust downward causing a pressure gradient within the crust thus causing the flow. Our S-P images show signal from the lithosphere-asthenosphere boundary (LAB) with an average LAB thickness of 100 km but with a sharp transition at about 1050 W from 75 km to over 100 km. The region with abnormally thick crust overlies a region where the lithosphere appears to have a break. We interpret our results as showing that lower lithosphere has and is delaminating near the edge of the Great Plains accompanied by lower crustal flow in some places determined by lower crustal viscosity. / text

Le travail réalisé pendant cette thèse a permis d'apporter de nouvelles contraintes sur le développement et la distribution de la déformation dans le manteau supérieur et plus particulièrement au niveau des grandes limites de plaques décrochantes. Grâce à l'apport de l'expérience USArray et d'une dizaine d'années d'enregistrements sismologiques supplémentaires, nous avons pu étudier, de manière précise, les variations d'anisotropie dans le voisinage de la Faille de San Andreas. Nous avons confirmé et étendu l'observation de deux couches anisotropes sous cette limite de plaque. On y observe une première couche localisée dans la lithosphère marquant la déformation induite à la limite de plaque, et une autre, asthénosphérique, cohérente avec l'anisotropie observée loin de la faille et d'origine plus discutée. Nous avons montré que la zone de déformation associée aux failles de San Andreas, Calaveras et d'Hayward a, vraisemblablement, une largeur d'au moins 40 kilomètres en base de lithosphère, sous chacune de ces failles. Nous avons ensuite procédé à la modélisation thermomécanique (ADELI) de la migration d'une limite de plaques décrochante couplée à une modélisation du développement de fabriques cristallographiques par une approche viscoplastique auto-cohérente (VPSC). Ceci nous a permis d'y observer le développement de la déformation et les conséquences des possibles interactions entre la déformation décrochante en surface et le cisaillement en base de lithosphère dû au déplacement horizontal des plaques. Les propriétés élastiques déduites des fabriques cristallographiques modélisées montrent que de telles interactions existent et provoquent, sous la limite de plaques, une rotation des orientations cristallographiques avec la profondeur. Le signal associé à ces rotations progressives n'est toutefois pas cohérent avec la présence de deux couches d'anisotropie comme proposée sous la faille de San Andreas. Nous pensons par conséquent qu'il existe, sous la Californie, une zone de découplage entre la lithosphère et l'asthénosphère, permettant d'individualiser une déformation lithosphérique d'une déformation asthénosphérique. Nous estimons, en outre, que l'anisotropie observée dans l'asthénosphère sous la Californie ne peut être expliquée seulement par le cisaillement induit par le déplacement de la lithosphère Nord Amérique. En effet, les propriétés anisotropes obtenues par modélisation à partir d'une plaque se déplaçant dans une direction et une vitesse proche de celle de la plaque Amérique du Nord montrent qu'on ne peut espérer guère plus que quelques dixièmes de seconde de délai au bout de 10 Ma de déplacement. Les déphasages mesurés en Californie étant de l'ordre de 1,5 s, il est donc nécessaire d'invoquer la présence d'écoulements mantelliques actifs sous cette région / This work provides new constraints on the development and on the distribution of the deformation in the upper mantle and particularly beneath transform plate boundaries. USArray experiment and the remarkable increase of the dataset in California for the past ten years allowed us to scrutinize the lateral variations of the anisotropy in the vicinity of the San Andreas Fault zone. We have confirmed and increased the detection of two layers of anisotropy beneath this plate boundary. The first layer, located in the lithosphere, is related to the deformation induced at the fault, and the other one, located in the asthenosphere, is coherent with the anisotropy observed far from it, its origin is however less clear. We show that the deformation zone associated both to the San Andreas, Calaveras and Hayward Faults, is likely 40 km wide at 70 km depth. We then performed numerical thermomechanical modeling (ADELI) of the displacement of a transform plate boundary associated with the computation of the development of crystallographic fabrics using a viscoplastic self-consistent approach (VPSC). We analyzed the distribution of the deformation in the model ant looked after the possible interactions at depth between deformation caused at surface by the strike-slip dynamic of the fault and the shearing at the base of the lithosphere caused by the horizontal displacement of the plates. Elastic properties derived from the crystallographic fabrics modeled, show that such interactions exist and induce, beneath the fault zone, a progressive rotation of the crystallographic fabrics with depth. Seismological signature of these smooth rotations is however not relevant with the presence of two anisotropic layers as proposed beneath California. We thus consider that a decoupling zone exists between the lithosphere and the asthenosphere beneath the California to account for the sharp separation between a lithospheric and an asthenospheric deformation. We furthermore estimate that anisotropy observed far form the San Andreas Fault in California cannot be explained only by the drag of the asthenosphere by the North America lithosphere as proposed in our article. Indeed, we can only expect few tenths of second of splitting delay from the anisotropic properties derived from the numerical modeling of a plate moving in the same direction and in the same velocity than the North American lithosphere only for 10 Ma of displacement. As delays observed in California rather reach 1.5 s, anisotropy in this region thus requires the existence of an active asthenospheric flow to be explained.

The Earth's climate system is driven by varying insolation from the Sun. The dominant variations in insolation are at 23 and 40 thousand year periods, yet for the past million years the Earth's climate has glacial cycles at approximately 100 kyr periodicity. These cycles are a coupled variation in temperature, ice volume, and atmospheric CO₂. Somehow the Earth system's collective response to 23 and 40 kyr insolation forcing produces 100 kyr glacial-interglacial cycles. Generally it has been assumed that the causative mechanisms are a combination of ice dynamics (high ice reflectivity controlling temperature) and ocean circulation (changing carbon partitioning between the deep ocean and the atmosphere, and heat transport to the poles). However, these proposed mechanisms have not yet resulted in a compelling theory for all three variations, particularly CO₂. This thesis explores the role of volcanic CO₂ emissions in glacial cycles. I calculate that glacial-driven sea level change alters the pressure on mid-ocean ridges (MORs), changing their CO₂ emissions by approximately 10%. This occurs because pressure affects the thermodynamics of melt generation. The delay between sea level change and the consequent change in MOR CO₂ emissions is several tens-of-thousands-of-years, conceptually consistent with a coupled non-linear oscillation that could disrupt glacial cycles from a 40 kyr mode to a multiple of that period. I develop an Earth system model to investigate this possibility, running for approximately one million years and explicitly calculating global temperatures, ice sheet configuration, and CO₂ concentration in the atmosphere. The model is driven by insolation, with all other components varying in response (and according to their own interactions). This model calculates that volcanism is capable of causing a transition to ̃100 kyr glacial cycles, however the required average volcanic CO₂ emissions are barely within the 95% confidence interval. Therefore it is possible for volcanic systems and glacial cycles to form a 100 kyr coupled oscillation.

Cette thèse vise à proposer de nouvelles méthodes permettant de reconstruire la circulation dans le manteau terrestre et l'évolution de la tectonique de surface pour les deux cents derniers millions d'années. Nous utilisons des modèles numériques de convection mantellique dans lesquels la dynamique de surface est comparable à la tectonique terrestre. En combinant ces modèles avec des reconstructions de la tectonique des plaques il est possible d'estimerla structure et l'évolution du champ de température dans le manteau. Jusqu'à présent, l'inclusion des reconstructions de la tectonique des plaques se faisait en imposant des conditions aux limites du modèle (équilibre des forces, vitesses imposées...). Ces techniques, bien que permettant de tester la validité de différents scénarios tectoniques alternatifs, n'autorisent pas de rétroaction dynamique de la convection mantellique sur la tectonique de surface.Dans ce travail, nous avons développé des techniques d'assimilation de données permettant d'intégrer les reconstructions de la tectonique des plaques dans un modèle numérique tout en laissant se développer de manière auto-cohérente cette rétroaction. Les techniques développées permettent également de prendre en compte les incertitudes associées aux reconstructions de la tectonique des plaques et de calculer les erreurs sur l'estimation finale de la circulationmantellique.Dans un premier temps, nous avons développé un filtre de Kalman suboptimal qui permet d'estimer la structure et l'évolution de la circulation mantellique la plus probable à partir d'un modèle numérique de convection et d'une sérietemporelle d'observations de surface, ainsi que de leurs incertitudes respectives.Ce filtre a été testé sur des expériences synthétiques. Celles-ci consistent à tenter de retrouver une évolution témoin à partir d'une série temporelle de données issues de cette évolution. Ces expériences ont montré qu'il était possible, enprincipe, de reconstruire la structure et l'évolution de l'ensemble du manteau à partir d'observations de vitesses et de flux de chaleur à la surface.Dans un second temps, nous avons développé un filtre de Kalman d'ensemble. Ce filtre permet non seulement d'estimer de manière plus précise la géométrie des structures mantelliques, mais aussi les incertitudes sur cette estimation. / This dissertation focuses on the developpement of data assimilation methods to reconstruct the circulation of the Earth's mantle and the evolution of its surface tectonics for the last 200~Myrs. We use numerical models of mantle convection in which the surface dynamics is similar to the Earth's. By combining these models with plate tectonics reconstructions, it is possible to estimate the structure and evolution of the temperature field of the mantle. So far, the assimilation of plate tectonics reconstructions was done by imposing specific boundary conditions in the model (force balance, imposed velocities...). These techniques, although insightful to test the likeliness of alternative tectonic scenarios, do not allow the full expression of the dynamical feedback between mantle convection and surface tectonics. We develop sequential data assimilation techniques able to assimilate plate tectonics reconstructions in a numerical model while simultaneously letting this dynamicalfeedback develop self-consistently. Moreover, these techniques take into account errors in plate tectonics reconstructions, and compute the error on the final estimation of mantle circulation.First, we develop a suboptimal Kalman filter. This filter estimates the most likely structure and evolution of mantle circulation from a numerical model of mantle convection, a time series of surface observations and the uncertainty on both. This filter was tested on synthetic experiments. The principle of a synthetic experiment is to apply the data assimilation algorithm to a set of synthetic observations obtained from a reference run, and to then compare the obtained estimation of the evolution with the reference evolution. The synthetic experiments we conducted showed that it was possible, in principle, to reconstruct the structure and evolution of the whole mantle from surface velocities and heat flux observations.Second, we develop an Ensemble Kalman Filter. Instead of estimating the most likely evolution, an ensemble of possible evolutions are computed. This technique leads to a better estimation of the geometry of mantle structures and a more complete estimation of the uncertainties associated.