Crustal shear-wave velocity structure beneath northeast India from teleseismic receiver function analysis

We investigated the seismic shear-wave velocity structure of the crust beneath nine broadband seismological stations of the Shillong–Mikir plateau and its adjoining region using teleseismic P-wave receiver function analysis. The inverted shear wave velocity models show ∼34–38 km thick crust beneath the Shillong Plateau which increases to ∼37–38 km beneath the Brahmaputra valley and ∼46–48 km beneath the Himalayan foredeep region. The gradual increase of crustal thickness from the Shillong Plateau to Himalayan foredeep region is consistent with the underthrusting of Indian Plate beyond the surface collision boundary. A strong azimuthal variation is observed beneath SHL station. The modeling of receiver functions of teleseismic earthquakes arriving the SHL station from NE backazimuth (BAZ) shows a high velocity zone within depth range 2–8 km along with a low velocity zone within ∼8–13 km. In contrast, inversion of receiver functions from SE BAZ shows high velocity zone in the upper crust within depth range ∼10–18 km and low velocity zone within ∼18–36 km. The critical examination of ray piercing points at the depth of Moho shows that the rays from SE BAZ pierce mostly the southeast part of the plateau near Dauki fault zone. This observation suggests the effect of underthrusting Bengal sediments and the underlying oceanic crust in the south of the plateau facilitated by the EW-NE striking Dauki fault dipping 30⁰ toward northwest.     

Crustal thickness and velocity structure beneath Singapore’s seismic network

We estimated the crustal thickness and velocity structure beneath the five stations comprising the Republic of Singapore’s seismic network. Our data set was composed of 697 teleseismic receiver functions and 7 months of broad-band data that was cross-correlated to produce inter-station Green’s functions. Surface wave group velocities were extracted from the Green’s functions to obtain dispersion data for a path from central Sumatra to Singapore in order to provide a complimentary data set to the receiver functions. Crustal thickness was estimated via an H − k stacking technique, and high-resolution 1D P-wave velocity profiles were generated beneath each station by jointly inverting receiver function stacks and the group velocity data using a linearised time-domain inversion scheme. Crustal thickness beneath four stations was found to be between 28.0 km and 32.0 km, while one station in the northeast of Singapore indicates 24.0 km thick crust. This implies a significant crustal thinning beneath Singapore over the lateral extent of 50.0 km. Inversion results exhibit several crustal features that are observable in the derived models at all five stations, indicating that they exist across Singapore as a whole. There appears to be an upper-crustal high-velocity zone beneath Singapore, underlain by a velocity inversion. Station NTU shows slower near-surface velocities than the other stations, consistent with its situation above the sedimentary Jurong formation. These results expand the available global velocity data set, as well as being useful for assessing the seismic hazard in Singapore.     

Moho depth variation beneath southwestern Japan revealed from the velocity structure based on receiver function inversion

The Philippine Sea plate is subducting under the Eurasian plate beneath the Chugoku-Shikoku region, southwestern Japan. We have constructed depth contours for the continental and oceanic Mohos derived from the velocity structure based on receiver function inversion. Receiver functions were calculated using teleseismic waveforms recorded by the high-density seismograph network in southwestern Japan. In order to determine crustal velocity structure, we first improved the linearized time-domain receiver function inversion method. The continental Moho is relatively shallow ( 30 km) at the coastline of the Sea of Japan and at the Seto Inland Sea, and becomes deeper–greater than 40 km–around 35°N and 133.8°E. Near the Seto Inland Sea, a low-velocity layer of thickness 10 km lies under the continental Moho. This low-velocity layer corresponds to the subducting oceanic crust of the Philippine Sea plate. The oceanic Moho continues to descend from south to northwest and exhibits complicated ridge and valley features. The oceanic Moho runs around 25 km beneath the Pacific coast and 45 km beneath the Seto Inland Sea, and it extends to at least to 34.5°N. The depth variation of the Moho discontinuities is in good qualitative agreement with the concept of isostasy. From the configurations of both the continental and oceanic Mohos, we demonstrate that the continental lower crust and the subducting oceanic crust overlap beneath the southern and central part of Shikoku and that a mantle wedge may exist beneath the western and eastern part of Shikoku. The southern edge of the overlapping region coincides with the downdip limit of the slip area of a megathrust earthquake.     

