海洋岩石圈板块有效弹性厚度研究

本文在前人研究大陆岩石圈板块有效弹性厚度的基础上,建立研究海洋岩石圈板块有效弹性厚度的理论模型,推导出与大陆岩石圈不同的海洋岩石圈板块响应函数 Z(k,Te) 理论计算公式.并分析海洋岩石圈板块响应函数 Z(k,Te) 的特点.文中对实际的海洋测量数据的响应函数 Z(k,Te) 进行计算和分析,估算我国南海南沙海域和南海中央海盆岩石圈板块有效弹性厚度分别约为10 km和6~7 km.     

Apparent velocity measurements on an oceanic lithosphere

A direct measurement of apparent velocities for oceanic paths was made with an array of sensitive ocean bottom seismographs. The measurement was performed by recording waves from shallow earthquakes which occurred in the area close to trench axes and which were accurately located by the land seismological network in Japan. The range of epicentral distances is from 500 to 1,800 km.The observed P travel times are less than those in the Jeffreys-Bullen tables by 6–10 s for the range of distances.Since the dimension of the OBS array is about 400 km, the apparent velocities are determined quite precisely and show little dependence on the epicentral distances. The average value of the apparent velocities for the range 500–1,700 km is 8.64 ± 0.13 km/s.An offset of travel times, which is thought to be associated with a low-velocity layer underneath the oceanic lithosphere, has been observed.These results indicate that a high-velocity layer with a velocity of 8.6 km/s exists in the lower part of the oceanic lithosphere. Beneath the 8.6-km/s layer there is a thin low-velocity layer under which the velocity of the P wave is again 8.6 km/s.     

Petrological models of the oceanic lithosphere: Geophysical and geochemical tests

Petrological models of the oceanic lithosphere are tested to satisfy geophysical and geochemical constraints within the framework of plate tectonics. Quartz eclogite, olivine eclogite, peridotite and dunite are considered as the material of the lithosphere. The temperature at the base of the lithosphere is assumed to be the solidus temperature. This temperature, the thermal conductivity, and the heat flow and topography changes with age are used as the geophysical constraints. The compressional wave velocity-depth profile is used to select preferred models. Among geophysically successful models, high-temperature models are preferred to wet low-temperature models, because the low-temperature models have difficulties in explaining the mechanism of generation of oceanic basalt magmas. A preferred model is a two-layer model 70 km thick consisting of peridotite at the upper lithosphere and olivine eclogite at the lower lithosphere bounded at the base by the dry solidus.     

A note on the lithosphere thickness and heat flow density of the Indian Craton from MAGSAT data

The Curie depth map of India compiled from MAGSAT data has been used for preparing the lithosphere thickness and the surface heat flow density maps of the Indian Craton, utilizing the concept of magnetothermometry. The lithosphere thicknesses of the major Indian geological units/provinces, as obtained from the prepared map, are found to be in reasonably good agreement with the previously published values for these regions. Also, the surface heat flow density values obtained from the prepared maps closely follow the previously published results. The maps are useful in providing first order estimates of lithosphere thickness and surface heat flow density of the important geological units/provinces of India.     

Age-dependent subduction of oceanic lithosphere beneath western South America

A recently established relation between the penetration depth of oceanic lithosphere and the lithospheric age appears to be of special interest to the understanding of the South American subduction zone. The main characteristics of this complicated zone, such as the absence of deep-focus earthquakes south of 30°S, the variations in the dip angle of the descending Nazca plate and the gap in seismic activity between depths of approximately 300 and 525 km, can be understood if the spatial and temporal variations in the age of the descending oceanic lithosphere are taken into account. In view of the significance of local aspects of the subduction process the South American-Nazca plate interaction cannot simply be considered as a type-example of the interaction between a continental and an oceanic plate.     

One thermal model of the young oceanic lithosphere: Direct problem

Summary The two-dimensional solution of the direct plate model of the oceanic lithosphere is computed for t ^1/2 dependences of its thickness using the multi-grid technique for the bilinear finite elements. If we assume the dependence of the lithospheric thickness on its age in the form h(t)=6·6t ^1/2 +7·5, then this model can roughly approximate the measured heat flow. In this model the horizontal conductive term is not negligible, the heat flow depends both on time and on the velocity of the sea-floor spreading which should be taken into account for the handling of heat flow measurements.

