Isostasy of aseismic tectonic units in the South Atlantic Ocean and geodynamic implications

Mechanical isostasy of the lithosphere in the South Atlantic Ocean was studied using information on gravity anomalies and bathymetry with additional constraints imposed by the altimetric geoid. The isostatic responses (admittances) over the Walvis Ridge, Rio Grande Rise and Trindade Seamount Chain were computed using a three-dimensional algorithm. The eastern Walvis Ridge and the Rio Grande Rise have the same response, which is well explained by an Airy model of isostasy. The other features are regionally supported. A variation in the thickness of the elastic plate was found along the western Walvis Ridge. A high value of the elastic plate thickness (20 km) was found under the Trindade Chain. Geodynamic implications are discussed in the light of these results.     

Isostasy of the Eastern Anatolia (Turkey) and Discontinuities of its Crust

This paper presents a probable isostatic model of the East Anatolian Region, which lies in a belt of significant plate movements. Probable locations of the horizontal and vertical discontinuities in the crust structure were determined using the normalized full gradient (NFG) method. For the purpose of explaining the mechanism that supports topography corresponding to the crust thickness in the region, calculations of effective elastic thickness (T _e) were carried out initially by utilizing admittance and misfit functions. According to these results, the effective elastic thickness value obtained was less than the crust thickness, even though the isostatic model does not conform with the Airy model. Consequently, it was assumed that there could be problems beneath the crust. Hence, the NFG method was applied on the Bouguer gravity data of the region in order to investigate probable discontinuities in the crust structure. According to the NFG results, vertical structural transitions were observed at a depth ranging between 10 and 30 km, which begin immediately north of the Bitlis Zagros Suture Zone (BZSZ) and continue in a northerly direction. The relationship between the effective elastic thickness (T _e; 13 km in average as determined in the last stage), and the seismogenic zone in the region was investigated. If the T _e value happens to be less then the crustal thickness, then one can say that there are problems in the crustal structure of the region similar to Eastern Anatolia. Indeed, when NFG results of the study area are examined, numerous vertical and horizontal discontinuities in the crust can be observed. These discontinuities, which correspond to low Bouguer gravity anomalies and shallow focal depth-earthquakes, are probably the source of the factors which rule the tectonic mechanism of the region.     

Origin and compensation of Chagos-Laccadive ridge, Indian ocean, from admittance analysis of gravity and bathymetry data

The Chagos-Laccadive ridge (CLR) is a prominent aseismic, volcanic ridge in the northern Indian ocean. The ridge, together with the Southern Mascarene plateau (SMP), to which it is genetically related, is considered as a volcanic trace of the Reunion hotspot. We have examined the isostatic compensation of the CLR through transfer function analysis of gravity and bathymetry data along seven profiles. The analysis suggests that the CLR is compensated locally, with an Airy crustal thickness (T_c) of 20 km. The rather low elastic plate thickness (T_e) of about 4 km implies that the volcanism of the ridge took place very near a spreading centre. The proximity of the Chagos fracture zone indicates that the emplacement was probably near a spreading centre-transform junction.     

Denudational isostatic rebound of intraplate highlands: The lachlan river valley,Australia

A range of evidence from the Lachlan valley in the southeast Australian highlands is consistent with Neogene isostatic rebound in response to denudational unloading. This evidence is found along the inland edge of the highlands in the transition zone between the highlands proper and the Lachlan's inland alluviated valley and the intracratonic Murray Basin. The amounts and rates of uplift indicated by offsets of suballuvium bedrock profiles and the long profiles of Tertiary valley-filling basalts are consistent with modelling of denudational rebound using known rates of highland denudation and basinal sedimentation, and reasonable crustal properties. The modelling shows that weak to moderately strong strong lithosphere (effective elastic thickness, T_e = 1-25 km) and strong lithosphere (T_e = 100 km) are all consistent with the observed amounts of rebound. Strong lithosphere must be broken, however, to be consistent with the field data. Even in the Australian setting, which is characterized by very low rates of denudation, isostatic rebound in response to denudational unloading must be a significant factor in maintaining highland elevation and must be incorporated in models of long-term landscape evolution. It would be expected that denudational isostatic rebound would be an even more significant component of long-term landscape evolution in areas of higher denudation rates.     