The three-dimensional shear velocity structure of lithosphere in the southern Baikal rift system and its surroundings

The three-dimensional shear velocity lithospheric structure at depths from 0 to 70 km beneath the southern Baikal rift system and its surroundings has been imaged by inversion of P-to-SV receiver functions from 46 digital stations operated in two teleseismic international projects in southern Siberia and Mongolia. The receiver functions were determined from teleseismic P waveforms and inverted to obtain depth dependences of S velocities at each station which were related to tectonic structures. The computed vertical and horizontal sections of the 3D shear velocity model imaged a transition from relatively thin crust of the southern Siberian craton to thicker crust in the folded area south and southeast of Lake Baikal, with a local zone of thin crust right underneath the South Baikal basin. The velocity structure beneath the Baikal rift, the mountains of Transbaikalia, Mongolia, and the southern craton margin includes several low-velocity zones at different depths in the crust. Some of these zones may record seismic anisotropy associated with mylonite alignment along large thrusts.     

Lithospheric structure of the Tien Shan region constrained by joint inversion of receiver function and Rayleigh wave dispersion

We obtain a lithospheric shear‐wave velocity model across the Tien Shan orogenic belt by jointly inverting Rayleigh wave group velocities and teleseismic P‐wave receiver functions at 61 broadband seismic stations deployed in this region. Our new model reveals prominent lateral variations of shear‐wave velocity in both the crust and uppermost mantle. This model reveals different structures in the upper and middle crust across the Talas Fergana Fault, which may suggest the presence of a tectonic boundary between the western and central Tien Shan beneath the fault. According to the velocity images, the depth extent of the fault is ~40 km and this is confined to the crust. Pronounced low‐velocity anomalies are imaged in the middle crust and uppermost mantle beneath the southern and middle Tien Shan, implying that the upwelling of the materials from the upper mantle could have played an important role in the mountain building.     

Lower-crust S-wave velocity beneath western Yunnan Province from waveform inversion of dense seismic observations

We use teleseismic body waveforms to explore S-wave layered velocity structures beneath 30 portable digital seismic stations deployed around western Yunnan Province. Results show that the Moho depth in this region is ∼40 km and decreases in general from north to south, consistent with previous geophysical studies. Associated with this lateral variation of the Moho depth, the lower crust above the Moho discontinuity has a 15–25 km thick zone with an S-wave velocity lower than that of the upper crust. This lower velocity zone might be interpreted as a lower crust weak channel, which may mechanically partially decouple the upper-crust deformation from the underlying mantle. Thus, the inverted S-wave velocity structure could provide new evidence for the lateral flow of lower crust in the build-up of the south-eastern Tibetan plateau.     

The Seismic Structure of Wilkes Land/Terre Adelie, East Antarctica and Comparison with Australia: First Steps in Reconstructing the Deep Lithosphere of Gondwana

A determination of the seismic structure of the crust and uppermost mantle of East Antarctica, in the region of Casey station, Wilkes Land and Dumont DUrville station, Terre Adelie is presented. High-fidelity waveforms from teleseismic earthquakes recorded at stations CASY and DRV (1996-2001) are used to calculate the seismic receiver function, the signal produced as energy passes through layers in the seismic velocity structure under the receiving station. The receiver functions are stacked to improve the signal-to-noise ratio and then modelled using an inverse algorithm to find the structure that best fits the observed waveform at each station. Inferences are made regarding the tectonic structure, in particular, the crustal thickness and character of the seismic Moho.The crustal thickness under Casey Station is found to be 30 km (+/- 2 km) with a fairly sharp Moho, considerably less than Dumont D'Urville Station, where the crustal thickness is 42 km, and there IS a significant low velocity region the deep crust. The structure of the Wilkes Land lithosphere is comparable to that of the Albany-Fraser Orogen, Western Australia, part of its conjugate margin. This places a new constraint on the relative position of East Antarctica and Australia in the reconstruction of Gondwana, and earlier, supercontinents. A recent reinterpretation of Antarctic geology proposes tectonic province boundaries trending perpendicular to the coast with counterparts in southern Australia. Seismic techniques, determining structure beneath regions with no surface exposure, are vital tools in testing such tectonic hypotheses, towards the reconstruction of Gondwana to full lithospheric depth.     

