The Gondwana supercontinent was an accretion of several cratons from different landmasses, namely South America, Antarctica, Africa, Madagascar, Australia, and the Indian subcontinent. The splitting of Gondwanaland during the Mesozoic led to the gradual rifting of these different cratons over geologic time. In this study, crustal structures are imaged by modeling receiver functions to understand the differences in the nature of the crust that was once part of Precambrian Gondwanaland. On comparing the overall crustal thickness with the age of the cratons, it was found that average bounds of crustal thickness varied from ~33 to 45 km in the Precambrian cratons of different ages, and composition varied from felsic to intermediate (Vp/Vs1.651.78). Observations of gradational Moho beneath few stations could indicate the possibility of mafic underplating at some point in their history of formation, growth, or evolution. Even if plate tectonics were dominant in the middle to late Archean, difference in spreading, drifting velocity, and distance travelled by the continents after Gondwana separation possibly led to crustal delamination, and destruction of thick crustal roots of cratons. Other than delamination, the role of episodic cycles of crustal growth is also observed in the pattern of crustal thickness across each division of the Precambrian. This could be due to alternative high and low crustal regeneration process, repeated episodic recycling, and reworking of early crust, supported by previous geodynamical models.

Cratons are stable portions of the Earth’s crust that have undergone very little deformation for prolonged time periods [1]. The Precambrian cratons have survived the cycles of merging and rifting of continents over geologic time, making them unique. In order to understand how the cratons formed and evolved, investigation of the lithospheric structure is essential. During the Gondwana period, the landmasses of South Africa, Antarctica, Australia, Madagascar, Sri Lanka, and the Indian subcontinent were all joined near the South Pole, making the Gondwana supercontinent an accretion of several cratons (Figure 1). The breakup of Gondwanaland during the Mesozoic led to the gradual rifting of these different cratons over geologic time at different rates. Sri Lanka had also moved along with the Dharwar craton and shared similar crustal properties with its Gondwana neighbours [2]. Cratonic roots are supposed to have a very strong thick lithosphere with thick crust [3]. But after the Gondwana breakup, the cratons have migrated along with the continents at different speeds, different directions, and some travelling very large distances. One of the important concerns is to what extent the crustal properties were affected during the continental evolution history after the Gondwana breakup. Another concern is whether there is a link between evolution of cratons and their seismic structure.

Figure 1

Archean cratons (brown colored) of South Africa, Antarctica, Australia, Madagascar, and the Indian subcontinent during the Gondwana period, used in the study.

Figure 1

Archean cratons (brown colored) of South Africa, Antarctica, Australia, Madagascar, and the Indian subcontinent during the Gondwana period, used in the study.

Models of crustal growth can be broadly categorized as crust based and mantle based. The crust-based models focus on cumulative crustal material at the present day, but the mantle-based models also shed light on crust that has been lost to mantle by reworking, subduction, or delamination, hence focusing on net crustal growth [4]. These studies are also aimed at estimating the formation age and how crustal processes lead to deviations in surface age, formation age, and average age of the crust. Cratons are the oldest regions on the earth preserving old crustal material, and because of this fact, investigating the crustal architecture and crustal properties with respect to the age of the cratons is crucial in understanding crustal evolution. Opinion remains divided in the scientific community with respect to crustal properties and ages of crustal material. One of the models by Smithson et al. [5] indicates that the crustal thickness and mean seismic velocity are proportional to age (except young orogens), implying greater crustal thickness should be observed in Archean crust compared to Proterozoic aged crust. An alternate view by Durrheim and Mooney [6] suggested that the velocity at the base of the Proterozoic crust is high (Vp>7.0) and is thicker (40-50 km), when compared to the Archean crust (27-40 km). Nelson [7] suggested that repeated episodes of the basaltic underplating cause crustal thickening. Rudnick and Fountain [8] found similar crustal thickness of Archean and Proterozoic. Arndt and Goldstein [9] had suggested that transportation processes between the lower continental crust and the upper mantle are potential contributors in building of the continental crust. It is still uncertain how crustal formation, growth, and composition have varied during the Precambrian; these open questions are addressed in this study.

The Moho or the Mohorovicic discontinuity defines the base of the crust and is the outermost seismic discontinuity in the Earth [10]. It is generally marked by an increase in velocity of P wave from <7.6 to >8 km/s [11]. Petrologically, it is called the crust-mantle transition zone, where the composition changes from mafic rock to ultramafic rock. Even though both terms are used simultaneously, studies suggest that they might vary [12]. According to Griffin and O’reilly [13], thick mafic roots may convert to the high-density mafic rock eclogite at base of crust. Eclogite has similar seismic velocities to the mantle, and hence, the crust-mantle boundary may not be observable as an abrupt increase in seismic velocity. Therefore, it is crucial to image the crustal properties and study the Moho nature. If, in a region, the lower crust is felsic and is underlain by a mafic upper mantle, then the nature of the Moho will be seismically sharp. On the other hand, if a mafic layer with high velocity is present at the base of the lower crust (as can be caused by intrusion from basaltic underplating) and is underlain by mafic upper mantle, then the Moho is seismically diffusive or gradational in nature.

Thus, we see that there is a link between evolution and seismic structure. Study of crustal structure helps in narrowing down on the competing hypothesis for crustal evolution. Repeated episodes of basaltic underplating have the possibility to manifest as a high velocity mafic layer at the base of the crust and could lead to crustal thickening. Study of crustal structure of different ages is also one of the best methods to test the different models of crustal evolution. The best areas to test this would be old crustal material that has not been deformed for a prolonged period of time and have survived the cycles of merging and continent rifting, the Precambrian cratons.

In this study, crustal properties have been imaged, comprising crustal thickness, crustal composition (Vp/Vs ratio), and nature of the Moho of the Precambrian cratons of Gondwanaland, using receiver function modeling. The technique of a single method (receiver function) used in this study for determination of crustal structure beneath the different areas of Precambrian cratons of Gondwanaland contributes to the subject of crustal evolution. It is worth investigating the crustal properties across the Precambrian of different ages and investigates the link between seismic structure and evolution. In order to obtain the detailed shear velocity structure along with Vp/Vs ratio variation with depth beneath the stations in different cratons, inversion of receiver function is needed. In this study, we obtain the shear velocity structure (Vs and Vp/Vs with depth) from shallow to the Moho depth using two different methods, the Hκ stacking and neighbourhood algorithm. Some of the previous studies in cratons just used the Hκ stacking approach, which is dependent on a priori P wave velocity (Vp) information, which may not be correct every time. Wrong selection of Vp may give erroneous H and κ values. Also, Hκ stacking provides crustal thickness and average crustal Vp/Vs ratio. Variation of Vs and Vp/Vs has important implications in terms of knowing any intracrustal discontinuities, or deep rocks in the crust, presence of fluids, etc.

Study area consists of Dharwar craton in India; Pilbara, Yilgarn, and Gawler cratons in Australia; Tanzania, Kaapvaal, and Zimbabwe cratons in Africa; Antananarivo craton in Madagascar; and Napier complex in Antarctica. The brief geology of all the cratons are discussed below.

2.1. Dharwar Craton

It is located in the southern part of India (Figure 2(a)). The crust of the Dharwar craton dates back to 3.4 Ga. The different lithological units are the Eastern Dharwar craton (EDC) and the Western Dharwar craton (WDC) that are separated by the Chitradurga schist belt (CSB) [14]. The Kaladagi basin of Proterozoic age borders the WDC in the north. Granulites and low-grade gneiss are present in the South of WDC, and greenstone is present in the North. To the west, the Western Ghats surrounds the WDC. There is a signature of thrusted contact between EDC and the Proterozoic Cuddapah Basin. On the eastern side, the Eastern Ghats surround the EDC. The Dharwar craton slowly transitions into the Neoarchean metamorphic terrane of the Southern Granulite Terrane in the south.

Figure 2

Maps of the different cratons showing the stations: (a) Dharwar, (b) Australia, (c) Tanzania, (d) Zimbabwe and Kaapvaal, (e) Antananarivo, and (f) Napier.

Figure 2

Maps of the different cratons showing the stations: (a) Dharwar, (b) Australia, (c) Tanzania, (d) Zimbabwe and Kaapvaal, (e) Antananarivo, and (f) Napier.

2.2. Pilbara Craton

The craton lies in Western Australia (Figure 2(b)). The Proterozoic Capricorn Orogen in the south separates it from the Yilgarn craton. Zircon crystals having an age of 3.8 Ga are exposed in the early Archean craton. Lithologically, it is divided into Western Pilbara terrane, Eastern Pilbara terrane, Kurrana terrane, and the Hamerseley Basin [15]. The Elizabeth Hill belt was created when the eastern and western part of the Pilbara craton collided. Granitic complexes are present in the southeastern part of the Kurrana terrane, and granite greenstone units are present in the western and eastern part. The Neoarchean Hamersley basin is present in the south and consists of volcanic sedimentary rocks. Tonalite-trondhjemite-granodiorite (TTG) rocks, greenstone belts, and dome and keel structures are also present in the Pilbara region that makes it unique. It has been suggested that vertical tectonics operated in the Pilbara, and due to partial convective overturn, the dome and keel structures were formed [16].

2.3. Yilgarn Craton

This craton (Figure 2(b)) and the Pilbara craton are separated by the Capricorn orogen, which is formed during Paleoproterozoic collision between the two cratons. The age of the oldest zircons dates back to 4.2 Ga [17] during the Hadean period. It consists of many terranes—Murchison province, Western gneiss terrane, Eastern Goldfield province, and the Southern cross province [1820].

2.4. Gawler Craton

The Gawler craton is situated in the southern part of Australia (Figure 2(b)). It comprises dominantly volcanic and sedimentary rocks. The age of the oldest granites in the craton dates back to 3.0 Ga. Tectonic activities in the Late Archean and early Proterozoic have led to the development of the craton. Felsic magmatism was dominant during the late Archean, and rifting was dominant initially during Proterozoic.

2.5. Tanzania Craton

Tanzania craton is located in the Eastern part of Africa (Figure 2(c)). The lithology consists of different units made of granites in the north and gneiss and granites in the South separated by the Dodoma Greenstone belt. The Mozambique belt runs along the eastern side of the craton and the Usagaran belt lies to the southeast. The Mozambique belt was created during Pan-African orogeny as a suture zone during the formation of Gondwanaland. The Tanzanian craton, which is formed due to the overriding plate of the subduction zone, created the Usagaran belt [21, 22]. Greenstone belts, TTG (Tonalite-trondhjemite granodiorite rocks), and gneisses are exposed on the craton [23]. Xenolith and geophysical studies indicate that a rigid and deep lithospheric keel is present beneath the craton [2426]. Presence of a superplume beneath the craton has led to the dynamic topography of the craton [27].

