乍得西南部中新生代火山作用与盆地演化

Located in the Central and West African, Chad, which is not well geological explored, is characterized by Mesozoic- Cenozoic intra-continental rift basins. The boreholes exposed that, during Mesozoic-Cenozoic times, volcanic activities were intense in these basins, but study on volcanic rocks is very weak, especially on those embedded in rift basins, and so far systematic and detailed work has still no carried out. Based on the project of China National Oil and Gas Exploration and Development Corporation, “The analysis of reservoir condition and the evaluation of exploration targets of seven basins in block H in Chad”, and the cooperative project between Institute of Geology and Geophysics, CAS and CNPC International (Chad) Co. Ltd., “Chronology and geochemistry studies on Mesozoic-Cenozoic volcanic rocks from southwestern Chad Basins”, systematic geochronology, geochemistry and Sr-Nd-Pb isotopic geochemistry studies on volcanic rocks from southwestern Chad basins have been done in the thesis for the first time. Detailed geochronological study using whole-rock K-Ar and Ar-Ar methods shows the mainly eruption ages of these volcanic rocks are Late Cretaceous- Paleogene. Volcanic rocks in the well Nere-1 and Figuier-1 from Doba basin are products of the Late Cretaceous which majority of the K-Ar (Ar-Ar) ages fall in the interval 95-75 Ma, whereas volcanic rocks in the well Ronier-1 from Bongor Basin and the Well Acacia-1 from Lake Chad Basin formed in the Paleogene which the ages concentrated in 66-52Ma. Two main periods of volcanic activity can be recognized in the study area, namely, the Late Cretaceous period and the Paleogene period. Volcanic activities have a general trend of south to north migration, but this may be only a local expression, and farther future studies should be carried on. Petrology study exhibits these volcanic rocks from southwestern Chad basins are mainly tholeiitic basalt. Major- and trace elements as well as Sr-Nd-Pb isotopic geochemistry studies show that the late Cretaceous and the Paleogene basalts have a definitely genetic relationship, and magmas which the basalts in southwestern Chad basins derived from were produced by fractional crystallization of olivine and clinopyroxene and had not do suffered from crustal contamination. These basalts are prominently enriched light rare earth elements (LREE), large-ion lithophile elements (LILE) and high field strength elements (HFSE) and depleted compatible elements. They have positive Ba, Pb, Sr, Nb, Ta, Zr, Hf anomalies and negative Th, U, P,Y anomalies. It is possible that the basalts from southwestern Chad basins mainly formed by mixing of depleted mantle (DM) and enriched mantle (EMⅡ) sources. The late Cretaceous basalts have higher (87Sr/86Sr)i ratios than the Paleogene basalts’, whereas have lower (143Nd/144Nd)i ratios than the latter, showing a significant temporal evolution. The mantle sources of the Late Cretaceous basalts may have more enriched mantle(EMⅡ) compositions, whereas those of the Paleogene basalts are relatively more asthenospheric mantle (DM) components. The mantle components with temporal change observed in basalts from Chad basins were probably correlated with the asthenospheric mantle upwelling and lithospheric thinning in Central and Western Africa since Mesozoic. Mesozoic- Cenozoic Volcanism in Chad basins probably is a product of intra- plate extensional stress regime, corresponded to the tectonic setting of the whole West and Central African during Cretaceous. Volcanism is closely correlated with rifting. As time passed from early period to late, the basaltic magma of Chad basins, characterized with shallower genetic depth, higher density and smaller viscosity, probably indicates the gradual strengthening evolution of the rifting. In the initial rife stage, volcanic activities are absent in the study area. Volcanic activities are basiccally corresponded with the strong extensional period of Chad basins, and the eruption of basalts was slightly lagged behind the extensional period. In the post-rift stage (30-0Ma), these basins shifted to the thermal sag phase, volcanic activities in the study area significantly decreased and then terminated.

南海张裂大陆边缘数值模拟研究