Crustal and lithospheric thinning beneath the seismogenic Kachchh rift zone,Gujarat (India): Its implications toward the generation of the 2001 Bhuj earthquake sequence

A high-resolution passive seismic experiment in the Kachchh rift zone of the western India has produced an excellent dataset of several thousands teleseismic events. From this network, 500 good teleseismic events recorded at 14 mobile broadband sites are used to estimate receiver functions (for the 30–310° back-azimuth ranges), which show a positive phase at 4.5–6.1 s delay time and a strong negative phase at 8.0–11.0 s. These phases have been modeled by a velocity increase at Moho (i.e. 34–43 km) and a velocity decrease at 62–92 km depth. The estimation of crustal and lithospheric thicknesses using the inversion of stacked radial receiver functions led to the delineation of a marked thinning of 3–7 km in crustal thickness and 6–14 km in lithospheric thickness beneath the central rift zone relative to the surrounding un-rifted parts of the Kachchh rift zone. On an average, the Kachchh region is characterized by a thin lithosphere of 75.9 ± 5.9 km. The marked velocity decrease associated with the lithosphere–asthenoshere boundary (LAB), observed over an area of 120 km × 80 km, and the isotropic study of xenoliths from Kachchh provides evidence for local asthenospheric updoming with pockets of partial melts of CO₂ rich lherzolite beneath the Kachchh seismic zone that might have caused by rifting episode (at 88 Ma) and the associated Deccan thermal-plume interaction (at 65 Ma) episodes. Thus, the coincidence of the area of the major aftershock activity and the Moho as well as asthenospheric upwarping beneath the central Kachchh rift zone suggests that these pockets of CO₂-rich lherzolite partial melts could perhaps provide a high input of volatiles containing CO₂ into the lower crust, which might contribute significantly in the seismo-genesis of continued aftershock activity in the region. It is also inferred that large stresses in the denser and stronger lower crust (at 14–34 km depths) induced by ongoing Banni upliftment, crustal intrusive, marked lateral variation in crustal thickness and related sub-crustal thermal anomaly play a key role in nucleating the lower crustal earthquakes beneath the Kachchh seismic zone.     

Crustal structure beneath Kottamiya broadband station,northern Egypt,from analysis of teleseismic receiver functions

We analyzed a total of 206 receiver functions beneath Kottamiya broadband station in northern Egypt to study the crustal structure and any azimuthal variations in the crustal thickness. The computed receiver functions are subdivided according to their azimuth into eight subgroups and analyzed separately using a genetic algorithm. The genetic algorithm is more appropriate than conventional linearized inversion schemes in regions where there is little a priori information about local crustal structures such as northern Egypt because it does not strongly depend on an initial model. The study region is located on the unstable shelf of Egypt in the northeastern corner of Africa. Little information about the deep structure of the crust beneath this region is available. For this reason, we have adopted the genetic algorithm to seismic waveform data recorded by Kottamiya broadband station. The crustal thickness varies slightly from 32 to 34 km with an average of 32.25 km, which is consistent with previous studies in the region. The crustal thickness shows a tendency of decrease toward the east and northeast being consistent with the general tectonic setting of the region including the opening of the Red Sea in the Tertiary times. Nonetheless, more teleseismic receiver functions from earthquakes recorded at denser seismic stations in northern Egypt and the southeastern Mediterranean combined with surface wave dispersion data as well as other geophysical investigations are necessary for more detailed imaging of the crustal structure which will deepen our understanding of the current tectonic and seismic activities of the region.     