Insights into the Crustal Structure and Geodynamic Evolution of the Southern Granulite Terrain,India, from Isostatic Considerations

The Southern Granulite Terrain of India, formed through an ancient continental collision and uplift of the earth’s surface, was accompanied by thickening of the crust. Once the active tectonism ceased, the buoyancy of these deep crustal roots must have supported the Nilgiri and Palani-Cardamom hills. Here, the gravity field has been utilized to provide new constraints on how the force of buoyancy maintains the state of isostasy in the Southern Granulite Terrain. Isostatic calculations show that the seismically derived crustal thickness of 43–44 km in the Southern Granulite Terrain is on average 7–8 km more than that required to isostatically balance the present-day topography. This difference cannot be solely explained applying a constant shift in the mean sea level crustal thickness of 32 km. The isostatic analysis thus indicates that the current topography of the Southern Granulite Terrain is overcompensated, and about 1.0 km of the topographic load must have been eroded from this region without any isostatic readjustment. The observed gravity anomaly, an order of magnitude lower than that expected (−125 mGal), however, shows that there is no such overcompensation. Thermal perturbations up to Pan-African, present-day high mantle heat flow and low Te together negate the possible resistance of the lithosphere to rebound in response to erosional unloading. To isostatically compensate the crustal root, compatible to seismic Moho, a band of high density (2,930 kg m⁻³) in the lower crust and low density (3,210 kg m⁻³) in the lithospheric mantle below the Southern Granulite Terrain is needed. A relatively denser crust due to two distinct episodes of metamorphic phase transitions at 2.5 Ga and 550 Ma and highly mobilized upper mantle during Pan-African thermal perturbation reduced significantly the root buoyancy that kept the crust pulled downward in response to the eroded topography.     

The effective elastic lithosphere under the Cook-Austral and Society islands

We have determined the elastic thicknessT_e of the oceanic lithosphere along two volcanic chains of the South Central Pacific: Cook-Austral and Society islands. We used a three-dimensional spatial method to model the lithospheric flexure assuming a continuous elastic plate. The model was constrained by geoid height data from the SEASAT satellite.Along the Cook-Austral chain the elastic thickness increases westward, from 2–4 km at McDonald hot spot to 14 km at Rarotonga. At McDonald seamount, however, the data are better explained by a local compensation model. The observed trend shows an increase ofT_e with age of plate at loading time. However, the elastic layer under the Cook-Austral appears systematically thinner by several kilometers than expected for “normal” seafloor, suggesting that substantial thermal thinning has taken place in this region. Considering the apparent thermal age of the plate instead of crustal age improves noticeably the results. Along the Society chainT_e varies from 20 km under Tahiti to 13 km under Maupiti which is located 500 km westward. When plotting together the Society and Cook-AustralT_e results versus age of load, we notice that within the first five million years after loading,T_e decreases significantly while tending rapidly to an equilibrium value. This may be interpreted as the effect of initial stress relaxation which occurs just after loading inside the lower lithosphere and suggests that the presently measured elastic thickness under the very young Tahiti load ( 0.8 Ma) is not yet the equilibrium thickness.     

地幔岩石圈热结构差异与中国大陆岩石圈均衡分析

