Éditeur(s) : HAL CCSD
Résumé : Joint inversions are now commonly used in Earth Sciences. They have been developed to better understand Earth structure since they provide more constraints on the inverted parameters. We propose a new process to simultaneously invert several data sets in order to better image 3-D crustal and upper mantle structures. Our inversion uses 3 kind of data which present a good complementarity : (1) P receiver functions to provide Moho depth variations, (2) teleseismic delay times of P waves to retrieve velocity anomalies in the crust and upper mantle, (3) gravity anomalies to image density variations at lithospheric scale. In our scheme, receiver functions are first inverted. The resulting Moho depths are then interpolated and incorporated as a priori information into the joint inversion of teleseismic delay times and gravity anomalies process. In our new approach, we perform a model space search for Moho variations, P-velocity and density structure to find acceptable fit to the three data sets. In order to preferentially sample the good data fit region, we chose the Neighborhood algorithm of Sambridge to optimistically survey the model space. We model the delay times with a 3-D raytracing using velocity nodes evenly spaced and linked with density nodes via a linear relationship. This relationship can evolve through depth, simulating its pressure dependance. We test our scheme on synthetic examples, and apply the method in a continental rifting region where deformation processes are complex and badly known. This joint inversion thus provides more constraints on the inverted parameters, and allows us to better image crust - mantle interactions in the case of the Hangai dome (Central Mongolia). Indeed, in this area, we expect an asthenospheric upwelling to dynamically support the topography, but we image a density - velocity anomaly in the lower crust. The effects on topography and the origin (composition vs. thermal) of this anomaly are still under debate.
Les inversions conjointes sont maintenant communément utilisées en Sciences de la Terre. Elles ont été développées afin d'améliorer notre connaissance et compréhension de la structure interne de la Terre puisqu'elles apportent de plus en plus de contraintes sur les paramètres inversés. Dans cette thèse, nous proposons un nouveau processus d'inversion conjointe qui prend en compte trois paramètres différents et qui mène à l'obtention d'un modèle lithosphérique tridimensionnel. Notre méthode utilise trois types de données gravimétriques et sismologiques qui présentent une bonne complémentarité : (1) Les fonctions récepteur P afin d'obtenir les variations de profondeur du Moho, (2) les délais de temps d'arrivées P des téléséismes pour retrouver les anomalies de vitesse dans la croûte et le manteau supérieur, (3) les anomalies gravimétriques qui donnent accès aux variations de densité à l'échelle lithosphérique. Dans le schéma d'inversion proposé, nous inversons d'abord les fonctions récepteur. Les variations de profondeur du Moho qui en résultent sont alors interpolées puis incorporées comme information a priori dans le processus d'inversion conjointe gravimétrie - tomographie télésismique. Que ce soit pour les fonctions récepteur ou pour la partie conjointe, notre méthode d'inversion se base sur un algorithme de recherche afin de trouver la structure qui minimise l'écart aux trois types de données. Nous avons opté pour l'algorithme de voisinage qui présente l'avantage de concentrer sa recherche dans les régions qui minimisent l'écart aux données sans toutefois abandonner les autres régions afin de limiter les risques de convergence vers un minimum local. Nous utilisons un modèle 3-D constitué de noeuds de vitesse afin de modéliser les délais de temps. Ces noeuds de vitesse sont liés aux noeuds de densité grâce à une relation linéaire entre vitesse et densité. Cette relation peut évoluer avec la profondeur de manière à simuler les effets dus aux variations de pression. Nous avons réalisé plusieurs tests synthétiques afin d'évaluer le comportement de notre nouvelle méthode. Ils ont permis de montrer que les variations du Moho ainsi que les anomalies de vitesse et densité sont bien retrouvées. Ils ont également montré l'importance de la prise en compte de ces trois types de données dans un même schéma d'inversion, notamment en ce qui concerne la faculté de discerner les variations d'interfaces des anomalies volumiques. Nous avons ensuite appliqué la méthode dans une région continentale en phase de rifting où les processus de déformation sont complexes et mal connus. Notre nouvelle méthode a apporté davantage de contraintes sur les paramètres inversés et a permis d'obtenir un modèle correctement résolu de la croûte et du manteau supérieur sous le dôme de Hangai (Mongolie Centrale). Une anomalie asthénosphérique sous le dôme est souvent invoquée pour expliquer sa topographie élevée mais c'est une anomalie (en vitesse et densité) crustale (croûte inférieure) qui s'exprime dans nos modèles. Ses effets sur la topographie et son origine (compositionnelle et/ou thermique) sont encore débattus et nécessiteraient des études plus poussées (paramétriques, VP =VS, ...).
https://tel.archives-ouvertes.fr/tel-00783524