Local earthquake tomography of the Andes at 20°S: Implications for the structure and building of the mountain range

Arrival-times of local events recorded in northern Chile and southern Bolivia were used to determine the P velocity structure above the subducted Nazca plate. The data were recorded between June and November 1994 by the French “Lithoscope” network: 41 vertical and 14 three-component short-period seismic stations were installed along a 700 km long profile crossing the main structures of the Andean chain, from the Coastal Cordillera to the Subandean Zone. The inversion method used is a modified version of Thurber’s 3D iterative simultaneous inversion code. The results were compared with a model obtained from previous German nearby refraction seismic studies and supplemented by field geological observations.The relocated seismicity is consistent with an ∼30° dipping slab between 0 and 170 km depth. We found a variation of about 30 km of the Moho depth along the profile. The crustal thickness is about 47 km under the Coastal Cordillera, 70 km under the Western Cordillera and the western part of the Eastern Cordillera, and 60–65 km beneath the Altiplano. Close to the surface, a good agreement between the velocity model and the geological structures is observed. Generally, in the upper crust, high velocities coincide with zones where basement is present near the surface. Low velocities are well correlated with the presence of very thick sedimentary basins or volcanic material. At greater depth, the trend of the velocity model is consistent with the existence of asymmetrical west-dipping imbricated blocks, overthrusting toward the east, which explain the asymmetrical pattern of the sedimentary basins. Beneath the Western Cordillera, the active volcanic arc, a large zone of low velocity is observed and interpreted to be due to partially molten material. A clear velocity contrast appears between the western and eastern parts of the upper mantle beneath the Andes; this geometry suggests the existence of a low velocity wedge in the mantle above the slab and the presence of a thick old lithosphere in the eastern part of the Andes.     

Crustal structure beneath two seismic broadband stations revealed from teleseismic P-wave receiver function analysis in the Virunga volcanic area,Western Rift Valley of Africa

The shear velocity structure beneath the Virunga volcanic area was estimated by using an average solution in the time domain inversion of stacked teleseismic receiver functions provided by two seismic broadband stations KUNENE (KNN) and KIBUMBA (KBB). These two stations are 29 km apart and located at the eastern and western escarpment of the Western Rift Valley of Africa in the Virunga area, respectively. The velocity model was presented as P-wave velocity models. From these models, the crust mantle transition zone beneath the area sampled by KNN and KBB in the Virunga area was determined at depth from about 36 to 39 km and 30 to 41 km, respectively. A low velocity zone was observed below stations KNN and KBB at depths between 20–30 km and 18–28 km, respectively, and with average velocity 5.9 km/s and 6.0 km/s. This low velocity zone may probably related to a magma chamber or a melt-rich sill. The models show also high velocity material (6.8–7.4 km/s) lying beneath stations KNN and KBB at depths 3–20 km and 3–10 km, respectively, which is indicative of magma cumulates within the volcanic edifice. The result obtained in this study was applied to the determination of epicentres during the period prior to the 27 November 2006 Nyamuragira eruption. This eruption was preceded by a swarm of hybrid volcanic earthquakes with clear P-waves onset. Using the receiver function model was found to improve the location of events. The located events correlate well with the location of the eruptive site and data provided by the INSAR observations of surface deformation associated with eruption.     

Crustal structure of the Central Precordillera of San Juan,Argentina (31°S) using teleseismic receiver functions