观测表明,大陆地区存在很多传统均衡模型无法解释的现象,其根本原因在于传统均衡理论中没有考虑地幔岩石圈部分由于热结构差异导致密度差异的影响.本文基于岩石圈尺度的质量平衡模型研究了中国大陆20个构造单元地壳及地幔岩石圈对地形海拔的贡献,以及各块体的均衡状态.计算结果表明,在一些地区,如塔里木盆地、北山和柴达木盆地,尽管岩石圈均衡模型和Airy模型得到了一致的海拔值,但岩石圈均衡模型更能体现均衡过程的物理本质;除青藏高原造山带外的多数块体,岩石圈均衡模型的计算结果更接近观测海拔和地表垂直运动状态;总体上,考虑地幔岩石圈热结构影响后,中国大陆各地区的均衡结果普遍优于传统的均衡模型.通过对均衡状态分析,我们得到以下主要结论：(1)构造稳定地区均衡程度较高;(2) 青藏高原及周边造山带现今地壳运动主要为区域构造过程及深部动力学过程所控制,均衡调整过程不是主要控制因素;(3) 现今地壳垂直运动比较明显的块体处于均衡调整阶段,地表垂直运动的大小反映了该区所受的均衡力作用的程度;(4) 构造稳定地区基于岩石层均衡的理论计算海拔与观测海拔之差值和现今地壳垂直运动速率有较好的相关性,据此我们可以通过均衡分析研究构造块体的运动趋势和动力学性质;(5) 地幔热结构对现今地形、海拔及地壳垂直运动有显著影响,在处理均衡问题时,地幔岩石圈热结构是我们必须考虑的重要因素.     

Some consequences of the differences in seismic properties in old and young continental plates for isostatic compensation

We present a short review of seismological data and show that a difference of upper mantle structure is clearly resolved between old and young continental plates. An Archean plate is approximately 1 second faster than a Paleozoic plate for P waves. Such a change in plate structure implies that the density of an Archean lithosphere should be greater than the density of a Paleozoic lithosphere, if they have the same chemical composition. Applying isostasy, we would expect a Precambrian plate to sink and its surface to stand at a level 8 km beneath that of a Paleozoic plate. We first check that the weight of a crustal column is not smaller for a Precambrian crust than for a Paleozoic crust. Explosion seismology data show that isostatic equilibrium is nearly achieved in the crust. Therefore a compensating mechanism has to be acting in the subcrustal lithosphere. Two hypotheses are presented. Provided that the chemical composition of the lower lithosphere does not vary with time, a simple condition for a continental plate not to sink when aging, would be that the density of the lower lithosphere be ?0.1 g/cm³ smaller than the density of the asthenosphere. Another condition that would prevent plates to sink would be a lowering of the density of the lower lithosphere with age due for example to a variation of its chemical composition and to a reduction in the mean atomic weight. Both hypotheses imply that isostatic equilibrium of old continental plates is maintained by a geochemical mechanism. Although thermal evolution is a key factor for explaining the evolution of continents, it has to be complemented by a hypothesis on the change of the chemistry of the lower continental lithosphere with respect to the asthenosphere and/or with respect to age.     

Crustal structure and origin of the Cape Verde Rise

The Cape Verde Islands are located on a mid-plate topographic swell and are thought to have formed above a deep mantle plume. Wide-angle seismic data have been used to determine the crustal and uppermost mantle structure along a ~ 440 km long transect of the archipelago. Modelling shows that ‘normal’ oceanic crust, ~ 7 km in thickness, exists between the islands and is gently flexed due to volcano loading. There is no direct evidence for high density bodies in the lower crust or for an anomalously low density upper mantle. The observed flexure and free-air gravity anomaly can be explained by volcano loading of a plate with an effective elastic thickness of 30 km and a load and infill density of 2600 kg m^− 3. The origin of the Cape Verde swell is poorly understood. An elastic thickness of 30 km is expected for the ~ 125 Ma old oceanic lithosphere beneath the islands, suggesting that the observed height of the swell and the elevated heat flow cannot be attributed to thermal reheating of the lithosphere. The lack of evidence for high densities and velocities in the lower crust and low densities and velocities in the upper mantle, suggests that neither a crustal underplate or a depleted swell root are the cause of the shallower than expected bathymetry and that, instead, the swell is supported by dynamic uplift associated with the underlying plume.     

Earthquake ground motion in Mexico City: An analysis of data recorded at Roma array

We analysed in detail three earthquakes recorded in a small-aperture accelerometric array in Mexico City, using the correlation of the records as a function of time along the accelerogram and frequency. Ground response is strongly conditioned by the fundamental period of the soft soils at the site of the array (T₀). Energy at periods longer than 2T₀ is guided by the crustal structure (with a thickness of 45 km). The wave field at periods between T₀ and 2T₀ also consists of surface waves but guided by the upper 2–3 km of volcanic sediments in central Mexico. For periods smaller than T₀, ground motion is uncorrelated among the stations. Our results indicate that seismic response of Mexico City, including its very long duration, results from deeply guided surface waves (between 2 and 45 km depth) interacting with the very local response of the soft surficial clay layer.     