---

¶rt; u ¶rt; n a¶rt;au ¶rt; mu ¶rt;u au um auum mu m t ^1/2 . um mua muu¶rt;a ¶rt; uu m. u n¶rt;num auum mu m u u¶rt; h(t)=6.6t ^1/2 +7.5, m ma ¶rt; nuuum annuum a¶rt;a mn nm. m ¶rt;u numa ma mnn¶rt;mu; mn nm auum m u u m mu auu au ¶rt;a u, ¶rt;am, nu am uu mn nma a¶rt; umam uu mu auu.

     

Volcano spacing,fractures, and thickness of the lithosphere

Both “hot-spot” type and possibly island-arc volcanoes may form at the intersections of fractures whose spacing is near the thickness of the lithosphere and increases with increasing thickness. An approximate equality between layer thickness and spacing of major fractures observed in some sedimentary rocks and clay cake models may thus extend to the “mega-joints” that have fractured the lithosphere and controlled volcano spacing on the earth, and possibly on Mars. If the hot-spot fractures are interpreted as due to shear, many hot-spot fracture systems suggest roughly north-south least principal stress, or, alternatively in some instances, a 90° rotation of this pattern.     

Rheological properties of deep subducted oceanic lithosphere and their geodynamic implications

According to the experimental studies on the rheology of two important mantle rocks (eclogite and harzburgite), the rheological properties of the deep subducted oceanic lithosphere are investigated by assuming a simplified harzburgite type slab model with moderate thickness of basaltic layer. When the mantle convergence rate is small or the subducting slab has been trapped in the mantle for an enough long time, the strength profile of the slab is characterized by a strong subducting crustal component lying on a weak subducting upper mantle. However, if the convergence rate is large enough, the subducting slab will be featured only by a rigid cold center. Our study suggests that the detachment of the subducting crust component from the underlying upper mantle is only likely to happen in hot slow subducting slabs, but not the cold fast subducting lithosphere. Rheological properties of the harzburgitic and the eclogitic upper mantle vary with depths. The eclogitic upper mantle is stronger than the peridotitic upper mantle across the upper mantle. Transition zone is the high strength and high viscosity layer in the upper mantle except the lithosphere.     

Geoid anomalies across fracture zones and the thickness of the lithosphere

Recent advances in the measurement and interpretation of geoid height anomalies provide a new way to estimate the thickness of the oceanic lithosphere as a function of crustal age. GEOS-III satellite altimetry measurements show abrupt changes in sea level across fracture zones which separate areas of lithosphere with different ages. These changes have the correct location, amplitude, and wavelength to be caused by the combined gravitational attraction of the relief across the fracture zone and the isostatic support of this relief. Eight profiles of geoid height and bathymetry across the Mendocino fracture zone are inverted to determine the depth of the isostatic compensation, assuming that the compensation occurs in a single layer. These depths are then interpreted with a thermal boundary layer model of lithospheric growth. To explain satisfactorily the geoid measurements, the thermal diffusivity of the upper mantle must be 3.3 × 10^?3 cm² s^?1 and the thickness of the lithosphere, defined as the depth at which the geotherm reaches 95% of its maximum value, must be9.1km m.y.^?1/2 × t^1/2, where t is lithospheric age.     

Thermally induced brittle deformation in oceanic lithosphere and the spacing of fracture zones

Brittle deformation of oceanic lithosphere due to thermal stress is explored with a numerical model, with an emphasis on the spacing of fracture zones. Brittle deformation is represented by localized plastic strain within a material having an elasto-visco-plastic rheology with strain softening. We show that crustal thickness, creep strength, and the rule governing plastic flow control the formation of cracks. The spacing of primary crack decreases with crustal thickness as long as it is smaller than a threshold value. Creep strength shifts the threshold such that crust with strong creep strength develops primary cracks regardless of crustal thicknesses, while only a thin crust can have primary cracks if its creep strength is low. For a thin crust, the spacing of primary cracks is inversely proportional to the creep strength, suggesting that creep strength might independently contribute to the degree of brittle deformation. Through finite versus zero dilatation in plastic strain, associated and non-associated flow rule results in nearly vertical and V-shaped cracks, respectively. Changes in the tectonic environment of a ridge system can be reflected in variation in crustal thickness, and thus related to brittle deformation. The fracture zone-free Reykjanes ridge is known to have a uniformly thick crust. The Australian-Antarctic Discordance has multiple fracture zones and thin crust. These syntheses are consistent with enhanced brittle deformation of oceanic lithosphere when the crust is thin and vice versa.     