The subduction of the Nazca plate under the South American plate around 31°S is characterized by flat slab geometry. The (Chilean) Pampean flat slab of Argentina associated with the subduction of the Juan Fernandez ridge lies in a region of a series of foreland uplifts corresponding to the thin-skinned Precordillera and basement cored Sierras Pampeanas ranges. The SIEMBRA project deployed 40 broadband stations in 2008–2009 in both the Precordillera and the Sierras Pampeanas with the aim to foster the understanding of the entire central Andean flat slab region. One of the SIEMBRA station (DOCA) located on the western flank of Sierra de la Invernada in the Central Precordillera appears particularly appropriate to study the crustal structure and eventually detect discontinuities related to terranes establishment. We thus performed a receiver function analysis using teleseismic data recorded at the DOCA station during the SIEMBRA project and from October 2011 to June 2012 using a broadband UNSJ (National University of San Juan) seismic station with the purpose to obtain crustal images with details of the intracrustal structure consistent with a mechanism that could explains both the observed earthquake depths and the uplift pattern in the Central Precordillera. Our results show that the Moho beneath the Precordillera lies at a depth of about 66 km. The Moho signal appears diminished and behaves irregularly as a function of azimuthal orientations. Although this observation could be the result of an irregular geometry it also correlates with the hypothesis of partial eclogitisation in the lower crust. Two mid-crustal discontinuities have also been revealed. The shallower one could correspond to a décollement level between the Precordilleran strata and the Cuyania basement at 21 km depth. The deeper one which the presence has been matched with a sharp decrease of the crustal seismic activity drove us to the hypothesis of a major change in crustal composition at 36 km. Finally the flat portion of the subducted slab has been imaged lying at about 100 km depth.     

The crustal structure beneath the Central Andean forearc and magmatic arc as derived from seismic studies — the PISCO 94 experiment in northern Chile (21°–23°S)

In this study, we present an interpretation of seismic refraction profiles from the PISCO 94 experiment in northern Chile. As the PISCO experiment was a combined active and passive seismological study, we also discuss results of the passive part in the context of the seismic refraction model. Previous seismic refraction and gravimetric studies indicate a maximum crustal thickness of about 70 km beneath the Pre- and Western Cordillera. The new seismic refraction data lead to a differentiated image of the Andean crust which shows strong varying characteristics. The crustal discontinuities (up to five are detected) dip from W to E. The upper crust has a thickness of 18 km (Precordillera) to 23 km (magmatic arc) underlain by the recent middle crust down to 35–45 km where the velocity increases to about 7 km/s at its base. This crustal level is interpreted as old continental lower crust and its base as blurred continental (paleo) Moho. Beneath the Precordillera, a strong discontinuity at 70 km depth with a velocity increase to about 8 km/s was detected, interpreted as the recent geophysical Moho. For the magmatic arc, this deep discontinuity could not be found by active seismic measurements. The tomographic models of the seismological studies, in general, confirm the seismic refraction results. Anomalously high v_p/v_s ratios in the deeper part of the forearc indicate a hydrated mantle wedge consisting of serpentine and amphibole-bearing peridotite and the 70 km discontinuity is interpreted as the boundary between these two different stages of the hydrated mantle wedge. A zone of high attenuation (Qp) and high v_p/v_s ratios beneath the magmatic arc coincides with the low velocity zones and indicates partially molten rocks from a depth of 20 km down to the asthenospheric wedge.     

Estimation of crustal discontinuities from reflected seismic waves recorded at Shillong and Mikir Hills Plateau,Northeast India

In this study, an attempt is made to determine seismic velocity structure of the crust and upper mantle beneath the Shillong-Mikir Hills Plateau in northeast India region. The principle of the technique is to relate seismic travel times with crustal thickness above the Conrad and Moho discontinuities. Broadband digital waveforms of the local earthquakes make a precise detection of the seismic phases possible that are reflected at these discontinuities. The results show that the Conrad discontinuity is at 18–20 (±0.5) km beneath the Shillong-Mikir Hills Plateau and the Moho discontinuity is at 30 ± 1.0 km beneath the Shillong Plateau and at 35 ± 1.0 km beneath the Mikir Hills.     

利用天然地震震相探讨阿尔金地区地壳结构

本文利用阿尔金地区的宽频地震数据,对布设在该区的１０个宽频地震台站用接收函数方法进行了速度结构反演,反演的初步结果发现,若至花土沟剖面在２０ｋｍ深度处有一条厚度达５～１０ｋｍ的低速带断续出现,莫霍界面呈台阶状展布,北部浅,南部深;塔里木盆地南缘的地壳厚度为４０～４２ｋｍ左右;在阿尔金南,北缘断裂两侧台站下方莫霍深度的错断约６．５～８ｋｍ,在柴达木盆地北缘,莫霍面的深度达５０ｋｍ以上,Ｓ波速为４．５     