Crustal structure variation along 30°N in the eastern Tibetan Plateau and its tectonic implications

We determined crustal structure along the latitude 30°N through the eastern Tibetan Plateau using a teleseismic receiver function analysis. The data came mostly from seismic stations deployed in eastern Tibet and western Sichuan region from 2004 to 2006. Crustal thickness and Vp/Vs ratio at each station were estimated by the H–k stacking method. On the profile, the mean crustal thickness and Vp/Vs ratio were found to be 62.3 km and 1.74 in the Lhasa block, 71.2 km and 1.79 near the Bangong–Nujiang suture, 66.3 km and 1.80 in the Qiangtang block, 59.8 km and 1.81 in the Songpan–Garze block, and 42.9 km and 1.76 in the Yangtze block, respectively. The estimated crustal thicknesses are consistent with predictions based on the topography and the Airy isostasy, except near the Bangong–Nujiang suture and in the Qiangtang block where the crust is 5–10 km thicker than predicted, indicating that the crust may be denser, possibly due to mafic underplating. We also inverted receiver functions for crustal velocity structure along the profile, which reveals a low S-wave velocity zone in the lower crust beneath the eastern Tibetan Plateau, although the extent of the low-velocity zone varies considerably. The low-velocity zone, together with previous results, suggests limited partial melting and localized crustal flow in the lower crust of the eastern Tibetan Plateau.     

The variation of crustal thickness across the Swiss Alps based on gravity and explosion seismic data

Summary Recently determined gravity anomalies along the NW-SE oriented Swiss Geotraverse from Basel to Bellinzona are used in combination with seismic refraction data to deduce a crustal section across the Swiss Alps. Topographic, Bouguer, free air, isostatic and geological corrections were applied to the data. Geological features considered in the corrections are the Swiss Molasse basin filled with sediments and the Ivrea body of high-density material. The resultant Bouguer anomaly over the Gotthard massif is 130 mgal lower than the Bouguer anomaly at the northern end of the profile near Basel. The Alpine region is associated with negative isostatic anomalies down to –20 mgal. The crustal thickness is found to increase gradually from the northern end of the profile (thicknessH=30 km) towards the Helvetic nappes at the northern margin of the Alps (H=38 km) and more rapidly towards the Gotthard massif (H=50 km) and further south to Biasca down to a depth of 58 km. From Biasca southward the crustal thickness thins quite rapidly to reach a depth of 30 km at the southern end of the profile near Bellinzona. Thus the Alps have a distinct asymmetric crustal root whose maximum thickness is almost twice the average crustal thickness in Central Europe. With the Mohorovii-discontinuity deduced from seismic observations an average constant density contrast of –0.33 gcm^–3 is found between the lower crust and upper mantle underneath the Alps.Institut für Geophysik, ETH Zürich, Contribution No. 130.     

Erosional unloading,hillslope geometry,and the height of the Cascade Range,Washington State,USA

Peaks in the Cascade Range in northern Washington State are on average ～800 m higher than in southern Washington. The influences of differential valley excavation and variations in hillslope length and average slope on these altitudinal trends were tested using a 3‐dimensional model for isostatic rock uplift and calculations of hillslope length and slope respectively. The magnitude of isostatic peak uplift calculated by the model is highly dependent on the flexural rigidity (D) and the related effective elastic thickness (T_e) of the crust of this region. Crustal rigidity was constrained using published estimates and by estimating the depth of the seismogenic zone in the area (D > 1 × 10²³ Nm and T_e > 24 km). With these constraints, isostatic compensation due to differential erosion added < 700 m and 300 m, or < 25% overall, of height to peaks in the northern and southern Cascades, respectively. Deeper valley incision in the northern Cascades accounts for < 300 m of the 800 m difference in peak altitudes between north and south. Similarly, variation in valley spacing and slope account for < 350 m of the difference in mean altitude between northern and southern regions. Hence, at least several hundred m difference in altitude between the northern and southern regions of the Cascades in Washington must be due to tectonic, geologic, or geophysical factors rather than surficial and geomorphic effects like isostatic response to valley incision and hillslope geometry. Copyright © 2009 John Wiley & Sons, Ltd.     