2.6. Kaapvaal Craton

It is located in the Southern part of Africa (Figure 2(d)). Similar rocks from this craton and the Pilbara suggest that they were part of the Vaalbara supercontinent [28]. Detrital zircons from the Gneiss complex in Swaziland date back to 3.7 Ga [29]. The lithology consists of granites, greenstone belts, and gneisses. The Bushfield igneous province is present in the North, along with the Limpopo belt that joined it with the Zimbabwe craton during late Archean metamorphism. Komatiites are present in the Barberton greenstone belt in the east.

2.7. Zimbabwe Craton

The Limpopo orogenic belt towards the Southern part of this craton separates it from the Kaapvaal craton (Figure 2(d)). It also consists of granites, greenstone belts, and gneisses. Lithologically, the Zimbabwe craton joined with the Kaapvaal craton in the Late Archean, and they are together called the Kalahari craton. The combination of sutures of Tokwe segment to the south and the Rhodesdale fragment in the north formed this cratonic block. Smith et al. [30] has found 3.7 Ga gneiss in the Tokwe segment. It has been suggested that the craton formed by arc accretion method [31]. In the northwest, it is surrounded by the Magondi belt and the Zimbabwe belt borders to the north of this craton.

2.8. Antananarivo Craton

It is located in central Madagascar (Figure 2(e)). The oldest rocks in the craton consist of migmatite gneiss that dates back to the Neoarchean age. The northern Sofia group and the southern Vondrozo group make up the supracrustal units of the craton. Intrusive origin migmatite gneisses are also present [32]. Three large synformal belts consisting of mafic schist and gneiss constitute the Tsaratanana complex in this craton. The rocks of the craton are of high metamorphic grade and complex in structure, and because it was formed in Neoarchean, it was subjected to magmatism, deformation, and metamorphism during the Proterozoic [33].

2.9. Napier Complex

The Napier complex lies in Enderby Land in Eastern Antarctica (Figure 2(f)). Before the breakup of Gondwanaland, the Dharwar craton of India shared a common boundary with this craton [34]. Various tectonothermal events, high-grade metamorphism, and plate tectonics led to Napier orogeny and formed the cratonic nucleus. Subsequent years of erosion and tectonic deformation have exposed the lower metamorphic core of the ancient mountains. The oldest tonalite and granitic gneiss record an age of 3.8 Ga [35]. Younger granulite facies rocks derived from younger granites and granodiorites and deformed gneisses from volcanic and sedimentary protoliths are also present [36]. In the late Archean, ultrahigh-temperature metamorphic conditions have led to their deformation. High-temperature metamorphism in the early Archean at 24-40 km crustal depth at 1050-1120°C was followed by sustained cooling periods at 20-30 km depth [37]. A subglacial depression separates the remaining Antarctica from this block.

3.1. Data Acquisition

The teleseismic events (epicentral distance in between 30° and 95°) for the different cratons were chosen such that all have a magnitude greater than or equal to 6. All the data was acquired from the Data Management Centre (DMC) of the Incorporated Research Institutions for Seismology (IRIS). A total of 52 seismic stations are used in our analysis for different cratons. In our H and Vp/Vs ratio calculation, we have used all freely available data from all the stations of Gondwanaland. If for a craton, two or more stations have similar H and κ estimates, we have taken representative station(s) to estimate the crustal parameters. The period of earthquake data in Dharwar was from 2006 to 2007 using instrument Streckeisen STS1 with sample rate of 20 samples/s. The data was from 2015 to 2017 with a sample rate of 40 samples/s in the Australian cratons using the instrument KS2000 in Pilbara, Guralp CMG3T in Yilgarn, and Guralp CMG3ESP in Gawler. In Tanzania craton, the earthquake data was extracted for the year 1994 with a rate of 20 samples/s that used instrument STS2 in stations MTOR, SING, and PUGE and with instrument cmg3esp in station MBWE. In the GETA station of Tanzania, the earthquake data was from 2011 to 2015 with a rate of 40 samples/s using the instrument trillium compact. In the DODT station of Tanzania, the data was from 2009 to 2015 with a rate of 40 samples/s using the instrument Streckeisen STS-2. In Kaapvaal craton, the period of earthquake data was from 1997 to 1999 using the instrument STS2 with a rate of 20 samples/s. In Zimbabwe craton, the period was from 1997 to 1998 in stations SA72, SA74, SA75, and SA77 and from 1997 to 1999 in station SA76, all using instrument STS2 with a sample rate of 20 samples/s. In Antananarivo craton, the period of earthquake data was from 2007 to 2010 in station ABPO with instrument Streckeisen with a rate of 20 samples/s and from 2012 to 2013 in stations BARY, BATG, BITY, and VINA using instrument Streckeisen with a rate of 40 samples/s. In Napier complex, the earthquake data was from 2013 to 2016 using the instrument Streckeisen with 40 samples/s. The coordinates of the stations in the cratons along with the networks are given in Table 1.

Table 1

Details of the stations operated in the study region with crustal thickness (H), average Vp/Vs ratio, average crustal shear wave velocity (avg. Vs) from Hκ stacking (1) and NA modeling (2).

StationNetworkLatitude (°)Longitude (°)Elevation (km)Hκ results (1)NA results (2)
H (1) (km)Vp/Vs (1)H (2) (km)Vp/Vs (2)Avg. Vs (2) (km/s)
Dharwar craton, India
HYBG17.41878.55251033.351.72932.71.7763.706
Pilbara craton, Australia
PSA00AU-21.5725119.84530529.151.71829.41.7073.86
PSAA1AU-21.576119.83632529.351.70229.41.7113.88
PSAA2AU-21.562119.84633029.351.71229.81.6913.8
PSAA3AU-21.578119.85431529.21.71430.31.6863.79
PSAB1AU-21.568119.82131729.71.69632.11.7293.88
PSAB2AU-21.554119.86130229.251.72228.21.7173.79
PSAC1AU-21.533119.80930330.051.66632.11.7123.87
PSAC3AU-21.621119.82735231.051.62229.01.663.71
PSAD2AU-21.617119.00236629.551.68832.11.7433.83
PSAD3AU-21.642119.77932729.41.71830.31.733.71
Yilgarn craton, Australia
BLDUAU-30.6147116.70928036.31.70237.11.763.8
MUNAU-31.9783116.20825334.61.71436.41.7543.774
NWAOAU-32.9277117.23926537.61.76641.61.7483.76
MEEKAU-26.638118.61553036.01.70636.21.7043.77
KMBLAU-31.3669121.88237036.851.73437.91.7443.8
Gawler craton, Australia
MULGAU30.2818134.058254.745.751.7346.21.7433.83
Tanzania craton, Africa
MBWEXD-4.9634.346110036.41.73736.21.7283.56
MTORXD-5.2535.4110038.71.70941.41.7513.76
PUGEXD-4.71433.184135036.41.73938.51.7473.64
SINGXD-4.6434.73146236.41.79938.11.7723.75
GETAAF-2.8633.24126735.91.71932.11.7573.44
DODTAF-6.18635.74111441.11.70440.21.753.68
Kaapvaal craton, Africa
SA13XA-25.8526.26160042.11.68937.71.7413.51
SA18XA-28.6334.305110236.21.70933.61.733.44
SA23XA-27.9323.4161240.11.74942.21.743.63
SA24XA-27.8824.236124438.31.69337.51.7513.66
SA25XA-27.8425.125125338.11.70435.21.7393.46
SA30XA-27.0724.165136136.51.68335.61.7293.65
SA32XA-26.8626.28138038.91.71535.21.7553.58
SA37XA-25.9723.72114233.91.71332.11.7423.47
SA38XA-25.9325.084121939.31.65936.11.7133.51
SA39XA-25.89526.151154541.51.70540.41.7313.66
SA40XA-25.89827.149151243.01.74542.41.7433.61
SA50XA-23.8727.16697638.41.74342.21.7413.75
SA55XA-22.97928.29891844.41.65142.21.7063.62
SA56XA-2329.07490942.641.69738.11.733.46
SA62XA-24.8525.135121440.651.76942.21.7563.72
SA139XA-25.8526.26160041.741.70539.11.7533.48
Zimbabwe craton
SA72XA-20.1428.61133738.01.71437.31.7473.72
SA74XA-21.9230.9348739.81.74539.61.7433.66
SA75XA-20.8628.99997139.11.69336.41.763.65
SA76XA-20.6329.84697836.21.7336.01.7113.6
SA77XAx-20.7530.9157640.11.67138.91.723.73
SA78XA-19.4630.772140136.51.71235.21.7713.55
Antananarivo craton, Madagascar
ABPOII-19.01847.229152843.51.62740.11.653.66
BARYXV-17.18446.857265.741.61.71940.11.7563.82
BATGXV-18.87846.1871552.240.01.8936.21.843.71
BITYXV-20.0647156641.91.73942.61.7523.75
VINAXV-18.17645.224293.836.11.73736.71.7493.74
Napier complex, Antarctica
MAWAU-67.662.8711246.841.63542.61.653.59
StationNetworkLatitude (°)Longitude (°)Elevation (km)Hκ results (1)NA results (2)
H (1) (km)Vp/Vs (1)H (2) (km)Vp/Vs (2)Avg. Vs (2) (km/s)
Dharwar craton, India
HYBG17.41878.55251033.351.72932.71.7763.706
Pilbara craton, Australia
PSA00AU-21.5725119.84530529.151.71829.41.7073.86
PSAA1AU-21.576119.83632529.351.70229.41.7113.88
PSAA2AU-21.562119.84633029.351.71229.81.6913.8
PSAA3AU-21.578119.85431529.21.71430.31.6863.79
PSAB1AU-21.568119.82131729.71.69632.11.7293.88
PSAB2AU-21.554119.86130229.251.72228.21.7173.79
PSAC1AU-21.533119.80930330.051.66632.11.7123.87
PSAC3AU-21.621119.82735231.051.62229.01.663.71
PSAD2AU-21.617119.00236629.551.68832.11.7433.83
PSAD3AU-21.642119.77932729.41.71830.31.733.71
Yilgarn craton, Australia
BLDUAU-30.6147116.70928036.31.70237.11.763.8
MUNAU-31.9783116.20825334.61.71436.41.7543.774
NWAOAU-32.9277117.23926537.61.76641.61.7483.76
MEEKAU-26.638118.61553036.01.70636.21.7043.77
KMBLAU-31.3669121.88237036.851.73437.91.7443.8
Gawler craton, Australia
MULGAU30.2818134.058254.745.751.7346.21.7433.83
Tanzania craton, Africa
MBWEXD-4.9634.346110036.41.73736.21.7283.56
MTORXD-5.2535.4110038.71.70941.41.7513.76
PUGEXD-4.71433.184135036.41.73938.51.7473.64
SINGXD-4.6434.73146236.41.79938.11.7723.75
GETAAF-2.8633.24126735.91.71932.11.7573.44
DODTAF-6.18635.74111441.11.70440.21.753.68
Kaapvaal craton, Africa
SA13XA-25.8526.26160042.11.68937.71.7413.51
SA18XA-28.6334.305110236.21.70933.61.733.44
SA23XA-27.9323.4161240.11.74942.21.743.63
SA24XA-27.8824.236124438.31.69337.51.7513.66
SA25XA-27.8425.125125338.11.70435.21.7393.46
SA30XA-27.0724.165136136.51.68335.61.7293.65
SA32XA-26.8626.28138038.91.71535.21.7553.58
SA37XA-25.9723.72114233.91.71332.11.7423.47
SA38XA-25.9325.084121939.31.65936.11.7133.51
SA39XA-25.89526.151154541.51.70540.41.7313.66
SA40XA-25.89827.149151243.01.74542.41.7433.61
SA50XA-23.8727.16697638.41.74342.21.7413.75
SA55XA-22.97928.29891844.41.65142.21.7063.62
SA56XA-2329.07490942.641.69738.11.733.46
SA62XA-24.8525.135121440.651.76942.21.7563.72
SA139XA-25.8526.26160041.741.70539.11.7533.48
Zimbabwe craton
SA72XA-20.1428.61133738.01.71437.31.7473.72
SA74XA-21.9230.9348739.81.74539.61.7433.66
SA75XA-20.8628.99997139.11.69336.41.763.65
SA76XA-20.6329.84697836.21.7336.01.7113.6
SA77XAx-20.7530.9157640.11.67138.91.723.73
SA78XA-19.4630.772140136.51.71235.21.7713.55
Antananarivo craton, Madagascar
ABPOII-19.01847.229152843.51.62740.11.653.66
BARYXV-17.18446.857265.741.61.71940.11.7563.82
BATGXV-18.87846.1871552.240.01.8936.21.843.71
BITYXV-20.0647156641.91.73942.61.7523.75
VINAXV-18.17645.224293.836.11.73736.71.7493.74
Napier complex, Antarctica
MAWAU-67.662.8711246.841.63542.61.653.59