珠江口盆地大地热流特征及其与热岩石圈厚度的关系

沉积盆地现今热流特征是岩石圈构造-热演化过程的综合反映和盆地热史恢复的必要约束条件,其总体变化趋势与热岩石圈厚度密切相关.本文根据新收集的珠江口盆地19口钻井温度数据,新增计算了19个大地热流数据,其中12个数据位于深水区（水深大于300 m）,丰富了该盆地深水区钻井地热数据.结合前人研究成果,绘制了该盆地的大地热流图,并分析了其热流分布特征.在此基础上通过求解一维热传导方程,计算得到36口井位处的热岩石圈厚度,量化了盆地大地热流与热岩石圈厚度间的关系.结果显示,珠江口盆地大地热流值介于24.2~121.0 mW·m-2,平均71.8±13.6 mW·m-2,新增盆地深水区钻井平均热流值高达84.5±4.4 mW·m-2.大地热流分布整体上从陆架区到陆坡区升高,而热岩石圈厚度整体分布趋势与大地热流相反.大地热流与热岩石圈厚度间存在良好的指数相关性.     

Stability of the oceanic lithosphere with variable viscosity: an initial-value approach

We have studied the problem concerning the onset of convective instabilities below the oceanic lithosphere. A system of linear partial differential equations, in which the background temperature field is time-dependent, is integrated in time to monitor the evolution of incipient disturbances. Two types of rheologies have been examined. One depends strongly on temperature. The other involves a viscosity which is both temperature- and pressure-dependent. The results from this initial-value approach, in which the viscosity profiles migrate downward with time, reveal the importance of considering temperature- and pressure-dependent rheology in issues regarding the development of local instabilities in upper mantle convection. For temperature-dependent viscosity, viscosities of 0(10²⁰P) are required to produce instabilities with growth-rates of 0(.1/Ma). In contrast, these same growth rates can be attained for a temperature- and pressure-dependent viscosity profile with a mean value close to 0(10²⁰P) in the upper mantle, owing to the presence of a low viscosity zone, 0(10²⁰P), existing right below the lithosphere. Unlike the results of temperature-dependent viscosity, whose growth-rates increase with time, the amplification of disturbances in a fluid medium with temperature- and pressure-dependent rheology reaches a maximum at an early age, < 50 Ma, and decreases thereafter with time. This suggests the potential importance played by initial disturbances in the evolution of the oceanic lithosphere.     

The elastic thickness of the lithosphere in the Pacific Ocean

In this study, we present determinations of the effective elastic thicknessT_e of the oceanic lithosphere along Pacific chains or archipelagoes.T_e is determined by computing the deflection of a continuous elastic plate under the load of volcanoes, and constrained by geoid heights provided by SEASAT. In the South Central Pacific, estimates of 14 km for the Marquesas and 6 km or less for the Pitcairn-Mururoa-Gloucester chain are in good agreement with a previous work in this region (Cook-Austral and Society chains). Around the Line Islands chain, SEASAT data reveal that the bathymetry is poorly known, preventing fine analysis. Meanwhile,T_e looks globally very low ( 6 km), except for three volcanoes but these results may be unreliable. The Easter chain features lowT_e values ( 6 km), with no noticeable variation along the chain. Higher values are found for a Samoan island, Manuae (24 km), and along the Hawaiian-Emperor seamounts chain (from 32 km at the eastern end of the chain to 21.5 km for the Hawaiian volcanoes, and from 25.5 to 15 km for the Emperor seamounts). The large number ofT_e estimates obtained in this study points out a noticeable difference between North and South Pacific results. Those from the North Pacific agree with the general trend (increase with the square root of age plate at loading time), while those from the South Central Pacific are much lower, according to their plate age. These lowT_e results from the South Pacific are only partly explained by taking account of thermal perturbations using the rejuvenation model. Therefore, these results then point out a regional difference in oceanic lithosphere.     