Lithospheric model with thick oceanic crust at the continental boundary: A mechanism for shallow spreading ridges in young oceans

In a general lithospheric model of a simple divergent ocean and continental margin that satisfies the constraints of isostasy and gravity anomalies, the free-air gravity anomaly at the margin is modelled by an oceanic crust that thickens exponentially toward the margin from its common value of 6.4 km about 600 km from the margin to 17.7 km at the margin; this postulated thickening is supported empirically by seismic refraction measurements made near continental margins. The thickness of the oceanic crust matches that of the continental lithosphere at breakup, as observed today in Afar and East Africa, and is interpreted as the initial oceanic surface layer chilled against the continental lithosphere. With continued plate accretion, the chilled oceanic crust thins exponentially to a steadystate thickness, which is achieved about 40 m.y. after breakup. These findings contrast with the generally held view that the oceanic crust has a uniform thickness.During the first 40 m.y. of spreading, the thicker oceanic crust, of density 2.86 g/cm³, displaces the denser (3.32 g/cm³) subjacent material; by isostasy, the spreading ridge and the rest of the seafloor thus stand higher in younger( <40m.y.) oceans than they do in older(>40m.y.) oceans. This is postulated to be the cause of the empirical relationship between the crestal depth of spreading ridges and the age (or half-width) of ocean basins.     

Intraplate seismicity in the western Bohemian Massif (central Europe): A possible correlation with a paleoplate junction

Locations of the Eger Rift, Cheb Basin, Quaternary volcanoes, crustal earthquake swarms and exhalation centers of CO₂ and ³He of mantle origin correlate with the tectonic fabric of the mantle lithosphere modelled from seismic anisotropy. We suggest that positions of the seismic and volcanic phenomena, as well as of the Cenozoic sedimentary basins, correlate with a “triple junction” of three mantle lithospheres distinguished by different orientations of their tectonic fabric consistent within each unit. The three mantle domains most probably belong to the originally separated microcontinents – the Saxothuringian, Teplá-Barrandian and Moldanubian – assembled during the Variscan orogeny. Cenozoic extension reactivated the junction and locally thinned the crust and mantle lithosphere. The rigid part of the crust, characterized by the presence of earthquake foci, decoupled near the junction from the mantle probably during the Variscan. The boundaries (transitions) of three mantle domains provided open pathways for Quaternary volcanism and the ascent of ³He- and CO₂-rich fluids released from the asthenosphere. The deepest earthquakes, interpreted as an upper limit of the brittle–ductile transition in the crust, are shallower above the junction of the mantle blocks (at about 12 km) than above the more stable Saxothuringian mantle lithosphere (at about 20 km), probably due to a higher heat flow and presence of fluids.     

Isostatic status in North China (1) Method and local compensation

Following Airy and Pratt principles, five kinds of local-compensation models are analysed and a rapid 3-D gravity formula is utilized to calculate isostatic anomalies for 66 models with different parameters. It is noted that isostatic gravity maps appear nearly identical in their main patterns and features. The optimum compensation model in North China is one of modified Airy models in which the different density distribution in the surface, upper crust and lower crust is taken into account and the standard crustal thickness is about 50km. The position of the complete compensation interface is located in the upper mantle. The North China platform as a whole is under sub-isostatic equilibrium status with an isostatic anomaly of about 18·10⁻⁵ m/s² on an average. The distribution of isostatic gravity anomaly shows an obvious blockwise pattern. Most positive anomaly areas occur over the eastern part, the Jiao-Liao Block, Mt. Yan block and northern margin of the Hebei-Shandong block, whereas a negative area occurs in the Shanxi graben. The comparison of models indicates that the Moho discontinuity is not suitable to be taken as a compensation interface, and the compensation effects in Airy model are better than that in Pratt model, which is consistent with the feature of dominant layered structure and less lateral inhomogeneity in crust. Some results about composite compensation, the basic characteristics of isostatic anomaly and deep stucture will be published later in the second part of this paper. Wang Bowen took part in some work in this paper.     