3.2. Receiver Function Calculation

The internal structure of the Earth near the receiver was estimated by modeling the time domain receiver function, which is calculated by deconvolving the horizontal component with the vertical component [3840]. When plane waves interact with a discontinuity (say, the Moho), a part of P wave gets converted to an S wave (denoted as Ps), and other parts undergo multiple reflections (multiples) in between the Earth’s surface and that discontinuity (denoted by PpPs and PpSs). To compute the receiver function, ,first, the TauP-Toolkit developed by Crotwell et al. [41] has been used to mark the P arrivals in the waveforms with reference to the IASP91 velocity model. The mean and trend are removed from the waveforms, and 150 s of the waveforms (from 30 s before P arrival time to 120 s after P arrival time) is selected for further processing. The radial and transverse components are then obtained by rotating the north-south and east-west components with respect to a great circular arc. The method of iterative time domain deconvolution [42] has been used to obtain radial and transverse receiver function by deconvolving the radial and transverse components, respectively, with the vertical component. Lowpass Gaussian filter of width (Gw) 2.0 (corner frequency ~1 Hz) was applied to the receiver functions. The receiver functions with more than 80% fitting are used in this study. The Ps arrival times are increased, and their multiples are decreased as ray parameters increase. So, move-out correction of the receiver functions has been done for individual stations using the IASP91 velocity model and reference epicentral distance of 67°. Each station consists of many good quality receiver functions; the minimum number is kept to 5.

Figure 3 shows the receiver functions of one representative station from each of the Archean cratons. Variations of Ps, PpPs, and PpSs are also shown in the figure with different colors, which can provide initial information of crustal thickness beneath each station. The Ps arrival time for the station HYB of Eastern Dharwar craton, India (Figure 3(a)), and PSA00 of Pilbara, Australia (Figure 3(b)), is observed at ~3.6 s; however, their multiples (PpPs and PpSs) have different arrivals. For the other two cratonic blocks of Australia, Yilgarn (PSA00) and Gawler craton (MULG), the Ps arrival times are observed at ~4.5 s and ~5.5 s, respectively. This suggests a thicker crust beneath the Gawler craton, compared to Pilbara and Yilgarn. Pilbara may have the thinnest crust in Australia. Based on the Ps arrivals of receiver functions estimated for DODT of Tanzania (Figure 3(e)), SA38 of Kaapvaal, and SA74 of Zimbabwe, which is ~4.5 s, with PpPs arrival ~15-16 s and PpSs arrival ~19-21 s; we can say crustal thickness may be similar for these cratons. In Antananarivo craton (ABPO) and Napier complex (MAW), the Ps, PpPs, and PpSs arrival times are observed to be ~4.5 s, ~17 s, and~22-23 s, respectively, suggesting similar crustal thickness in both the cratons.

Figure 3

Receiver functions plot with increasing back azimuths for one representative station from each of the Archean cratons, namely, (a) the Dharwar, (b) the Pilbara, (c) the Yilgarn, (d) the Gawler, (e) the Tanzania, (f) the Kaapvaal, (g) the Zimbabwe, (h) the Antananarivo, and (i) the Napier. The red, blue, and green lines are probable markers of the Ps, PpPs, and PpSs timings.

Figure 3

Receiver functions plot with increasing back azimuths for one representative station from each of the Archean cratons, namely, (a) the Dharwar, (b) the Pilbara, (c) the Yilgarn, (d) the Gawler, (e) the Tanzania, (f) the Kaapvaal, (g) the Zimbabwe, (h) the Antananarivo, and (i) the Napier. The red, blue, and green lines are probable markers of the Ps, PpPs, and PpSs timings.

3.3. Calculation of Thickness of Crust (H) and Vp/Vs Ratio (κ)

The arrival time of the phases Ps, PpPs, and PpSs provides only the preliminary information of the crustal thickness. The average Moho depth and ratio of Vp/Vs in the stations in different cratons is calculated by the Hκ method [43] considering the arrival times of the receiver functions recorded in different stations. The original receiver functions have been used for the method. In this method, we optimize the objective function SH,κ that is illustrated in Equation (1). The variation of H is in between 20 and 60 km (increment of 0.05 km) and κ between 1.6 and 2.0 (with 0.002 increment), and each time, S is calculated for the different pairs. The corresponding H and κ values having the highest S value give the actual value of best fit and average Vp/Vs ratio.
where rjt is the receiver function amplitude corresponding to the jth event, N is the total number of receiver functions, and t1, t2, and t3 are anticipated Ps, PpPs, and PpSs+PsPs arrival times representing different H and κ values. The weighting factors w1, w2, and w3 are defined such that wi=1 and contributions from all 3 phases in receiver function are weighted to assign wi. Out of all the phases, the amplitude of Ps is greater than PpPs and PpSs, and accordingly, it has been assigned higher weight than the other two. Various combinations of wis have been used (w1=0.6, w2=0.3, and w3=0.1 or w1=0.7, w2=0.3, and w3=0) depending on the amplitude of the converted and multiple phases. Also, tis are calculated using the following equation.

The average P wave velocity of the crust is considered 6.4-6.5 km/s for the calculation of t1, t2, and t3. Vp was varied between 6.1 and 6.6 km/s to calculate the uncertainties in H and κ for the error analysis, along with changing the combination of weighting factors. The uncertainties varied from 3.0 km for H and 0.5 for Vp/Vs. The Hκ result of station PSAD3 in Pilbara craton is shown in Figure 4.

Figure 4

The Hκ result for station PSAD3 in Pilbara craton, Australia. (a) Receiver function plotted with increasing back azimuth values and (b) the Hκ result for the station PSAD3 along with their respective values shown in the top right corner. The black circle is the actual value of the crustal thickness (H) and average Vp/Vs ratio (κ) for the optimal value of SH,κ.

Figure 4

The Hκ result for station PSAD3 in Pilbara craton, Australia. (a) Receiver function plotted with increasing back azimuth values and (b) the Hκ result for the station PSAD3 along with their respective values shown in the top right corner. The black circle is the actual value of the crustal thickness (H) and average Vp/Vs ratio (κ) for the optimal value of SH,κ.

3.4. Inversion Modeling Using Neighbourhood Algorithm (NA)

The Hκ method calculates the crustal properties by assuming an average Vp, which may not be the same for all the stations. Hence, only calculation of Moho depth and an average Vp/Vs ratio is possible. To find out Vs and Vp/Vs ratio variation with depths, inversion modeling is used. Some other previous studies used tomography or gravity inversion for large-scale regions; hence, the resolution for the crust could be low. Inversion of receiver function using the neighbourhood algorithm method is a nonlinear method free from a priori information, which makes the estimations more robust in the study. To reduce the effect of nonuniqueness and nonlinearity, which are common problems in receiver function inversion (Ammon et al., 1990), Sambridge [44] developed a fully nonlinear derivative free direct search approach, the neighbourhood algorithm (NA) for finding models of acceptable data fit. Through the modeling, we can get an ensemble of models that preferentially sample the good data-fitting regions of parameter space, rather than seeking a single optimal model. Linearized inversion is strongly dependent on the initial model. However, NA is less sensitive to the initial model, and in most cases, it converges to the global minimum. We also observed that for different initial models, NA produces a similar final model by only changing either the number of iterations and/or number of samples/resamples.

The final model is selected from the set of 1000 best models, based on the minimum misfit between the synthetic and observed receiver functions. Misfit is the measure of the discrepancy between the true and predicted RF from an Earth model. Minimum misfit does not guarantee the best model, and in many cases, it provides an unrealistic Earth model. We have given preference to the model providing minimum misfit (may not be the least) with a realistic Earth model. Detailed velocity modeling at depth is needed in order to explore the crustal nature. In this method, the lithosphere consists of 6 horizontal layers, namely, sediment, basement, upper, middle, and lower crusts, and upper mantle, with 4 model parameters, the Vp/Vs ratio, the thickness of layer, and the shear wave velocities at the uppermost and lowermost of the layer. This makes it a 24- (6×4) dimensional parameter space. In this method, synthetic receiver functions are constructed for ns models that were created for every iteration. An objective function is defined based on the L2 norm of the difference between observed and synthetic receiver function, which is considered the misfit. Then nr models were selected based on lower misfit values. Then again, ns new models are created based on random walk through the Voronoi cell of nr models. In this way, new models are produced in areas where the fitting is good. The fitting and inversion result of station PSAD3 in Pilbara craton, Australia, is shown in Figures 5(a) and 5(b), respectively. 40,000 velocity models that have been generated by iterating the neighbourhood algorithm 400 times are shown in Figure 5(b), and the original stacked receiver function is compared with the synthetic receiver functions corresponding to each model. Best fitted (colored section) 1000 models are selected based on the minimum misfit between the synthetic and observed receiver functions. The red line represents the best fitted model. The observed stacked receiver function (black line) along with the best fitted synthetic receiver function (red line) bound by ±1σ error bounds (black dashed line) is shown in Figure 6 from one station in each craton. The Moho could be observed as a sharp step jump in the shear wave velocity inversion which is easy to identify, or as a gradational change in velocity. In the latter case, it is picked up where the velocity crosses 4.2-4.3 km/s (Figure 7). The same process is followed for all other stations (included in supplementary available (here)). The uncertainties were calculated from the best fitting models. The uncertainties for κ and H are shown in Figure 8.