The South China Sea (SCS) is one of the largest marginal seas in the western Pacific, which is located at the junction of Eurasian plate, Pacific plate and Indian-Australian plate. It was formed by continent breakup and sea-floor spreading in Cenozoic. The complicated interaction among the three major plates made tectonic movement complex and geological phenomena very rich in this area. The SCS is an ideal place to study the formation and evolution of rifted continental margin and sea-floor spreading since it is old enough to have experienced the major stages of the basin evolution but still young enough to have preserved its original nature. As the demand for energy grows day by day in our country, the deep water region of the northern continental margin in the SCS has become a focus of oil and gas exploration because of its huge hydrocarbon potential. Therefore, to study the rifted continental margin of the SCS not only can improve our understanding of the formation and evolution processes of rifted continental margin, but also can provide theoretical support for hydrocarbon exploration in rifted continental margin. This dissertation mainly includes five topics as follows: (1) Various classic lithosphere stretching models are reviewed, and the continuous non-uniform stretching model is modified to make it suitable for the case where the extension of lithopheric mantle exceeds that of the crust. Then simple/pure shear flexural cantilever model is applied to model the basement geometries of SO49-18 profile in the northern continental margin of the SCS. By fitting the basements obtained by using 2DMove software with modeling results, it is found that the reasonable effective elastic thickness is less than 5km in this region. According to this result, it is assumed that there is weak lower crust in the northern continental margin in the SCS. (2) We research on the methods for stretching factor estimation based on various lithosphere stretching models, and apply the method based on multiple finite rifting model to estimate the stretching factors of several wells and profiles in the northern continental margin of the SCS. (3) We improve one-dimension strain rate inversion method with conjugate gradient method, and apply it to invert the strain rate of several wells in the northern continental margin of the SCS. Two-dimension strain rate forward modeling is carried out, and the modeling results show that effective elastic thickness is a key parameter to control basin’s geometry. (4) We simulate divergent upwelling mantle flow model using finite difference method, and apply this newly developed model to examine the formation mechanism of the northwest and central sub-basin in the SCS. (5) We inverse plate thickness and basal temperature of oceanic lithosphere using sea-floor ages and bathymetries of the North Pacific and the North Atlantic based on varied-parameters plate model, in which the heat conductivity, heat capacity and coefficient of thermal expansion depend on temperature or depth. A new empirical formula is put forward based the inversed parameters, which depicts the relation among sea-floor age, bathymetry and heat flow. Then various similar empirical formulae, including the newly developed one, are applied to examine the sea-floor spread issue in the SCS based on the heat flow and bathymetry data of the abyssal sub-basin.

新增生岩石圈与软流圈的相互作用：来自地幔橄榄岩捕虏体的证据