Thermal cooling of the oceanic lithosphere: new constraints from geoid height data

Up to now, tests of thermal models of the oceanic lithosphere as it cools and moves away from the ridge crest have been based mainly on topography and heat flow data. However, large areas of the ocean floor deviate from the normal subsidence due to thermal contraction and heat flow data are not very sensitive to the form of the model.

Cooling of the lithosphere causes a short-wavelength step in the geoid across fracture zones that can also be used to constrain thermal models. We have analyzed geoid data at fracture zones from the SEASAT altimeter measurements in the entire Pacific Ocean and redetermined parameters of the cooling models. We find that the data reveal two distinct regimes of cooling; one for seafloor ages in the range 0–30 Ma, the other beyond 30 Ma; this does not appear to be correlated with particular fracture zones but rather it is representative of the whole area studied, i.e., the entire south Pacific and northeast Pacific Ocean. These two trends may be interpreted in terms of two different (asymptotic) thermal thicknesses of the plate model. The smaller thermal thickness ( 65 km) found for ages <30 Ma—compared to 90 km in the age range 30–50 Ma—calls for some kind of thermal perturbation in the vicinity of the ridge crest.

From the results obtained in this study, we conclude that the half-space cooling model is unable to explain the data, that beyond 30 Ma, a simple plate model gives a satisfactory fit to the data but in the younger plate portion (ages < 30 Ma) the cooling history of the oceanic lithosphere is much more complex than predicted by the usual cooling models. Furthermore, the depth-age relationship obtained from the geoid-derived thermal parameters departs significantly beyond 30 Ma from the widely used Parsons and Sclater's depth-age curve, predicting a lesser subsidence.     

Longshot experiments to study velocity anisotropy in the oceanic lithosphere of the northwestern Pacific

Several long-range explosion seismology experiments have been conducted in the northwestern Pacific basin, where one of the oldest oceanic lithospheres is postulated to exist. The experiments were conducted from 1974 to 1980. Highly sensitive ocean-bottom seismographs which had been developed for longshot experiments were used. The lengths of the profiles ranged from 1000 to 1800 km, and the directions were chosen to provide wide azimuthal coverage. One of the aims of this series of experiments was to test the existence of velocity anisotropy on a large, regional scale.The results show that the oceanic lithosphere has anisotropy wherein the velocity changes by 4–7%. The anisotropy extends from a depth of at least 40 to 140 km beneath the sea bottom; however, the magnitude of the anisotropy may vary with depth. The azimuth of the maximum velocity is 150–160° clockwise from north, and coincides with the “fossil” direction of spreading of the Pacific plate, whereas it differs from the present direction of plate motion by ～ 30°. The azimuth does not seem to depend on depth. In the direction of maximum velocity, the lithosphere is basically two-layered: 8.0–8.2 and 8.6 km s^?1. The depth of the interface is 50–60 km beneath the sea floor.     

A model of the lithosphere thickness in the region of the Carpathians

Summary Directionally independent average P residuals computed for waves of teleseismic events arriving under various azimuths and incidence angles provided the basis for estimating the lithosphere thickness beneath the Carpathians and their surroundings. A thin lithosphere (60–80 km) was determined for the Pannonian Basin and the Transylvanian Basin, the thickest lithosphere(about 180 km) beneath the South Carpathians at the contact with the Moesian Platform. In other parts of the Carpathian belt the lithosphere thickens beneath the outer parts towards the SW margin of the East-European and the Moldavian Platforms. The lithosphere thicknesses derived from P residuals correlate well with the magnetotelluric determinations of a layer of increased electrical conductivity in the upper mantle.

---