Figure 5

Inversion result for station PDAD3 in Pilbara craton, Australia. (a) Comparison of the original receiver function and synthetic receiver function obtained from NA in black and blue lines, respectively. (b) Variation of the Vp/Vs ratio and shear wave velocity with depth. The gray region shows ~40,000 models generated during inversion of which the best 1000 models are shown in yellow and green regions that have the least difference between observed and synthetic receiver function. The white and red lines in the model density plot represent the shear velocity average model and best fitting model, respectively. The red line on the left indicates the variation of Vp/Vs with depth. Red arrows mark the Moho depth.

Figure 5

Inversion result for station PDAD3 in Pilbara craton, Australia. (a) Comparison of the original receiver function and synthetic receiver function obtained from NA in black and blue lines, respectively. (b) Variation of the Vp/Vs ratio and shear wave velocity with depth. The gray region shows ~40,000 models generated during inversion of which the best 1000 models are shown in yellow and green regions that have the least difference between observed and synthetic receiver function. The white and red lines in the model density plot represent the shear velocity average model and best fitting model, respectively. The red line on the left indicates the variation of Vp/Vs with depth. Red arrows mark the Moho depth.

Figure 6

Comparison of synthetic receiver function (red line) with observed stacked receiver function (black line) with ±1σ error bounds (dashed line) from one station in each craton.

Figure 6

Comparison of synthetic receiver function (red line) with observed stacked receiver function (black line) with ±1σ error bounds (dashed line) from one station in each craton.

Figure 7

Inversion results of one station from each craton. The craton name is on the top left corner of each plot. The station name is shown on the left top in each figure. Other details are the same as in Figure 5.

Figure 7

Inversion results of one station from each craton. The craton name is on the top left corner of each plot. The station name is shown on the left top in each figure. Other details are the same as in Figure 5.

Figure 8

The estimated Vp/Vs ratio versus Moho depth along with the error bounds for the different cratons in this study.

Figure 8

The estimated Vp/Vs ratio versus Moho depth along with the error bounds for the different cratons in this study.

Table 1 incorporates the Hκ and NA results from the stations used in this study. The Hκ results, fit of observed and synthetic receiver functions, and NA results from one station in each craton are shown in Figures 6, 7, and 9, respectively, and the remaining results from the other stations have been illustrated in supplementary and are documented in Table 1. The crustal thickness and Vp/Vs ratio estimated using both the methods (Hκ and NA) are very similar, with few exceptional cases. Slight variations of those values could be a result of selection of a priori value of average Vp used in the Hκ method. Since, no a priori information is used in NA, the results estimated by NA are used for the results and discussion. Use of both Hκ and NA makes the estimation more robust in this study. Detailed discussion is focused mainly on those cratons having sufficient stations (Pilbara, Yilgarn, Gawler, Tanzania, Kaapvaal, Zimbabwe, and Antananarivo cratons). But previous studies in these cratons have been included for comparison in crustal properties. Baranov and Morelli [45] suggested that Napier crustal thickness to be uniform of around 40 km. Here, we obtain a higher limit of ~43 km using NA (Figure 7) and ~47 km using the Hκ method (Figure 9). The Moho depth of the station in Gawler craton (~46 km) is slightly higher than the highest Moho depth (~44 km) indicated by Salmon et al. [46]. They have also been considered in this study for comparative study in crustal properties. We will at first discuss the crustal thickness and Vp/Vs ratio results of the six cratons (Pilbara, Yilgarn, Gawler, Tanzania, Kaapvaal, Zimbabwe, and Antananarivo cratons) obtained in this study and comparing our results with previous models of crustal studies. This is followed by discussion of crustal nature and rocks in the areas using Vp/Vs and Vs values. Crustal parameters have also been used to infer the degree of crustal reworking and relative density. Our results also support the link between seismic structure, evolution, and Gondwana connection.

Figure 9

Crustal thickness (H) and Vp/Vs (κ) results from Hκ stacking of one station from each craton. The craton name is on the top left corner of each plot. The top right corner shows the station name along with H and Vp/Vs values. The colored part indicates higher S values, with red being the highest. Other details are the same as in Figure 4.

Figure 9

Crustal thickness (H) and Vp/Vs (κ) results from Hκ stacking of one station from each craton. The craton name is on the top left corner of each plot. The top right corner shows the station name along with H and Vp/Vs values. The colored part indicates higher S values, with red being the highest. Other details are the same as in Figure 4.

Results estimated in our study are in agreement with the previous studies done on those areas (for example in Australia by Kennett et al. [47] and Salmon et al. [46], in Africa by Youssof et al. [48], and in Madagascar by Andriampenomanana et al. [49]). Results obtained from our study show significant variation in crustal thicknesses (~28-43 km) within and across the different Archean cratons ~28-32 km in Pilbara, ~36-42 km in Yilgarn, 36-41 km in Tanzania, 32-42 km in Kaapvaal, 35-40 km in Zimbabwe, and 36-43 km in Antananarivo (Table 1; Figure 8; Figure 10). The crustal thickness in single stations in EDC, Gawler, and Napier are ~33 km, ~46 km, and~43 km, respectively (Figure 7). The Eastern Dharwar Craton was studied extensively by Borah et al. [50] and Sarkar et al. [51]. In Napier and Gawler, along with our results, we refer to the results of Baranov and Morelli [45] and Salmon et al. [46]. These variations of crustal thicknesses are very much similar to values calculated by Durrheim and Mooney [6] in different Archean cratons (27-40 km), with a higher limit obtained in our study. The average crustal Vs varies between ~3.4 and 3.8 km/s in the stations. The average crust Vp/Vs ratio of all the Archean cratons varies between 1.65 and 1.78, with the exception of station BATG that has a Vp/Vs value of 1.84, which indicates that the crust beneath the Archean cratons is mostly of felsic to intermediate in nature. Crustal delamination could possibly have resulted in the stabilization of the Archean cratons, leading to dominant felsic to intermediate composition of the Archean cratons [52].

Figure 10

Variation of Moho depth with respect to age of craton. The dashed lines in black are the boundaries of the Paleoarchean, Mesoarchean, Neoarchean, and the Hadean. The colored circles represent the estimated Moho depth and age of the stations in the different cratons in this study. The colored horizontal lines are from the previous studies of Moho depths. The red dashed lines represent the average Moho depth in each craton considering this study and previous studies. The faint blue vertical lines on the left and right show the bounds for the average range of crustal thickness.

Figure 10

Variation of Moho depth with respect to age of craton. The dashed lines in black are the boundaries of the Paleoarchean, Mesoarchean, Neoarchean, and the Hadean. The colored circles represent the estimated Moho depth and age of the stations in the different cratons in this study. The colored horizontal lines are from the previous studies of Moho depths. The red dashed lines represent the average Moho depth in each craton considering this study and previous studies. The faint blue vertical lines on the left and right show the bounds for the average range of crustal thickness.

4.1. Pilbara Craton