The petrology and geochemistry of peridotites entrained in Beiyan Cenozoic alkaline basalts within the middle segment of Tan-Lu fault zone and clinopyroxene megacrysts in the late Mesozoic and Cenozoic alkaline basaltic rocks from the North China Craton, have been systematically investigated. The main conclusions are obtained as follows. The peridotites entrained in alkaline basalts at Beiyan, Shandong Province, China are comprised of dominantly spinel lherzolites and spinel wehrlites with porphyroclastic, granuloblastic textures to resorption textures. The xenoliths are fertile in major element compositions (High CaO, TiO2, Low MgO, Cr2O3). The olivine Fo (= 100×Mg / (Mg+Fe) possesses a low and very large range of 81.0 to 91.0. The peridotites contain high percentages (Lherzolites: 10 - 19% in volume; Wehrlites: 24 - 28% in volume) of clinopyroxene with spongy textures. The Sr and Nd isotopic ratios of clinopyroxene separates from peridotites and pyroxenite xenoliths have a depleted and small range fall within the area of MORB, similar to newly-accreted lithospheric mantle. However, the appearance of many wehrlites and highly enriched LREE pattern suggest that this newly-accreted lithospheric mantle was considerably modified and reconstructed recently through the peridotite-asthenospheric melt interaction. The upwelling of asthenosphere from late Cretaceous to Eogene and upper mantle shearing of the Tan-Lu fault played an important role in the modification and reconstruction of the newly-accreted lithospheric mantle. The clinopyroxene megacrysts in the late Mesozoic and Cenozoic alkaline basaltic rocks from the eatern North China Craton are different in aspects of major elements, trace elements and isotopic composition. The characteristics of texture, mineral compositions and geochemistry as well as the Fe-Mg partitioning between the crystal and the melt indicates that the Al-augites in the Cenozoic basalts represent high-pressure crystallization products of alkaline basaltic melts. Thus, both of clinopyroxene megacrysts and host basalts could be derived from a same primitive magma. However, the Al-augites in the late Mesozoic basaltic rocks represent accidentally-included xenocrysts of basaltic components which had crystallized in the depth from a previously melting episode. The more depleted Sr-Nd isotopic compositions of Cenozoic megacrysts compared with those of host alkaline basalts and tholeiites demonstrate that even the alkali basalts could not completely represent primitive magma initiating in asthenosphere. That is to say, the Cenozoic alkaline basalts were more or less modified by some enriched Sr-Nd isotopic components during their eruption. Meanwhile, the tholeiites were not the products formed only by fractional crystallization of alkaline basaltic magma or different degrees of partial melting. It may result from the contribution of lithospheric mantle materials or crust contamination in magma chamber to alkali basaltic magmas.

藏北可可西里碰撞后钾质火山作用和岩石成因

Post-collisional, potassic igneous rocks are widely distributed in the Hoh Xil area of the northern Tibetan Plateau. Based on the field work, petrography, mineral chemistry, K-Ar geochronology, element and Sr-Nd-Pb isotope geochemistry, this thesis systematically studied the spatial and temporal distribution of the volcanic rocks, chemical characteristics, formation mechanism and partial melting mechanism of the magma source region, geodynamic setting of magmatism, as well as crustal assimilation and fractional crystallization (AFC). The results show that: 1. The Miocene (7.77-17.82 Ma) volcanic products dominantly are trachandesite and trachy, and subordinate rhyolites, associated with stike-slip faults and thrust faults, formed morphology of small lava platforms and cinder cones. 2. Phenocrysts in the lavas are augite, andesine, sanidine, calcic amphibole and subordinate orthopyroxene, biotite and Ti-Fe oxides, displaying typical quench texture. Equilibrium temperatures and pressures of clinopyroxene phenocrysts indicate the magma chamber is located in upper-middle crust. 3. Rhyolites are the products of crustal melting and fractionation of shoshonitic magmas. The source region of intermediate magmas is enriched continental lithospheric mantle, which contains residual minerals such as phlogopite, rutile and spinel, and enriched by subducted sediments during earlier multi-episodes of subduction. 4. Upwelling of asthenosphere provides heat for source region melting, and faults provide channels for magma eruption. 5. Northward underthrusting of Indian continental lithosphere and southward of backstop of Asian continental lithosphere resulted in upwelling of hot asthenosphere. Geochemical characteristics of the potassic magmatism in North Tibet are dominantly controlled by source region composition, partial melting, and crustal assimilation and fractional crystallization (AFC).

大别－苏鲁造山带与华北陆块东部岩石层三维速度结构

The 3-D velocity images of the crest and upper mantle beneath the region of 112° -124°E, 28°-39°N including the Dabie-Sulu orogenic belt are reconstructed by using 36405 P-wave arrivals of 3437 regional and 670 distant earthquakes during the period from 1981 to 1996, and gridding the area of 0.5° * 0.5°. The results of tomography demonstrate that: 1. The results of tomographic imaging show a broad heterogeneity in P wave velocity structure for the lithosphere beneath the Dabie-Sulu orogenic belt. 2. In the Dabie orogenic belt, the velocity patterns in the crust are different among various tectonic units. The Dabie and Qinling orogenic belts are remarkable in the tomographic images, and in mm the Hongan and Dabie blocks in the Dabie orogenic belt are also imaged very distinguishably. 3. A velocity (about 5.9～6.0 km/s) layer exists in the Dabie block at depth between 15～25 km, which is coincident with the low-resistance layer at the depth of 12-23 km, being inferred to be the tectonic detachment zone and suggesting that the extension detachment structure was formed in the middle crust. Beneath the southern and northerm Dabie tectonic units, the north-dipping high-velocity (at level of 6.5 ～ 6.6 km/s) block was developed in the crust, which might be correlated with the UHP rockswith low content of the meta-ultramafic rocks. This result is in agreement with the geological observation on the surface. 4. The velocity image at 40 km depth reveals the features at the top of mantle and the configuration of the Moho discontinuity. The depth of the Moho changes slightly along the trend of the orogenic belt. It in Hongan block is less than 40 km, but it is different in the western and eastern parts of the Dabie block, the former is more than 40 km, and the latter less than or equal to 40 km. The remnant of the mountain root exists between the Shangcheng-Macheng fault and the line of Huoshan-Yuexi-Yingshan in the Dabie orogenic belt, and beneath the southern and northern Dabie tectonic units. However, the thickness of the Moho is about 40 km and there is no obvious changes, which suggest that the Dabie orogenic belt has been experienced quite in the gravity equilibration. The Moho's depth in the Sulu is less than 40 km. 5. There is a dipping slab-like high-velocity body in the uppermost mantle. It is sandwiched by slow velocities and exists beneath the Dabie-Sulu orogenic belt in the range of depths between the Moho discontinuity and 110 km at least. This high-velocity body outlines a picture of the slab interpreted as the remnant of the Triassic subducted YZ. 6. The Sulu orogenic belt displays "crocodilian" velocity structure, the upper crust of the Yangtze thrusted over the Huabei crest, and the Huabei crust indented into the Yangtze crust, where the ancient subduction zone of the Yangtze lithosphere located. Based on the previous geological data, this structure is not related with the collision between the Yangtze and Sino-Korean Blocks, but caused by the sinistral offset of the Tan-Lu Fault. Studied on the velocity structure of the eastern Huabei lithosphere indicates: 1. The 'present-day' lithosphere of the eastern Huabei is between 40-100 km thick with greatly thinned lithosphere around the Bohai Sea. Generally, thickness of the lithosphere in this region decreased eastwards. 2. The attenuation of the lithosphere is attributed to the strongly uplift of the asthenosphere. In the area between the Taihang Mountains and the Tan-Lu Fault, there is a 'lever' with red low velocity belt, it is clearly defined, transverse continuity, depth between 100-150 km, local variations visible, and an upwards trend towards the Bohai Sea. Generally, the velocity structure in the mantle beneath the lithosphere displays irregular column-shape consisting of alternating high and low velocities, and when cold high velocity ancient lithosphere connects with the hot low velocity mantle materials forming precipitous compact structure. More heat pathways from the mantle occur towards the Tan-Lu Fault. 3. The strongly irregular characteristics of the contact between the asthenosphere and the lithosphere is induced by the long-term hot, chemical erosion and alteration on the contact. 4. There are still preserved high velocity lithosphedc root beneath Huabei with 'block-shape' distribution and surrounded by hot materials. Results of our studies indicate that the evolution models of the eastern China mantle are characterized by the direct contact between the uplifted lithosphere and the Huabei Craton accompanying the upwelling of the deep mantle materials. At the contact betwen the lithosphere and the asthenosphere, the upwelled mantle materials replaced and altered the lower lithosphere forming the metasome through the hot and chemical modifications impacted on the Craton lithosphere, and changed it into the lithosphere gradually, resulting in the lithospheric thinning. Thus, the lithospheric thinning is the result of the upwelling of the asthenosphere.

扬子板块西部大陆地幔地球化学特征及壳幔演化－有关基性岩超基性岩的地球化学研究

The continental mantle geochemical characteristics and crust-mantle evolution in the west of Yangtze Plate was discussed through the study of some within-plate basic-ultrabasic rocks from Lower Proterozoic to Later Paleozoic in this paper. In the Lower Proterozoic, the plate subduction between the pre-Tethys Proterozoic Ocean Plate and paleo-Yangtze Plate induced some basic volcanic formed in the island arc-back arc surrounding, which were represented by Ailaoshan Group-Dibadu Formation-Dahongshan Group, and there existed EM I component in the mantle source. The Middle Proterozoic Caiziyuan peridotite was formed in the epicontinental basin at the ocean-land boundary or within-continent rift basin. Its mantle source could be metasomatized by the dehydration fluid of subducted plate, and much initial radioactive ~(143)Nd was added to the source. In the Later Proterozoic, some rifts at the epicontinent or within-continent was formed due to the pre-Tethys oceanic plate subduction, and within-plate hot-spot Dahongshan diabase came into being. The whole-rock isochronal age of diabase is 1066±110Ma, and its mantle source was enriched Nd isotope and trace element which was related to the primary volatile component from asthenosphere and mantle plume. Its mantle source was included "FOZO" component representing mantle plume. The layer ultramafic rocks located at the Panxi Rift in the Middle-Later Paleozoic were resulted from different period and source. The early ultramafic indicated the incipient action of Panxi Rift, which is residue of continental lithospheric partial melting. Its mantle source involved subducted material and had distinct EM II component. The Emeishan basalt in the Later Paleozoic was typical continental flood basalt and its source also contained EM II component. The subduction of paleo-Tethys Ocean Plate provided essential dynamic condition for the large-scale opening of Panxi Rift, while the mantle plume supplied much material for Emeishan basalt. However, the plume was contaminated by the metasomatized continental mantle lithosphere in its upwelling process, which resulted in the Sr isotopic and incompatible elemental enrichment, and the Nd isotope kept down the weak-depleted character of mantle plume. The magmatic history in the west of Yangtze Plate is the tectonic process between pre-Tethys, paleo-Tethys Oceanic Plate and Yangtze Plate in a long history. Due to the subduction of oceanic plate, the crustal source material took part in the crust-mantle evolution widely. the continental mantle lithosphere in the west of Yangtze Plate was metasomatized by the fluid released by the subducted plate and the primary volatile from deeper mantle, and the mantle source include obvious enriched component.

华南震旦纪-寒武纪两期成磷事件及其地球动力学意义

Phosphorus is an important biological and ecological element that to a certain degree constrains ecological environment and nutrient (including carbon) cycling. Marine sedimentary phosphorites are the principal phosphorus supply of the mankind. In the eastern to southern margins of the Yangtze Craton, South China, there are two phosphogenetic events at the Doushantuo stage of the Late Sinian and the Meishucun stage of the Early Cambrian respectively, corresponding two explosion events of life across the Precambrian\Cambrian boundary. Phosphorus ores from the Sinian and Cambrian phosphate in South China can be classified roughly into two categories, namely, grained and non-grained phosphorites. Grained phosphorites, hosted in dolostone type of phosphogenetic sequences and with larger industrial values, occur mainly in margins of the Upper Yangtze Platform, formed in shallow-water environments with high hydraulic energy and influenced by frequent sea-level change. Non-grained phosphorites, hosted principally in black-shale type of phosphogenetic sequences and with smaller industrial values, are distributed mainly in the Jiangnan region where deeper-water sub-basins with low hydraulic energy were prevailing at the time of phosphogenesis. Secular change ofδ~(13)C, δ~(18) O, ~(86)Sr/~(87)Sr values of carbonates from Sinian and Cambrian sequences were determined. A negative abnormal ofδ~(13)C, δ~(18)O values and positive abnormal of 86Sr/87Sr values from the fossiliferous section of the Lowest Cambrian Meishucun Formation implies life depopulation and following explosion of life across the PrecambriamCambrian boundary. Based on a lot of observations, this paper put forward a six-stage genetic model describing the whole formational process of industrial phosphorites: 1) Phosphorus was transported from continental weathering products and stored in the ocean; 2) dissolved phosphates in the seawater were enriched in specific deep seawater layer; 3) coastal upwelling currents took this phosphorus-rich seawater to a specific coastal area where phosphorus was captured by oceanic microbes; 4) clastic sediments in this upwelling area were enriched in phosphorus because of abundant phosphorus-rich organic matters and because of phosphorus absorption on grain surfaces; 5) during early diagenesis, the phosphorus enriched in the clastic sediments was released into interstitial water by decomposition and desorption, and then transported to the oxidation-reduction interface where authigenic phosphates were deposited and enriched; 6) such authigenic phosphate-rich layers were scoured, broken up, and winnowed in shallow-water environments resulting in phosphate enrichment. The Sinian-Cambrian phosphorites in South China are in many aspects comparable with coastal-upwelling phosphorites of younger geological ages, especially with phosphorites from modern coastal upwelling areas. That implies the similarities between the Sinian-Cambrian ocean and the modern ocean. Although Sinian-Cambrian oceanic life was much simpler than modern one, but similar oceanic planktons prevail, because oceanic planktons (particularly phytoplanktons) are crucial for phosphate enrichment related to coastal upwelling. It implies also a similar seawater-layering pattern between the Sinian-Cambrian ocean and the modern ocean. The two global phosphate-forming events and corresponding life-explosion events at the Sinian and Cambrian time probably resulted from dissolved-phosphate accumulation in seawater over a critical concentration during the Earth's evolution. Such an oceanic system with seawater phosphorus supersaturation is evidently unstable, and trends to return to normal state through phosphate deposition. Accordingly, this paper put forward a new conception of "normal state <=> phosphorus-supersaturation state" cycling of oceanic system. Such "normal state <=> phosphorus-supersaturation state" cycling was not only important for the three well-known global phosphate-forming events, also related to the critical moments of life evolution on the Earth. It might be of special significance. The favorable paleo-oceanic orientation in regard to coastal-upwelling phosphorite formation suggests a different orientation of the Yangtze Craton between the Sinian time and the present time (with a 135° clockwise difference), and a 25° anti-clockwise rotation of the Yangtze Craton from late Sinian to early Cambrian. During the Sinian-Cambrian time, the Yangtze Craton might be separated from the Cathaysia Block, but might be still associated with the North China Craton.