Резюме Осре?rt;ненные временные невязкu Р волн, незaвuсuмые оm нanрaвленuя u рaссчumaнные ?rt;ля у?rt;aленных землеmрясенuŭ uз волн, nрuхо?rt;ящuх nо?rt; рaзнымu aзuмуmaмu u у лaмu na?rt;енuя, ?rt;aюm основaнuе ?rt;ля оnре?rt;еленuя мощносmu лumосферы Кaрnam u uх окресmносmеŭ. Тонкaя лumосферa(60–80 км) оnре?rt;еленa ?rt;ля Пaннонско о u Трaнсuльвaнско о бaссеŭнов, a сaмaя мощнaя лumосферa(около 180 км) нaхо?rt;umся nо?rt; Южнымu Кaрnamaмu. В ?rt;ру uх чaсmях Кaрnamско о nоясa мощносmь лumосферы нaрaсmaеm во внешнuх зонaх в нanрaвленuu к ю о-зana?rt;ноŭ окрauне Восmочно-Евроnеŭскоŭ u Мол?rt;aвскоŭ nлamформ. Нaблю?rt;aеmся хорошaя корреляцuя мощносmu лumосферы, оnре?rt;еленноŭ нa основе невязок временu Р волн, с мa нumоmеллурuческuмu ?rt;aннымu nо лубuне слоя nовышенноŭ элекmроnрово?rt;uмосmu в верхнеŭ мaнmuu.

     

Interarc spreading and Cordilleran tectonics as alternates related to the age of subducted oceanic lithosphere

Old, cold oceanic lithosphere is denser and therefore gravitationally more unstable than younger, hotter oceanic lithosphere. Hence, whereas old lithosphere will sink under its own weight, subduction of young lithosphere may require an additional force. Interarc spreading occurs or occurred recently in the western Pacific, in the southern Atlantic, and possibly in the Mediterranean, where the subducted sea floor appears to be more than 50 m.y. old, and in many cases, is more than 100 m.y. old. In most of these regions, the ease with which the old dense lithosphere sinks may have contributed to a seaward migration of the trenches, which led to interarc spreading. Cordilleran tectonics, including high mountains and broad zones of deformation, are present on the margins of the eastern Pacific where the subducted oceanic lithosphere is younger than about 50 m.y. An extra force, which we presume to be necessary to cause subduction of the young lithosphere, may be responsible for the deformation and mountains just as an extra force seems necessary to drive continental collision in Asia. The extensive early Tertiary deformation across a broad zone of western North America may be related to the long-term, continuous subduction of young lithosphere of the Farallon and Kula plates.     

Effective elastic thickness of island arc lithosphere under Japan

Abstract Using topography and observed gravity anomalies, we have estimated the effective elastic thickness as a measure of strength of Japanese island arc lithosphere. The thickness is found to range from about 3 km to >20 km. The thickness seems to be controlled primarily by the thermal state of the lithosphere. The higher the heat flow, the thinner is the elastic plate. However, several areas show significant deviations. The smaller effective elastic thickness in the northern Ryukyu arc than that inferred from heat flow may be attributed to the stress regime. In Japan, extensional tectonics are going on only in the Ryukyu arc region. Shallow subducting slab under the south-western Japan frontal arc probably increases the effective thickness by several kilometers. The determined effective elastic thickness suggests that when we consider vertical movements in the volcanic arc, we should take account of topographic and subsurface loading over a few hundred kilometers. However, if the dip of the slab is shallow, the flexural responses of the underlying slab, not only that of the island arc lithosphere, should be taken into account for the compensation, as is the case of the south-western Japan frontal arc.     

Effect of the oceanic lithosphere velocity on free convection in the asthenosphere beneath mid-ocean ridges

The paper presents results obtained in experiments on a horizontal layer heated from below in its central part and cooled from above; the layer models the oceanic asthenosphere. Flow velocity and temperature profiles are measured and the flow structure under boundary layer conditions is determined (at Rayleigh numbers Ra > 5 × 10⁵). The flow in the core of a plane horizontal layer heated laterally and cooled from above develops under conditions of a constant temperature gradient averaged over the layer thickness. The flow core is modeled by a horizontal layer with a moving upper boundary and with adiabatic bounding surfaces under conditions of a constant horizontal gradient of temperature. Exact solutions of free convection equations are found for this model in the Boussinesq approximation. Model results are compared with experimental data. Temperature and flow velocity ranges are determined for the boundary layer regime. Based on the experimental flow velocity profiles, an expression is found for the flow velocity profile in a horizontal layer with a mobile upper boundary heated laterally and cooled from above. Free convection velocity profiles are obtained for the asthenosphere beneath a mid-ocean ridge (MOR) with a mobile lithosphere. An expression is obtained for the tangential stress at the top of the asthenosphere beneath an MOR and the total friction force produced by the asthenospheric flow at the asthenosphere-lithosphere boundary is determined.