The crustal thickness beneath the stations varies from 28 to 32 km with an average towards the lower limit. This suggests that the crust beneath Pilbara is thin. Our results show good correlation with Moho depth models (28-33 km) estimated in this region [46, 47, 53]. Our observations show that the Vp/Vs ratio varies from 1.66 to 1.74, with an average of ~1.71. This estimation of crustal properties beneath the Pilbara in this study is consistent with the average crustal thickness of 30 km and a Vp/Vs ratio of 1.70 in earlier models [54]. The average crustal shear wave velocity varies from 3.71 to 3.85 km/s. Gradational Moho is also observed beneath a few stations (Figure 7 and Supplementary) Mafic underplating was reported in previous studies ([55]; Drummond and Collins [56], which could be a possible reason for the existence of transitional Moho beneath few stations.

4.2. Yilgarn Craton

Previous studies have estimated the Yilgarn crust to be relatively thicker (>35 km) than Pilbara (Yuan, 2015). The crustal thickness observed in this study beneath the stations in this region is much higher than the Pilbara and varies from ~36 to 42 km. The Vp/Vs ratio varies from ~1.70 to 1.76, with an average of ~1.74. This is consistent with the previous Vp/Vs average model (~1.73 beneath Yilgarn) revealed by Yuan [57]. The average crustal Vs varies from 3.76 to 3.8 km/s.

4.3. Tanzania Craton

The crustal thickness beneath the stations in our study varies from 36 to 41 km in Tanzania. This is comparable with the average Moho depth of 39 km in the Archean terranes of Africa [58] using the gravity field model. Our results also show a good correlation with the crustal thickness values (34-40 km) estimated beneath the Tanzania craton in the Kenya Rift International seismic project [59, 60]. The crustal average Vp/Vs ratio in our study varies from 1.73 to 1.77. This is consistent with the Vp/Vs values (1.75-1.77) reported by Hodgson et al. [61] with a slightly lower limit. The crustal shear wave velocity beneath the stations in our study varies from 3.56 to 3.78 km/s. Gradational Moho is also observed beneath few stations (Figure 7 and Supplementary 3), indicating the presence of mafic underplating. Magmatic events such as intrusion of mafic dyke swarms were revealed in studies in the Tanzania craton [62], thus confirming the reason for the gradational Moho obtained in few stations.

4.4. Kalahari Craton (Kaapvaal and Zimbabwe)

The Kalahari cratons show approximately similar variations in crustal thickness (32-42 km in Kaapvaal and 35-40 km in Zimbabwe). The Vp/Vs ratio varies from 1.73 to 1.75 in Kaapvaal and 1.72 to 1.77 in Zimbabwe. These values are in good correlation with crustal studies done in Africa [48]. Moho is observed clearly beneath the stations in good correlation with studies by Youssof et al. [48], except for station SA55.

4.5. Antananarivo Craton

The crustal thickness beneath the stations in our study varies from 36 to 42 km, with a Vp/Vs ratio variation of 1.65-1.84, with an average of 1.75. This is consistent with the average Vp/Vs ratio of 1.75 in Antananarivo suggested by Rindraharisaona et al. [63]. The presence of mafic lower crust in this craton [63] leads to higher value of Vp/Vs (>1.8) beneath some stations. We also observe the Moho transition to be gradational in a few stations (Figure 7 and Supplementary 3). The gradational nature of Moho obtained in this study in some stations confirms the mafic underplating beneath this craton. Similar range of Moho depth (36-43 km) was also reported by Andriampenomanana et al. [49]. The crustal shear wave velocity beneath the stations in this region varies from 3.66 to 3.84 km/s. An average crustal velocity of 3.7 km/s was suggested by Andriampenomanana et al. [49]. This is consistent with our obtained value of average crustal velocity of 3.74 km/s. Pacquette et al. [33] suggested that the Neoarchean rocks in this craton were subjected to magmatism, deformation, and metamorphism that resulted in complex structured higher metamorphic grade rocks. Mafic intrusive rocks subjected to extreme conditions of metamorphosis had resulted in the formation of mafic gneisses in Madagascar cratons (Boger et al., [64]; [65, 66]), thus having a high Vp/Vs ratio. This geological study is coherent with the geophysical Vp/Vs ratios obtained in this study beneath the stations.

4.6. Extrapolation of Vp/Vs to Infer the Possible Rock Types

Laboratory experiments conducted by Christensen [67] indicated that rock silica percentage affects the Vp/Vs ratios. The Vp/Vs ratio for felsic crust having SiO2%>66% is less than 1.75 (e.g., Vp/Vs of felsic rock granite ~1.71); the ratio is 1.75-1.80 for intermediate crust having SiO2%~52-66% (e.g., Vp/Vs of intermediate rock diorite ~1.78), and more than 1.80 for mafic crust having SiO2%~45-52% (e.g., Vp/Vs of mafic rock gabbro is ~1.87). The presence of fluids can also affect the Vp/Vs ratios, increasing it even further. This typically occurs in subducted oceanic crust [68, 69]. The Vp/Vs ratio results estimated in the different cratons in this study range from 1.65 to 1.78 indicating that the nature of the crustal composition on average is felsic to intermediate. In order to check the possible rock types in the Precambrian cratons, we have averaged the laboratory data of the Vp/Vs ratio and Vs for the crustal layers (upper crust at pressure of 200 MPa, middle crust at pressure of 600 MPa, and lower crust at pressure of 800 MPa), with reference to Christensen [67], illustrated in Figure 11. We have then superimposed on Figure 11 the average crustal Vp/Vs ratio and Vs values for the stations denoted in colored circles that we have estimated using the NA modeling. The rocks inferred from the Vp/Vs ratios and Vs values in our study were found to be similar to the geological studies of rocks in those areas. Granitic complexes are present in the southeastern part of the Kurrana terrane, and granite greenstone units are present in the western and eastern part of the Pilbara region. The lithology of the Kalahari cratons consists of granites, greenstone belts, and gneisses, consistent with the Vp/Vs ratios (~1.72). The high obtained Vs values in Madagascar craton (as high as 3.8 km/s) are similar to the Vs of mafic granulites (~3.8 km/s). Large belts of mafic schist and gneiss constitute the Tsaratanana complex in Antananarivo, which might result in high Vs values. The high-grade structurally complex metamorphic rocks of the Neoarchean age were subjected to magmatism, deformation, and metamorphism during the Proterozoic [33]. Beneath Kaapvaal craton, biotite (tonalite) gneiss has been observed [70], leading to the observed Vp/Vs and Vs values as in Figure 11. Braun and Kriegsman [71] suggested the presence of felsic paragranulites in the Dharwar. Our results indicate that in Dharwar, the Vp/Vs and Vs values are close to paragranulite, felsic granulite, and diorite (Figure 11).

Figure 11

Comparison of the average crust Vs and Vp/Vs ratio estimated for different stations in the study with the different Vs and Vp/Vs ratio of the rocks calculated from Christensen [67]. The boundaries between felsic to intermediate and intermediate to mafic are marked by black dashed line.

Figure 11

Comparison of the average crust Vs and Vp/Vs ratio estimated for different stations in the study with the different Vs and Vp/Vs ratio of the rocks calculated from Christensen [67]. The boundaries between felsic to intermediate and intermediate to mafic are marked by black dashed line.

4.7. Crustal Parameters to Infer about the Degree of Crustal Reworking and Relative Density

Although some undeformed ancient crust can preserve the original crustal properties, denudations alter the surface topography and the geological history of the craton. Tectonics can change the crustal properties of the ancient cratonic crust during the long evolutionary processes, such as magmatism, underplating, or delamination due to rifting or compression. Chevrot and van der Hilst [72] suggested that negative regression line between Vp/Vs and Moho has the possibility to provide evidence for past tectonic activities and crustal growth such as orogenic processes or erosion. If a crustal block undergoes collision, orogeny pushes the upper crust upwards and thickens it, hence increasing the felsic component. This leads to lowering of the Vp/Vs ratio. The opposite occurs in erosional process, i.e., the upper felsic part is removed, resulting in thinner mafic crust (high Vp/Vs ratio). So, two opposite types of tectonic processes can result in negative correlation of Vp/Vs ratio with Moho. After comparing Vp/Vs with Moho (Figure 8), it was found that no concrete relation in trend of regression was observed between the Vp/Vs and Moho in our study. We think that crustal reworking has a great contribution in differences in expected and observed Vp/Vs values when compared to the Moho. Also, to compile the relation between H and Vp/Vs, we have also prepared Figure 12 that shows the crustal thickness (H) versus Vp/Vs from previous literature for the cratons in this study superimposed with our estimations to check the correlation. We see that from having many data points of H and Vp/Vs from our study region, there exists no correlation between the crustal thickness and Vp/Vs ratio. Relation between crustal thickness and Vp/Vs is not observed globally and depends on the tectonics of the region. Cratonic regions are highly stable regions. The relations may hold for active orogenic regions [72].

Figure 12

Moho vs. Vp/Vs plots from the Precambrian cratons from previous literature (colored circles) along with our estimations (plus sign).

Figure 12

Moho vs. Vp/Vs plots from the Precambrian cratons from previous literature (colored circles) along with our estimations (plus sign).

Additionally, the relative crustal density can be inferred by plotting elevation versus crustal thickness, and the Airy isostatic equilibrium line assuming 3300 kg/m3 mantle density and 2760 kg/m3 average crustal density as illustrated in Figure 13. The points lying above the red line are relatively denser than isostatic equilibrium. The green triangles show the dynamic topography of the Tanzania craton. It is possibly due to the presence of a superplume beneath the craton that has led to the dynamic topography of the craton [27]. The relative crustal density is also high in most of the Tanzanian stations. For Australian cratons, the relative density of Pilbara craton is high (Figure 13), compared to Yilgarn. This also supports the vertical tectonics that operated in the region [16].

Figure 13

Elevation versus crustal thickness for the different stations in the study. The red line is the Airy isostatic equilibrium line, assuming average crustal density of 2760 kg/m3 and a mantle density of 3300 kg/m3.

Figure 13

Elevation versus crustal thickness for the different stations in the study. The red line is the Airy isostatic equilibrium line, assuming average crustal density of 2760 kg/m3 and a mantle density of 3300 kg/m3.

4.8. Comparison of Crustal Properties with Age and Episodic Evolution of Crust

The Hadean and Archean periods were important periods during Earth’s history, and the crustal properties of the cratons belonging to the rocks of the Precambrian are important as they can give valuable clues in understanding continental structure and evolution. Current global crustal studies suggest the average crustal thickness of continents in the range of ~35-40 km. As the Precambrian cratons have survived the cycles of continental merging and rifting, it is worth investigating to see if the continental crusts of the cratonic regions follow similar range and also whether there are Moho variations with respect to the age of the Precambrian cratons as illustrated in Figure 10. We have also compiled previous estimates and included them in the plot. The rocks of the Yilgarn date back to the oldest period, Hadean. The Moho estimates are almost similar to the present-day continental estimates. The Archean eon can be divided into 4 geologic eras—the Eoarchean (4-3.6 Ga), the Paleoarchean (3.6-3.2 Ga), the Mesoarchean (3.2-2.8 Ga), and the Neoarchean (2.8-2.4 Ga). The Napier and Pilbara cratons belong to the Eoarchean, and the Moho depths fall outside the global continental average estimates. It is interesting to note the difference in Moho depth from two of the oldest Archean cratons: very shallow in Pilbara stations (~28 km-31 km) and very deep in the station of Napier (~43 km), even though both belong to the same age of 3.8 Ga. The difference in the crustal formation process in the Pilbara and the Antarctic craton could possibly have an effect on the crustal nature. It has been suggested that due to vertical tectonics operated in the Pilbara and due to partial convective overturn, the dome and keel structures were formed [16]. The Napier complex is a unique ultrahigh-temperature metamorphic complex that has ever been recorded with the highest grade of metamorphism. Various tectonothermal events, high-grade metamorphism, and plate tectonics led to Napier orogeny and formed the thick cratonic nucleus. The Paleoarchean includes the Kaapvaal and Zimbabwe cratons. Tanzania and Gawler belong to the Mesoarchean, and EDC and Antananarivo belong to the Neoarchean. Careful observation of Figure 10 shows that overall (excluding Pilbara and Gawler), there is not much variation in range of average crustal thickness in this study (~33-45 km) with respect to age of craton during the Precambrian. This range (~33-45 km) is slightly wider than the current global average ~35-40 km Moho depth. Previous studies indicate that the conditions of crustal growth have varied during the Archean and plate tectonics were dominant in the middle to late Archean. Plate tectonics control the crustal thickness to a large extent [73]. In the early Earth, the temperature was very high, and the plates were less dense and buoyant, so plate tectonics was probably not that effective [74]. Heat flow in the earth controls the rate of cooling of rocks. It takes some time for the lithosphere to cool down and the plates to become denser [75]. Neoarchean rocks were subjected to magmatism, deformation, and metamorphism during the Proterozoic [33], and underplating in the lower crust may lead to accumulation of more crustal material. But crustal melting processes like delamination can lead to melting and destruction of the lower crust. Hence, an increase in crustal thickness could not be preserved during the later part of Archean. Durrheim and Mooney [6] suggested that the mafic layer in the lower crust started appearing in the Proterozoic, with thicker crust appearing in the Proterozoic. The observations of few gradational Moho in our study beneath few stations indicate that the Archean crusts were subjected to the possibility of mafic underplating at some point, or thick crusts were present even during the Archean. In fact, Baranov and Morelli [45] reported that the central part of Antarctic Grunehogna craton shows an exceptionally thicker crustal root, which could be added later or formed initially.

Other than delamination, the role of episodic cycles of crustal evolution is also observed in the pattern of crustal thickness across each division of the Precambrian to some extent. Figure 10 shows that the range in the variation of observed present crustal thicknesses follows an interesting relatively high and low pattern alternatively with respect to age. The Hadean crust shows ~7 km variation in crustal thickness, followed by ~16-17 km in the Eoarchean and Paleoarchean. The Mesoarchean shows ~11 km Moho depth variation, increasing to ~15 km in the Neoarchean. Such a situation could take place due to alternative high and low crustal regeneration process, repeated episodic recycling, and reworking of early crust. Repeated episodes of basaltic underplating as indicated by Nelson [7] could also account for such differences in evolution of crust through time. Geodynamic models indicate that episodic cycles of large high and low melt production and crustal growth were prevalent in the early Earth [3]. Not only of the formation process, there is also an equal possibility of the other situation of later processes affecting the crustal thicknesses that we see today. The early crust after formation can go through many later processes that impact its original thickness. But the strong cratonic roots and strong lithosphere have allowed them to survive the later processes.

4.9. Link between Evolution, Seismic Structure, and Gondwana

The cratons have survived the cycles of merging and rifting of continents. Plate tectonics have controlled the conditions of crustal formation, growth, and evolution. During Gondwana amalgamation, the collision of Archean and Proterozoic continental fragments led to Pan-African orogeny that comprises the bedrock geology of the African cratons, and this is similar to the Vp/Vs and Vs of the rocks in those areas. The Mozambique belt was formed as a suture between plates during the orogeny. Stern [76] suggested that it resulted from subduction and collision of east and west Gondwana. Different types of collision and also collision rate could have been an important parameter in crustal thickness variation of the cratonic areas. Crustal modification alters the composition. Geological studies indicate the presence of dominantly felsic to intermediate rocks in the cratonic areas, which is consistent with the Vp/Vs ratios obtained in our study. The breakup of the supercontinent Gondwanaland led to the separation of the continents and their gradual movement in different directions. This may lead to crustal evolution and change in crustal properties. Difference in spreading and migration rate of the continents also leads to crustal delamination, and destruction of crustal roots. The present-day crustal thickness and composition were affected by plate tectonics to a large extent. Some parts of East Gondwanaland, Australia, and India moved at a faster rate than other continents. Rifting histories of Madagascar, Antarctica, and Africa also played an important role in their crustal evolution.

The Antananarivo craton in Madagascar and the African cratons—Zimbabwe and Tanzania—show a similar range of crustal thickness. Sarkar et al. [51] had studied the Dharwar craton in detail and revealed the average Moho depths of ~37 km in Dharwar craton (~34 km in Eastern Dharwar and~41 km in Western Dharwar). Sri Lanka also had drifted along with India, and an average Moho depth of ~37 km was recorded in the study by Mukherjee et al. [2]. Although the current position of Madagascar is close to Africa, Madagascar and India initially were part of a single tectonic plate during Gondwana configuration and had started rifting post Gondwana breakup. Then India separated from Madagascar and started drifting at a high speed towards Asia in the current position. The similar range of crustal thicknesses in Antananarivo, Tanzania, Zimbabwe, Dharwar, and Sri Lanka supports their Gondwana connection. The rifting history of Africa also possibly led to its crustal modification. At first, it separated along with South America from East Gondwanaland, and then later, it separated from South America giving birth to the Atlantic Ocean. The African cratons also, especially Kaapvaal, show a huge variation of thickness. The eastern part of Africa served as the connecting region between West and East Gondwana, and the cratons in this region went through crustal modifications to a large extent, resulting in varied Moho depths. Antarctica was also closer to India initially, and Southern Australia was attached to east Antarctica. Then Antarctica moved Southwards, and Australia moved eastwards at a high rate. Geochemical evidence indicates that Gawler craton was part of the Antarctic Shield prior to breakup of Gondwana [77], and the closer Moho depth in this study (~46 km in Gawler and~43 km in Napier) also supports this conclusion. Zegers et al. [28] had observed large shear zones in the Pilbara region. The Pilbara was exposed to vertical convective overturns, and the extensive extensional environment also led to crustal modifications [78]. All of this suggests a reduction in its crustal thickness, in agreement to our estimations. The crusts of the Precambrian cratons evolved to a great extent after the separation from Gondwana configuration. The present crustal properties seemed to have been affected by plate tectonics, spreading velocities and large distances covered by the continents. Australia and India covered a large distance after splitting up completely from Gondwana.

This work analyses crustal properties across the ancient Precambrian cratonic regions that previously made up the Gondwana supercontinent. Using the teleseismic analysis method of the receiver function analysis, we have analysed and compiled the crustal thickness, average Vp/Vs, ratios and crustal Vs. These are then compared across cratonic locations, and we attempt to infer details of ancient crustal formation processes and subsequent crustal remodeling.

  • (i)

    On comparing the crustal thickness with the age of the cratons, it was found that average bounds of crustal thickness were similar between Hadean, early Archean, and late Archean cratons (~33-45 km). The cratons have a dynamic history of migration and separation after Gondwana breakup with crustal modifications. Even if plate tectonics was dominant in the middle to late Archean, difference in spreading and migration rate of the continents led to crustal delamination, and destruction of crustal roots. Comparison of the crustal properties was also found to support their Gondwana connection

  • (ii)

    There are intracratonic variations of Archean crustal thicknesses beneath the studied stations of eastern Gondwana: ~28-32 km in Pilbara; ~36-42 km in Yilgarn; ~46 km in Gawler; 36-41 km in Tanzania; 32-42 km in Kaapvaal; ~35-40 km in Zimbabwe; and ~36-43 km in Antananarivo. It also shows nonuniform crustal thickness variation of the Archean crusts. The stations in this study show the thickness ranges from less than 30 km in parts of Pilbara to more than 40 km beneath parts of Yilgarn, Tanzania, Kaapvaal, Zimbabwe, and Antananarivo. Thick crusts of early Archean age suggest the possibility of high rate of crustal formation or crustal modification, due to which materials could be added later leading to crustal thickening. Vp/Vs ratios mostly in the range ~1.65-1.78 indicate felsic to intermediate nature of the Archean cratons, with evidence of gradational Moho beneath few stations. The rocks inferred from the Vp/Vs ratios and Vs values in our study and silica percentages were found to be similar to the global laboratory experiments and geological studies of rocks in those areas

  • (iii)

    Other than delamination, the role of episodic cycles of crustal growth and evolution is also observed in the pattern of crustal thickness observed at present stage in each division of the Precambrian. The Hadean crust shows ~7 km variation in crustal thickness, followed by ~16-17 km in the Eoarchean and Paleoarchean. The Mesoarchean shows ~11 km Moho depth variation, increasing to ~15 km in the Neoarchean. A number of factors could account for such oscillating differences: operation of plate tectonics through time, alternative high and low crustal regeneration process, repeated episodic recycling, and reworking of crust. Investigation of a new dense network of stations in the Precambrian cratons and integrating numerous geological and geophysical detailed observations in future would help in understanding early Earth in more detail.

Seismic waveform data included in this study can be obtained from the following: Incorporated Research Institutions for Seismology (IRIS) Data Management Centre (DMC), Institut de Physique du Globe de Paris (IPGP), and Réseau sismologique et géodésique français (RESIF) data centers. The networks are as follows: (i) G (doi:10.18715/GEOSCOPE.G), (ii) AU (https://www.fdsn.org/networks/detail/AU/), (iii) XD (doi:10.7914/SN/XD_1994), (iv) XA (doi:10.7914/SN/XA_1997), (v) II (doi:10.7914/SN/II), and (vi) XV (doi:10.7914/SN/XV_2011). Data citations for the above networks are given as follows: (i) Network G (network name: GEOSCOPE). Institut De Physique Du Globe De Paris (IPGP), & Ecole Et Observatoire Des Sciences De La Terre De Strasbourg (EOST) (1982). GEOSCOPE, French Global Network of broadband seismic stations. Institut de Physique du Globe de Paris (IPGP), Université de Paris (doi:10.18715/GEOSCOPE.G). (ii) Network AU (network name: Australian National Seismograph Network ANSN). Operated by Geoscience Australia (https://www.fdsn.org/networks/detail/AU/). (iii) Network XD (network name: Seismic investigations of the Lithospheric Structure of the Tanzanian Craton (Tanzania)). Tom Owens and Andy Nyblade (1994). Seismic investigations of the Lithospheric Structure of the Tanzanian Craton (data set). International Federation of Digital Seismograph Networks (doi:10.7914/SN/XD_1994). (iv) Network XA (network name: Anatomy of an Archean Craton, South Africa, a Multidisciplinary Experiment across the Kaapvaal Craton (SASEK)). Paul Silver (1997). Anatomy of an Archean Craton, South Africa, a Multidisciplinary Experiment across the Kaapvaal Craton (data set). International Federation of Digital Seismograph Networks (doi:10.7914/SN/XA_1997). (v) Network II (network name: Global Seismograph Network- IRIS/IDA (GSN)). Scripps Institution of Oceanography (1986). Global Seismograph Network-IRIS/IDA (data set). International Federation of Digital Seismograph Networks (doi:10.7914/SN/II). (vi) Network XV (network name: INVESTIGATION OF SOURCES OF INTRAPLATE VOLCANISM USING PASSCAL BROADBAND INSTRUMENTS IN MADAGASCAR, THE COMORES, AND MOZAMBIQUE (MACOMO)). Michael Wysession, Douglas Wiens, & Andy Nyblade (2011). INVESTIGATION OF SOURCES OF INTRAPLATE VOLCANISM USING PASSCAL BROADBAND INSTRUMENTS IN MADAGASCAR, THE COMORES, AND MOZAMBIQUE (data set). International Federation of Digital Seismograph Networks (doi:10.7914/SN/XV_2011).

The authors declare that they have no conflicts of interest.

Seismic waveform data has been obtained from the following: Incorporated Research Institutions for Seismology (IRIS) Data Management Centre (DMC), Institut de Physique du Globe de Paris (IPGP), and Réseau sismologique et géodésique français (RESIF) data centers. The networks are as follows: (i) G (doi:10.18715/GEOSCOPE.G), (ii) AU (https://www.fdsn.org/networks/detail/AU/), (iii) XD (doi:10.7914/SN/XD_1994), (iv) XA (doi:10.7914/SN/XA_1997), (v) II (doi:10.7914/SN/II), and (vi) XV (doi:10.7914/SN/XV_2011). Data citation is included after this section. Most of the data processing was done using Seismic Analysis Code (SAC), and figures were made using Generic Mapping Tools (GMT). Author P. Mukherjee is grateful to the INSPIRE fellowship, Department of Science and Technology (DST), India, and MEXT fellowship, Japan Science and Technology (JST), for supporting the study. K. Borah gratefully acknowledges the financial support by the Science and Engineering Research Board (SERB), India (Project No. MTR/2019/001260 and CRG/2020/003726), and Academic Research Funding (ARF) from IISER Kolkata. P. Mukherjee thanks others from IISER Kolkata Earth Science Department who have helped in this study. P. Mukherjee thanks Dr. Y. Ito for guidance and research budget.

1.
Neuendorf
K. K. E.
,
Mehl
J. P.
,
Jackson
J. A.
.
Glossary of Geology
 ,
2005
, 5th,
Alexandria, Virginia
,
American Geological Institute
2.
Mukherjee
P.
Borah
K.
Bora
D. K.
Nature of crust beneath Sri Lanka using teleseismic receiver function
Journal of Asian Earth Sciences
 , 
2020
, vol. 
187, article 104096
 
3.
Capitanio
F. A.
Nebel
O.
Cawood
P. A.
Thermochemical lithosphere differentiation and the origin of cratonic mantle
Nature
 , 
2020
, vol. 
588
 
7836
(pg. 
89
-
94
)
[PubMed]
4.
Korenaga
J.
Crustal evolution and mantle dynamics through Earth history
Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences
 , 
2018
, vol. 
376
 
2132, article 20170408
[PubMed]
[PubMed]
5.
Smithson
S. B.
Johnson
R. A.
Wong
Y. K.
Mean crustal velocity: a critical parameter for interpreting crustal structure and crustal growth
Earth and Planetary Science Letters
 , 
1981
, vol. 
53
 
3
(pg. 
323
-
332
)
[PubMed]
6.
Durrheim
R. J.
Mooney
W. D.
Evolution of the Precambrian lithosphere; seismological and geochemical constraints
Journal of Geophysical Research
 , 
1994
, vol. 
99
 
B8
(pg. 
15359
-
15374
)
7.
Nelson
K. D.
A unified view of craton evolution motivated by recent deep seismic reflection and refraction results
Geophysical Journal International
 , 
1991
, vol. 
105
 
1
(pg. 
25
-
35
)
[PubMed]
8.
Rudnick and Fountain
Nature and composition of the continental crust: a lower crustal perspective
Reviews of Geophysics
 , 
1995
, vol. 
33
 
3
(pg. 
267
-
309
)
[PubMed]
9.
Arndt
N. T.
Goldstein
S. L.
An open boundary between lower continental crust and mantle: its role in crust formation and crustal recycling
Tectonophysics
 , 
1989
, vol. 
161
 
3-4
(pg. 
201
-
212
)
[PubMed]
10.
Jarchow
C. M.
Thompson
G. A.
The nature of the Mohorovicic discontinuity
Annual Review of Earth and Planetary Sciences
 , 
1989
, vol. 
17
 
1
(pg. 
475
-
506
)
11.
Steinhart
J. S.
,
Meyer
R. P.
.
Explosion studies of continental structure
 ,
1961
,
602
,
Carnegie Institute of Washington
12.
Cook
F. A.
White
D. J.
Jones
A. G.
Eaton
D. W. S.
Hall
J.
Clowes
R. M.
How the crust meets the mantle: Lithoprobe perspectives on the Mohorovičić discontinuity and crust–mantle transitionThis article is one of a series of papers published in this Special Issue on the theme Lithoprobe — parameters, processes, and the evolution of a continent
Canadian Journal of Earth Sciences
 , 
2010
, vol. 
47
 
4
(pg. 
315
-
351
)
[PubMed]
13.
Griffin
W. L.
O’reilly
S. Y.
Nixon
P. H.
The composition of the lower crust and the nature of the continental Moho-xenolith evidence
Mantle Xenoliths
 , 
1987
(pg. 
413
-
430
)
14.
Drury
S. A.
Harris
N. B.
Holt
R. W.
Reeves-Smith
G. J.
Wightman
R. T.
Precambrian tectonics and crustal evolution in south India
Journal of Geology
 , 
1984
, vol. 
92
 
1
(pg. 
3
-
20
)
[PubMed]
15.
Hickman
A. H.
,
van Kranendonk
M. J.
.
Archean crustal evolution and mineralization of the northern Pilbara Craton—a field guide
 ,
2008
,
Geological Survey of Western Australia, Record 2008/13
16.
Hickman
A. H.
.
Pilbara Supergroup of the East Pilbara Terrane, Pilbara Craton: Updated Lithostratigraphy and Comments on the Influence of Vertical Tectonics
 ,
2011
,
Geological Survey of Western Australia
,
Annual Review
17.
Menneken
M.
Nemchin
A. A.
Geisler
T.
Pidgeon
R. T.
Wilde
S. A.
Hadean diamonds in zircon from Jack Hills, Western Australia
Nature
 , 
2007
, vol. 
448
 
7156
(pg. 
917
-
920
)
[PubMed]
[PubMed]
18.
Myers
J. S.
Precambrian history of the West Australian craton and adjacent orogens
Annual Review of Earth and Planetary Science
 , 
1993
, vol. 
21
 
1
(pg. 
453
-
485
)
[PubMed]
19.
Myers
J. S.
Preface: Archaean geology of the Eastern Goldfields of Western Australia—regional overview
Precambrian Research
 , 
1997
, vol. 
83
 (pg. 
1
-
10
)
20.
Myers
J. S.
,
Hocking
R. M.
.
Geological Map of Western Australia, 1:2,500,000
 ,
1998
, 13th,
Geological Survey of Western Australian, Perth
21.
Collins
A.
Reddy
S.
Buchan
C.
Mruma
A.
Temporal constraints on Palaeoproterozoic eclogite formation and exhumation (Usagaran Orogen, Tanzania)
Earth and Planetary Science Letters
 , 
2004
, vol. 
224
 
1-2
(pg. 
175
-
192
)
[PubMed]
22.
Möller
A.
Appel
P.
Mezger
K.
Schenk
V.
Evidence for a 2 Ga subduction zone eclogites in the Usagaran belt of Tanzania
Geology
 , 
1995
, vol. 
23
 
12
(pg. 
1067
-
1070
)
23.
Schulter
T.
.
Geology of East Africa
 ,
1997
,
Stuttgart
,
Borntraeger Verlagsbuchhandlung, Science Publishers
24.
Chesley
J. T.
Rudnick
R. L.
Lee
C. T.
Re-Os systematics of mantle xenoliths from the East African Rift: age, structure, and history of the Tanzanian craton
Geochima et Cosmochimica Acta
 , 
1999
, vol. 
63
 
7–8
(pg. 
1203
-
1217
)
[PubMed]
25.
Lee
C. T.
Rudnick
R. L.
Compositionally stratified cratonic lithosphere: petrology and geochemistry of peridotite xenoliths from the Labait tuff cone, Tanzania
Proceedings of the Seventh International Kimberlite Conference
1999
Cape Town, South Africa
(pg. 
503
-
521
)
26.
Ritsema
J.
Ni
S.
Helmberger
D.
Crotwell
H.
Evidence for strong shear velocity reductions and velocity gradients in the lower mantle beneath Africa
Geophysical Research Letters
 , 
1998
, vol. 
25
 
23
(pg. 
4245
-
4248
)
[PubMed]
27.
Mulibo
D. G.
Nyblade
A. A.
Mantle transition zone thinning beneath eastern Africa: evidence for a whole-mantle superplume structure
Geophysical Research Letters
 , 
2013
, vol. 
40
 
14
(pg. 
3562
-
3566
)
[PubMed]
28.
Zegers
T. E.
de Wit
M. J.
Dann
J.
White
S. H.
Vaalbara, Earth's oldest assembled continent? A combined structural, geochronological, and palaeomagnetic test
Terra Nova
 , 
1998
, vol. 
10
 
5
(pg. 
250
-
259
)
[PubMed]
29.
Kröner
A.
Tegtmeyer
A.
Gneiss-greenstone relationships in the Ancient Gneiss Complex of southwestern Swaziland, southern Africa, and implications for early crustal evolution
Precambrian Research
 , 
1994
, vol. 
67
 
1-2
(pg. 
109
-
139
)
[PubMed]
30.
Smith
C. B.
Pearson
D. G.
Bulanova
G. P.
Beard
A. D.
Carlson
R. W.
Wittig
N.
Sims
K.
Chimuka
L.
Muchemwa
E.
Extremely depleted lithospheric mantle and diamonds beneath the southern Zimbabwe craton
Lithos
 , 
2009
, vol. 
112
 
2
(pg. 
1120
-
1132
)
[PubMed]
31.
Blenkinsop
T. G.
Martin
A.
Jelsma
H. A.
Vinyu
M. L.
Wit
M. J.
de
Ashwal
L. D.
The Zimbabwe Craton
Greenstone Belts
 , 
1997
Oxford, England
Clarendon Press
(pg. 
567
-
580
)
32.
BGS (British geological Survey), USGS (US Geological Survey), and GLW (Green Left Weekly)
.
Revision de la cartographie ge ´ologique et minie `re des zones Nord et Centre de Madagascar (Zones A, B et D)
 ,
2008
,
Republique de Madagascar
,
Ministe `rede L’energie et des Mines (MEM/SG/DG/UCP/PGRM)
33.
Paquette
J. L.
Goncalves
P.
Devouard
B.
Nicollet
C.
Micro-drilling ID–TIMS U–Pb dating of single monazites: a new method to unravel complex polymetamorphic evolutions. Application to the UHT granulites of Andriamena (northcentral Madagascar)
Contributions to Mineralogy and Petrology
 , 
2004
, vol. 
147
 
1
(pg. 
110
-
122
)
34.
Ishikawa
M.
Kanao
M.
Structure and collision tectonics of Pan-African orogenic belt—scientific significance of the geotransect for a supercontinent: Gondwanaland
Bulletin of Earthquake Research Institute, University of Tokyo
 , 
2002
, vol. 
77
 (pg. 
287
-
302
)
35.
Harley
S. L.
Black
L. P.
A revised Archaean chronology for the Napier Complex, Enderby Land, from SHRIMP ion-microprobe studies
Antarctic Science
 , 
1997
, vol. 
9
 
1
(pg. 
74
-
91
)
36.
Sheraton
J. W.
Thomson
J. W.
Collerson
K. D.
Fahrig
W. F.
Mafic dyke swarms of Antarctica Halls H.C.
Mafic Dyke Swarms
 , 
1987
Geological Association of Canada
(pg. 
419
-
432
)
37.
Harley
S. L.
Motoyoshi
Y.
Al zoning in orthopyroxene in a sapphirine quartzite: evidence for >1120°C UHT metamorphism in the Napier Complex, Antarctica, and implications for the entropy of sapphirine
Contributions to Mineralogy and Petrology
 , 
2000
, vol. 
138
 
4
(pg. 
293
-
307
)
38.
Ammon
J. C.
The isolation of receiver effects from teleseismicPwaveforms
Bulletin of the Seismological Society of America
 , 
1991
, vol. 
81
 
6
(pg. 
2504
-
2510
)
39.
Langston
C. A.
Corvallis, oregon, crustal and upper mantle receiver structure from teleseismic P and S waves
Bulletin of the Seismological Society of America
 , 
1977
, vol. 
67
 
3
(pg. 
713
-
724
)
40.
Langston
C. A.
Structure under Mount Rainer, Washington, inferred from teleseismic body waves
Journal of Geophysical Research
 , 
1979
, vol. 
84
 (pg. 
4749
-
4762
)
41.
Crotwell
P.
Owens
H. J.
Ritsema
J.
The TauP Toolkit: flexible seismic travel-time and raypath utilities
Seismological Research Letters
 , 
1999
, vol. 
70
 
2
(pg. 
154
-
160
)
42.
Ligorría
J. P.
Ammon
C. J.
Iterative deconvolution and receiver function estimation
Bulletin of the Seismological Society of America
 , 
1999
, vol. 
89
 
5
(pg. 
1395
-
1400
)
43.
Zhu
L.
Kanamori
H.
Moho depth variation in southern California from teleseismic receiver functions
Journal of Geophysical Research
 , 
2000
, vol. 
105
 
B2
(pg. 
2969
-
2980
)
[PubMed]
44.
Sambridge
M. S.
Geophysical inversion with a neighbourhood algorithm: I. searching a parameter space
Geophysical Journal International
 , 
1999
, vol. 
138
 
2
(pg. 
479
-
494
)
[PubMed]
45.
Baranov
A.
Morelli
A.
The Moho depth map of the Antarctica region
Tectonophysics
 , 
2013
, vol. 
609
 (pg. 
299
-
313
)
[PubMed]
46.
Salmon
M.
Kennett
B. L. N.
Stern
T.
Aitken
A. R. A.
The Moho in Australia and New Zealand
Tectonophysics
 , 
2013
, vol. 
609
 (pg. 
288
-
298
)
[PubMed]
47.
Kennett
B. L. N.
Salmon
M.
Saygin
E.
Group AMW
AusMoho: the variation of Moho depth in Australia
Geophysical Journal International
 , 
2011
, vol. 
187
 
2
(pg. 
946
-
958
)
[PubMed]
48.
Youssof
M.
Thybo
H.
Artemieva
I. M.
Levander
A.
Moho depth and crustal composition in Southern Africa
Tectonophysics
 , 
2013
, vol. 
609
 (pg. 
267
-
287
)
49.
Andriampenomanana
F.
Nyblade
A. A.
Wysession
M. E.
Durrheim
R. J.
Tilmann
F.
Julià
J.
Pratt
M. J.
Rambolamanana
G.
Aleqabi
G.
Shore
P. J.
Rakotondraibe
T.
The structure of the crust and uppermost mantle beneath Madagascar
Geophysical Journal International
 , 
2017
, vol. 
210
 
3
(pg. 
1525
-
1544
)
[PubMed]
50.
Borah
K.
Rai
S. S.
Gupta
S.
Prakasam
K. S.
Kumar
S.
Sivaram
K.
Preserved and modified mid-Archean crustal blocks in Dharwar craton: seismological evidence
Precambrian Research
 , 
2014
, vol. 
246
 (pg. 
16
-
34
)
[PubMed]
51.
Sarkar
D.
Kumar
M. R.
Saul
J.
Kind
R.
Raju
P. S.
Chadha
R. K.
Shukla
A. K.
A receiver function perspective of the Dharwar craton (India) crustal structure
Geophysical Journal International
 , 
2003
, vol. 
154
 
1
(pg. 
205
-
211
)
[PubMed]
52.
Kay
R. W.
Kay
M. S.
Delamination and Delamination magmatism
Tectonophysics
 , 
1983
, vol. 
219
 (pg. 
177
-
189
)
53.
Drummond
B. J.
A crustal profile across the Archaean Pilbara and Northern Yilgarn cratons
Journal of Australian Geology and Geophysics
 , 
1979
, vol. 
4
 
1979
(pg. 
171
-
180
)
54.
Yuan
H.
Secular change in Archaean crust formation recorded in Western Australia
Nature Geoscience
 , 
2015
, vol. 
8
 
10
(pg. 
808
-
813
)
[PubMed]
55.
Beintema
K.
White
S.
Wijbrans
J. R.
Compression-extension tectonics in the West Pilbara Craton, Western Australia
Journal of the Virtual Explorer
 , 
2003
, vol. 
13
 (pg. 
1
-
42
)
56.
Drummond
D. J.
Collins
C. D. N.
Seismic evidence for underplating of the lower continental crust of Australia
Earth and Planetary Science Letters
 , 
1986
, vol. 
79
 
3-4
(pg. 
361
-
372
)
[PubMed]
57.
Yuan
H.
Crustal composition of the Western Australian craton and implications for craton formation tectonics
AGU 2014 Fall Meeting abstracts, id. V43C-4898
2015b
58.
Tedla
E. G.
van der Meijde
M.
Nyblade
A. A.
van der Meer
F. D.
A crustal thickness map of Africa derived from a global gravity field model using Euler deconvolution
Geophysical Journal International
 , 
2011
, vol. 
187
 
1
(pg. 
1
-
9
)
[PubMed]
59.
Prodehl
C.
Keller
G. R.
Khan
M. A.
Crustal and upper mantle structure of the Kenya Rift
Tectonophysics
 , 
1994
, vol. 
236
 
483
(pg. 
437
-
452
)
60.
Fuchs
K.
Altherr
B.
Muller
B.
Prodehl
C.
Structure and dynamic processes in the lithosphere of the Afro-Arabian rift system
Tectonophysics
 , 
1997
, vol. 
278
 (pg. 
1
-
352
)
61.
Hodgson
I.
Finnigan
K.
Gallacher
R.
Keir
D.
Ebinger
C.
Mtelela
K.
Crustal structure at a young continental rift: a receiver function study from the Tanganyika rift
Tectonics
 , 
2017
, vol. 
36
 
12
(pg. 
2806
-
2822
)
[PubMed]
62.
Halls
H. C.
Burns
K. G.
Bullock
S. J.
Batterham
P. M.
Halls
H. C.
Fahrig
W. F.
Mafic dyke swarms of Tanzania interpreted from areomagnetic data
Mafic dyke swarms
 , 
1987
, vol. 
34
 
Geol. Soc. Canada Spec. Pap.
(pg. 
173
-
186
)
63.
Rindraharisaona
E. J.
Tilmann
F.
Yuan
X.
Rümpker
G.
Giese
J.
Rambolamanana
G.
Barruol
G.
Crustal structure of southern Madagascar from receiver functions and ambient noise correlation: implications for crustal evolution
Journal of Geophysical Research
 , 
2017
, vol. 
122
 
2
(pg. 
1179
-
1197
)
[PubMed]
64.
Boger
S. D.
,
Ferreira
C. A. M.
,
Hirdes
W.
,
Schulte
B.
,
Jenett
T.
,
Dallwig
R.
.
Carte géologique de Madagascar 1/500,000, feuille N°10 – Fianarantsoa
 ,
2009
,
Antananarivo
,
Ministere de L’energie et des mines, Project de Gouvernance des Ressources Minér- ales (PGRM)
65.
Handke
M. J.
Tucker
R. D.
Ashwal
L. D.
Neoproterozoic continental arc magmatism in west–central Madagascar
Geology
 , 
1999
, vol. 
27
 (pg. 
351
-
354
)
66.
Randriamananjara
T.
.
Carte Géologique de Madagascar 1/100,000, Feuille N°R51 – Vohitrandriana
 ,
2008
,
Antananarivo
,
Ministere de L’energie et des mines, Project de Gouvernance des Ressour- ces Minérales
67.
Christensen
N. I.
Poisson’s ratio and crustal seismology
Journal of Geophysical Research: Solid Earth
 , 
1996
, vol. 
101
 
B2
(pg. 
3139
-
3156
)
[PubMed]
68.
Audet
P.
Burgman
R.
Possible control of subduction zone slow-earthquake periodicity by silica enrichment
Nature
 , 
2014
, vol. 
510
 
7505
(pg. 
389
-
392
)
[PubMed]
[PubMed]
69.
Mukherjee
P.
Ito
Y.
Garcia
E. S.
Martinez
R. P.
Shibutani
T.
The Chile subduction zone: nature of the lithosphere between alternating regions of slow earthquakes versus non slow earthquakes
EGU General Assembly
 , 
2020
, vol. 
2020, article 4035
 
70.
Agangi
A.
Hoffmann
A.
Elburg
M.
A review of Palaeoarchaean felsic volcanism in the eastern Kaapvaal craton: linking plutonic and volcanic records
Geoscience Frontiers
 , 
2018
, vol. 
9
 
3
(pg. 
667
-
688
)
[PubMed]
71.
Braun
I.
Kriegsman
L. M.
Proterozoic crustal evolution of southernmost India and Sri Lanka
Geological Society London. Special Publications
 , 
2003
, vol. 
206
 
1
(pg. 
169
-
202
)
[PubMed]
72.
Chevrot
S.
Van der Hilst
R. D.
The Poisson ratio of the Australian crust: geological and geophysical implications
Earth and Planetary Science Letters
 , 
2000
, vol. 
183
 
1-2
(pg. 
121
-
132
)
[PubMed]
73.
Artemieva
I. M.
Meissner
R.
Crustal thickness controlled by plate tectonics: a review of crust-mantle interaction processes illustrated by European examples
Tectonophysics
 , 
2012
, vol. 
530-531
 (pg. 
18
-
49
)
[PubMed]
74.
Brown
M.
Johnson
T.
Gardiner
N. J.
Plate tectonics and the Archean earth
Annual Review of Earth and Planetary Sciences
 , 
2020
, vol. 
48
 
1
(pg. 
291
-
320
)
75.
Sleep
N.
Evolution of the continental lithosphere
Annual Review of Earth and Planetary Sciences
 , 
2005
, vol. 
33
 
1
(pg. 
369
-
393
)
[PubMed]
76.
Stern
R. J.
Arc assembly and continental collision in the Neoproterozoic East African Orogen: Implications for the consolidation of Gondwanaland
Annual Review of Earth and Planetary Sciences
 , 
1994
, vol. 
22
 (pg. 
319
-
351
)
77.
Fitzsimons
I. C. W.
Yoshida
M.
Windley
B. F.
Dasgupta
S.
Proterozoic basement provinces of southern and southwestern Australia, and their correlation with Antarctica
Proterozoic East Gondwana: Supercontinent Assembly and Breakup
 , 
2002
Geological Society of London Special Publication 206
(pg. 
93
-
130
)
78.
Betts
P. G.
Giles
D.
Lister
G. S.
Frick
L.
Evolution of the Australian lithosphere
Australian Journal of Earth Sciences
 , 
2002
, vol. 
49
 
4
(pg. 
661
-
695
)
Exclusive Licensee GeoScienceWorld. Distributed under a Creative Commons Attribution License (CC BY 4.0).