地幔对流理论及其在解释地球物理观测的应用

With the improving of mantle convection theory, the developing of computing method and increasing of the measurement data, we can numerically simulate more clearly about the effects on some geophysical observed phenomenons such as the global heat flow and global lithospheric stress field in the Earth's surface caused by mantle convection, which is the primary mechanism for the transport of heat from the Earth's deep interior to its surface and the underlying force mechanism of dynamics in the Earth.Chapter 1 reviews the historical background and present research state of mantle convection theory.In Chapter 2, the basic conception of thermal convection and the basic theory about mantle flow.The effects on generation and distribution of global lithospheric stres s field induced by mantle flow are the subject of Chapter 3. Mantle convection causes normal stress and tangential stresses at the bottom of the lithosphere, and then the sublithospheric stress field induces the lithospheric deformation as sixrface force and results in the stress field within the lithosphere. The simulation shows that the agreement between predictions and observations is good in most regions. Most of subduction zones and continental collisions are under compressive. While ocean ridges, such as the east Pacific ridge, the Atlantic ridge and the east African rift valley, are under tensile. And most of the hotspots preferentially occur in regions where calculated stress is tensile. The calculated directions of the most compressive principal horizontal stress are largely in accord with that of the observation except for some regions such as the NW-Pacifie subduction zone and Qinghai-Tibet Plateau, in which the directions of the most compressive principal horizontal stress are different. It shows that the mantel flow plays an important role in causing or affecting the large-scale stress field within the lithosphere.The global heat flow simulation based on a kinematic model of mantle convection is given in Chapter 4. Mantle convection velocities are calculated based on the internal loading theory at first, the velocity field is used as the input to solve the thermal problem. Results show that calculated depth derivatives of the near surface temperature are closely correlated to the observed surface heat flow pattern. Higher heat flow values around midocean ridge systems can be reproduced very well. The predicted average temperature as a function of function of depth reveals that there are two thermal boundary layers, one is close to the surface and another is close to the core-mantle boundary, the rest of the mantle is nearly isothermal. Although, in most of the mantle, advection dominates the heat transfer, the conductive heat transfer is still locally important in the boundary layers and plays an important role for the surface heat flow pattern. The existence of surface plates is responsible for the long wavelength surface heat flow pattern.In Chapter 5, the effects on present-day crustal movement in the China Mainland resulted from the mantle convection are introduced. Using a dynamic method, we present a quantitative model for the present-day crustal movement in China. We consider not only the effect of the India-Eurasia collision, the gravitational potential energy difference of the Tibet Plateau, but also the contribution of the shear traction on the bottom of the lithosphere induced by the global mantle convection. The comparison between our results and the velocity field obtained from the GPS observation shows that our model satisfactorily reproduces the general picture of crustal deformation in China. Numerical modeling results reveal that the stress field on the base of the lithosphere induced by the mantle flow is probably a considerable factor that causes the movement and deformation of the lithosphere in continental China with its eflfcet focuing on the Eastern China A numerical research on the small-scale convection with variable viscosity in the upper mantle is introduced in Chapter 6. Based on a two-dimensional model, small-scale convection in the mantle-lithosphere system with variable viscosity is researched by using of finite element method. Variation of viscosity in exponential form with temperature is considered in this paper The results show that if viscosity is strongly temperature-dependent, the upper part of the system does not take a share in the convection and a stagnant lid, which is identified as lithosphere, is formed on the top of system because of low temperature and high viscosity. The calculated surface heat flow, topography and gravity anomaly are associated well with the convection pattern, namely, the regions with high heat flow and uplift correspond to the upwelling flow, and vice versa.In Chapter 7, we give a brief of future research subject: The inversion of lateral density heterogeneity in the mantle by minimizing the viscous dissipation.