长江中下游成矿带及邻区地壳结构——MASH成矿过程的P波接收函数成像证据？

利用长江中下游成矿带多学科深部探测剖面于2009年11月至2011年3月间采集的天然地震数据,通过天然地震接收函数成像等分析研究,得到了研究区地壳和上地幔结构的清晰图像。接收函数成像结果显示研究区内Moho面深度存在着明显的起伏变化,在长江中下游成矿带(指剖面穿过的长江中下游成矿带宁芜矿集区,下同)下方存在着"幔隆构造"。在剖面东南端(即扬子克拉通北缘),Moho面相对稳定,深度约为30km;在茅山和江南断裂附近,Moho面存在上下起伏现象;在剖面中部或宁芜矿集区下方,Moho面存在明显隆起,深度只有28km;在郯庐断裂带下方,Moho面明显加深,深度达到36km;进一步向北到华北地台南缘,Moho面深度逐渐恢复到了32km左右的平均深度水平。其次,我们在接收函数成像结果中发现,长江中下游成矿带与其周边下地壳结构存在着明显的差异,成矿带的下地壳具有显著的地震波方位各向异性。扬子克拉通北缘的下地壳呈高速的近水平状结构,地震波各向异性特征不明显;与此相比,长江中下游成矿带的下地壳虽然也呈近水平状结构特征,但是,对于沿成矿带走向方向传播的地震波,其下地壳具有高速特征,而对于垂直于成矿带走向方向上传播的地震波,其下地壳却又表现为低速特征,这意味着成矿带的下地壳存在着平行于成矿带走向(即近北东—南西)方向的地震波各向异性,我们解释其是下地壳熔融并沿成矿带走向水平流动导致矿物晶体定向排列的结果。最后,在郯庐断裂以西的华北地台南缘观测到一条从上地壳延伸到中下地壳的南南东向倾斜的转换震相,我们推测它可能是合肥盆地内地壳伸展构造的反映。此外,我们发现接收函数成像结果中观测到的"幔隆构造"与远震P波层析成像结果在成矿带下方150km深度上显示的上地幔低速异常(江国明等,另文发表)存在着良好的对应关系,我们解释它们是软流圈物质上涌的遗迹。综合天然地震接收函数成像、远震P波层析成像和前人关于岩浆岩等方面的研究成果,我们认为长江中下游成矿带现今的下地壳可能是中生代发生成矿作用的多级岩浆房系统的一部分,成矿带的形成可能是类似MASH过程的产物。首先,软流圈物质上涌导致了长江中下游成矿带及其周边拉张环境的形成,在其上部地壳中形成了一系列伸展构造;然后,软流圈物质通过底侵进入长江中下游成矿带的原下地壳并与原下地壳物质发生同化作用,形成类埃达克质岩浆;接着,类埃达克质岩浆沿着伸展、拆离构造上升到地壳浅部形成不同层次的岩浆房和侵入岩体,并与围岩作用形成矿床。     

Evidence of a partial melt zone beneath Stromboli volcano (Italy) from inversion of teleseismic receiver functions

Teleseismic body waves from broadband seismic stations are used to investigate the crustal and uppermost mantle structure of Stromboli volcano through inversion of the receiver functions (RFs). First, we computed RFs in the frequency domain using a multiple-taper spectral correlation technique. Then, the non-linear neighbourhood algorithm was applied to estimate the seismic shear wave velocity of the crust and uppermost mantle and to define the main seismic velocity discontinuities. The stability of the inversion solution was tested using a range of initial random seeds and model parameterizations. A shallow Moho, present at depth of 14.8 km, is evidence of a thinned crust beneath Stromboli volcano. However, the most intriguing and innovative result is a low S velocity layer in the uppermost mantle, below 32 km. The low S velocity layer suggests a possible partial melt region associated with the volcanism, as also recently supported by tomographic studies and petrological estimations.     