喜马拉雅地区与天山地区地壳强地震的孕震因素探讨

利用哈佛大学GCMT数据中心和前人积累的历史地震资料(1962~2016年M W>4.0地震)以及Crust2.0地壳结构统计分析了喜马拉雅地区、天山地区的地壳区域构造与地震活动间的相关性。此外,利用GFZ地学研究中心提供的静态卫星重力模型GGM03S/EGM2008和地形模型Topo计算了2个地区的各类重力异常场,同时还模拟了不同地壳弹性参数下的重力异常场,结果表明喜马拉雅地区重力异常场在水平、垂直方向的梯度特征远大于天山地区的异常特征,且喜马拉雅地区的有效弹性板厚度Te(6~15km)小于天山地区的有效弹性板厚度Te(20~30km)。最后,利用喜马拉雅地区与天山地区的GPS震间三维形变场约束了断层运动模型,结果显示两者主前缘断裂的断层闭锁深度及应力积累状态存在较大的差异。因此认为,造成青藏高原及邻区的边界地壳区域地震活动性差异的动力学因素,与地壳有效弹性板厚度、孕震断层参数及区域应力积累状态等密切相关。     

Gravity anomaly, lithospheric structure and seismicity of Western Himalayan Syntaxis

A compiled gravity anomaly map of the Western Himalayan Syntaxis is analysed to understand the tectonics of the region around the epicentre of Kashmir earthquake of October 8, 2005 (Mw = 7.6). Isostatic gravity anomalies and effective elastic thickness (EET) of lithosphere are assessed from coherence analysis between Bouguer anomaly and topography. The isostatic residual gravity high and gravity low correspond to the two main seismic zones in this region, viz. Indus–Kohistan Seismic Zone (IKSZ) and Hindu Kush Seismic Zones (HKSZ), respectively, suggesting a connection between siesmicity and gravity anomalies. The gravity high originates from the high-density thrusted rocks along the syntaxial bend of the Main Boundary Thrust and coincides with the region of the crustal thrust earthquakes, including the Kashmir earthquake of 2005. The gravity low of HKSZ coincides with the region of intermediate–deep-focus earthquakes, where crustal rocks are underthrusting with a higher speed to create low density cold mantle. Comparable EET (∼55 km) to the focal depth of crustal earthquakes suggests that whole crust is seismogenic and brittle. An integrated lithospheric model along a profile provides the crustal structure of the boundary zones with crustal thickness of about 60 km under the Karakoram–Pamir regions and suggests continental subduction from either sides (Indian and Eurasian) leading to a complex compressional environment for large earthquakes.     

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.     

Body-wave wave-form modelling and source parameters for the indochina border earthquakes

Wave-form modelling of body waves has been done to study the seismic source parameters of three earthquakes which occurred on October 21, 1964 (M _b =5.9), September 26, 1966 (M _b =5.8) and March 14, 1967 (M _b =5.8). These events occurred in the Indochina border region where a low-angle thrust fault accommodates motion between the underthrusting Indian plate and overlying Himalaya. The focal depths of all these earthquakes are between 12–37 km. The total range in dip for the three events is 5°–20°. TheT axes are NE-SW directed whereas the strikes of the northward dipping nodal planes are generally parallel to the local structural trend. The total source durations have been found to vary between 5–6 seconds. The average values of seismic moment, fault radius and dislocation are 1.0–11.0×10²⁵ dyne-cm, 7.7–8.4km and 9.4–47.4 cm, respectively whereas stress drop, apparent stress and strain energy are found to be 16–76 bars, 8.2–37.9 bars and 0.1–1.7×10²¹ ergs, respectively. These earthquakes possibly resulted due to the tension caused by the bending of the lithospheric plate into a region of former subduction which is now a zone of thrusting and crustal shortening.