tel-00783524
https://tel.archives-ouvertes.fr/tel-00783524
https://tel.archives-ouvertes.fr/tel-00783524/document
https://tel.archives-ouvertes.fr/tel-00783524/file/Basuyau2011.pdf

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
Mapping forest aboveground biomass (AGB) has become an important task, particularly for the reporting of carbon stocks and changes. AGB can be mapped using synthetic aperture radar data (SAR) or passive optical data. However, these data are insensitive to high AGB levels (>150 Mg/ha, and >300 Mg/ha for P-band), which are commonly found in tropical forests. Studies have mapped the rough variations in AGB by combining optical and environmental data at regional and global scales. Nevertheless, these maps cannot represent local variations in AGB in tropical forests. In this paper, we hypothesize that the problem of misrepresenting local variations in AGB and AGB estimation with good precision occurs because of both methodological limits (signal saturation or dilution bias) and a lack of adequate calibration data in this range of AGB values. We test this hypothesis by developing a calibrated regression model to predict variations in high AGB values (mean >300 Mg/ha) in French Guiana by a methodological approach for spatial extrapolation with data from the optical geoscience laser altimeter system (GLAS), forest inventories, radar, optics, and environmental variables for spatial inter-and extrapolation. Given their higher point count, GLAS data allow a wider coverage of AGB values. We find that the metrics from GLAS footprints are correlated with field AGB estimations (R 2 =0.54, RMSE=48.3 Mg/ha) with no bias for high values. First, predictive models, including remote-sensing, environmental variables and spatial correlation functions, allow us to obtain "wall-to-wall" AGB maps over French Guiana with an RMSE for the in situ AGB estimates of ~51 Mg/ha and R²=0.48 at a 1-km grid size. We conclude that a calibrated regression model based on GLAS with dependent environmental data can produce good AGB predictions even for high AGB values if the calibration data fit the AGB range. We also demonstrate that small temporal and spatial mismatches between field data and GLAS footprints are not a problem for regional and global calibrated regression models because field data aim to predict large and deep tendencies in AGB variations from environmental gradients and do not aim to represent high but stochastic and temporally limited variations from forest dynamics. Thus, we advocate including a greater variety of data, even if less precise and shifted, to better represent high AGB values in global models and to improve the fitting of these models for high values.
ISSN: 0303-2434

Droits : http://creativecommons.org/licenses/by/
hal-01381116
https://hal.archives-ouvertes.fr/hal-01381116
https://hal.archives-ouvertes.fr/hal-01381116/document
https://hal.archives-ouvertes.fr/hal-01381116/file/Fayad2016a.pdf
DOI : 10.1016/j.jag.2016.07.015

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
Joint inversions are now commonly used in the earth sciences. They have been developed to better understand the structure of the earth, since they provide more constraints on the inverted parameters. We propose a new process to simultaneously invert several data sets in order to better image 3D crustal and upper mantle structures. Our inversion uses three kinds of data that present good complementarity: (1) P-wave receiver functions to provide Moho depth variations, (2) teleseismic delay times of P-waves to retrieve velocity anomalies in the crust and the upper mantle, and (3) gravity anomalies to image density variations at the lithospheric scale. We use a stochastic scheme, where receiver functions are first inverted. The interpolated resulting Moho depths are incorporated as a priori information into the joint inversion of teleseismic delay times and gravity anomalies process. Moreover, velocity and density can be linked by empirical relationships, which justifies the joint inversion of those parameters. In our stochastic approach, we perform a model space search for Moho variations, P-velocity, and density structure to find an acceptable fit to the three data sets. In order to preferentially sample the good data fit regions, we chose the neighborhood algorithm of Sambridge to optimistically survey the model space. We model the delay times with 3D raytracing using evenly spaced velocity-density nodes. We present here the first results given by this method on synthetic tests. (C) 2011 Elsevier Ltd. All rights reserved.
ISSN: 0098-3004