Lateral growth of NE Tibetan Plateau restricted by the Asian lithosphere: Results from a dense seismic profile

We present high-resolution receiver function images along a 700-km long dense seismic array extending from northern Tibetan Plateau to the Alxa block, crossing the entire Qilian thrust belt (QTB). The dense stations, with less than ~2 km station intervals, allow the receiver functions to unveil unprecedented details of crustal structures across the northern frontier of the growing Tibetan Plateau. The migration image shows a thickened and strongly deformed QTB crust, with an uneven Moho and complex internal structures that are indicative of pure-shear shortening. The Alxa block, in contrast, has a thinner crust, a flat Moho, and little internal crustal deformation. These results suggest that the lateral growth of NE Tibetan Plateau is restricted by the strong Asian lithosphere, which shows no visible subduction beneath the Tibetan Plateau as previously suggested.     

Crustal thickening processes in the Central Andes and the different natures of the Moho-discontinuity

New seismic data from the Central Andes allow us to clarify the crustal structure of this mountain chain and to address the problem of crustal thickening. Evidence for the deep crustal root can be observed in both gravimetric and seismological data. Crustal structure and composition change significantly from east to west. In the eastern part of the backarc the Moho discontinuity is clearly recognisable. However only poor Moho arrivals are observed by active seismic measurements beneath the Altiplano and the Western Cordillera where broad-band seismology data indicate such a discontinuity. In the Precordillera, a pronounced discontinuity is detected at a depth of 70 km. Along the coast, the oceanic Moho is developed at a depth of 40 km. There are several processes which can change the petrological and petrophysical properties of the rocks forming the crust. Variations of the classical Moho discontinuity are presented which do not correspond to the petrological crust/mantle boundary. Tectonic shortening in the backarc is the dominant process contributing to at least 50–55% to the root formation along 21°S. In the forearc and arc, hydration of the mantle wedge produced ≈15–20% of crustal thickening. Magmatic thickening and tectonic erosion contributed only ≈5%. The other ≈25% is not yet explained.     

Crustal complexity in the Lachlan Orogen revealed from teleseismic receiver functions

There is an ongoing debate about the tectonic evolution of southeast Australia, particularly about the causes and nature of its accretion to a much older Precambrian core to the west. Seismic imaging of the crust can provide useful clues to address this issue. Seismic tomography imaging is a powerful tool often employed to map elastic properties of the Earth's lithosphere, but in most cases does not constrain well the depth of discontinuities such as the Mohorovi?i? (Moho). In this study, an alternative imaging technique known as receiver function (RF) has been employed for seismic stations near Canberra in the Lachlan Orogen to investigate: (i) the shear-wave-velocity profile in the crust and uppermost mantle, (ii) variations in the Moho depth beneath the Lachlan Orogen, and (iii) the nature of the transition between the crust and mantle. A number of styles of RF analyses were conducted: H-K stacking to obtain the best compressional–shear velocity (V _P /V _S) ratio and crustal thickness; nonlinear inversion for the shear-wave-velocity structure and inversion of the observed variations in RFs with back-azimuth to investigate potential dipping of the crustal layers and anisotropy. The thick crust (up to 48 km) and the mostly intermediate nature of the crust?mantle transition in the Lachlan Orogen could be due to the presence of underplating at the base of the crust, and possibly to the existing thick piles of Ordovician mafic rocks present in the mid and lower crust. Results from numerical modelling of RFs at three seismic stations (CAN, CNB and YNG) suggest that the observed variations with back-azimuth could be related to a complex structure beneath these stations with the likelihood of both a dipping Moho and crustal anisotropy. Our analysis reveals crustal thickening to the west beneath CAN station which could be due to slab convergence. The crustal thickening may also be related to the broad Macquarie volcanic arc, which is rooted to the Moho. The crustal anisotropy may arise from a strong N–S structural trend in the eastern Lachlan Orogen and to the preferred crystallographic orientation of seismically anisotropic minerals in the lower and middle crust related to the paleo-Pacific plate convergence.