hal-00643158
https://hal.archives-ouvertes.fr/hal-00643158
DOI : 10.1016/j.cageo.2010.11.017

Éditeur(s) : HAL CCSD European Geosciences Union
Résumé : International audience
We develop and implement a new method to take into account the impact of waves into the 3-D circulation model SYMPHONIE (Marsaleix et al., 2008, 2009a), following the simplified equations of Bennis et al. (2011) which use glm2z-RANS theory (Ardhuin et al., 2008b). These adiabatic equations are completed by additional parameterizations of wave breaking, bottom friction and wave-enhanced vertical mixing, making the forcing valid from the surf zone through to the open ocean. The wave forcing is performed by wave generation and propagation models WAVEWATCH III® (Tolman, 2008, 2009; Ardhuin et al., 2010) and SWAN (Booij et al., 1999). The model is tested and compared with other models for a plane beach test case, previously tested by Haas and Warner (2009) and Uchiyama et al. (2010). A comparison is also made with the laboratory measurements of Haller et al. (2002) of a barred beach with channels. Results fit with previous simulations performed by other models and with available observational data. Finally, a realistic case of energetic waves travelling over a coast of the Gulf of Lion (in the northwest of the Mediterranean Sea) for which currents are available at different depths as well as an accurate bathymetric database of the 0-10 m depth range, is then simulated. A grid nesting approach is used to account for the different forcings acting at different spatial scales. The simulation coupling the effects of waves and currents is successful to reproduce the powerful northward littoral drift in the 0-15 m depth zone. More precisely, two distinct cases are identified: when waves have a normal angle of incidence with the coast, they are responsible for complex circulation cells and rip currents in the surf zone, and when they travel obliquely, they generate a northward littoral drift. These features are more complicated than in the test cases, due to the complex bathymetry and the consideration of wind and non-stationary processes. Wave impacts in the inner shelf are less visible since wind and regional circulation seem to be the predominant forcings. Besides, a discrepancy between model and observations is noted at that scale, possibly linked to an underestimation of the wind stress. Lastly, this three-dimensional method allows a good representation of vertical current profiles and permits to calculate the shear stress associated with wave and current. Future work will focus on the combination with a sediment transport model.
ISSN: 1812-0806

hal-00766548
https://hal.archives-ouvertes.fr/hal-00766548
https://hal.archives-ouvertes.fr/hal-00766548/document
https://hal.archives-ouvertes.fr/hal-00766548/file/os-8-657-2012.pdf
DOI : 10.5194/osd-8-2417-2011

Éditeur(s) : HAL CCSD EDP Sciences
Résumé : International audience
Aggregation patterns are often visually detected in sets of location data. These clusters may be the result of interesting dynamics or the effect of pure randomness. We build an asymptotically Gaussian test for the hypothesis of randomness corresponding to a homogeneous Poisson point process. We first compute the exact first and second moment of the Ripley K-statistic under the homogeneous Poisson point process model. Then we prove the asymptotic normality of a vector of such statistics for different scales and compute its covariance matrix. From these results, we derive a test statistic that is chi-square distributed. By a Monte-Carlo study, we check that the test is numerically tractable even for large data sets and also correct when only a hundred of points are observed
ISSN: 1292-8100

hal-01502640
https://hal-agroparistech.archives-ouvertes.fr/hal-01502640
DOI : 10.1051/ps/2012027