Hello Guest User, You are visiting this website from a computer with an IP address of 172.70.126.84 with the name of '?' since Tue Apr 23, 2024 at 3:04:39 PM PT for approx. 0 minutes now.
The Sheahan Diamond Literature Reference Compilation - Scientific and Media Articles based on Major Keyword - Geophysics - Seismic
The Sheahan Diamond Literature Reference Compilation is compiled by Patricia Sheahan who publishes on a monthly basis a list of new scientific articles related to diamonds as well as media coverage and corporate announcements called the Sheahan Diamond Literature Service that is distributed as a free pdf to a list of followers. Pat has kindly agreed to allow her work to be made available as an online digital resource at Kaiser Research Online so that a broader community interested in diamonds and related geology can benefit. The references are for personal use information purposes only; when available a link is provided to an online location where the full article can be accessed or purchased directly. Reproduction of this compilation in part or in whole without permission from the Sheahan Diamond Literature Service is strictly prohibited. Return to Diamond Keyword Index
Sheahan Diamond Literature Reference Compilation - Scientific Articles by Author for all years
Each article reference in the SDLRC is tagged with one or more key words assigned by Pat Sheahan to highlight the main topics of the article. In an effort to make it easier for users to track down articles related to a specific topic, KRO has extracted these key words and developed a list of major key words presented in this Key Word Index to which individual key words used in the article reference have been assigned. In most of the individual Key Word Reports the references are in crhonological order, though in some such as Deposits the order is first by key word and then chronological. Only articles classified as "technical" (mainly scientific journal articles) and "media" (independent media articles) are included in the Key Word Index. References that were added in the most recent monthly update are highlighted in yellow.
While "gravity", "magnetic" and "induced polarity" geophysical surveys focus on lateral differences in the composition of the earth at different scales, the "seismic surveys" referred to by the key word Geophysics - Seismic attempt to describe differences in the earth's composition in the vertical dimension. These articles are all of a scientific nature and have little to do with diamonds other than describing the inner structure of the earth and highlighting those regions conducive to diamond formation and kimberlite eruptions.
Integration of COCORP deep relfection and magnetic anomaly analysis in the southeast USA:implications for origin of the Brunswick and East Coast magneticanomalies
Geological Society of America (GSA) Bulletin, Vol. 102, No. 2, February pp. 271-279
Layered Proterozoic rocks and a Proterozoic angular unconformity beneath the U.S. midcontinent on Cocorp and reprocessed industry seismic reflectiondata
Geological Society of America Annual Meeting Abstract Volume, Vol. 23, No. 5, San Diego, p. A 315
Anisotropy of magnetic remanence- a brief review of mineralogical physical origins and geological applications, and comparison with susceptibilityanisotrop
Pure and Applied Geophysics, Vol. 136, No. 1, May pp. 1-28
Constraints on the structure and tectonic development of the early Mesozoic south Georgia Rift, southeastern United States; seismic reflection data processing &int
Preliminary results of new Cocorp deep seismic reflection profiling across the intracratonic Williston Basin and underlying early Proterozoic Trans-HudsonOrogen
Geological Society of America Annual Meeting Abstract Volume, Vol. 23, No. 5, San Diego, p. A 315
Hutchinson, D.R., Lee, M.W., Behrendt, J., Cannon, W.F., Green
Variations in the reflectivity of the Moho transition zone beneath The midcontinent Rift System of North America. Results from true amplitude Glimpcedata
Journal of Geophysical Research, Vol. 97, No. B4, April 10, pp. 4721-4738
A Tertiary asthenospheric flow beneath the southern French Massif Central indicated by upper mantle seismic anisotropy and related to west mediterranean extension.
Earth and Planetary Science Letters, Vol. 202, 1, pp.31-47.
Products of 2.65 - 2.58 Ga orogenesis in the Slave Province correlated with Slave - Northern Cordillera Lithospheric Evolution ( SNORCLE) seismic reflection
Canadian Journal of Earth Science, Vol. 39,8,August pp. 1189-1200.
Northwest Territories
Geophysics - seismic, Slave Province - tectonics, geodynamics
Contrasting lower crustal evolution across an Archean Proterozoic suture: physical, chemical and geochronologic studies of lower crustal xenoliths in southern Wyoming and northern Colorado.
American Geophysical Union, Geophysical Monograph, No. 154, pp. 139-162.
Pilidou, SA., Priestley, K., Gudmundsson, O., Debayle, E.
Upper mantle S-wave speed heterogeneity and anisotropy beneath the North Atlantic from regional surface wave tomography: the Iceland and Azores plumes.
Geophysical Journal International, Vol. 159, 3, pp. 1057-1076.
High stress tests show surprising change in a mantle mineral's behaviour. Compressing common iron bearing minerals to the pressures found deep within Earth.
Skeleton, A., Whitmarsh, R., Arghe, F., Crill, P., Koyi, H.
Constraining the rate and extent of mantle serpentinization from seismic and petrological data: implications for chemosynthesis and tectonic processes.
Skelton, A., Whitmarsh, R., Arghe, F., Crill, P., Koyi, H.
Constraining the rate and extent of mantle serpentinization from seismic and petrological data: implications for chemosynthesis and tectonic processes.
Upper mantle tomographic Vp and Vs images of the Rocky Mountains in Wyoming, Colorado, New Mexico: evidence for a thick heterogeneous chemical lithosphere.
American Geophysical Union, Geophysical Monograph, No. 154, pp. 329-346.
Using shear wave splitting measurements to investigate the upper mantle anisotropy beneath the North Chin a Craton: distinct variation from east to west.
Geophysical Research Letters, Vol. 32, 10, May 28, DOI 10.1029/2005 GLO22585
An investigation of upper mantle heterogeneity beneath the Archean and Proterozoic crust of western Canada from lithoprobe controlled source seismic experiments.
Tectonophysics, Vol. 416, 1-4, April 5, pp. 187-207.
Stability and compressibility of the high pressure phases of AL203 up to 200 GPa: implications for the electrical conductivity of the base of the lower mantle.
Earth and Planetary Science Letters, Vol. 246, 3-4, pp. 326-335.
The Tornquist Zone, north east inclining lithospheric transition at the south western margin of the Baltic Shield: revealed through a nonlinear teleseismic tomographic inversion.
Tectonophysics, Vol. 416, 1-4, April 5, pp. 151-166.
A compositional anomaly at the Earth's core mantle boundary as an anchor to the relatively slowly moving surface hotspots and as source to the DUPAL anomaly.
Earth and Planetary Science Letters, Vol. 246, 1-2, June 15, pp. 138-148.
Crust mantle boundary in the central Fennoscandian shield: constraints from wide angle P and S wave velocity models and new results of reflection profiling in Finland
Journal of Geophysical Research, Vol. 112, B4, B04302.
Using xenoliths to explore variations in upper mantle composition and the relation of composition to seismic velocity structure beneath the Colorado Plateau.
Geological Society of America Annual Meeting 2007, Denver Oct. 28, 1p. Abstract
Constraints on the current rate of deformation and surface uplift of the Australian continent from a new seismic database and low T thermochronological data.
Australian Journal of Earth Sciences, Vol. 56, 2, pp. 99-110.
Cobden, L., Goes, S., Ravenna, M., Styles, E., Cammarano, F., Gallagher, K., Connolly, J.
Thermochemical interpretation of 1-D seismic dat a for the lower mantle: the significance of nonadiabiatic thermal gradients and compositional heterogeneity.
Journal of Geophysical Research, Vol. 114, B 11, B11309
Janik, T., Kozlovskaya, E., Helikkinen, P., Tliniemi, J.
Evidence for preservation of crustal root beneath the Proterozoic Lapland-Kola orogen ( northern Fennoscandian shield) derived from P and S wave models.
Journal of Geophysical Research, Vol. 114. B 6, B06308.
Yoneda, A., Chen, G., Spetzler, H.A., Getting, I.G.
The effect of composition, temperature and pressure on the elasticity of olivine and garnet: implications for interpreting seismic velocity variations in mantle.
The lithospheric root beneath Hudson Bay, Canada from Rayleigh wave dispersion: no clear seismological distinction between Archean and Proterozoic mantle.
Joint seismic geodynamic mineral physical modeling of African geodynamics: a reconciliation of deep mantle convection with surface geophysical constraints.
Earth and Planetary Science Letters, Vol. 295, 3-4, pp. 329-341.
Crustal structure of the Zimbabwe craton and the Limpopo belt of southern Africa: new constraints from seismic dat a and implications for its evolution.
South African Journal of Geology, Vol. 112, pp. 213-228.
Jones, A.G., Plomerova, J., Korja, T., Sodoudi, F., Spakman, W.
Europe from the bottom up: a statistical examination of the central and northern European lithosphere asthenosphere boundary comparing seismological & EM
Craig, T.J., Jackson, J.A., Priestley, K., McKenzie, D.
Earthquake distribution patterns in Africa: their relationship to variations in lithospheric and geological structure, and their rheological implication
Geophysical Journal International, Vol. 185, 1, pp. 403-404.
Obrebski, M., Allen, R.M., Pollitz, F., Hung, S-H.
Lithosphere asthenosphere interaction beneath the western United States from the joint inversion of body-wave traveltimes and surface wave phase veolocities.
Geophysical Journal International, March 25, In press available
Mantle flow and lithosphere asthenosphere coupling beneath the southwestern edge of the North American craton: constraints from shear wave splitting measurements.
Earth and Planetary Science Letters, Vol. 402, pp. 209-220.
Intricate heterogeneous structures of the top 300 km of the Earth's inner core inferred from global array data.II. Frequency dependence of inner core attentuation and its implication.
Earth and Planetary Science Letters, Vol. 405, Nov. pp. 231-243.
Abstract: The problem of the existence of the asthenosphere for old Precambrian cratons is still discussed. In order to study the seismic lithosphere-asthenosphere boundary (LAB) beneath the Baltic Shield, we used records of nine local earthquakes with magnitudes ranging from 2.7 to 5.9. To model the LAB, original data were corrected for topography and Moho depth using a reference model with a 46-km-thick crust. For two northern events at Spitsbergen and Novaya Zemlya, we observe a low-velocity layer, 60-70-km-thick asthenosphere, and the LAB beneath Barents Sea was found at depth of c. 200 km. Sections for other events show continuous first arrivals of P-waves with no evidence for "shadow zone" in the whole range of registration, which could either be interpreted as the absence of the asthenosphere beneath the central part of the Baltic Shield, or that the LAB in this area occurs deeper (>200 km). The relatively thin low-velocity layer found beneath southern Sweden, 15 km below the Moho, could be interpreted as small-scale lithospheric heterogeneities, rather than asthenosphere. Differentiation of the lower lithosphere velocities beneath the Baltic Shield could be interpreted as regional heterogeneity or as anisotropy of the Baltic Shield lithosphere, with high velocities approximately in the east-west direction, and slow velocities approximately in the south-north direction.
Geophysical Journal International, Vol. 201, 1, pp. 61-69.
South America, Brazil
Geophysics - seismic
Abstract: Results of P-wave traveltime seismic tomography in central Brazil unravel the upper-mantle velocity structure and its relationship with the tectonic framework. Data were recorded between 2008 and 2012 at 16 stations distributed over the study area, and were added to the database used by Rocha et al. to improve the resolution of anomalies, and to image the surrounding regions. The main objective was to observe the upper-mantle boundary zone between the Amazonian and São Francisco cratons, represented by mobile belts, inside the Tocantins Province, and to study the lithosphere related to the collision between these two cratons during the Neoproterozoic. A set of low-velocity anomalies was observed crossing the study area in the NE–SW direction, in agreement with the main trend of the Transbrasiliano lineament. The region where the anomalies are located was interpreted as the zone separating the Amazonian and São Franciscan palaeoplates. There is a good correlation between the low-velocity anomalies and the high seismicity of this region, suggesting that it is a region of weakness, probably related to lithospheric thinning. High velocities were observed under the Amazonian and São Francisco cratons. A model is proposed for the lithospheric subsurface in central Brazil, emphasizing the boundary zone between the main palaeoplates in the study area. After merging both databases, the low-velocity anomalies in the central part of the study area suggest tectonic partitioning of the lithosphere. Synthetic tests show that the tomography results are robust.
Density, temperature and composition of the North American lithosphere - new insights from a joint analysis of seismic, gravity and mineral physics data: 1. density structure of the crust and upper mantle.
Geochemistry, Geophysics, Geosystems: G3, Vol. 15, 12, pp. 4781-4807.
Earth and Planetary Science Letters, Vol. 429, pp. 101-109.
Mantle
Seismic -anisotropy
Abstract: The seismic structure of Earth’s inner core is highly complex, displaying strong anisotropy and further regional variations. However, few seismic waves are sensitive to the inner core and fundamental questions regarding the origin of the observed seismic features remain unanswered. Thus, newtechniques to observe different types of inner corewaves are imperative to improve data coverage. Here, we detail our method for detecting exotic inner core phases such as PKJKP and PKIIKP, using inner core compressional waves as proof of concept.We use phase weighted stacking on long period data from a global distribution of stations, and employ several synthetic methods, including normal mode summation and SPECFEM, to identify and confirm the inner core phases. We present evidence for two observations of exotic inner core compressional waves, and apply the technique to a previously detected inner core shear wave. A possible new inner core shear wave remains unconfirmed. Additionally, we show how our method is important for rejecting potential observations, and distinguishing between waves with similar traveltime and slowness. The method is most successful for detecting exotic inner core compressional waves, and will provide a new approach for studying the compressional wave structures in the upper inner core.
Abstract: North America provides an important test for assessing the coupling of large continents with heterogeneous Archean- to Cenozoic-aged lithospheric provinces to the mantle flow. We use the unprecedented spatial coverage of the USArray seismic network to obtain an extensive and consistent data set of shear wave splitting intensity measurements at 1436 stations. Overall, the measurements are consistent with simple shear deformation in the asthenosphere due to viscous coupling to the overriding lithosphere. The fast directions agree with the absolute plate motion direction with a mean difference of 2° with 27° standard deviation. There are, however, deviations from this simple pattern, including a band along the Rocky Mountain front, indicative of flow complication due to gradients in lithospheric thickness, and variations in amplitude through the central United States, which can be explained through varying contributions of lithospheric anisotropy. Thus, seismic anisotropy may be sourced in both the asthenosphere and lithosphere, and variations in splitting intensity are due to lithospheric anisotropy developed during deformation over long time scales.
Geophysical Journal International, Vol. 203, pp. 910-926.
India
Geophysics - seismics
Abstract: We present a comprehensive study of thickness and composition of the crust; and the nature of crust-mantle boundary beneath Southern India using P-wave receiver function from 119 seismic stations. Data from distributed network of seismograph location encompass geological domains like mid to late Archean Dharwar craton, Archean and Proterozoic metamorphic terrains, Proterozoic basin, rifted margins and escarpments, and Deccan volcanics. Except for the mid to lower crust exhumed Archean terrains (of West Dharwar and Southern Granulite) all other geological domains have crustal thickness in the range 33-40 km. In the western Dharwar, crustal thickness increases from ?40 km in the north to over 50 km in the south. The Archean domain of granulite terrain is thicker (40-45 km) and more mafic compared to its counterpart in south deformed at 550 Ma. Most of the crustal blocks have low to moderate Vp/Vs (1.72-1.76) representing a felsic to intermediate composition. Exception to the above include Archean granulite terrain with high Vp/Vs (1.76–1.81) suggestive of more mafic crust beneath them. When accounted for the paleo burial depth of 15-25 km, the study suggests a possible Himalaya-Tibet like scenario beneath the mid-late Archean in southwestern Dharwar and north granulite terrain whose deeper crust has progressively densified. This led to a gradational crust-mantle transition that is otherwise sharp elsewhere. The study suggests a more homogenized and felsic nature of the Precambrian crust beneath the terrains formed after 2.6 Ga, possibly due to delamination of the mafic lower crust. Our study does not suggest any distinction between late Archean and Proterozoic crust. The Deccan volcanism at 65 Ma does not appear to have altered the crustal character beneath it and is similar to the adjoining late Archean east Dharwar craton. The western Ghat escarpment and the coastal plain formed due to separation of India from Madagascar are underlain by mafic lower crust.
South African Journal of Geology, Vol. 118, 3, pp. 225-248.
Africa, South Africa
Geophysics - seismics
Abstract: Here we present a comprehensive depth and thickness map of the main Karoo and Cape Basins using borehole and reflection seismic data. The depth to the Whitehill Formation, which is the focus of current shale gas interest within the Karoo, is also mapped. Change: The deepest part of the basin is in the south, along the northern boundary of the Cape Fold Belt (~4000 m in the southwest Karoo and ~5000 m in the southeast; ~5500 to 6000 m sediment thickness). The Whitehill Formation along this boundary reaches a depth of ~3000 m in the southwest and ~4000 m in the southeast. Limited borehole data in the southeastern Karoo show a broad deepening of the basin here compared to the southwestern Karoo. In the southeast near East London faulting has resulted in deepening of the basin close to the coast, with the Whitehill Formation deepening to over ~5000 km. Seismic and borehole data show that the Cape Supergroup pinches out below the Karoo Basin around Beaufort West and Graaff-Reinet in the southern Karoo (32.6°S for the Bokkeveld and 32.4°S for the Table Mountain Group). The Cape Supergroup reaches thicknesses of around 4 km in the south. The gravity effect of these sediments does not account for the Cape Isostatic Anomaly (CIA) in the southern part of the Karoo Basin near Willowmore and Steytlerville, i.e., an ~45 mGal Bouguer gravity low. A refraction seismic profile over the anomaly shows this region is associated with a large volume of low velocity/density shallow sediments (4.5 m/s2, 2500 kg/m3), as well as a low velocity/density anomaly associated with a normal fault and the Klein Winterhoek Thrust Fault (5.5 m/s2, 2650 kg/m3). These low density shallow sediments are explained by uplift of Karoo and Cape sediments of ~2 km or greater that is evident on Soekor reflection seismic data. This deformation has brought lower density shales (1800 to 2650 kg/m3) of the Ecca Group closer to the surface. These shallower features along with a deeper lower crust in this region (6.5 m/s2, 2900 kg/m3) are interpreted to account for the CIA.
Journal of Geophysical Research, Vol. 120, 8, pp. 5782-5803.
United States, Arkansas
Geophysics - seismics
Abstract: Global Navigation Satellite System data across the New Madrid seismic zone (NMSZ) in the central United States over the period from 2000 through 2014 are analyzed and modeled with several deformation mechanisms including the following: (1) creep on subsurface dislocations, (2) postseismic frictional afterslip and viscoelastic relaxation from the 1811–1812 and 1450 earthquakes in the NMSZ, and (3) regional strain. In agreement with previous studies, a dislocation creeping at about 4 mm/yr between 12 and 20 km depth along the downdip extension of the Reelfoot fault reproduces the observations well. We find that a dynamic model of postseismic frictional afterslip from the 1450 and February 1812 Reelfoot fault events can explain this creep. Kinematic and dynamic models involving the Cottonwood Grove fault provide minimal predictive power. This is likely due to the smaller size of the December 1811 event on the Cottonwood Grove fault and a distribution of stations better suited to constrain localized strain across the Reelfoot fault. Regional compressive strain across the NMSZ is found to be less than 3?×?10?9/yr. If much of the present-day surface deformation results from afterslip, it is likely that many of the earthquakes we see today in the NMSZ are aftershocks from the 1811–1812 New Madrid earthquakes. Despite this conclusion, our results are consistent with observations and models of intraplate earthquake clustering. Given this and the recent paleoseismic history of the region, we suggest that seismic hazard is likely to remain significant.
Geological Society of America Special Paper, No. 514, pp. SPE514-13.
Mantle
Geophysics - seismics
Abstract: We use the fact that geoid anomalies are directly related to the local dipole moment of the density-depth distribution to help constrain density variations within the lithosphere and the associated tectonic stresses. The main challenge with this approach is isolating the upper mantle geoid contribution from the full geoid (which is dominated by sources in the lower mantle). We address this issue by using a high-pass spherical harmonic filtering of the EGM2008-WGS84 geoid to produce an "upper mantle" geoid. The tectonic implications of the upper mantle are discussed in terms of plate tectonics and intraplate stresses. We find that globally there is about a 9 meter geoid step associated with the cooling oceanic lithosphere that imparts a net force of ~2.5x1012 N/m in the form of "ridge push" - a magnitude that is consistent with 1-d models based on first-order density profiles. Furthermore, we ind a consistent 6 meter geoid step across passive a continental margin which has the net effect of educing the compressive stresses in the continents due to the ridge force. Furthermore, we use the pper mantle geoid to reevaluate the tectonic reference state which previously studies estimated using n assumption of Airy-based isostasy. Our evaluation of the upper mantle geoid confirms the near quivalence of the gravitational potential energy of continental lithosphere with an elevation of about 750 meters and the mid-ocean ridges. This result substantiates early conclusions about the tectonic reference state and further supports the prediction that continental regions are expected to be in a slightly extensional state of stress.
Geophysical Journal International, Vol. 202, pp. 142-162.
Canada, Alberta
Geophysics - seismics
Abstract: We analyse continuous recordings from 23 broadband seismic stations near Alberta, the southwestern sector of the Western Canada Sedimentary Basin. Noise correlation tomo-graphy based on vertical-component seismograms reveals below-average shear velocities at shallow and middle crustal depths in central Alberta, spanning across Proterozoic accreted terranes and Archean microcontinents. This observation likely results from extensive plate convergence and crustal melting during the Proterozoic eon. The overall correlation between the crustal velocities and presumed basement domains is lower than expected, however. In the lower crust, the main pattern of shear velocities is relatively concordant with the reported domain boundaries and key Precambrian structures appear to be intact. The shear velocities beneath the Loverna Block, the largest constituent of the Hearne craton, are 10?per?cent higher than the regional average. This prominent northeast striking seismic anomaly is moderately correlated with the regional heat flow and potentially represents the remnant core of the Archean Hearne province. The associated high velocities extend into the western part of the Medicine Hat Block, a possible Archean microcontinent with a debatable origin, and contribute to a strong east-west structural gradient in the lower crust. The presence and the continuity of this anomalous structure imply extensive communications among the various basement domains in southern Alberta during the assembly of the North American continent.
Abstract: Stable continental cratons are the oldest geologic features on the planet. They have survived 3.8 to 2.5 billion years of Earth’s evolution1, 2. The key to the preservation of cratons lies in their strong and thick lithospheric roots, which are neutrally or positively buoyant with respect to surrounding mantle3, 4. Most of these Archaean-aged cratonic roots are thought to have remained stable since their formation and to be too viscous to be affected by mantle convection2, 3, 5. Here we use a combination of gravity, topography, crustal structure and seismic tomography data to show that the deepest part of the craton root beneath the North American Superior Province has shifted about 850?km to the west-southwest relative to the centre of the craton. We use numerical model simulations to show that this shift could have been caused by basal drag induced by mantle flow, implying that mantle flow can alter craton structure. Our observations contradict the conventional view of cratons as static, non-evolving geologic features. We conclude that there could be significant interaction between deep continental roots and the convecting mantle.
Earth and Planetary Science Letters, Vol. 431, pp. 225-235.
Australia
Geophysics - seismics - cratons
Abstract: A direct image of P-wave reflectivity in the lithosphere and asthenosphere beneath seismic stations is extracted from stacked autocorrelograms of continuous component records. The autocorrelograms emphasise near vertically travelling waves, so that multiples are more muted than in receiver function studies and it is possible to work at higher frequencies than for receiver functions. Across a wide range of geological environments in Australia, in the 0.5-4.0 Hz frequency band, distinct reflections are seen in the crust underlain by weaker reflectivity in the lithosphere and asthenosphere. The base of crustal reflectivity fits well with Moho estimates from other classes of information. Few mantle reflectors have been seen in conventional reflection profiling at frequencies above 10 Hz; the presence of reflections in the 0.5-4.0 Hz band suggests variations on vertical scales of a few hundred metres with amplitudes of the order of 1%. There are slight indications of a change of reflection character in the lower part of the lithosphere in the transition to the asthenosphere. At a few stations there is a very clear lamination at asthenospheric depth, as well as reflections from the base of the S wave low velocity zone. Reflection bands often occur at depths where discontinuities have been inferred from S wave receiver function work at the same station, but would not by themselves be distinctive of a mid-lithosphere discontinuity.
Abstract: We analyzed satellite gravity and geoid anomaly and topography data to determine the 3D lithospheric density structure of the Singhbhum Protocontinent. Our density model shows that distinct vertical density heterogeneities exist throughout the lithosphere beneath the Singhbhum Protocontinent. The crustal structure identified includes a lateral average crustal density variation from 2800 to 2890 kg/m3 as well as a relatively flat Moho at 35-40 km depth in Singhbhum Protocontinent and Bastar Craton. A similar Moho depth range is found for the Mahanadi, Damodar, and Bengal basins. In the northern part of the area, Moho undulates between more than 40 km under the confluence of Mahanadi-Damodar Gondwana basins and the Ganga foreland basin, and 36-32 km under the Eastern Ghats Mobile belt and finally reaches 24 km in the Bay of Bengal. The lithosphere-asthenosphere boundary (LAB) across the Singhbhum Protocontinent is at a depth of about 130-140 km. In the regions of Bastar Craton and Bengal Basin, the LAB dips to about 155 ± 5 km depth. The confluence of Mahanadi and Damodar Gondwana basins toward the north-west and the foreland Ganga Basin toward the north are characterized by a deeper LAB lying at a depth of over 170 and 200 km, respectively. In the Bay of Bengal, the LAB is at a shallower depth of about 100-130 km except over the 85 0E ridge (150 km), and off the Kolkata coast (155 km). Significant density variation as well as an almost flat crust-mantle boundary indicates the effect of significant crustal reworking. The thin (135-140 km) lithosphere provides compelling evidence of lithospheric modification in the Singhbhum Protocontinent. Similarities between the lithospheric structures of the Singhbhum Craton, Chhotanagpur Gneiss Complex, and Northern Singhbhum Mobile Belt confirm that the repeated thermal perturbation controlled continental lithospheric modification in the Singhbhum Protocontinent.
Abstract: The formation mechanisms for early Archaean continental crust are controversial. Continental crust may have accumulated via horizontal accretion in modern-style subduction zones or via vertical accretion above upper mantle upwelling zones. However, the characteristics of the continental crust changes at the transition between the Archaean and Proterozoic eons, suggesting that continental crust did not form in subduction zones until at least the late Archaean. Here I use seismic receiver function data to analyse the bulk properties of continental crust in Western Australia, which formed and stabilized over a billion years in the Archaean. I find that the bulk seismic properties of the crust cluster spatially, with similar clusters confined within the boundaries of tectonic terranes. I use local Archaean crustal growth models to show that both plume and subduction processes may have had a role in creating crust throughout the Archaean. A correlation between crustal age and the bulk seismic properties of the crust reveals a trend: from about 3.5?Gyr ago (Ga) to the end of the Archaean, the crust gradually thickened and simultaneously became more evolved in composition. I propose that this trend reflects the transition between crust dominantly formed above mantle plumes, to crust formed in subduction zones-a transition that may reflect secular cooling of Earth’s mantl
Geochemistry, Geophysics, Geosystems: G3, Vol. 16, 10, pp. 3436-3455.
Mantle
Geophysics - seismics
Abstract: Relating seismic anisotropy to mantle flow requires detailed understanding of the development and evolution of olivine crystallographic preferred orientation (CPO). Recent experimental and field studies have shown that olivine CPO evolution depends strongly on the integrated deformation history, which may lead to differences in how the corresponding seismic anisotropy should be interpreted. In this study, two widely used numerical models for CPO evolution—D-Rex and VPSC—are evaluated to further examine the effect of deformation history on olivine texture and seismic anisotropy. Building on previous experimental work, models are initiated with several different CPOs to simulate unique deformation histories. Significantly, models initiated with a preexisting CPO evolve differently than the CPOs generated without preexisting texture. Moreover, the CPO in each model evolves differently as a function of strain. Numerical simulations are compared to laboratory experiments by Boneh and Skemer (2014). In general, the D-Rex and VPSC models are able to reproduce the experimentally observed CPOs, although the models significantly over-estimate the strength of the CPO and in some instances produce different CPO from what is observed experimentally. Based on comparison with experiments, recommended parameters for D-Rex are: M*?=?10, ?*?=?5, and ??=?0.3, and for VPSC: ??=?10-100. Numerical modeling confirms that CPO evolution in olivine is highly sensitive to the details of the initial CPO, even at strains greater than 2. These observations imply that there is a long transient interval of CPO realignment which must be considered carefully in the modeling or interpretation of seismic anisotropy in complex tectonic settings.
Geochemistry, Geophysics, Geosystems: G3, Vol. 16, 10, pp. 3343-3361.
United States
Geophysics - seismics
Abstract: Teleseismic P wave amplitude spectra from deep earthquakes recorded by USArray are inverted for maps of upper mantle ?t* for multiple frequency bands within 0.08-2 Hz. All frequency bands show high ?t* regions in the southwestern U.S., southern Rocky Mountains, and Appalachian margin. Low ?t* is more common across the cratonic interior. Inversions with narrower frequency bands yield similar patterns, but greater ?t* magnitudes. Even the two standard deviation ?t* magnitude for the widest band is ?2-7 times greater than predicted by global QS tomography or an anelastic olivine thermal model, suggesting that much of the ?t* signal is nonthermal in origin. Nonthermal contributions are further indicated by only a moderate correlation between ?t* and P travel times. Some geographic variations, such as high ?t* in parts of the cratonic interior with high mantle velocities and low heat flow, demonstrate that the influence of temperature is regionally overwhelmed. Transverse spectra are used to investigate the importance of scattering because they would receive no P energy in the absence of 3-D heterogeneity or anisotropy. Transverse to vertical (T/Z) spectral ratios for stations with high ?t* are higher and exhibit steeper increases with frequency compared to T/Z spectra for low ?t* stations. The large magnitude of ?t* estimates and the T/Z spectra are consistent with major contributions to ?t* from scattering. A weak positive correlation between intrinsic attenuation and apparent attenuation due to scattering may contribute to ?t* magnitude and the moderate correlation of ?t* with travel times.
Abstract: The viscosity structure of Earth’s deep mantle affects the thermal evolution of Earth, the ascent of mantle plumes, settling of subducted oceanic lithosphere, and the mixing of compositional heterogeneities in the mantle. Based on a reanalysis of the long-wavelength nonhydrostatic geoid, we infer viscous layering of the mantle using a method that allows us to avoid a priori assumptions about its variation with depth. We detect an increase in viscosity at 800- to 1200-kilometers depth, far greater than the depth of the mineral phase transformations that define the mantle transition zone. The viscosity increase is coincident in depth with regions where seismic tomography has imaged slab stagnation, plume deflection, and changes in large-scale structure and offers a simple explanation of these phenomena.
Journal of Geophysical Research, Vol. 120, 10, pp. 6993-7012.
Canada, Alberta
Geophysics - seismics LVZs
Abstract: Crustal low-velocity zones (LVZs) have been reported in active orogens such as the Himalayas and the Andes but rarely in stable cratonic regions. In this study, we provide compelling evidence for a significant midcrustal LVZ beneath eastern-central Alberta, an integral part of the Precambrian Canadian Shield covered by thick Phanerozoic sedimentary deposits. This 200?km wide, over 10?km thick midcrustal LVZ is well resolved by shear velocity inversions using P-to-S receiver functions from more than 4600 earthquakes. It is generally overlain by a high-velocity upper crust in the depth range of 8-15?km, especially in western-central Alberta, which coincides with the previously documented Winagami reflection sequence. We interpret the LVZ to be of granitic composition, potentially in connection with the crystallization of partially molten crust during the Paleoproterozoic eon. In addition to the Precambrian tectonic history of western Laurentia, which featured plate convergence conducive to crustal melting, our crustal model is further supported by (1) a moderate spatial correlation between the LVZ and heat flow, and (2) shear velocities consistent with that of granite. The well preserved Winagami reflection sequence and the LVZ are potential evidence of distinct episodes of magmatism and crust modification in the Precambrian basement of the Western Canada Sedimentary Basin. The existence of a broad crustal LVZ suggests extensive subduction, orogenesis, and crustal melting during the Precambrian assembly of the North American craton.
Geochemistry, Geophysics, Geosystems: G3, Vol. 16, 11, Nov. pp. 3982-4005.
Mantle
Geophysics - seismic
Abstract: Seismic models provide constraints on the thermal and chemical properties of the cratonic upper mantle. Depth profiles of shear velocity from global and regional studies contain positive velocity gradients in the uppermost mantle and often lack a low-velocity zone, features that are difficult to reconcile with the temperature structures inferred from surface heat flow data and mantle-xenolith thermobarometry. Furthermore, the magnitude and shape of the velocity profiles vary between different studies, impacting the inferences drawn about mantle temperature and composition. In this study, forward modeling is used to identify the suite of one-dimensional shear-velocity profiles that are consistent with phase-velocity observations made for Rayleigh waves traversing Precambrian cratons. Two approaches to the generation of 1-D models are considered. First, depth profiles of shear velocity are predicted from thermal models of the cratonic upper mantle that correspond to a range of assumed values of mantle potential temperature, surface heat flow, and radiogenic heat production in the lithosphere. Second, shear velocity-depth profiles are randomly generated. In both cases, Rayleigh wave phase velocity is calculated from the Earth models, and acceptable models are identified on the basis of comparison to observed phase velocity. The results show that it is difficult but not impossible to find acceptable Earth models that contain a low-velocity zone in the upper mantle and that temperature structures that are consistent with constraints from mantle xenoliths yield phase-velocity predictions lower than observed. For most acceptable randomly generated Earth models, shear velocity merges with the global average at approximately 300 km.
Earth and Planetary Science Letters, Vol. 433, pp. 360-369.
Mantle
Geophysics - seismics
Abstract: In regions of the mantle where multi-phases coexist like at the olivine-wadsleyite-ringwoodite transitions, the stress induced by the seismic waves may drive a mineralogical reaction between the low to high pressure phases, a possible source of dissipation. In such a situation, the amount of attenuation critically depends on the timescale for the phase transformations to reach equilibrium relative to the period of the seismic wave. Here we report synchrotron-based measurements of the kinetics of the olivine to ringwoodite transformation at pressure-temperature conditions of the co-stability loop, for iron-rich olivine compositions. Both microstructural and kinetic data suggest that the transformation rates are controlled by growth processes after the early saturation of nucleation sites along olivine grain boundaries. Transformation-time data show an increase of reaction rates with temperature and iron content, and have been fitted to a rate equation for interface-controlled transformation: G=k0?T?exp?[n?XFa]?exp?[?(?Ha+PV?)/RT]×[1?exp?(?Gr/RT)]G=k0?T?exp?[n?XFa]?exp?[?(?Ha+PV?)/RT]×[1?exp?(?Gr/RT)], where XFaXFa is the fayalite fraction, the exponential factor n=9.7n=9.7, View the MathML sourceln?k0=?9.1 ms?1. View the MathML sourceXFa?1 and ?Ha=199 kJ/mol?Ha=199 kJ/mol, assuming V?=0 cm3/molV?=0 cm3/mol. Including these new kinetic results in a micro-mechanical model of a two-phase loop (Ricard et al., 2009), we predict View the MathML sourceQK?1 and View the MathML sourceQ??1 significantly higher than the PREM values for both body waves and normal modes. This attests that the olivine-wadsleyite transition can significantly contribute to the attenuation of the Earth's mantle transition zone.
Abstract: Using free-board modeling, we examine a vertically-averaged mantle density beneath the Archean-Proterozoic Siberian Craton in the layer from the Moho down to base of the chemical boundary layer (CBL). Two models are tested: in Model 1 the base of the CBL coincides with the LAB, whereas in Model 2 the base of the CBL is at a 180 km depth. The uncertainty of density model is < 0.02 t/m3 or < 0.6% with respect to primitive mantle. The results, calculated at in situ and at room temperature (SPT) conditions, indicate a heterogeneous density structure of the Siberian lithospheric mantle with a strong correlation between mantle density variations and the tectonic setting. Three types of cratonic mantle are recognized from mantle density anomalies. ‘Pristine’ cratonic regions not sampled by kimberlites have the strongest depletion with density deficit of 1.8-3.0% (and SPT density of 3.29-3.33 t/m3 as compared to 3.39 t/m3 of primitive mantle). Cratonic mantle affected by magmatism (including the kimberlite provinces) has a typical density deficit of 1.0-1.5%, indicative of a metasomatic melt-enrichment. Intracratonic sedimentary basins have a high density mantle (3.38-3.40 t/m3 at SPT) which suggests, at least partial, eclogitization. Moderate density anomalies beneath the Tunguska Basin imply that the source of the Siberian LIP lies outside of the Craton. In situ mantle density is used to test the isopycnic condition of the Siberian Craton. Both CBL thickness models indicate significant lateral variations in the isopycnic state, correlated with mantle depletion and best achieved for the Anabar Shield region and other intracratonic domains with a strongly depleted mantle. A comparison of synthetic Mg# for the bulk lithospheric mantle calculated from density with Mg# from petrological studies of peridotite xenoliths from the Siberian kimberlites suggests that melt migration may produce local patches of metasomatic material in the overall depleted mantle.
Geochemistry, Geophysics, Geosystems: G3, Vol. 17, 1, pp. 2-15.
United States
Geophysics - seismics
Abstract: Seismic anisotropy in the upper mantle beneath continental interiors is generally complicated, with contributions from both the lithosphere and the asthenosphere. Previous studies of SKS splitting beneath the eastern United States have yielded evidence for complex and laterally variable anisotropy, but until the recent arrival of the USArray Transportable Array (TA) the station coverage has been sparse. Here we present SKS splitting measurements at TA stations in eastern North America and compare the measured fast directions with indicators such as absolute plate motion, surface geology, and magnetic lineations. We find few correlations between fast directions and absolute plate motion, except in the northeastern U.S. and southern Canada, where some stations exhibit variations in apparent splitting with backazimuth that would suggest multiple layers of anisotropy. A region of the southeastern U.S. is dominated by null SKS arrivals over a range of backazimuths, consistent with previous work. We document a pattern of fast directions parallel to the Appalachian mountain chain, suggesting a contribution from lithospheric deformation associated with Appalachian orogenesis. Overall, our measurements suggest that upper mantle anisotropy beneath the eastern United States is complex, with likely contributions from both asthenospheric and lithospheric anisotropy in many regions.
Physics of the Earth and Planetary Interiors, Vol. 254, pp. 12-24.
Africa, Europe
Geophysics - seismics, anisotropy, boundary
Abstract: The inner core under Africa is thought to be a region where the nature of inner core texture changes: from the strongly anisotropic ‘western’ part of the inner core to the weakly anisotropic, or isotropic ‘eastern’ part of the inner core. Additionally, observations of a difference in isotropic velocity between the two hemispheres have been made. A very large new dataset of simultaneous PKPdf and PKPbc observations, on which differential travel times have been measured, is used to examine the upper 360 km of the inner core under Europe, Africa and the surrounding oceans. Inversion of the differential travel time data for laterally varying inner core anisotropy reveals that inner core anisotropy is stronger under central Africa and the Atlantic Ocean than under the western Indian Ocean. No hemispherical pattern is present in Voigt isotropic velocities, indicating that the variation in anisotropy is due to differing degrees of crystal alignment in the inner core, not material differences. When anisotropy is permitted to change with depth, the upper east-most part of the study region shows weaker anisotropy than the central and western regions. When depth dependence in the inner core is neglected the hemisphere boundary is better represented as a line at 40°E than one at 10°E, however, it is apparent that the variation of anisotropy as a function of depth means that one line of longitude cannot truly separate the more and less anisotropic regions of the inner core. The anisotropy observed in the part of the inner core under Africa which lies in the ‘western’ hemisphere is much weaker than that under central America, showing that the western hemisphere is not uniformly anisotropic. As the region of low anisotropy spans a significant depth extent, it is likely that heterogeneous heat fluxes in the core, which may cause variations in inner core anisotropy, have persisted for several hundred million years.
Physics of the Earth and Planetary Interiors, Vol. 254, pp. 12-24.
Africa, Europe
Geophysics - seismics, anisotropy, boundary
Abstract: The inner core under Africa is thought to be a region where the nature of inner core texture changes: from the strongly anisotropic ‘western’ part of the inner core to the weakly anisotropic, or isotropic ‘eastern’ part of the inner core. Additionally, observations of a difference in isotropic velocity between the two hemispheres have been made. A very large new dataset of simultaneous PKPdf and PKPbc observations, on which differential travel times have been measured, is used to examine the upper 360 km of the inner core under Europe, Africa and the surrounding oceans. Inversion of the differential travel time data for laterally varying inner core anisotropy reveals that inner core anisotropy is stronger under central Africa and the Atlantic Ocean than under the western Indian Ocean. No hemispherical pattern is present in Voigt isotropic velocities, indicating that the variation in anisotropy is due to differing degrees of crystal alignment in the inner core, not material differences. When anisotropy is permitted to change with depth, the upper east-most part of the study region shows weaker anisotropy than the central and western regions. When depth dependence in the inner core is neglected the hemisphere boundary is better represented as a line at 40°E than one at 10°E, however, it is apparent that the variation of anisotropy as a function of depth means that one line of longitude cannot truly separate the more and less anisotropic regions of the inner core. The anisotropy observed in the part of the inner core under Africa which lies in the ‘western’ hemisphere is much weaker than that under central America, showing that the western hemisphere is not uniformly anisotropic. As the region of low anisotropy spans a significant depth extent, it is likely that heterogeneous heat fluxes in the core, which may cause variations in inner core anisotropy, have persisted for several hundred million years.
Geophysical Research Letters, Vol. 43, 5, pp. 1928-1933.
Canada, Ontario
Geophysics - seismics
Abstract: The Superior Province of North America has not experienced major internal deformation for nearly 2.8?Gyr, preserving the Archean crust in its likely original state. We present seismological evidence for a sharp (less than 1?km) crust-mantle boundary beneath three distinct Archean terranes and for a more vertically extensive boundary at sites likely affected by the 1.2-0.9?Ga Grenville orogeny. At all sites crustal thickness is smaller than expected for the primary crust produced by melting under higher mantle potential temperature conditions of Archean time. Reduced thickness and an abrupt contrast in seismic properties at the base of the undisturbed Archean crust are consistent with density sorting and loss of the residues through gravitational instability facilitated by higher temperatures in the upper mantle at the time of formation. Similar sharpness of crust-mantle boundary in disparate Archean terranes suggests that it is a universal feature of the Archean crustal evolution.
Earth and Planetary Science Letters, Vol. 449, pp. 164-175.
United States, Canada
Geophysics - seismics
Abstract: We estimate crustal and uppermost mantle shear velocity structure beneath 30 stations in North America by jointly inverting the high frequency scattered wavefield observed in the P wave coda, together with long period surface wave phase and group dispersion data. Several features distinguish our approach from previous such joint inversions. 1) We apply a cross-convolution method, rather than more standard deconvolution approaches used in receiver function studies, and consider both Love and Rayleigh wave dispersion, allowing us to infer profiles of radial anisotropy. 2) We generate probabilistic 1D radially anisotropic depth profiles across the whole uppermost mantle, down to ?350 km depth. 3) The inverse problem is cast in a trans-dimensional Bayesian formalism, where the number of isotropic and anisotropic layers is treated as unknown, allowing us to obtain models described with the least number of parameters. Results show that the tectonically active region west of the Rocky Mountain Front is marked by a Lithospheric Asthenosphere Boundary and a Lehmann Discontinuity occurring at relatively shallow depths (60-150 km and 100-200 km, respectively), whereas further east, in the stable craton, these discontinuities are deeper (170-200 km and 200-250 km, respectively). In addition, in the stable part of the continent, at least two Mid-Lithospheric Discontinuities are present at intermediate depths, suggesting the existence of strong lithospheric layering, and a mechanism for lithospheric thickening by underplating of additional layers as cratonic age increases. The Moho across the continent as well as mid-crustal discontinuities in the craton are also imaged, in agreement with independent studies.
Abstract: With convergent plate boundaries at some distance, the sources of the lithospheric stress field of the North Atlantic Realm are mainly mantle tractions at the base of the lithosphere, lithospheric density structure and topography. Given this, we estimate horizontal deviatoric stresses using a well-established thin sheet model in a global finite element representation. We adjust the lithospheric thickness and the sub-lithospheric pressure iteratively, comparing modelled in plane stress with the observations of the World Stress Map. We find that an anomalous mantle pressure associated with the Iceland and Azores melt anomalies, as well as topography are able to explain the general pattern of the principle horizontal stress directions. The Iceland melt anomaly overprints the classic ridge push perpendicular to the Mid Atlantic ridge and affects the conjugate passive margins in East Greenland more than in western Scandinavia. The dynamic support of topography shows a distinct maximum of c. 1000 m in Iceland and amounts <150 m along the coast of south-western Norway and 250 -350 m along the coast of East Greenland. Considering that large areas of the North Atlantic Realm have been estimated to be sub-aerial during the time of break-up, two components of dynamic topography seem to have affected the area: a short-lived, which affected a wider area along the rift system and quickly dissipated after break-up, and a more durable in the close vicinity of Iceland. This is consistent with the appearance of a buoyancy anomaly at the base of the North Atlantic lithosphere at or slightly before continental breakup, relatively fast dissipation of the fringes of this, and continued melt generation below Iceland.
Abstract: Seismic images of Earth's interior reveal two massive anomalous zones at the base of the mantle, above the core, where seismic waves travel slowly. The mantle materials that surround these anomalous regions are thought to be composed of cooler rocks associated with downward advection of former oceanic tectonic plates. However, the origin and composition of the anomalous provinces is uncertain. These zones have long been depicted as warmer-than-average mantle materials related to convective upwelling. Yet, they may also be chemically distinct from the surrounding mantle, and potentially partly composed of subducted or primordial material, and have therefore been termed thermochemical piles. From seismic, geochemical and mineral physics data, the emerging view is that these thermochemical piles appear denser than the surrounding mantle materials, are dynamically stable and long-lived, and are shaped by larger-scale mantle flow. Whether remnants of a primordial layer or later accumulations of more-dense materials, the composition of the piles is modified over time by stirring and by chemical reactions with material from the surrounding mantle, underlying core and potentially from volatile elements transported into the deep Earth by subducted plates. Upwelling mantle plumes may originate from the thermochemical piles, so the unusual chemical composition of the piles could be the source of distinct trace-element signatures observed in hotspot lavas.
Physics of the Earth and Planetary Interiors, Vol. 254, pp. 12-24.
Africa, Europe
Anisotropy
Abstract: The inner core under Africa is thought to be a region where the nature of inner core texture changes: from the strongly anisotropic ‘western’ part of the inner core to the weakly anisotropic, or isotropic ‘eastern’ part of the inner core. Additionally, observations of a difference in isotropic velocity between the two hemispheres have been made. A very large new dataset of simultaneous PKPdf and PKPbc observations, on which differential travel times have been measured, is used to examine the upper 360 km of the inner core under Europe, Africa and the surrounding oceans. Inversion of the differential travel time data for laterally varying inner core anisotropy reveals that inner core anisotropy is stronger under central Africa and the Atlantic Ocean than under the western Indian Ocean. No hemispherical pattern is present in Voigt isotropic velocities, indicating that the variation in anisotropy is due to differing degrees of crystal alignment in the inner core, not material differences. When anisotropy is permitted to change with depth, the upper east-most part of the study region shows weaker anisotropy than the central and western regions. When depth dependence in the inner core is neglected the hemisphere boundary is better represented as a line at 40°E than one at 10°E, however, it is apparent that the variation of anisotropy as a function of depth means that one line of longitude cannot truly separate the more and less anisotropic regions of the inner core. The anisotropy observed in the part of the inner core under Africa which lies in the ‘western’ hemisphere is much weaker than that under central America, showing that the western hemisphere is not uniformly anisotropic. As the region of low anisotropy spans a significant depth extent, it is likely that heterogeneous heat fluxes in the core, which may cause variations in inner core anisotropy, have persisted for several hundred million years.
Geochemistry, Geophysics, Geosystems: G3, Vol. 17, 2, pp. 5056-5077.
Mantle
Geophysics - seismics
Abstract: The large low shear-wave velocity provinces (LLSVP) are thermochemical anomalies in the deep Earth's mantle, thousands of km wide and ?1800 km high. This study explores the hypothesis that the LLSVPs are compositionally subdivided into two domains: a primordial bottom domain near the core-mantle boundary and a basaltic shallow domain that extends from 1100 to 2300 km depth. This hypothesis reconciles published observations in that it predicts that the two domains have different physical properties (bulk-sound versus shear-wave speed versus density anomalies), the transition in seismic velocities separating them is abrupt, and both domains remain seismically distinct from the ambient mantle. We here report underside reflections from the top of the LLSVP shallow domain, supporting a compositional origin. By exploring a suite of two-dimensional geodynamic models, we constrain the conditions under which well-separated "double-layered" piles with realistic geometry can persist for billions of years. Results show that long-term separation requires density differences of ?100 kg/m3 between LLSVP materials, providing a constraint for origin and composition. The models further predict short-lived "secondary" plumelets to rise from LLSVP roofs and to entrain basaltic material that has evolved in the lower mantle. Long-lived, vigorous "primary" plumes instead rise from LLSVP margins and entrain a mix of materials, including small fractions of primordial material. These predictions are consistent with the locations of hot spots relative to LLSVPs, and address the geochemical and geochronological record of (oceanic) hot spot volcanism. The study of large-scale heterogeneity within LLSVPs has important implications for our understanding of the evolution and composition of the mantle.
Abstract: Continental lithosphere houses the oldest and thickest regions of the Earth's surface. Locked within this deep and ancient rock record lies invaluable information about the dynamics that has shaped and continue to shape the planet. Much of that history has been dominated by the forces of plate tectonics which has repeatedly assembled super continents together and torn them apart - the Wilson Cycle. While the younger regions of continental lithosphere have been subject to deformation driven by plate tectonics, it is less clear whether the ancient, stable cores formed and evolved from similar processes. New insight into continental formation and evolution has come from remarkable views of deeper lithospheric structure using enhanced seismic imaging techniques and the increase in large volumes of broadband data. Some of the most compelling observations are that the continental lithosphere has a broad range in thicknesses (< 100 to > 300 km), has complex internal structure, and that the thickest portion appears to be riddled with seismic discontinuities at depths between ~ 80 and ~ 130 km. These internal structural features have been interpreted as remnants of lithospheric formation during Earth's early history. If they are remnants, then we can attempt to investigate the structure present in the deep lithosphere to piece together information about early Earth dynamics much as is done closer to the surface. This would help delineate between the differing models describing the dynamics of craton formation, particularly whether they formed in the era of modern plate tectonics, a transitional mobile-lid tectonic regime, or are the last fragments of an early, stagnant-lid planet. Our review paper (re)introduces readers to the conceptual definitions of the lithosphere and the complex nature of the upper boundary layer, then moves on to discuss techniques and recent seismological observations of the continental lithosphere. We then review geodynamic models and hypotheses for the formation of the continental lithosphere through time and implications for the formation and preservation of deep structure. These are contrasted with the dynamical picture of modern day continental growth during lateral accretion of juvenile crust with reference to examples from the Australian Tasmanides and the Alaskan accretionary margin.
European Geoscience Union Conference, April, Vienna, Vol. 19, EGU2017-2344. 1p. Abstract
Mantle
Geophysics - seismics
Abstract: The piezoelectric and seismo-electrokinetic phenomena are manifested by electrical and electromagnetic processes that occur in rocks under the influence of elastic oscillations triggered by shots or mechanical impacts.
Geophysical Research Letters, Vol. 43, 23, pp. 12,054-62.
United States
Geophysics - seismics
Abstract: To constrain crustal structures in the southern Appalachians and the suture zone with the Gondwanan-affinity Suwannee terrane, we applied the 2-D generalized Radon transform wavefield migration method to the scattered incident P wavefield recorded by the EarthScope Southeastern Suture of the Appalachian Margin Experiment and adjacent Transportable Array stations. We resolve the root of thickened crust beneath the high topography of the Blue Ridge Mountains and estimate its density contrast with the mantle to be only 104?±?20?kg/m3. A weak velocity contrast across the crustal root Moho is observed and may be related to an ongoing crustal delamination event, possibly contributing to local tectonic rejuvenation. Beneath the Suwannee terrane, we confirm prior observations of a gently south-southeastward dipping crustal suture, indicating the terminal collision of Laurentia and Gondwana involved several hundred kilometers of overthrusting.
Abstract: Progress in seismic methodology and ambitious large-scale seismic projects are enabling high-resolution imaging of the continental crust. The ability to constrain interpretations of crustal seismic data is based on laboratory measurements on rock samples and calculations of seismic properties. Seismic velocity calculations and their directional dependence are based on the rock micro fabric, which consists of mineral aggregate properties including crystallographic preferred orientation (CPO), grain shape and distribution, grain boundary distribution, and misorientation within grains. Single mineral elastic constants and density are crucial for predicting seismic velocities, preferably at conditions that span the crust. However, high temperature and pressure properties are not as common as elastic constants at standard temperature and pressure (STP) at atmospheric conditions. Continental crust has a very diverse mineralogy, however a select number appear to dominate seismic properties because of their high volume fraction contribution. Calculations of micro fabric-based seismic properties and anisotropy are performed with averaging methods that in their simplest form takes into account the CPO and modal mineral composition. More complex methods can take into account other microstructural characteristics, including the grain shape and distribution of mineral grains, and cracks and pores. A challenge for the geophysics and rock physics communities is the separation of intrinsic factors affecting seismic anisotropy, due to properties of crystals within a rock and apparent sources due to extrinsic factors like cracks, fractures and alteration. This is of particular importance when trying to deduce the state of crustal composition and deformation from seismic parameters.
Abstract: The composition of the lower mantle—comprising 56% of Earth’s volume—remains poorly constrained. Among the major elements, Mg/Si ratios ranging from ~0.9-1.1, such as in rocky Solar-System building blocks (or chondrites), to ~1.2-1.3, such as in upper-mantle rocks (or pyrolite), have been proposed. Geophysical evidence for subducted lithosphere deep in the mantle has been interpreted in terms of efficient mixing, and thus homogenous Mg/Si across most of the mantle. However, previous models did not consider the effects of variable Mg/Si on the viscosity and mixing efficiency of lower-mantle rocks. Here, we use geodynamic models to show that large-scale heterogeneity associated with a 20-fold change in viscosity, such as due to the dominance of intrinsically strong (Mg, Fe)SiO3-bridgmanite in low-Mg/Si domains, is sufficient to prevent efficient mantle mixing, even on large scales. Models predict that intrinsically strong domains stabilize mantle convection patterns, and coherently persist at depths of about 1,000-2,200?km up to the present-day, separated by relatively narrow up-/downwelling conduits of pyrolitic material. The stable manifestation of such bridgmanite-enriched ancient mantle structures (BEAMS) may reconcile the geographical fixity of deep-rooted mantle upwelling centres, and geophysical changes in seismic-tomography patterns, radial viscosity, rising plumes and sinking slabs near 1,000?km depth. Moreover, these ancient structures may provide a reservoir to host primordial geochemical signatures.
Abstract: Hints from seismic tomography and geochemistry indicate that Earth's mantle is heterogeneous at large scale. Numerical simulations of mantle convection show that, if it started enriched in silicates, the lower mantle may remain unmixed today.
Earth and Planetary Science Letters, Vol. 470, pp. 54-63.
Mantle
geophysics - seismic
Abstract: Seismic tomography resolves anomalies interpreted as oceanic lithosphere subducted deep into Earth's lower mantle. However, the fate of the compositionally distinct oceanic crust that is part of the lithosphere is poorly constrained but provides important constraints on mixing processes and the recycling process in the deep Earth. We present high-resolution seismic array analyses of anomalous P-waves sampling the deep mantle, and deterministically locate heterogeneities in the lowermost 300 km of the mantle. Spectral analysis indicates that the dominant scale length of the heterogeneity is 4 to 7 km. The heterogeneity distribution varies laterally and radially and heterogeneities are more abundant near the margins of the lowermost mantle Large Low Velocity Provinces (LLVPs), consistent with mantle convection simulations that show elevated accumulations of deeply advected crustal material near the boundaries of thermo-chemical piles. The size and distribution of the observed heterogeneities is consistent with that expected for subducted oceanic crust. These results thus suggest the deep mantle contains an imprint of continued subduction of oceanic crust, stirred by mantle convection and modulated by long lasting thermo-chemical structures. The preferred location of the heterogeneity in the lowermost mantle is consistent with a thermo-chemical origin of the LLVPs. Our observations relate to the mixing behaviour of small length-scale heterogeneity in the deep Earth and indicate that compositional heterogeneities from the subduction process can survive for extended times in the lowermost mantle.
Abstract: Progress in seismic methodology and ambitious large-scale seismic projects are enabling high-resolution imaging of the continental crust. The ability to constrain interpretations of crustal seismic data is based on laboratory measurements on rock samples and calculations of seismic properties. Seismic velocity calculations and their directional dependence are based on the rock micro fabric, which consists of mineral aggregate properties including crystallographic preferred orientation (CPO), grain shape and distribution, grain boundary distribution, and misorientation within grains. Single mineral elastic constants and density are crucial for predicting seismic velocities, preferably at conditions that span the crust. However, high temperature and pressure properties are not as common as elastic constants at standard temperature and pressure (STP) at atmospheric conditions. Continental crust has a very diverse mineralogy, however a select number appear to dominate seismic properties because of their high volume fraction contribution. Calculations of micro fabric-based seismic properties and anisotropy are performed with averaging methods that in their simplest form takes into account the CPO and modal mineral composition. More complex methods can take into account other microstructural characteristics, including the grain shape and distribution of mineral grains, and cracks and pores. A challenge for the geophysics and rock physics communities is the separation of intrinsic factors affecting seismic anisotropy, due to properties of crystals within a rock and apparent sources due to extrinsic factors like cracks, fractures and alteration. This is of particular importance when trying to deduce the state of crustal composition and deformation from seismic parameters.
Earth and Planetary Science Letters, Vol. 472, pp. 164-173.
Mantle
geophysics, seismic
Abstract: Subduction zones are likely a major source of compositional heterogeneities in the mantle, which may preserve a record of the subduction history and mantle convection processes. The fine-scale structure associated with mantle heterogeneities can be studied using the scattered seismic wavefield that arrives as coda to or as energy preceding many body wave arrivals. In this study we analyse precursors to PP by creating stacks recorded at globally distributed stations. We create stacks aligned on the PP arrival in 5° distance bins (with range 70–120°) from 600 earthquakes recorded at 193 stations stacking a total of 7320 seismic records. As the energy trailing the direct P arrival, the P coda, interferes with the PP precursors, we suppress the P coda by subtracting a best fitting exponential curve to this energy. The resultant stacks show that PP precursors related to scattering from heterogeneities in the mantle are present for all distances. Lateral variations are explored by producing two regional stacks across the Atlantic and Pacific hemispheres, but we find only negligible differences in the precursory signature between these two regions. The similarity of these two regions suggests that well mixed subducted material can survive at upper and mid-mantle depth. To describe the scattered wavefield in the mantle, we compare the global stacks to synthetic seismograms generated using a Monte Carlo phonon scattering technique. We propose a best-fitting layered heterogeneity model, BRT2017, characterised by a three layer mantle with a background heterogeneity strength (?=0.8%) and a depth-interval of increased heterogeneity strength (?=1%) between 1000 km and 1800 km. The scalelength of heterogeneity is found to be 8 km throughout the mantle. Since mantle heterogeneity of 8 km scale may be linked to subducted oceanic crust, the detection of increased heterogeneity at mid-mantle depths could be associated with stalled slabs due to increases in viscosity, supporting recent observations of mantle viscosity increases due to the iron spin transition at depths of ?1000 km.
Geophysical Research Letters, Vol. 42, 10. pp. 3829-3846.
Australia
geophysics, seismic
Abstract: A new 3-D S wave speed model for the Australian region derived from multimode surface waves allows us to examine the nature of the lithosphere-asthenosphere transition (LAT) and its relation to radial anisotropy. In eastern Phanerozoic Australia the estimated depths of the LAT tie well with those from receiver functions. However, in the Archean and Proterozoic lithosphere in western and central Australia, the LAT derived from the surface wave model is generally much deeper than the discontinuities recognized from receiver functions and shows a smooth transition. There is significant radial anisotropy (SH?>?SV) in the upper lithosphere as well as in the LAT and the underlying asthenosphere. Strong anisotropy in the asthenosphere reflects the effects of present shear flow in the mantle beneath the continent. The lateral variation of lithospheric anisotropy correlates well with the suture zones between cratonic blocks, representing frozen anisotropy associated with the ancient tectonics of Australia.
Abstract: It is now accepted that significant volumes of pyroxenites are generated in the subduction factory and remain trapped in the mantle. In ophiolites and orogenic massifs the geometry of pyroxenite layers and their relationships with the host peridotite can be observed directly. Since a large part of what is known about the upper mantle structure is derived from the analysis of seismic waves, it is crucial to integrate pyroxenites in the interpretations. We modeled the seismic properties of a peridotitic mantle rich in pyroxenite layers in order to determine the impact of layering on the seimsic properties. To do so, EBSD data on deformed and undeformed pyroxenites from the Cabo Ortegal complex (Spain) and the Trinity ophiolite (California, USA) respectively are combined with either A or B-type olivine fabrics in order to model a realistic pyroxenite-rich upper mantle. Consideration of pyroxeniterich domains within the host mantle wall rock is incorporated in the calculations using the Schoenberg and Muir group theory [1]. This quantification reveals the complex dependence of the seismic signal on the deformational state and relative abundance of each mineral phase. The incorporation of pyroxenites properties into geophysical interpretations in understanding the lithospheric structure of subduction zones will lead to more geologically realistic models.
Abstract: The Earth’s lowermost mantle large low velocity provinces are accompanied by small-scale ultralow velocity zones in localized regions on the core-mantle boundary. Large low velocity provinces are hypothesized to be caused by large-scale compositional heterogeneity (i.e., thermochemical piles). The origin of ultralow velocity zones, however, remains elusive. Here we perform three-dimensional geodynamical calculations to show that the current locations and shapes of ultralow velocity zones are related to their cause. We find that the hottest lowermost mantle regions are commonly located well within the interiors of thermochemical piles. In contrast, accumulations of ultradense compositionally distinct material occur as discontinuous patches along the margins of thermochemical piles and have asymmetrical cross-sectional shape. Furthermore, the lateral morphology of these patches provides insight into mantle flow directions and long-term stability. The global distribution and large variations of morphology of ultralow velocity zones validate a compositionally distinct origin for most ultralow velocity zones.
Geophysics Research Letters, Vol. 44, 5, pp. 2148-2158.
Mantle
geophysics - seismics
Abstract: We explore the effect of Mg/Fe substitution on the sound velocities of iron-rich (Mg1 ? xFex)O, where x = 0.84, 0.94, and 1.0. Sound velocities were determined using nuclear resonance inelastic X-ray scattering as a function of pressure, approaching those of the lowermost mantle. The systematics of cation substitution in the Fe-rich limit has the potential to play an important role in the interpretation of seismic observations of the core-mantle boundary. By determining a relationship between sound velocity, density, and composition of (Mg,Fe)O, this study explores the potential constraints on ultralow-velocity zones at the core-mantle boundary.
Earth and Planetary Science Letters, Vol. 477, pp.
Mantle
geophysics - seismics
Abstract: We exploit conversions between P and S waves for large-scale, high-resolution imaging of the mantle transition zone beneath Northwest Pacific and the margin of Eastern Asia. We find pervasive reflectivity concentrated in two bands with apparent wave-speed reduction of ?2% to ?4% about 50 km thick at the top of the transition zone and 100 km thick at the bottom. This negative reflectivity associated with the scattered-waves at depth is interpreted jointly with larger-scale mantle tomographic images, and is shown to delineate the stagnant portions of the subducted Pacific plate in the transition zone, with largely positive shear-wave velocity contrasts. The upper reflectivity zone connects to broad low-velocity regions below major intra-plate volcanoes, whereas the lower zone coincides locally with the occurrence of deep-focus earthquakes along the East Asia margin. Similar reflectivity is found in Pacific Northwest of the USA. We demonstrate that the thermal signature of plates alone is not sufficient to explain such features. Alternative explanations for these reflective zones include kinetic effects on olivine phase transitions (meta-stability), compositional heterogeneities within and above stagnant plates, complex wave-propagation effects in the heterogeneous slab structure, or a combination of such factors. We speculate that part of the negative reflectivity is the signature of compositional heterogeneities, as revealed by numerous other studies of seismic scattering throughout the mantle, and that such features could be widespread across the globe.
Earth and Planetary Science Letters, Vol. 479, pp. 170-178.
Canada, Nunavut, Baffin Island
geophysics - seismics
Abstract: The northern Hudson Bay region in Canada comprises several Archean cratonic nuclei, assembled by a number of Paleoproterozoic orogenies including the Trans-Hudson Orogen (THO) and the Rinkian-Nagssugtoqidian Orogen. Recent debate has focused on the extent to which these orogens have modern analogues such as the Himalayan-Karakoram-Tibet Orogen. Further, the structure of the lithospheric mantle beneath the Hudson Strait and southern Baffin Island is potentially indicative of Paleoproterozoic underthrusting of the Superior plate beneath the Churchill collage. Also in question is whether the Laurentian cratonic root is stratified, with a fast, depleted, Archean core underlain by a slower, younger, thermally-accreted layer. Plate-scale process that create structures such as these are expected to manifest as measurable fossil seismic anisotropic fabrics. We investigate these problems via shear wave splitting, and present the most comprehensive study to date of mantle seismic anisotropy in northern Laurentia. Strong evidence is presented for multiple layers of anisotropy beneath Archean zones, consistent with the episodic development model of stratified cratonic keels. We also show that southern Baffin Island is underlain by dipping anisotropic fabric, where underthrusting of the Superior plate beneath the Churchill has previously been interpreted. This provides direct evidence of subduction-related deformation at 1.8 Ga, implying that the THO developed with modern plate-tectonic style interactions.
Earth and Planetary Science Letters, Vol. 481, pp. 223-235.
United States, Illinois, Indiana, Kentucky
geophysics - seismics Reelfoot Rift
Abstract: Seismic discontinuities between the Moho and the inferred lithosphere-asthenosphere boundary (LAB) are known as mid-lithospheric discontinuities (MLDs) and have been ascribed to a variety of phenomena that are critical to understanding lithospheric growth and evolution. In this study, we used S-to-P converted waves recorded by the USArray Transportable Array and the OIINK (Ozarks-Illinois-Indiana-Kentucky) Flexible Array to investigate lithospheric structure beneath the central U.S. This region, a portion of North America's cratonic platform, provides an opportunity to explore how terrane accretion, cratonization, and subsequent rifting may have influenced lithospheric structure. The 3D common conversion point (CCP) volume produced by stacking back-projected Sp receiver functions reveals a general absence of negative converted phases at the depths of the LAB across much of the central U.S. This observation suggests a gradual velocity decrease between the lithosphere and asthenosphere. Within the lithosphere, the CCP stacks display negative arrivals at depths between 65 km and 125 km. We interpret these as MLDs resulting from the top of a layer of crystallized melts (sill-like igneous intrusions) or otherwise chemically modified lithosphere that is enriched in water and/or hydrous minerals. Chemical modification in this manner would cause a weak layer in the lithosphere that marks the MLDs. The depth and amplitude of negative MLD phases vary significantly both within and between the physiographic provinces of the midcontinent. Double, or overlapping, MLDs can be seen along Precambrian terrane boundaries and appear to result from stacked or imbricated lithospheric blocks. A prominent negative Sp phase can be clearly identified at 80 km depth within the Reelfoot Rift. This arrival aligns with the top of a zone of low shear-wave velocities, which suggests that it marks an unusually shallow seismic LAB for the midcontinent. This boundary would correspond to the top of a region of mechanically and chemically rejuvenated mantle that was likely emplaced during late Precambrian/early Cambrian rifting. These observations suggest that the lithospheric structure beneath the Reelfoot Rift may be an example of a global phenomenon in which MLDs act as weak zones that facilitate the removal of cratonic lithosphere that lies beneath.
Abstract: The existence of undulations of the geoid, gravity and bathymetry in ocean basins, as well as anomalies in heat flow, point to the existence of small scale convection beneath tectonic plates. The instabilities that could develop at the base of the lithosphere are sufficiently small scale (< 500 km) that they remain mostly elusive from seismic detection. We take advantage of 3D spherical numerical geodynamic models displaying plate-like behavior to study the interaction between large-scale flow and small-scale convection. We find that finger-shaped instabilities develop at seafloor ages > 60 Ma. They form networks that are shaped by the plate evolution, slabs, plumes and the geometry of continental boundaries. Plumes impacting the boundary layer from below have a particular influence through rejuvenating the thermal lithosphere. They create a wake in which new instabilities form downstream. These wakes form channels that are about 1000 km wide, and thus are possibly detectable by seismic tomography. Beneath fast plates, cold sinking instabilities are tilted in the direction opposite to plate motion, while they sink vertically for slow plates. These instabilities are too small to be detected by usual seismic methods, since they are about 200 km in lateral scale. However, this preferred orientation of instabilities below fast plates could produce a pattern of large-scale azimuthal anisotropy consistent with both plate motions and the large scale organisation of azimuthal anisotropy obtained from recent surface wave models.
Journal of Asian Earth Sciences, in press available, 11p.
India
geophysics - seismics
Abstract: The time-bound crustal evolution and subsequent deformation of the Cuddapah basin, Nellore Schist Belt and Eastern Ghats terrain of Eastern Dharwar Craton, which have undergone sustained geodynamic upheavals since almost 2.0 billion years, remain enigmatic. An attempt is made here to integrate newly available potential field data and other geophysical anomalies with deep seismic structure, to examine the generative mechanism of major crustal features, associated with this sector. Our study indicates that the initial extent of the Cuddapah basin sedimentation may have been much larger, extending by almost 50-60?km west of Tadipatri during Paleoproterozoic period, which subsequently shrank due to massive erosion following thermal uplift, caused by SW Cuddapah mantle plume. Below this region, crust is still quite warm with Moho temperatures exceeding 500?°C. Similarly, Nallamalai Fold Belt rocks, bounded by two major faults and extremely low gravity, may have occupied a large terrain in western Cuddapah basin also, before their abrasion. No geophysical signatures of thrusting are presently seen below this region, and thus it could not be an alien terrain either. In contrast, Nellore Schist Belt is associated with strikingly high positive gravity, possibly caused by a conspicuous horst structure and up dipping mafic crustal layers underneath, that resulted due to India-east Antarctica collision after the cessation of prolonged subduction (1.6-0.95?Ga). Further, the crustal seismic and gravity signatures would confirm presence of a totally distinct geological terrain east of the Cuddapah basin, but the trace of Eastern Ghats Belt is all together missing. Instead, all the geophysical signatures, point out to presence of a Proterozoic sedimentary terrain, east of Nellore Schist Belt. It is likely that the extent of Prorerozoic sedimentation was much larger than thought today. In addition, presence of a seismically detected Gondwana basin over Nellore Schist Belt, apart from some recently discovered similar subsurface Gondwana occurrences in intracratonic parts, would indicate that Dharwar Craton was rifting even during Gondwana period, thereby challenging the long held view of cratonic stability.
Earth and Planetary Science Letters, Vol. 458, 1, pp.405-417.
Africa, Madagascar
geophysics - seismics
Abstract: The crust and upper mantle of the Madagascar continental fragment remained largely unexplored until a series of recent broadband seismic experiments. An island-wide deployment of broadband seismic instruments has allowed the first study of phase velocity variations, derived from surface waves, across the entire island. Late Cenozoic alkaline intraplate volcanism has occurred in three separate regions of Madagascar (north, central and southwest), with the north and central volcanism active until <1 Ma, but the sources of which remains uncertain. Combined analysis of three complementary surface wave methods (ambient noise, Rayleigh wave cross-correlations, and two-plane-wave) illuminate the upper mantle down to depths of 150 km. The phase-velocity measurements from the three methods for periods of 8-182 s are combined at each node and interpolated to generate the first 3-D shear-velocity model for sub-Madagascar velocity structure. Shallow (upper 10 km) low-shear-velocity regions correlate well with sedimentary basins along the west coast. Upper mantle low-shear-velocity zones that extend to at least 150 km deep underlie the north and central regions of recent alkali magmatism. These anomalies appear distinct at depths <100 km, suggesting that any connection between the zones lies at depths greater than the resolution of surface-wave tomography. An additional low-shear velocity anomaly is also identified at depths 50-150 km beneath the southwest region of intraplate volcanism. We interpret these three low-velocity regions as upwelling asthenosphere beneath the island, producing high-elevation topography and relatively low-volume magmatism.
Abstract: Lateral variations of seismic wave speeds and attenuation (dissipation of strain energy) in the Earth’s upper mantle have the potential to map key characteristics such as temperature, major-element composition, melt fraction and water content1,2,3. The inversion of these data into meaningful representations of physical properties requires a robust understanding of the micromechanical processes that affect the propagation of seismic waves2,3. Structurally bound water (hydroxyl) is believed to affect seismic properties2,3 but this has yet to be experimentally quantified. Here we present a comprehensive low-frequency forced-oscillation assessment of the seismic properties of olivine as a function of water content within the under-saturated regime that is relevant to the Earth’s interior. Our results demonstrate that wave speeds and attenuation are in fact strikingly insensitive to water content. Rather, the redox conditions imposed by the choice of metal sleeving, and the associated defect chemistry, appear to have a substantial influence on the seismic properties. These findings suggest that elevated water contents are not responsible for low-velocity or high-attenuation structures in the upper mantle. Instead, the high attenuation observed in hydrous and oxidized regions of the upper mantle (such as above subduction zones) may reflect the prevailing oxygen fugacity. In addition, these data provide no support for the hypothesis whereby a sharp lithosphere-asthenosphere boundary is explained by enhanced grain boundary sliding in the presence of water.
Abstract: Seismic waves that propagate through a layer of Earth’s upper mantle are highly attenuated. Contrary to general thinking, this attenuation seems to be strongly affected by oxidation conditions, rather than by water content.
Journal of Asian Earth Sciences, Vol. 157, pp. 302-316.
India
geophysics - seismics
Abstract: The time-bound crustal evolution and subsequent deformation of the Cuddapah basin, Nellore Schist Belt and Eastern Ghats terrain of Eastern Dharwar Craton, which have undergone sustained geodynamic upheavals since almost 2.0 billion years, remain enigmatic. An attempt is made here to integrate newly available potential field data and other geophysical anomalies with deep seismic structure, to examine the generative mechanism of major crustal features, associated with this sector. Our study indicates that the initial extent of the Cuddapah basin sedimentation may have been much larger, extending by almost 50-60?km west of Tadipatri during Paleoproterozoic period, which subsequently shrank due to massive erosion following thermal uplift, caused by SW Cuddapah mantle plume. Below this region, crust is still quite warm with Moho temperatures exceeding 500?°C. Similarly, Nallamalai Fold Belt rocks, bounded by two major faults and extremely low gravity, may have occupied a large terrain in western Cuddapah basin also, before their abrasion. No geophysical signatures of thrusting are presently seen below this region, and thus it could not be an alien terrain either. In contrast, Nellore Schist Belt is associated with strikingly high positive gravity, possibly caused by a conspicuous horst structure and up dipping mafic crustal layers underneath, that resulted due to India-east Antarctica collision after the cessation of prolonged subduction (1.6-0.95?Ga). Further, the crustal seismic and gravity signatures would confirm presence of a totally distinct geological terrain east of the Cuddapah basin, but the trace of Eastern Ghats Belt is all together missing. Instead, all the geophysical signatures, point out to presence of a Proterozoic sedimentary terrain, east of Nellore Schist Belt. It is likely that the extent of Prorerozoic sedimentation was much larger than thought today. In addition, presence of a seismically detected Gondwana basin over Nellore Schist Belt, apart from some recently discovered similar subsurface Gondwana occurrences in intracratonic parts, would indicate that Dharwar Craton was rifting even during Gondwana period, thereby challenging the long held view of cratonic stability.
Journal of Geophysical Research, Vol. 123, 1. pp. 384-400.
Mantle
geophysics - seismic
Abstract: The strongest evidence to support the classical plume hypothesis comes from seismic imaging of the mantle beneath hot spots. However, imaging results are often ambiguous and it is questionable whether narrow plume tails can be detected by present?day seismological techniques. Here we carry out synthetic tomography experiments based on spectral element method simulations of seismic waves with period T > 10 s propagating through geodynamically derived plume structures. We vary the source?receiver geometry in order to explore the conditions under which lower mantle plume tails may be detected seismically. We determine that wide?aperture (4,000-6,000 km) networks with dense station coverage (<100-200 km station spacing) are necessary to image narrow (<500 km wide) thermal plume tails. We find that if uncertainties on traveltime measurements exceed delay times imparted by plume tails (typically <1 s), the plume tails are concealed in seismic images. Vertically propagating SKS waves enhance plume tail recovery but lack vertical resolution in regions that are not independently constrained by direct S paths. We demonstrate how vertical smearing of an upper mantle low?velocity anomaly can appear as a plume originating in the deep mantle. Our results are useful for interpreting previous plume imaging experiments and guide the design of future experiments.
Abstract: Investigation of large scale suture zones in old continental interiors offers insights into the evolution of continents. The Dharwar Craton (DC) and the Southern Granulite Terrain(SGT) of the Indian shield represent large segments of Precambrian middle to lower crust and preserve a geological record spanning from Mesoarchean to Cambrian. This study illuminates the deep structure of the Palghat-Cauvery Shear Zone System (PCSS) and the Palghat-Cauvery Suture Zone (PCSZ) that comprise crustal-scale structures related to multiple episodes of orogeny, crust formation and reworking. We utilize here 3202 high quality P-receiver functions computed using new data from a 23 station seismic network operated by us. Results show a thick (>38?km) mafic (Poisson's ratio >0.25) crust beneath the SGT. The change in crustal thickness is gradual, with a shallower Moho towards the south of PCSZ. We found little evidence for drastic changes in crustal thickness across prominent shear zones like the PCSZ and Moyar-Bhavani. Few seismic stations located along these boundaries have shown evidence for dipping reflectors around 8-20?km depth, with strikes matching well with the trends of surface geological sutures. We opine that these suture zones do not show indications of a terrane boundary. However, a drastic change in the crustal thickness is observed around the prograde metamorphic transition zone or broadly, the "Fermor line", which separates rocks of Chanockitic (Orthopyroxene bearing granitoid) and non-Charnockitic (Orthopyroxene-free granitoid) mineral assemblage, further north beneath the DC. We suggest that thicknening of crust north of Moyar-Attur Shear Zone (MASZ) and around Fermor line is related to subduction processes operative during the Precambrian.
Earth and Planetary Science Letters, Vol. 492, pp. 197-205.
United States, Colorado
geophysics - seismic
Abstract: We analyze teleseismic P-to-S conversions for high-resolution imaging of the mantle transition zone beneath the Colorado Rocky Mountains using data from a dense PASSCAL seismic broadband deployment. A total of 6,021 P-to-S converted receiver functions are constructed using a multi-channel minimum-phase deconvolution method and migrated using the common converted point technique with the 3-D teleseismic P- and S-wave tomography models of Schmandt and Humphreys (2010). The image finds that the average depths of the 410-km discontinuity (the 410) and 660-km discontinuity (the 660) at and respectively. The peak-to-peak topography of both discontinuities is 33 km and 27 km respectively. Additionally, prominent negative polarity phases are imaged both above and below the 410. To quantify the mean properties of the low-velocity layers about 410 km, we utilize double gradient layer models parameterization to fit the mean receiver function waveform. This waveform fitting is accomplished as a grid-search using anelastic synthetic seismograms. The best-fitting model reveals that the olivine-wadsleyite phase transformation width is 21 km, which is significantly larger than anhydrous mineral physics prediction (4-10 km) (Smyth and Frost, 2002). The findings of a wide olivine-wadsleyite phase transformation and the negative polarity phases above and below the 410, suggest that the mantle, at least in the 350-450 km depth range, is significantly hydrated. Furthermore, a conspicuous negative polarity phase below the 660 is imaged in high velocity region, we speculate the low velocity layer is due to dehydration flux melting in an area of convective downwelling. Our interpretation of these results, in tandem with the tomographic image of a Farallon slab segment at 800 km beneath the region (Schmandt and Humphreys, 2010), is that hydrous and upwelling mantle contributes to the high-standing Colorado Rocky Mountains.
Abstract: The Capricorn Orogen in central Western Australia played important roles in initializing and finalizing the West Australian craton. Surface geological mapping and isotopic studies show that the crust has recorded over a billion years of tectonic history spanning from its crustal formation in the Archean to episodes of tectono?thermal events during the Proterozoic cratonization processes. The region therefore provides us with an ideal laboratory to characterize the seismic signature associated with tectonic processes. We constructed a crustal shear?wave velocity model of the core region of the orogen, the Glenburgh Terrane and its north boundary, by inverting the array group velocity dispersion data measured from a high density temporary array. A modified Bayesian Transdimensional tomography technique, which incorporates a smooth?varying regional reference velocity model and Moho topography, was used to invert for the crustal velocity variations. The inverted velocity model adds great detail to the intra?crustal structure, and provides complementary seismic velocity information to refine the regional tectonic processes. Distinct patterns in the velocity structure support that the Glenburgh Terrane is a microcontinent originated in the Archean, and favor the operation of Paleoproterozoic subduction/accretion leading to the 2.2 Ga Ophthalmian orogeny that initiated the assembly of the West Australian craton.
Abstract: On the early Earth, oceanic plateaux similar to present-day Iceland are thought to have evolved into less dense microcontinents as they thickened by continued melt intrusion and crustal fractionation. These earliest continents may have been so weak on a hotter Earth that they collapsed laterally in response to thickening by further magmatic growth or tectonic imbrication. This continental spreading is likely to have resulted in the development of pervasive ductile strain fabrics in the deeper crust, which, if preserved, could generate seismic reflections. Here we present seismic images from the ancient core of the Archaean Yilgarn Craton of Australia that reveal shallowly dipping to horizontal reflections that pervade the middle and lower crust. We interpret these reflective fabrics as the result of widespread lateral crustal flow during the late stage of craton evolution approximately 2.66 to 2.61?billion years ago, which coincided with the widespread intrusion of high-temperature crustal melts, as thickened early continental crust collapsed. The consequent subsidence of large regions of the upper crust, including volcanic and sedimentary greenstone rocks, in the hanging walls of listric mid-lower crustal ductile flow fabrics caused these rocks to drop beneath the granitic melts rising towards the surface, and did not involve Rayleigh-Taylor instabilities within a mostly mobile crust.
Geochemistry, Geophysics, Geosystems, Vol. 19, 5, pp. 1464-1483.
Mantle
geophysics - seismic
Abstract: SubMachine is a collection of web-based tools for the interactive visualisation, analysis, and quantitative comparison of global-scale, volumetric (3-D) data sets of the subsurface, with supporting tools for interacting with other, complementary models and data sets as listed below. In short, SubMachine is a computational engine (Machine) to visualize models and datasets of the sub-surface (Sub).
Abstract: Investigation of large scale suture zones in old continental interiors offers insights into the evolution of continents. The Dharwar Craton (DC) and the Southern Granulite Terrain(SGT) of the Indian shield represent large segments of Precambrian middle to lower crust and preserve a geological record spanning from Mesoarchean to Cambrian. This study illuminates the deep structure of the Palghat-Cauvery Shear Zone System (PCSS) and the Palghat-Cauvery Suture Zone (PCSZ) that comprise crustal-scale structures related to multiple episodes of orogeny, crust formation and reworking. We utilize here 3202 high quality P-receiver functions computed using new data from a 23 station seismic network operated by us. Results show a thick (>38?km) mafic (Poisson's ratio >0.25) crust beneath the SGT. The change in crustal thickness is gradual, with a shallower Moho towards the south of PCSZ. We found little evidence for drastic changes in crustal thickness across prominent shear zones like the PCSZ and Moyar-Bhavani. Few seismic stations located along these boundaries have shown evidence for dipping reflectors around 8-20?km depth, with strikes matching well with the trends of surface geological sutures. We opine that these suture zones do not show indications of a terrane boundary. However, a drastic change in the crustal thickness is observed around the prograde metamorphic transition zone or broadly, the “Fermor line”, which separates rocks of Chanockitic (Orthopyroxene bearing granitoid) and non-Charnockitic (Orthopyroxene-free granitoid) mineral assemblage, further north beneath the DC. We suggest that thicknening of crust north of Moyar-Attur Shear Zone (MASZ) and around Fermor line is related to subduction processes operative during the Precambrian.
Earth and Planetary Science Letters, Vol. 498, pp. 77-87.
Mantle
geophysics - seismic
Abstract: Earth's transition zone at depths between 410 km and 660 km plays a key role in Earth's deep water cycle since large amounts of hydrogen can be stored in the nominally anhydrous minerals wadsleyite and ringwoodite, . Previous mineral physics experiments on iron-free wadsleyite proposed low seismic velocities as an indicative feature for hydration in the transition zone. Here we report simultaneous sound wave velocity and density measurements on iron-bearing wadsleyite single crystals with 0.24 wt-% . By comparison with earlier studies, we show that pressure suppresses the velocity reduction caused by higher degrees of hydration in iron-bearing wadsleyite, ultimately leading to a velocity cross-over for both P-waves and S-waves. Modeling based on our experimental results shows that wave speed variations within the transition zone as well as velocity jumps at the 410-km seismic discontinuity, both of which have been used in previous work to detect mantle hydration, are poor water sensors. Instead, the impedance contrast across the 410-km seismic discontinuity that is reduced in the presence of water can serve as a more robust indicator for hydrated parts of the transition zone.
Abstract: On the early Earth, oceanic plateaux similar to present-day Iceland are thought to have evolved into less dense microcontinents as they thickened by continued melt intrusion and crustal fractionation. These earliest continents may have been so weak on a hotter Earth that they collapsed laterally in response to thickening by further magmatic growth or tectonic imbrication. This continental spreading is likely to have resulted in the development of pervasive ductile strain fabrics in the deeper crust, which, if preserved, could generate seismic reflections. Here we present seismic images from the ancient core of the Archaean Yilgarn Craton of Australia that reveal shallowly dipping to horizontal reflections that pervade the middle and lower crust. We interpret these reflective fabrics as the result of widespread lateral crustal flow during the late stage of craton evolution approximately 2.66 to 2.61?billion years ago, which coincided with the widespread intrusion of high-temperature crustal melts, as thickened early continental crust collapsed. The consequent subsidence of large regions of the upper crust, including volcanic and sedimentary greenstone rocks, in the hanging walls of listric mid-lower crustal ductile flow fabrics caused these rocks to drop beneath the granitic melts rising towards the surface, and did not involve Rayleigh-Taylor instabilities within a mostly mobile crust.
Geochemistry, Geophysics, Geosystems, Vol. 19, 7, pp. 2062-2086. doi.org/10/1029/ 2018GC007534
Mantle
geophysics - seismics
Abstract: Some seismic models derived from tomographic studies indicate elevated shear?wave velocities (?4.7 km/s) around 120-150 km depth in cratonic lithospheric mantle. These velocities are higher than those of cratonic peridotites, even assuming a cold cratonic geotherm (i.e., 35 mW/m2 surface heat flux) and accounting for compositional heterogeneity in cratonic peridotite xenoliths and the effects of anelasticity. We reviewed various geophysical and petrologic constraints on the nature of cratonic roots (seismic velocities, lithology/mineralogy, electrical conductivity, and gravity) and explored a range of permissible rock and mineral assemblages that can explain the high seismic velocities. These constraints suggest that diamond and eclogite are the most likely high?Vs candidates to explain the observed velocities, but matching the high shear?wave velocities requires either a large proportion of eclogite (>50 vol.%) or the presence of up to 3 vol.% diamond, with the exact values depending on peridotite and eclogite compositions and the geotherm. Both of these estimates are higher than predicted by observations made on natural samples from kimberlites. However, a combination of ?20 vol.% eclogite and ~2 vol.% diamond may account for high shear?wave velocities, in proportions consistent with multiple geophysical observables, data from natural samples, and within mass balance constraints for global carbon. Our results further show that cratonic thermal structure need not be significantly cooler than determined from xenolith thermobarometry.
Abstract: We investigate the mantle of central-eastern Greenland by using recordings with data from 24 local broad-band seismograph stations. We apply P wave receiver function technique and evaluate the difference in the arrival times of seismic phases that are formed by P to SV mode conversion at the 410-km and 660-km seismic discontinuities. These boundaries mark the top and bottom of the mantle transition zone (MTZ). The difference in the arrival time of the phases from the 410-km and 660-km discontinuities is sensitive to the thickness of the MTZ and relatively insensitive to volumetric velocity anomalies above the 410-km discontinuity. Near the east coast of Greenland in the region of the Skaergaard basalt intrusions we find two regions where the differential time is reduced by more than 2 s. The 410-km discontinuity in these regions is depressed by more than 20 km. The depression may be explained by a temperature elevation of 150 °C. We hypothesize that the basaltic intrusions and the temperature anomalies at a depth of 400 km are, at least partly, effects of the passage of Greenland over the Iceland hotspot at about 55 Ma. This explanation is consistent with the concept of tectosphere and implies that the upper mantle to a depth of 400 km translates coherently with the Greenland plate.
Journal of Geophysical Research, Vol. 123, 7, pp. 5690-5709.
Canada, Nunavut
Geophysics - seismic
Abstract: The geology of northern Hudson Bay, Canada, documents more than 2 billion years of history including the assembly of Precambrian and Archean terranes during several Paleoproterozoic orogenies, culminating in the Trans?Hudson Orogen (THO) ?1.8 Ga. The THO has been hypothesized to be similar in scale and nature to the ongoing Himalaya?Karakoram?Tibetan orogen, but the nature of lithospheric terrane boundaries, including potential plate?scale underthrusting, is poorly understood. To address this problem, we present new P and S wave tomographic models of the mantle seismic structure using data from recent seismograph networks stretching from northern Ontario to Nunavut (60-100?W and 50-80?N). The large size of our network requires careful mitigation of the influence of source side structure that contaminates our relative arrival time residuals. Our tomographic models reveal a complicated internal structure in the Archean Churchill plate. However, no seismic wave speed distinction is observed across the Snowbird Tectonic Zone, which bisects the Churchill. The mantle lithosphere in the central region of Hudson Bay is distinct from the THO, indicating potential boundaries of microcontinents and lithospheric blocks between the principal colliders. Slow wave speeds underlie southern Baffin Island, the leading edge of the generally high wave speed Churchill plate. This is interpreted to be Paleoproterozoic material underthrust beneath Baffin Island in a modern?style subduction zone setting.
Journal of Geophysical Research, Vol. 123, 7, pp. 5644-5660.
Africa, east Africa
geophysics - seismic
Abstract: Although the East African rift system formed in cratonic lithosphere above a large?scale mantle upwelling, some sectors have voluminous magmatism, while others have isolated, small?volume eruptive centers. We conduct teleseismic shear wave splitting analyses on data from 5 lake?bottom seismometers and 67 land stations in the Tanganyika?Rukwa?Malawi rift zone, including the Rungwe Volcanic Province (RVP), and from 5 seismometers in the Kivu rift and Virunga Volcanic Province, to evaluate rift?perpendicular strain, rift?parallel melt intrusion, and regional flow models for seismic anisotropy patterns beneath the largely amagmatic Western rift. Observations from 684 SKS and 305 SKKS phases reveal consistent patterns. Within the Malawi rift south of the RVP, fast splitting directions are oriented northeast with average delays of ~1 s. Directions rotate to N?S and NNW north of the volcanic province within the reactivated Mesozoic Rukwa and southern Tanganyika rifts. Delay times are largest (~1.25 s) within the Virunga Volcanic Province. Our work combined with earlier studies shows that SKS?splitting is rift parallel within Western rift magmatic provinces, with a larger percentage of null measurements than in amagmatic areas. The spatial variations in direction and amount of splitting from our results and those of earlier Western rift studies suggest that mantle flow is deflected by the deeply rooted cratons. The resulting flow complexity, and likely stagnation beneath the Rungwe province, may explain the ca. 17 Myr of localized magmatism in the weakly stretched RVP, and it argues against interpretations of a uniform anisotropic layer caused by large?scale asthenospheric flow or passive rifting.
Geochemistry, Geophysics, Geosystems, Vol. 19, 8, pp. 2593-2614.
United States
geophysics - seismic
Abstract: Juxtaposed terranes of highly varied tectonic history make up the contiguous U.S.: the tectonically active western U.S., the largely quiescent Archean and Proterozoic cratons of the central U.S., and the Phanerozoic orogen and rifted margin of the eastern U.S. The transitions between these regions are clearly observed with Sp converted wave images of the uppermost mantle. We use common conversion point stacked Sp waves recorded by EarthScope's Transportable Array and other permanent and temporary broadband stations to image the transition from a strong velocity decrease at the lithosphere?asthenosphere boundary (or LAB) beneath the western U.S. to deeper, less continuous features moving east that largely lie within the lithosphere. Only sparse, localized, weak phases are seen at LAB depths beneath the cratonic interior. Instead, we observe structures within the cratonic lithosphere that are most prominent within the Archean lithosphere of the Superior Craton. The transition from west to east is clearly revealed by cluster analysis, which also shows eastern U.S. mantle velocity gradients as more similar to the western U.S. than the ancient interior, particularly beneath New England and Virginia. In the western U.S., the observed strong LAB indicates a large enough velocity gradient (an average velocity drop of 10?±?4.5% distributed over 30?±?15 km) to imply that melt has ponded beneath the lithosphere.
Abstract: Two-station surface-wave analysis was used to measure Rayleigh-wave phase velocities between 105 station pairs in western Canada, straddling the boundary between the tectonically active Cordillera and the adjacent stable craton. Major variations in phase velocity are seen across the boundary at periods from 15 to 200 s, periods primarily sensitive to upper mantle structure. Tomographic inversion of these phase velocities was used to generate phase velocity maps at these periods, indicating a sharp contrast between low-velocity Cordilleran upper mantle and high-velocity cratonic lithosphere. Depth inversion along selected transects indicates that the Cordillera-craton upper mantle contact varies in dip along the deformation front, with cratonic lithosphere of the Taltson province overthrusting Cordilleran asthenosphere in the northern Cordillera, and Cordilleran asthenosphere overthrusting Wopmay lithosphere further south. Localized high-velocity features at sub-lithospheric depths beneath the Cordillera are interpreted as Farallon slab fragments, with the gap between these features indicating a slab window. A high-velocity feature in the lower lithosphere of the Slave province may be related to Proterozic or Archean subduction.
Abstract: We present original results that contribute to the understanding of lithospheric structures modification of regions that have witnessed superimposition of multiple tectonic events throughout their geological history. We analyze satellite gravity data through two-dimensional radially-averaged power spectral analysis as well as passive seismic data through thermal modeling to image the depth to the Moho and the lithosphere - asthenosphere boundary (LAB beneath the Central Africa Orogenic Belt (CAOB). The CAOB is an ENE-trending deformation belt extending from Cameroon in the west to Sudan in the east. In Cameroon, it is found on the northern edge of the Congo craton represented by the Oubanguides orogenic belt (the Western Cameroon, the Adamawa - Yade, and the Yaoundé domains). It coincides with the Adamawa plateau and the Benue Trough, and it is spotted by the Cenozoic Cameroon Volcanic Line (CVL). The CAOB was formed during the Precambrian Greater Gondwana assembly but was reactivated during the Mesozoic as a result of Gondwana breakup. We find deeper Moho and LAB) beneath Congo craton and the Yaoundé domain reaching ~50?km and ~200?km, respectively. We map shallower Moho and LAB beneath the CAOB (together with the Adamawa plateau and the Benue trough) reaching ~25?km and ~70?km, respectively. We interpret the shallower LAB beneath the CAOB as due to zonal sub-continental lithospheric mantle (SCLM) delamination along the northern edge of the Congo craton that occurred in association with collisional assembly of Greater Gondwana. This allowed for channelization of mantle flow during the Cenozoic resulting in the formation of the CVL and the uplift of the Adamawa plateau. Our approach can be used to understand the modification of lithospheric structures beneath other terrains that have long tectonic history.
Earth and Planetary Science Letters, Vol. 503, pp. 108-117.
Africa, South Africa, Zimbabwe
geophysics - seismic
Abstract: The depths of the 410 km (d410) and 660 km (d660) discontinuities beneath southern Africa, which is presumably underlain by the lower-mantle African Superswell, are imaged in 1? radius consecutive circular bins using over 6400 P-to-S receiver functions (RFs) recorded by 130 seismic stations over a 27 yr period. When the IASP91 standard Earth model is utilized for moveout correction and time-depth conversion, a normal mantle transition zone (MTZ) thickness of 246 ± 7 km is observed, suggesting that the Superswell has no discernible effect on mantle transition zone temperature. Based on the negligible disparity of the mean MTZ thicknesses between on (246 ± 6 km) and off (246 ± 8 km) cratonic regions, we conclude that the deep Archean cratonic keels possess limited influence on MTZ thermal structure. The apparently shallower-than-normal MTZ discontinuities and the parallelism between the d410 and d660 are mostly the results of upper mantle high wave speed anomalies probably corresponding to a thick lithosphere with a mean thickness of about 245 km beneath the Kaapvaal and 215 km beneath the Zimbabwe cratons. In contradiction to conclusions from some of the previous studies, the resulting spatial distribution of the stacking amplitudes of the P-to-S converted phases at the discontinuities is inconsistent with the presence of an excessive amount of water in the MTZ and atop the d410.
Abstract: Across the entire mantle we interpret 94 positive seismic wave-speed anomalies as subducted lithosphere and associate these slabs with their geological record. We document this as the Atlas of the Underworld, also accessible online at www.atlas-of-the-underworld.org, a compilation comprising subduction systems active in the past ~ 300 Myr. Deeper slabs are correlated to older geological records, assuming no relative horizontal motions between adjacent slabs following break-off, using knowledge of global plate circuits, but without assuming a mantle reference frame. The longest actively subducting slabs identified reach the depth of ~ 2500 km and some slabs have impinged on Large Low Shear Velocity Provinces in the deepest mantle. Anomously fast sinking of some slabs occurs in regions affected by long-term plume rising. We conclude that slab remnants eventually sink from the upper mantle to the core-mantle boundary. The range in subduction-age versus - depth in the lower mantle is largely inherited from the upper mantle history of subduction. We find a significant depth variation in average sinking speed of slabs. At the top of the lower mantle average slab sinking speeds are between 10 and 40 mm/yr, followed by a deceleration to 10-15 mm/yr down to depths around 1600-1700 km. In this interval, in situ time-stationary sinking rates suggest deceleration from 20 to 30 mm/yr to 4-8 mm/yr, increasing to 12-15 mm/yr below 2000 km. This corroborates the existence of a slab deceleration zone but we do not observe long-term (> 60 My) slab stagnation, excluding long-term stagnation due to compositional effects. Conversion of slab sinking profiles to viscosity profiles shows the general trend that mantle viscosity increases in the slab deceleration zone below which viscosity slowly decreases in the deep mantle. This is at variance with most published viscosity profiles that are derived from different observations, but agrees qualitatively with recent viscosity profiles suggested from material experiments.
Abstract: We present original results that contribute to the understanding of lithospheric structures modification of regions that have witnessed superimposition of multiple tectonic events throughout their geological history. We analyze satellite gravity data through two-dimensional radially-averaged power spectral analysis as well as passive seismic data through thermal modeling to image the depth to the Moho and the lithosphere - asthenosphere boundary (LAB beneath the Central Africa Orogenic Belt (CAOB). The CAOB is an ENE-trending deformation belt extending from Cameroon in the west to Sudan in the east. In Cameroon, it is found on the northern edge of the Congo craton represented by the Oubanguides orogenic belt (the Western Cameroon, the Adamawa - Yade, and the Yaoundé domains). It coincides with the Adamawa plateau and the Benue Trough, and it is spotted by the Cenozoic Cameroon Volcanic Line (CVL). The CAOB was formed during the Precambrian Greater Gondwana assembly but was reactivated during the Mesozoic as a result of Gondwana breakup. We find deeper Moho and LAB) beneath Congo craton and the Yaoundé domain reaching ~50?km and ~200?km, respectively. We map shallower Moho and LAB beneath the CAOB (together with the Adamawa plateau and the Benue trough) reaching ~25?km and ~70?km, respectively. We interpret the shallower LAB beneath the CAOB as due to zonal sub-continental lithospheric mantle (SCLM) delamination along the northern edge of the Congo craton that occurred in association with collisional assembly of Greater Gondwana. This allowed for channelization of mantle flow during the Cenozoic resulting in the formation of the CVL and the uplift of the Adamawa plateau. Our approach can be used to understand the modification of lithospheric structures beneath other terrains that have long tectonic history.
McPeak, S., Samson, C., Lamontagne, M., Elliott, B.
Application of passive seismic methodologies to the determination of overburden thickness.
2018 Yellowknife Geoscience Forum , p. 111-112. abstract
Canada, Northwest Territories
geophysics - seismics
Abstract: Diamond mining is central to the economic development of the Canadian North. Innovative methods are needed to identify new prospective targets, as many of them are hidden beneath a thick overburden of glacial sediments. Passive seismics is an emerging method used to map the thickness of near-surface geological layers. Vibrations from distant earthquakes are used as a source of signal and data is processed to estimate the depth of the interface between the overburden and the underlying bedrock. In July 2018, four Tromino seismographs were taken to a study site located approximately ten minutes driving north of Yellowknife. A total of 146 Tromino measurements and associated GPS elevation measurements were taken at 6 m intervals along a dirt road. Elevation measurements were averaged over the course of four days and the survey line was approximately 740 m long. Results indicated that depth to bedrock decreases near outcrops and increases in valleys however; another geophysical dataset is needed to validate the passive seismic data.
Earth and Planetary Science Letters, Vol. 509, 1, pp. 88-99.
Mantle
geophysics - seismic
Abstract: An endothermic phase transition in mantle material at 660-km depth constitutes a barrier that in most cases prevents the direct penetration of subducted slabs. Seismic tomography shows that subducted material is in many subduction zones trapped at the bottom of the transition zone, just above the 660-km phase boundary. Recent tomographic models however also report subducted material that penetrates to the shallow lower mantle, and there it is observed to flatten at about 1000-km depth. Models of slab dynamics that generally assume sharp rheological transition at 660-km depth, however, mostly predict slab stagnation at the bottom of the transition zone. Multiple lines of evidence, including recent experiments, indicate that viscosity may gradually increase in the uppermost ?300 km of the lower mantle, rather than simply changing abruptly at the upper-lower mantle boundary. Here we present the results of a modeling study focused on the effects of rheological transition between upper and lower mantle material on slab deformation and stagnation. We test the effects of smoothing the viscosity increase over 300 km and shifting it to a depth of 1000 km or even deeper. We show that slab ability to penetrate to the lower mantle is mainly controlled by the trench migration rate, which in turn is affected by crustal viscosity. Coupling between the subducting and overriding plates thus plays a key role in controlling slab penetration to the lower mantle and stagnation in the deep transition zone or shallow lower mantle. Models with strong crust and consequently negligible rollback display penetration to the lower mantle without much hindrance and no stagnation above or below the 660-km interface, regardless of viscosity stratification in the shallow lower mantle. Models with weak crust are characterized by fast rollback, and penetration is very limited as slabs buckle horizontally and flatten above the 660-km boundary. Most interesting from the point of view of shallow lower mantle stagnation are models with intermediate crustal viscosity. Here rollback is efficient, though slower than in weak-crust cases. Horizontally lying slab segments are trapped in the transition zone if the sharp viscosity increase occurs at 660 km, but shifting the viscosity increase to 1000 km depth allows for efficient sinking of the flat-lying part and results in temporary stagnation below the upper-lower mantle boundary at about 1000 km depth.
Abstract: Measuring seismic anisotropy within the Earth is essential as it constitutes a proxy for inferring upper mantle deformation related to mantle flow, that develops preferred orientations of the minerals in response to tectonic strain. The most-used method to detect anisotropy beneath a seismic station is the measurement of teleseismic SKS wave splitting on two horizontal recordings, i.e. measuring the delay time (?t) between two fast- and slow- polarized shear-waves and the orientation of polarization (?). This technique allows a integrative measurement (SKS data, hereinafter) that estimates the average ? and ?t along the entire SKS ray-path. Despite its importance for large-scale anisotropy within the upper mantle, the analysis of SKS data suffers from several limitations : (1) SKS data become difficult to interpret in regions where several anisotropic layers occur; (2) SKS waves fail to provide robust information about anisotropy related to thin layers; and (3) SKS data can investigate rock volumes with an horizontal symmetry axis only. During the last decade a new method, called harmonic decomposition of teleseismic Receiver Functions (RFs) has been developed in order to detect more complex anisotropic layering. This methodology is based on the extraction of back-azimuth harmonics of the RF dataset. Briefly, it constitutes a tool to appreciate the value of ? and ?t at every depth-level affording a detailed study of the rock anisotropy with both plunging and horizontal symmetry axis. RFs studies are however commonly limited to the first 10-15s of the signal and do not sample the deepest anisotropy. In this work we investigate in details both SKS data and RFs harmonic decomposition for a pool of stations deployed in northeastern Brasil, in order to understand how results from the analysis of these two observables can be jointly interpreted. We focus our study on the permanent station RCBR and on temporary seismic stations deployed in the area. We show that comparison and/or joint interpretation is not straightforward as both results can vary according to the amount of data available and their distribution in back-azimuth, and filtering. However, tacking into account those issues, the integration of these two observables represent a great step-forward for robust detection of upper mantle anisotropy.
Journal of the Geological Society of India, Vol. 93, 2, pp. 163-170.
India
geophysics - seismic
Abstract: The crustal structures of the Narmada region in Central India bounded by fault system (Narmada- North and South faults : NNF and NSF) has been derived from deep seismic sounding (DSS) studies along the two profiles trending almost north-south direction. The wide-angle phases have been modeled kinematically and dynamically using the 2-D asymptotic ray tracing technique. The combined seismic and gravity modeling reveals a multilayer crust in the region. The crustal wide-angle reflection phases map the Moho discontinuity, where the P-wave velocity jumps from 7.2 km s-1 to 8.0-8.1 km s-1, at depth varying between 38 km and 44 km. A layer with velocity 7.2 km s-1, exists above the Moho in most parts of the profiles and is attributed to the magmatic underplating related to the Deccan volcanism (~65 Ma). The intriguing observation of the study is a zone characterized by anomalous high velocity (6.5-6.6 km s-1) within the upper crust. 2-D gravity modeling demonstrates that this anomalous layer has a density of ~2.9 gm cm-3, which is equivalent to the rocks metamorphosed to granulite/amphibolite facies. This high velocity layer probably represents the granulite enclaves within the Archaean granites/gneiss rocks and was formed during the cratonization of the Achaean crust. Importantly, this high velocity layer shows an average upward displacement of ~8.5 km within the region bounded by NNF and NSF as compared to the regions beyond it. The studies suggest that the observed displacement in the high velocity layer of the upper crust is a result of repeated reactivation of the Narmada fault system.
Abstract: Cratons are, by definition, the most tectonically stable and oldest parts of the continental lithosphere on Earth. The Archean Slave craton is located in the northwestern part of the Canadian Shield. The propensity of diamondiferous kimberlite pipes in the central Slave craton raises many questions regarding their structural environment and source. Here, we provide the most robust teleseismic P and S body wave tomography models over the Slave craton region based on 20,547 P-wave delay times, 6,140 direct S-wave delay times and 3,381 SKS delay times. The P-wave model reveals an alternating pattern of relative positive and negative anomalies over a fine broad scale region within the central Slave craton. Furthermore, the P-wave model revealed two fine structures located in the lithosphere beneath the Lac de Gras kimberlite cluster, with relatively slow anomalies (B - C) that extend from 75 km to 350 km depths with an apparent dip to the north. These relatively slow P- and S-wave anomalies are associated with metasomatised regions within the lithosphere. The S-wave model displays a slow S-wave anomaly lying from 300 km depth to the transition zone beneath the central Slave craton. This anomaly is located beneath the Lac de Gras kimberlite cluster. We suggest that this anomaly is not the cause of the actual kimberlites at the surface since last eruption occurred 75-45 Ma ago but may be related to a potential kimberlite magma ascent in the asthenosphere.
Abstract: Upper mantle anisotropy investigated using 172 core-refracted (SKS, SKKS) seismic phases along a ~660?km long profile at 10 to 20?km intervals from the west to the east coast of South India reveals significant lateral variations in its magnitude and direction. This profile, with 38 broadband seismic stations, covers mid-Archean Western Dharwar craton (WDC), late-Archean Eastern Dharwar Craton (EDC), Proterozoic Cuddapah Basin (CB) and the passive continental margins along the west and east coast. The observed fast polarization directions (FPDs) show lateral variability: NW50o to NW5o beneath the WDC, NW40o to NE30o beneath the EDC and N5o to N85o beneath the CB and further east. The delay time varies between 0.4 and 2.0?s with an average of 1?s. However, we are unable to fit a two layers anisotropy model for the region due to sparse azimuth coverage. Beneath the WDC, the direction of the fast axis follows trends of shear zones and faults, suggesting “frozen-in” anisotropy in the lithosphere, possibly established during the lithospheric evolution in mid-late Archean. In the EDC, the fast axis does not only follow the plate motion direction but it deviates, manifesting late Archean to Proterozoic deformation may still be present as fossil mantle anisotropy. The splitting trend beneath the CB and Eastern Ghat (EG) follows the strike of the rift along with plate motion direction, indicating anisotropy is influenced by the combination of “frozen” anisotropy due to continental rifting along the eastern margin of Indian plate and active asthenospheric flow.
Earth and Planetary Science Letters, Vol. 514, 1, pp. 143-155.
Mantle
geophysics - seismic
Abstract: The alignment of intrinsically anisotropic olivine crystals under convection is typically invoked as the cause of the bulk of seismic anisotropy inferred from shear-wave splitting (SWS). This provides a means of constraining the interplay between continental dynamics and the deep mantle, in particular for densely instrumented regions such as North America after USArray. There, a comparison of “fast orientations” from SWS with absolute plate motions (APM) suggests that anisotropy is mainly controlled by plate motions. However, large regional misfits and the limited realism of the APM model motivate us to further explore SWS based anisotropy. If SWS is estimated from olivine alignment in mantle circulation instead, plate-driven flow alone produces anisotropy that has large misfits with SWS. The addition of large-scale mantle density anomalies and lateral viscosity variations significantly improves models. Although a strong continental craton is essential, varying its geometry does, however, not improve the plate-scale misfit. Moreover, models based on higher resolution tomography degrade the fit, indicating issues with the flow model assumptions and/or a missing contributions to anisotropy. We thus compute a “lithospheric complement” to achieve a best-fit, joint representation of asthenospheric and frozen-in lithospheric anisotropy. The complement shows coherent structure and regional correlation with independently imaged crustal and upper mantle anisotropy. Dense SWS measurements therefore provide information on depth-dependent anisotropy with implications for tectonics, but much remains to be understood about continental anisotropy and its origin.
Earth and Planetary Science Letters, Vol. 519, pp. 297-307.
United States
geophysics - seismics
Abstract: The passive margin of the eastern coast of the United States is known to be geologically active, with recently rejuvenated topography, intraplate seismicity, and volcanism of Eocene age. This study uses seismic data from the Mid-Atlantic Geophysical Integrative Collaboration (MAGIC) experiment to constrain lateral variations in the attenuation of teleseismic P waves beneath the central Appalachian Mountains to shed light on the structure and dynamics of the upper mantle at this “active” passive margin. We use a Monte Carlo approach to estimate variations in attenuation along with both data and model uncertainties. The quality factor of the upper mantle dramatically decreases over a distance of less than 50 km on the western side of the central Appalachian Mountains, where a low-velocity anomaly has been previously inferred. Extrinsic factors such as scattering or focusing are rejected as explanations for the observations on the basis of finite-difference waveform modeling experiments. The peak in attenuation beneath the crest of the Appalachian Mountains requires that near- to super-solidus conditions occur in the upper mantle and is co-located with volcanism of Eocene age. Our preferred interpretation is that the attenuation reflects the removal of the mantle lithosphere via delamination beneath the mountains, followed by ongoing small-scale convection.
Geotreks.com.au/ work/giant-ring-structures/ north-america/ superior-craton, May CS#1
Canada, Ontario
geophysics - seismics
Abstract: This paper the first of a series of papers to describe the genesis and mineralisation of the North American Superior Greater Craton from the 60 to 300 km depth using detailed seismic tomography. Greater Superior Craton occupies the core of the North American continent.
Abstract: We studied a suite of mantle xenoliths carried by Cenozoic volcanism in the Borborema Province, NE Brazil. These xenoliths sample a subcontinental lithospheric mantle affected by multiple continental convergence and rifting events since the Archean. Equilibrium temperatures indicate a rather hot geotherm, implying a ca. 80?km thick lithosphere. Most xenoliths have coarse-granular and coarse-porphyroclastic microstructures, recording variable degrees of annealing following deformation. The high annealing degree and equilibrated pyroxene shapes in coarse-granular peridotites equilibrated at ~900?°C indicate that the last deformation event that affected these peridotites is several hundreds of Ma old. Coarse-porphyroclastic peridotites equilibrated at 950-1100?°C probably record younger (Cretaceous?) deformation in the deep lithospheric mantle. In addition, a few xenoliths show fine-porphyroclastic microstructures and equilibrium temperatures ?1200?°C, which imply recent deformation, probably related to the dykes that fed the Cenozoic volcanism. Chemical and microstructural evidence for reactive percolation of melts is widespread. Variation in textural and chemical equilibrium among samples implies multiple melt percolation events well spaced in time (from Neoproterozoic or older to Cenozoic). Crystal preferred orientations of olivine and pyroxenes point to deformation controlled by dislocation creep with dominant activation of the [100](010) and [001]{0kl} slip systems in olivine and pyroxenes, respectively, for all microstructures. Comparison of xenoliths' seismic properties to SKS splitting data in the nearby RCBR station together with the equilibrated microstructures in the low-temperature xenoliths point to coupled crust-mantle deformation in the Neoproterozoic (Brasiliano) continental-scale shear zones, which is still preserved in the shallow lithospheric mantle. This implies limited reworking of the lithospheric mantle in response to extension during the opening of the Equatorial Atlantic in the Cretaceous, which in the present sampling is restricted to the base of the lithosphere.
Abstract: A 5.6-km-long line of refraction and reflection seismic data spanning the Pliocene-Pleistocene fill of the Olduvai Basin, Tanzania is presented. The line is oriented along a northwest-southeast profile through the position of Olduvai Gorge Coring Project (OGCP) Borehole 2A. Our aims are to (1) delineate the geometry of the basin floor by tracing bedrock topography of the metaquartzitic and gneissic basement, (2) map synsedimentary normal faults and trace individual strata at depth, and (3) provide context for the sequence observed in OGCP cores. Results with refraction tomography and poststack migration show that the maximum basin depth is around 405?m (±25?m) in the deepest portion, which quadruples the thickness of the basin-fill previously known from outcrops. Variations in seismic velocities show the positions of lower density lake claystones and higher density well-cemented sedimentary sequences. The Bed I Basalt lava is a prominent marker in the refraction seismic results. Bottom-most sediments are dated to >2.2?Ma near where Borehole 2A bottoms out at the depth of 245?m. However, the seismic line shows that the basin-fill reaches a maximum stratigraphic thickness of around 380?m deep at Borehole 2A, in the western basin where the subsidence was greatest. This further suggests that potential hominin palaeoenvironments were available and preserved within the basin-fill possibly as far back as around 4?Ma, applying a temporal extrapolation using the average sediment accretion rate.
Progress in Earth and Planetary Science, doi.org/10.1186 /s40645-019-0299-5
Mantle
geophysics - seismic
Abstract: The depth extent of the crustal seismogenic zone is closely related to the size of earthquakes. The mechanisms that control the depth of the lower transition of the seismogenic zone are important issues in seismology and disaster mitigation. Laboratory studies have shown that the mechanism of earthquake nucleation is controlled by the frictional properties of fault materials. We measured the velocity dependences of the steady-state friction of quartz and feldspar, two major components of crustal rocks, under dry and wet conditions at temperatures up to 600?°C. In the presence of water, the temperature range over which the velocity dependence of steady-state friction was negative was wider for feldspar than for quartz, thus indicating that the temperature range of earthquake nucleation is wider for feldspar than for quartz. Considering that temperature increases with depth, our findings indicate that the material properties of feldspar likely play a dominant role in limiting the depth extent of the seismogenic zone.
Abstract: The southern U.S. continental margin records a history spanning ca. 1.2 Ga, including two Wilson cycles. However, due to a thick sediment cover, the paucity of significant local seismicity, and, until recently, sparse instrumentation, details of this passive margin’s tectonomagmatic evolution remain disputed. This paper compares recent S-wave tomography and crustal thickness models based on USArray data to help establish a framework for geodynamic interpretation. Large-scale patterns of crustal velocity anomalies, corresponding to major regional features such as the Ouachita orogenic front and the Precambrian margin, are generally consistent between the models. The spatial extent of smaller-scale tectonic features, such as the Sabine Uplift and Wiggins block, remains poorly resolved. An inverse relationship between crustal thickness and Bouguer gravity across the continental margin is observed. This model comparison highlights the need for additional P-wave tomography studies and targeted, higher density station deployments to better constrain tectonic features.
Journal of Geophysical Research: Solid Earth, Vol. 124, 7, pp. 7105-7119.
United States
geophysics - seismic
Abstract: Detailed P wave velocity and anisotropy structure of the uppermost mantle below the central United States is presented based on a tomographic inversion of Pn traveltimes for earthquakes in the range 2 to 14°. Dense raypath coverage throughout the northern Mississippi Embayment is obtained using the Northern Embayment Lithosphere Experiment and U.S. Transportable Array data sets. A detailed analysis of the trade?off between velocity and anisotropy variations demonstrates that both are well resolved over most of the study area. Anomalously fast Pn velocities are identified below the northern Mississippi Embayment, centered on the New Madrid seismic zone. A prominent region of low velocity coincides with the southwestern margin of the Illinois basin. Pn anisotropy displays complex patterns and differs from absolute plate motion directions and SKS splitting directions. A circular pattern of fast anisotropy directions is centered on the New Madrid seismic zone and may be related to the presence of the mafic “rift pillow.”
Abstract: As mineral exploration seeks deeper targets there will be a greater reliance on geophysical data and a better understanding of the geological meaning of the responses will be required, and this must be achieved with less geological control from drilling. Also, exploring based on the mineral system concept requires particular understanding of geophysical responses associated with altered rocks. Where petrophysical datasets of adequate sample size and measurement quality are available, physical properties show complex variations, reflecting the combined effects of various geological processes. Large datasets, analysed as populations, are required to understand the variations. We recommend the display of petrophysical data as frequency histograms as the nature of the data distribution is easily seen with this form of display. A petrophysical data set commonly contains a combination of overlapping sub?populations, influenced by different geological factors. To understand the geological controls on physical properties in hard rock environments it is necessary to analyse the petrophysical data not only in terms of the properties of different rock types. It is also necessary to consider the effects of processes such as alteration, weathering, metamorphism, and strain, and variables such as porosity and stratigraphy. To address this complexity requires that much more supporting geological information be acquired than is current practice. The widespread availability of field portable instruments means quantitative geochemical and mineralogical data can now be readily acquired, making it unnecessary to rely primarily on categorical rock classification schemes. The petrophysical data can be combined with geochemical, petrological and mineralogical data to derive explanations for observed physical property variations based not only on rigorous rock classification methods, but also in combination with quantitative estimates of alteration and weathering. To understand how geological processes will affect different physical properties it is useful to define three end?member forms of behaviour. Bulk behaviour depends on the physical properties of the dominant mineral components. Density and, to a lesser extent, seismic velocity show such behaviour. Grain and texture behaviour occur when minor components of the rock are the dominate controls on its physical properties. Grain size and shape control grain properties, and for texture properties the relative positions of these grains are also important. Magnetic and electrical properties behave in this fashion. Thinking in terms of how geological processes change the key characteristics of the major and minor mineralogical components allows the resulting changes in physical properties to be understood and anticipated.
Journal of Geophysical Research: Solid Earth, Vol. 174, 7, pp. 6625-6652.
United States
geophysics - seismic
Abstract: The Southeastern United States is an ideal location to understand the interactions between mountain building, rifting, and magmatism. Line 2 of the Suwannee suture and Georgia Rift basin refraction seismic experiment in eastern Georgia extends 420 km from the Inner Piedmont to the Georgia coast. We model crustal and upper mantle VP and upper crustal VS. The most dramatic model transition occurs at the Higgins?Zietz magnetic boundary, north of which we observe higher upper crustal VP and VS and lower VP/VS. These observations support the interpretation of the Higgins?Zietz boundary as the Alleghanian suture. North of this boundary, we observe a low?velocity zone less than 2 km thick at ~5?km depth, consistent with a layer of sheared metasedimentary rocks that forms the Appalachian detachment. To the southeast, we interpret synrift sediments and decreasing crustal thickness to represent crustal thinning associated with the South Georgia Rift Basin and subsequent continental breakup. The correspondence of the northern limit of thinning with the interpreted suture location suggests that the orogenic suture zone and/or the Gondwanan crust to the south of the suture helped localize subsequent extension. Lower crustal VP and VP/VS preclude volumetrically significant mafic magmatic addition during rifting or associated with the Central Atlantic Magmatic Province. Structures formed during orogenesis and/or extension appear to influence seismicity in Georgia today; earthquakes localize along a steeply dipping zone that coincides with the northern edge of the South Georgia Basin and the change in upper crustal velocities at the Higgins?Zietz boundary.
Geophysical Research Letters, Vol. 46, 16, pp. 9509-9518.
Africa, South Africa
geophysics - seismics
Abstract: P and S waves travel times from large, distant earthquakes recorded on seismic stations in Botswana and South Africa have been combined with existing data from the region to construct velocity models of the upper mantle beneath southern Africa. The models show a region of higher velocities beneath the Rehoboth Province and parts of the northern Okwa Terrane and the Magondi Belt, which can be attributed to thicker cratonic lithosphere, and a region of lower velocities beneath the Damara?Ghanzi?Chobe Belt and Okavango Rift, which can be attributed a region of thinner off?craton lithosphere. This finding suggests that the spatial extent of thick cratonic lithosphere in southern Africa is greater than previously known. In addition, within the cratonic lithosphere an area of lower velocities is imaged, revealing parts of the cratonic lithosphere that may have been modified by younger magmatic events.
Geochemistry, Geophysics, Geosystems, Vol. 20, 8, pp. 4130-4145.
Mantle
geophysics - seismics
Abstract: Geochemical variations in volcanic rocks erupted at Earth's surface indicate differences in mantle composition, but our understanding of the location, formation, and history of compositionally distinct mantle domains remains incomplete. In particular, some hotspot lavas contain signatures of primitive regions within the mantle that have remained relatively isolated and unprocessed throughout Earth's history. Here we use models of mantle flow to predict the locations within the mantle that are sampled by the mantle plumes associated with hotspot volcanism. Combining these models of mantle flow with state-of-the-art seismic images and a comprehensive catalog of hotspot lava geochemistry, we find that hotspots with a more primitive geochemical signature (as indicated by the isotopes of helium) sample the two large low shear-velocity provinces in the lowermost mantle. Complementary constraints from xenon and tungsten isotope ratios associated with primitive materials then require these continent-sized provinces in Earth's deep interior formed early in Earth's history, survived the violent Moon-forming giant impact, and remained relatively unmixed with the rest of the solid Earth over the past 4.5 billion years.
Geophysical Research Letters, Vol. 46, 22, pp. 12862-12869.
Australia
geophysics - seismic
Abstract: Despite decades of study, the mechanisms that lead to the localization of intracontinental seismicity remain vigorously debated. We find a very strong correlation between the attenuation of teleseismic P waves and the occurrence of intraplate seismicity in Australia. The regions with the highest attenuation host ~2 orders of magnitude more earthquakes per unit of area than the least attenuating regions. We argue that the attenuation we observe is produced by lateral variations in the thickness and/or viscosity of the lithospheric mantle and further suggest that the correlation we document implies that lithospheric mantle structure exerts first?order controls on the localization of intraplate seismicity.
Journal of Geophysical Research: Solid Earth, Vol. 124, 12, pp. 13124-131249. ( open access)
South America, Venezuela
geophysics - seismic
Abstract: Northwestern Venezuela is located in the complex deformation zone between the Caribbean and South American plates. Several models regarding the lithospheric structure of the Mérida Andes have been proposed. Nevertheless, they lack relevant structural information in order to support the interpretation of deeper structures. Therefore, a 560?km?long refraction profile across the northern part of Mérida Andes, oriented in a NNW direction, covering areas from the Proterozoic basement in the south, to both Paleozoic and Meso?Cenozoic terranes of northwestern Venezuela to the north, is analyzed in this contribution. Thirteen land shots were recorded by 545 short?deployment seismometers, constraining P wave velocity models from first?arrival seismic tomography and layer?based inversion covering the whole crust in detail, with some hints to upper mantle structures. The most prominent features imaged are absence of a crustal root associated to the Mérida Andes, as the Northern Andes profile is located marginal to the Andean crustal domain, and low?angle subduction of the Caribbean oceanic slab (~10-20°) beneath northwestern South America. Further crustal structures identified in the profile are (a) crustal thinning beneath the Falcón Basin along the western extension of the Oca?Ancón fault system interpreted as a back?arc basin; (b) suture zones between both the Proterozoic and Paleozoic provinces (Ouachita?Marathon?related suture?), and Paleozoic and Meso?Cenozoic terranes (peri?Caribbean suture) interpreted from lateral changes in seismic velocity; and (c) evidence of a deep Paleozoic(?) extensional basin, underlying thick Mesozoic and Cenozoic sequences (beneath the Guárico area).
Abstract: From 2009 to 2016, a drastic increase in seismic activity occurred in the Central and Eastern US (CEUS), particularly in the Oklahoma-Kansas region. The majority of hypocenters were focused in the crystalline basement rock. Information regarding the physical properties (elastic wave velocity, peak strength, etc.) of rocks in the CEUS basement to date is sparse. Forecasting future seismic hazard and predicting the in situ response of the crystalline basement requires their geomechanical parameters be adequately constrained. This work assesses the mechanical and petrophysical properties of several sets of basement rocks from Oklahoma to provide a better framework for understanding intraplate seismicity and overall basement deformation in the continental United States. Laboratory experiments were conducted with granite, rhyolite and diabase basement rock samples collected from southern Oklahoma. Evolution of compressional and shear wave velocity with increasing confinement was measured through a series of ultrasonic velocity tests. A suite of uniaxial and triaxial tests were conducted to measure the elastic and inelastic deformation behavior of the basement rocks. Deformation data was evaluated using the Mohr-Coulomb criterion and compared with additional preexisting deformation data of igneous basement rocks. Dynamic and static elastic properties compare favorably with available field measurements and demonstrate the role physical properties can play in varying mechanical behavior. Granitic samples demonstrate moderate variation of intrinsic physical properties can alter elastic properties and failure behavior significantly. Water-weakening in the basement rocks may indicate fluid-assisted processes such as stress corrosion cracking enhance deformation in the crystalline basement.
Journal of Geophysical Research: Solid Earth, in press available, 24p. Pdf
Canada
geophysics - seismics
Abstract: The geological record of SE Canada spans more than 2.5Ga, making it a natural laboratory for the study of crustal formation and evolution over time. We estimate the crustal thickness, Poisson's ratio, a proxy for bulk crustal composition, and shear velocity (Vs) structure from receiver functions at a network of seismograph stations recently deployed across the Archean Superior craton, the Proterozoic Grenville and the Phanerozoic Appalachian provinces. The bulk seismic crustal properties and shear velocity structure reveal a correlation with tectonic provinces of different ages: the post-Archean crust becomes thicker, faster, more heterogenous and more compositionally evolved. This secular variation pattern is consistent with a growing consensus that crustal growth efficiency increased at the end of the Archean. A lack of correlation among elevation, Moho topography, and gravity anomalies within the Proterozoic belt is better explained by buoyant mantle support rather than by compositional variations driven by lower crustal metamorphic reactions. A ubiquitous ?20km thick high-Vs lower-crustal layer is imaged beneath the Proterozoic belt. The strong discontinuity at 20km may represent the signature of extensional collapse of an orogenic plateau, accommodated by lateral crustal flow. Wide anorthosite massifs inferred to fractionate from a mafic mantle source are abundant in Proterozoic geology and are underlain by high Vs lower crust and a gradational Moho. Mafic underplating may have provided a source for these intrusions and could have been an important post-Archean process stimulating mafic crustal growth in a vertical sense.
Earth and Planetary Science Letters, Vol. 539, 116240 9p. Pdf
Mantle
geophysics - seismics
Abstract: Earth's thermo-chemical structure exerts a fundamental control on mantle convection, plate tectonics, and surface volcanism. There are indications that mantle convection occurs as an intermittent-stage process between layered and whole mantle convection in interaction with a compositional stratification at 660 km depth. However, the presence and possible role of any compositional layering in the mantle remains to be ascertained and understood. By interfacing inversion of a novel global seismic data set with petrologic phase equilibrium calculations, we show that a compositional boundary is not required to explain short- and long-period seismic data sensitive to the upper mantle and transition zone beneath stable continental regions; yet, radial enrichment in basaltic material reproduces part of the complexity present in the data recorded near subduction zones and volcanically active regions. Our findings further indicate that: 1) cratonic regions are characterized by low mantle potential temperatures and significant lateral variability in mantle composition; and 2) chemical equilibration seems more difficult to achieve beneath stable cratonic regions. These findings suggest that the lithologic integrity of the subducted basalt and harzburgite may be better preserved for geologically significant times underneath cratonic regions.
Annual Review of Earth and Planetary Sciences, Vol.48, 25p. Pdf
Mantle
geophysics - seismics
Abstract: Seismology provides important constraints on the structure and dynamics of the deep mantle. Computational and methodological advances in the past two decades improved tomographic imaging of the mantle and revealed the fine-scale structure of plumes ascending from the core-mantle boundary region and slabs of oceanic lithosphere sinking into the lower mantle. We discuss the modeling aspects of global tomography including theoretical approximations, data selection, and model fidelity and resolution. Using spectral, principal component, and cluster analyses, we highlight the robust patterns of seismic heterogeneity, which inform us of flow in the mantle, the history of plate motions, and potential compositionally distinct reservoirs. In closing, we emphasize that data mining of vast collections of seismic waveforms and new data from distributed acoustic sensing, autonomous hydrophones, ocean-bottom seismometers, and correlation-based techniques will boost the development of the next generation of global models of density, seismic velocity, and attenuation. Seismic tomography reveals the 100-km to 1,000-km scale variation of seismic velocity heterogeneity in the mantle. Tomographic images are the most important geophysical constraints on mantle circulation and evolution.
Progress in Earth and Planetary Science , Vol. 7, 16p. Pdf
Mantle
geophysics - seismics
Abstract: Mantle tomography reveals the existence of two large low-shear-velocity provinces (LLSVPs) at the base of the mantle. We examine here the hypothesis that they are piles of oceanic crust that have steadily accumulated and warmed over billions of years. We use existing global geodynamic models in which dense oceanic crust forms at divergent plate boundaries and subducts at convergent ones. The model suite covers the predicted density range for oceanic crust over lower mantle conditions. To meaningfully compare our geodynamic models to tomographic structures, we convert them into models of seismic wavespeed and explicitly account for the limited resolving power of tomography. Our results demonstrate that long-term recycling of dense oceanic crust naturally leads to the formation of thermochemical piles with seismic characteristics similar to the LLSVPs. The extent to which oceanic crust contributes to the LLSVPs depends upon its density in the lower mantle for which accurate data is lacking. We find that the LLSVPs are not composed solely of oceanic crust. Rather, they are basalt rich at their base (bottom 100-200 km) and grade into peridotite toward their sides and top with the strength of their seismic signature arising from the dominant role of temperature. We conclude that recycling of oceanic crust, if sufficiently dense, has a strong influence on the thermal and chemical evolution of Earth’s mantle.
Physics of the Earth and Planetary Interiors, in press available, 51p. Pdf
United States, Canada
geophysics - seismics
Abstract: We examined SKS-SKKS splitting intensity discrepancies for phases that sample the lowermost mantle beneath North America, which has previously been shown to exhibit seismic anisotropy using other analysis techniques. We examined data from 25 long-running seismic stations, along with 244 stations of the temporary USArray Transportable Array, located in the eastern, southeastern and western U.S. We identified 279 high-quality SKS-SKKS wave pairs that yielded well-constrained splitting intensity measurements for both phases. Of the 279 pairs, a relatively small number (15) exhibited discrepancies in splitting intensity of 0.4 s or greater, suggesting a contribution to the splitting of one or both phases from anisotropy in the lowermost mantle. Because only a small minority of SK(K)S phases examined in this study show evidence of being affected by lowermost mantle anisotropy, the traditional interpretation that splitting of these phases primarily reflects anisotropy in the upper mantle directly beneath the stations is appropriate. The discrepant pairs exhibited a striking geographic trend, sampling the lowermost mantle beneath the southern U.S. and northern Mexico, while other regions were dominated by non-discrepant pairs. We carried out ray theoretical modeling of simple anisotropy scenarios that have previously been suggested for the lowermost mantle beneath North America, invoking the alignment of post-perovskite due to flow induced by the impingement of the remnant Farallon slab on the core-mantle boundary. We found that our measurements are generally consistent with this model and with the idea of slab-driven flow, but relatively small-scale lateral variations in the strength and/or geometry of lowermost mantle anisotropy beneath North America are also likely present.
Abstract: The Natal Valley, offshore Mozambique, is a key area for understanding the evolution of East Gondwana. Within the scope of the integrated multidisciplinary PAMELA project, we present new wide?angle seismic data and interpretations, which considerably alter Geoscience paradigms. These data reveal the presence of a 30?km?thick crust that we argue to be of continental nature. This falsifies all the most recent palaeo?reconstructions of the Gondwana. This 30?km?thick continental crust 1,000 m below sea level implies a complex history with probable intrusions of mantle?derived melts in the lower crust, connected to several occurrences of magmatism, which seems to evidence the crucial role of the lower continental crust in passive margin genesis.
Abstract: The Canadian High Arctic preserves a long and complex tectonic history, including craton formation, multiple periods of orogenesis, extension and basin formation, and the development of a passive continental margin. We investigate the possible preservation of deformational structures throughout the High Arctic subcontinental lithosphere using measurements of seismic anisotropy from shear wave splitting at 11 seismograph stations across the region, including a N-S transect along Ellesmere Island. The majority of measurements indicate a fast-polarisation orientation that parallels tectonic trends and boundaries, suggesting that lithospheric deformation is the dominant source of seismic anisotropy in the High Arctic; however, a sub-lithospheric contribution cannot be ruled out. Beneath Resolute in the central Canadian Arctic, distinct back-azimuthal variations in splitting parameters can be explained by two anisotropic layers. The upper layer is oriented E-W and correlates with tectonic trends and the inferred lithospheric deformation history of the region. The lower layer has a ?NNE-SSW orientation and may arise from present-day convective mantle flow beneath locally-thinned continental lithosphere. In addition to inferences of anisotropic structure beneath the Canadian High Arctic, measurements from the far north of our study region suggest the presence of an anisotropic zone in the lowermost mantle beneath northwest Alaska.
Abstract: Deep within Earth's interior, at ?2900 km beneath the surface, lies the boundary between the solid silicate rock mantle and the liquid iron-nickel alloy core (the core-mantle boundary). Geophysicists have studied the complex thermal and chemical dynamics that take place in this boundary layer. In the early 20th century, Gutenberg investigated the structure of the lowermost region, or base, of the mantle by recording with only a few seismograms from a small number of large-magnitude earthquakes that occurred thousands of kilometers away (1). The structure of the rocks just above the core-mantle boundary—designated as D? by Jeffreys in 1939 (2)—forms a distinct layer with surprising complexity. Now, on page 1223 of this issue, Kim et al. (3) describe new structural heterogeneities in the lowermost mantle with the use of a learning algorithm that does not require any a priori knowledge of Earth.
Annual Review of Earth and Planetary Sciences, Vol. 48, 1, pp. w77-401.
Mantle
geophysics - seismics
Abstract: Seismology provides important constraints on the structure and dynamics of the deep mantle. Computational and methodological advances in the past two decades improved tomographic imaging of the mantle and revealed the fine-scale structure of plumes ascending from the core-mantle boundary region and slabs of oceanic lithosphere sinking into the lower mantle. We discuss the modeling aspects of global tomography including theoretical approximations, data selection, and model fidelity and resolution. Using spectral, principal component, and cluster analyses, we highlight the robust patterns of seismic heterogeneity, which inform us of flow in the mantle, the history of plate motions, and potential compositionally distinct reservoirs. In closing, we emphasize that datamining of vast collections of seismic waveforms and new data from distributed acoustic sensing, autonomous hydrophones, ocean-bottom seismometers, and correlation-based techniques will boost the development of the next generation of global models of density, seismic velocity, and attenuation.
Journal of Earth System Science, Vol. 128, 43, 8p. Pdf
India
geophysics, seismics
Abstract: An improved shear wave velocity (Vs) structure of the lithosphere of peninsular India using the surface wave tomography from the ambient noise and earthquake waveforms suggests its bipolar character. While most of the geological domains of India are characterised by a uniform lithospheric mantle of Vs?4.5 km/s, the three cratonic regions, eastern Dharwar, Bastar and Singhbhum, hosting most of the diamondiferous kimberlite fields, show significantly high Vs of 4.7 km/s and above in their lower lithosphere beyond ?90 km depth. The higher velocity could best be explained by the presence of diamond and/or eclogite along with peridotite in mantle. This unique relationship suggests the regional seismic image of lithosphere as a guide for exploration of diamonds.
Science Adavances, Vol. 6, eaba7118 May 27, 9p. Pdf
Asia, Java
geophysics -seismic
Abstract: New passive- and active-source seismic experiments reveal unusually high mantle P-wave speeds that extend beneath the remnants of the world’s largest known large igneous province, making up the 120-million-year-old Ontong-Java-Manihiki-Hikurangi Plateau. Sub-Moho Pn phases of ~8.8 ± 0.2 km/s are resolved with negligible azimuthal seismic anisotropy, but with strong radial anisotropy (~10%), characteristic of aggregates of olivine with an AG crystallographic fabric. These seismic results are the first in situ evidence for this fabric in the upper mantle. We show that its presence can be explained by isotropic horizontal dilation and vertical flattening due to late-stage gravitational collapse and spreading in the top 10 to 20 km of a depleted, mushroom-shaped, superplume head on a horizontal length scale of 1000 km or more. This way, it provides a seismic tool to track plumes long after the thermal effects have ceased.
Precambrian Research, in press available, 43p. Pdf
United States, Iowa
geophysics - seismics
Abstract: Large amplitude aeromagnetic and gravity anomalies over a ~9500 km2 area of northeast Iowa and southeast Minnesota have been interpreted to reflect the northeast Iowa intrusive complex (NEIIC), a buried intrusive igneous complex composed of mafic/ultramafic rocks in the Yavapai Province (1.8-1.7 Ga). Hundreds of meters of Paleozoic sedimentary cover and a paucity of basement drilling have prevented detailed studies of the NEIIC. Long considered, but not proven, to be related to the ~1.1 Ga Midcontinent Rift System (MRS), the NEIIC is comparable in areal extent to the richly mineralized Duluth Complex and is similarly located near the margin of the MRS. New geochronological and geophysical data together support an MRS affinity for the NEIIC. A dike swarm imaged in aeromagnetic data is cut by intrusions of the NEIIC, and a new apatite U-Pb date of ~1170 Ma on one of the dikes thus represents a maximum age for the NEIIC. A minimum age constraint is suggested by (1) large-volume magmatism associated with the MRS that was the last such event to affect the region; and (2) the presence of reversely magnetized dikes, similar in character to MRS-related dikes elsewhere, that cut several intrusions of the NEIIC. The NEIIC is largely characterized by the presence of multiple zoned intrusions, many of which contain large volumes of mafic-ultramafic rocks and have strong geophysical similarities to alkaline intrusive complexes elsewhere, including the MRS-related Coldwell Complex of Ontario. The largest of the zoned intrusions are ~40 km in diameter and are interpreted to have thicknesses of many kilometers. Suspected faults, alignments of intrusions, and intrusive margins tend to be aligned along northwest and northeast trends that match the trends of the Belle Plaine fault zone and Fayette structural zone, both previously interpreted as pre-MRS, possibly lithospheric-scale discontinuities that may have controlled NEIIC emplacement. These interpretations collectively imply notable potential for the NEIIC to host several different types of undiscovered base metal and critical mineral deposits.
Physics of the Earth and Planetary Interiors, Vol. 306, 106543, 9p. Pdf
India
geophysics -seismic
Abstract: The south Indian shield, primarily consisting of Archean cratons and Cretaceous-Tertiary Deccan Volcanic Province (DVP), has undergone several major tectonic episodes during its evolution. The Deccan volcanism at Cretaceous-Tertiary (~65 Ma) is the last major tectono-thermal event, which influenced a substantial part of the south Indian shield. To understand the influence of the Deccan volcanism on the evolution of the south Indian shield, we study the crustal seismic structure of the ~65 Ma Deccan Volcanic Province and the adjacent ~2.6 Ga Eastern Dharwar Craton (EDC), which forms the basement of the volcanic terrain. We calculate teleseismic receiver functions for 18 broadband seismic stations along a ~1000 km long seismological profile that cut across both the EDC and DVP. The analysis and modelling, using H-Vp/Vs stacking and generalized neighbourhood algorithm inversion of the receiver functions show distinct crustal structure (crustal thickness, average composition, shear wave velocity variation, nature of crust-mantle boundary, etc.) across the EDC and DVP. The results clearly indicate that the crustal structure is heterogeneous beneath the DVP compared to a relatively uniform structure below the EDC. Using results from this study along with earlier results, we infer that the present Eastern Dharwar Craton terrain is not affected by any tectono-thermal event for a long geological time, including the Deccan volcanism. Whereas, the present Deccan Volcanic Province is highly affected by the Reunion mantle plume-crust interaction.
Physics of the Earth and Planetary Interiors, Vol. 305, 106504, 15p. Pdf
United States
geophysics - seismics
Abstract: We examined SKS-SKKS splitting intensity discrepancies for phases that sample the lowermost mantle beneath North America, which has previously been shown to exhibit seismic anisotropy using other analysis techniques. We examined data from 25 long-running seismic stations, along with 244 stations of the temporary USArray Transportable Array, located in the eastern, southeastern and western U.S. We identified 279 high-quality SKS-SKKS wave pairs that yielded well-constrained splitting intensity measurements for both phases. Of the 279 pairs, a relatively small number (15) exhibited discrepancies in splitting intensity of 0.4 s or greater, suggesting a contribution to the splitting of one or both phases from anisotropy in the lowermost mantle. Because only a small minority of SK(K)S phases examined in this study show evidence of being affected by lowermost mantle anisotropy, the traditional interpretation that splitting of these phases primarily reflects anisotropy in the upper mantle directly beneath the stations is appropriate. The discrepant pairs exhibited a striking geographic trend, sampling the lowermost mantle beneath the southern U.S. and northern Mexico, while other regions were dominated by non-discrepant pairs. We carried out ray theoretical modeling of simple anisotropy scenarios that have previously been suggested for the lowermost mantle beneath North America, invoking the alignment of post-perovskite due to flow induced by the impingement of the remnant Farallon slab on the core-mantle boundary. We found that our measurements are generally consistent with this model and with the idea of slab-driven flow, but relatively small-scale lateral variations in the strength and/or geometry of lowermost mantle anisotropy beneath North America are also likely present.
Abstract: Knowledge of the crust and lithospheric structure of the Indian sub-continent primarily comes from several active and passive seismic experiments. These studies are i) controlled source, ii) surface wave studies, iii) receiver functions and v) tomographic studies. The results from these studies in the Indian shield have emanated several interesting features that were hitherto unknown. The peninsular, central and north-western part of the shields, Himalayan and Andaman-Nicobar regions have shown that continental collision and extension from the Proterozoic to Recent time has played an important role in formation and geodynamics of these features. The granulites, in the southern granulite terrain, are formed primarily due to the release of the carbonic fluids from the supracrustal rocks of the subduction zone and volcanic arc environment. These were later exhumed from the deep crust during the collision process. In the central Indian shield the Narmada-Son lineament and the central Indian suture are the main features of the crust. In the Narmada region, mafic intrusion in the upper crust appears to have played an important role in shaping the present structural trends. The Central Indian suture is a collision zone developed due to the interaction of the Bastar and Bundelkhand cratons. In the northwesternpart of the India, the Aravalli-Delhi trend is the controlling feature for the tectonics of the region. Demarcation of the various boundaries between different crustal units are marked across the trend, by changes in the dip direction and steeply dipping reflections, cutting across the nearly horizontal reflections at various depths in the crust. Plate tectonics appears to be responsible for generation of this belt. In the crustal block between the Delhi-Aravalli system and the Narmada-Son Lineament, which is running to the south of the Saurashtra peninsula the crust up uplifted by as much as 4 to 6 km as compared to the regions outside these trends. Apart from the deep crustal structure, lithospheric and upper mantle studies till 660km depth have also been conducted in the entire Indian plate using seismological tools e.g. P-to-s and S-to-p receiver function, surface waves dispersion and tomographic studies. The Himalayan region shows the architecture of the under thrusting Indian plate beneath the Tibetan plate in the north and north-west, while the subduction beneath the Burmese arc has been mapped in the eastern part. Further, a number of studies have been conducted in the Andaman-Nicobar Islands to image the subduction of Indian oceanic plate in order to understand the genesis of earthquakes in these regions.
Pure and Applied Geophysics, Vol. 177, 8, pp. 3775-3795. pdf
Europe, Finland
geophysics - seismics
Abstract: The Kokkola-Kymi Deep Seismic Sounding profile crosses the Fennoscandian Shield in northwest-southeast (NW-SE) direction from Bothnian belt to Wiborg rapakivi batholith through Central Finland granitoid complex (CFGC). The 490-km refraction seismic line is perpendicular to the orogenic strike in Central Finland and entirely based on data from quarry blasts and road construction sites in years 2012 and 2013. The campaign resulted in 63 usable seismic record sections. The average perpendicular distance between these and the profile was 14 km. Tomographic velocity models were computed with JIVE3D program. The velocity fields of the tomographic models were used as starting points in the ray tracing modelling. Based on collected seismic sections a layer-cake model was prepared with the ray tracing package SEIS83. Along the profile, upper crust has an average thickness of 22 km average, and P-wave velocities (Vp) of 5.9-6.2 km/s near the surface, increasing downward to 6.25-6.40 km/s. The thickness of middle crust is 14 km below CFGC, 20 km in SE and 25 km in NW, but Vp ranges from 6.6 to 6.9 km/s in all parts. Lower crust has Vp values of 7.35-7.4 km/s and lithospheric mantle 8.2-8.25 km/s. Moho depth is 54 km in NW part, 63 km in the middle and 43 km in SW, yet a 55-km long section in the middle does not reveal an obvious Moho reflection. S-wave velocities vary from 3.4 km/s near the surface to 4.85 km/s in upper mantle, consistently with P-wave velocity variations. Results confirm the previously assumed high-velocity lower crust and depression of Moho in central Finland.
White-Gaynor, A.L., Nyblade, A.A., Durrheim, R., Raveloson, R., van der Meijde, M., Fadel, I., Paulssen, H., Kwadiba, M., Ntibinyane, O., Titus, N., Sitali, M.
Geochemistry, Geophysics, Geosystems, 10.1029/GC008925 20p. Pdf
Africa, South Africa
Geophysics, seismic
Abstract: We report new P and S wave velocity models of the upper mantle beneath southern Africa using data recorded on seismic stations spanning the entire subcontinent. Beneath most of the Damara Belt, including the Okavango Rift, our models show lower than average velocities (?0.8% Vp; ?1.2% Vs) with an abrupt increase in velocities along the terrane's southern margin. We attribute the lower than average velocities to thinner lithosphere (~130 km thick) compared to thicker lithosphere (~200 km thick) immediately to the south under the Kalahari Craton. Beneath the Etendeka Flood Basalt Province, higher than average velocities (0.25% Vp; 0.75% Vs) indicate thicker and/or compositionally distinct lithosphere compared to other parts of the Damara Belt. In the Rehoboth Province, higher than average velocities (0.3% Vp; 0.5% Vs) suggest the presence of a microcraton, as do higher than average velocities (1.0% Vp; 1.5% Vs) under the Southern Irumide Belt. Lower than average velocities (?0.4% Vp; ?0.7% Vs) beneath the Bushveld Complex and parts of the Mgondi and Okwa terranes are consistent with previous studies, which attributed them to compositionally modified lithosphere resulting from Precambrian magmatic events. There is little evidence for thermally modified upper mantle beneath any of these terranes which could provide a source of uplift for the Southern African Plateau. In contrast, beneath parts of the Irumide Belt in southern and central Zambia and the Mozambique Belt in central Mozambique, deep?seated low velocity anomalies (?0.7% Vp; ?0.8% Vs) can be attributed to upper mantle extensions of the African superplume structure.
Abstract: The model, which will incorporate 227 million surface wave measurements, could help with everything from earthquake characterization to neutrino geosciences.
Physics of the Earth and Planetary Interiors, Vol. 308, 106552, 14p. Pdf
Mantle
geophysics - seismics
Abstract: Constraining Earth's bulk composition is fundamental to understanding our planet's formation and evolution. While the lower mantle accounts for a majority of the bulk silicate Earth, it is also the least accessible. As experimental and theoretical mineral physics constraints on mineral elasticity at lower mantle temperatures and pressures have improved, comparisons between predicted seismic velocity and density profiles for hypothesized bulk compositions and 1D seismic models have become commonplace. However, the degree to which a given composition is a better or worse fit than another composition is not always reported, nor are the influences of the assumed temperature profile and other uncertainties discussed. Here we compare seismic velocities and densities for perovskitite, pyrolite, and harzburgite bulk compositions calculated using advanced ab initio techniques to explore the extent to which the associated uncertainties affect our ability to distinguish between candidate compositions. We find that predicted differences between model compositions are often smaller than the influence of temperature uncertainties and therefore these comparisons lack discriminatory power. The inability to distinguish between compositions is largely due to the high sensitivity of seismic properties to temperature accompanied by uncertainties in the mantle geotherm, coupled with diminished sensitivity of seismic velocity to composition toward the base of the mantle. An important exception is the spin transition in (Mg,Fe)O-ferropericlase, which is predicted to give a distinct variation in compressional wave velocity that should distinguish between relatively ferro-magnesian and silica-rich compositions. However, the absence of an apparent spin transition signature in global 1D seismic profiles is a significant unresolved issue in geophysics, and it has important geochemical implications. The approach we present here for establishing discriminatory power for such comparisons can be applied to any estimate of seismic velocities and associated uncertainties, and offers a straightforward tool to evaluate the robustness of model comparisons.
Geophysical Journal International, Vol. 224, 2, pp. 961-972.
Mantle
geophysics - seismics
Abstract: Previous studies have shown that a low viscosity upper mantle can impact the wavelength of mantle flow and the balance of plate driving to resisting forces. Those studies assumed that mantle viscosity is independent of mantle flow. We explore the potential that mantle flow is not only influenced by viscosity but can also feedback and alter mantle viscosity structure owing to a non-Newtonian upper-mantle rheology. Our results indicate that the average viscosity of the upper mantle, and viscosity variations within it, are affected by the depth to which a non-Newtonian rheology holds. Changes in the wavelength of mantle flow, that occur when upper-mantle viscosity drops below a critical value, alter flow velocities which, in turn, alter mantle viscosity. Those changes also affect flow profiles in the mantle and the degree to which mantle flow drives the motion of a plate analogue above it. Enhanced upper-mantle flow, due to an increasing degree of non-Newtonian behaviour, decreases the ratio of upper- to lower-mantle viscosity. Whole layer mantle convection is maintained but upper- and lower-mantle flow take on different dynamic forms: fast and concentrated upper-mantle flow; slow and diffuse lower-mantle flow. Collectively, mantle viscosity, mantle flow wavelengths, upper- to lower-mantle velocities and the degree to which the mantle can drive plate motions become connected to one another through coupled feedback loops. Under this view of mantle dynamics, depth-variable mantle viscosity is an emergent flow feature that both affects and is affected by the configuration of mantle and plate flow.
Physics of the Earth and Planetary Interiors, Vol. 309, di.org/10.1016/ jpepi.2020.106600 20p. Pdf
Mantle
geophysics - seismics
Abstract: Over the last decade, seismological studies have shed new light on the properties of the mantle lithosphere and their physical and chemical origins. This paper synthesizes recent work to draw comparisons between oceanic and continental lithosphere, with a particular focus on isotropic velocity structure and its implications for mantle temperature and partial melt. In the oceans, many observations of scattered and reflected body waves indicate velocity contrasts whose depths follow an age-dependent trend. New modeling of fundamental mode Rayleigh waves from the Pacific ocean indicates that cooling plate models with asymptotic plate thicknesses of 85-95 km provide the best overall fits to phase velocities at periods of 25 s to 250 s. These thermal models are broadly consistent with the depths of scattered and reflected body wave observations, and with oceanic heat flow data. However, the lithosphere-asthenosphere velocity gradients for 85-95 km asymptotic plate thicknesses are too gradual to generate observable Sp phases, both at ages less than 30 Ma and at ages of 80 Ma or more. To jointly explain Rayleigh wave, scattered and reflected body waves and heat flow data, we propose that oceanic lithosphere can be characterized as a thermal boundary layer with an asymptotic thickness of 85-95 km, but that this layer contains other features, such as zones of partial melt from hydrated or carbonated asthenosphere, that enhance the lithosphere-asthenosphere velocity gradient. Beneath young continental lithosphere, surface wave constraints on lithospheric thickness are also compatible with the depths of lithosphere-asthenosphere velocity gradients implied by converted and scattered body waves. However, typical steady-state conductive models consistent with continental heat flow produce thermal and velocity gradients that are too gradual in depth to produce observed converted and scattered body waves. Unless lithospheric isotherms are concentrated in depth by mantle upwelling or convective removal, the presence of an additional factor, such as partial melt at the base of the thermal lithosphere, is needed to sharpen lithosphere-asthenosphere velocity gradients in many young continental regions. Beneath cratons, numerous body wave conversions and reflections are observed within the thick mantle lithosphere, but the velocity layering they imply appears to be laterally discontinuous. The nature of cratonic lithosphere-asthenosphere velocity gradients remains uncertain, with some studies indicating gradual transitions that are consistent with steady-state thermal models, and other studies inferring more vertically localized velocity gradients.
Physics of the Earth and Planetary Interiors, Vol. 308, di.org/10.1016 /jpepi.2020. 106552 14p. Pdf
Mantle
geophysics - seismics
Abstract: Constraining Earth's bulk composition is fundamental to understanding our planet's formation and evolution. While the lower mantle accounts for a majority of the bulk silicate Earth, it is also the least accessible. As experimental and theoretical mineral physics constraints on mineral elasticity at lower mantle temperatures and pressures have improved, comparisons between predicted seismic velocity and density profiles for hypothesized bulk compositions and 1D seismic models have become commonplace. However, the degree to which a given composition is a better or worse fit than another composition is not always reported, nor are the influences of the assumed temperature profile and other uncertainties discussed. Here we compare seismic velocities and densities for perovskitite, pyrolite, and harzburgite bulk compositions calculated using advanced ab initio techniques to explore the extent to which the associated uncertainties affect our ability to distinguish between candidate compositions. We find that predicted differences between model compositions are often smaller than the influence of temperature uncertainties and therefore these comparisons lack discriminatory power. The inability to distinguish between compositions is largely due to the high sensitivity of seismic properties to temperature accompanied by uncertainties in the mantle geotherm, coupled with diminished sensitivity of seismic velocity to composition toward the base of the mantle. An important exception is the spin transition in (Mg,Fe)O-ferropericlase, which is predicted to give a distinct variation in compressional wave velocity that should distinguish between relatively ferro-magnesian and silica-rich compositions. However, the absence of an apparent spin transition signature in global 1D seismic profiles is a significant unresolved issue in geophysics, and it has important geochemical implications. The approach we present here for establishing discriminatory power for such comparisons can be applied to any estimate of seismic velocities and associated uncertainties, and offers a straightforward tool to evaluate the robustness of model comparisons.
Earth and Planetary Science Letters, Vol. 554, doi:10.1016/j.epsl.2020.116685
Mantle
geophysics - seismic
Abstract: Here, we provide a reappraisal of potential LLSVPs compositions based on an improved mineralogical model including, for instance, the effects of alumina. We also systematically investigate the effects of six parameters: FeO and Al2O3 content, proportion of CaSiO3 and bridgmanite (so that the proportion of ferropericlase is implicitly investigated), Fe3+/?Fe and temperature contrast between far-field mantle and LLSVPs. From the 81 millions cases studied, only 79000 cases explain the seismic observations. Nevertheless, these successful cases involve a large range of parameters with, for instance, FeO content between 12--25~wt\% and Al2O3 content between 3--17~wt\%. We then apply a principal component analysis (PCA) to these cases and find two robust results: (i) the proportion of ferropericlase should be low (<6vol\%); (ii) the formation of Fe3+-bearing bridgmanite is much more favored than other iron-bearing phases. Following these results, we identify two end-member compositions, Bm-rich and CaPv-rich, and discuss their characteristics. Finally, we discuss different scenarios for the formation of LLSVPs and propose that investigating the mineral proportion produced by each scenario is the best way to evaluate their relevance. For instance, the solidification of a primitive magma ocean may produce FeO and Al2O3 content similar to those suggested by our analysis. However, the mineral proportion of such reservoirs is not well-constrained and may contain a larger proportion of ferropericlase than what is allowed by our results.
Journal of Geophysical Research, Solid Earth, https://eartharxiv.org /repository/ view/1694/ 21p. Pdf
Mantle
Geophysics - seismics
Abstract: The inner core is a Moon?sized ball of solid metal at the very center of the Earth. Vibrations from earthquakes move faster through the inner core if they travel parallel to Earth's axis (the line between the North and South Poles) than if they travel parallel to the Equator. This probably means that the grains of metal in the inner core are themselves aligned with Earth's axis. Previous studies determined that this alignment likely happened after the inner core had formed, but those experiments were done on pure iron, whereas the inner core is mostly iron but also contains other elements. We did experiments at high pressures and temperatures on a more realistic core metal containing iron, nickel, and silicon. We found that this metal would be much stronger than pure iron at inner core pressures and temperatures; it is still possible for it to produce a north-south alignment, but it is much more difficult for it to do so. This could mean that the alignment occurred while the inner core was forming (rather than afterward), which might change how we think about the forces present in the deep Earth today.
Abstract: New evidence from seismic tomography reveals a unique mineral fabric restricted to the thick mantle lithosphere beneath ancient continental cratons, providing an important clue to the formation of these prominent and influential features in Earth’s geological history. Olivine, the dominant mineral of Earth’s upper mantle, has elastic properties that differ along its three crystallographic axes, and preferential alignment of individual olivine grains during plastic deformation can affect the bulk nature of seismic-wave propagation. Surface-wave tomography has shown that over most of Earth, deformation of the mantle lithosphere has oriented olivine crystals with the fast axis in the horizontal plane, but at depths centered at ?150 km within cratonic continental-lithosphere roots, the fast crystallographic axis is preferentially aligned vertically. Because of the high viscosity of the cratonic roots, this fabric is likely to be a vestige from craton formation. Geochemical and petrological studies of upper-mantle garnet-peridotite nodules demonstrate that the cratonic mantle roots are stabilized by their reduced density, which was caused by melt removal at much shallower depths than those from which the nodules were subsequently extracted. The mineral fabric inferred from surface-wave tomography suggests that horizontal shortening carried the depleted zone downward after the melt-depletion event to form the thick continental roots, stretching the depleted material in the vertical dimension by pure shear and causing the fast crystallographic axis to be aligned vertically. This seismological fabric at ?150 km is evidence of the shortening event that created the cratonic roots.
Journal of South American Earth Sciences, Vol. 107, 103088, 17p. Pdf
South America, Venezuela
Geophysics - seismics
Abstract: As part of the lithosphere, the crust represents Earth's rigid outer layer. Some of the tools to study the crust and its thickness are wide-angle seismic studies. To date, a series of seismic studies have been carried out in Venezuela to determine in detail the crustal thickness in the southern Caribbean, in the region of the Caribbean Mountain System in northern Venezuela, as well as along the Mérida Andes and surrounding regions. In this study, a review of the wide-angle seismic data is given, incorporating new data from the GIAME project for western Venezuela, resulting in a map of Moho depth north of the Orinoco River, which serves as the basis for future integrated models. Differences in Moho depths from seismic data and receiver function analysis are discussed. From the Caribbean plate, Moho depth increases from 20 to 25 km in the Venezuela Basin to about 35 km along the coast (except for the Falcón area where a thinning to less than 30 km is observed) and 40-45 km in Barinas - Apure and Guárico Basins, and Guayana Shield, respectively. Values of more than 50 km are observed in the Maturín Basin and in the southern part of the Mérida Andes.
Abstract: We explore the shear-wave lithospheric velocity structure of the Eastern Caribbean and Northern Venezuela using ambient noise tomography with stations deployed around the study area. We construct cross-correlation functions from continuous seismic records, and measure phase velocities of fundamental-mode Rayleigh waves. These velocities are further projected onto 0.6°x0.6° phase velocity grids for each period between 5 s and 50 s. The pseudo-dispersion curve at each grid point is inverted for 1D shear velocity profiles by using a Markov Chain Monte Carlo scheme. The interpolated 3D velocity model shows that the mean shear velocity of the Eastern Caribbean lithospheric mantle is lower than the global average, which is in agreement with values reported in other large igneous provinces. We interpret that low velocities in the lithospheric keel are associated with an anomalous composition and/or an elevated thermal state; this gives the Caribbean plate a high buoyancy that determines the subduction polarities in the region. The results also indicate that: (a) the mantle beneath Northern Venezuela retains compositional anomalies related to extension processes of different ages; (b) the overriding of the Caribbean plate by the Great Antilles arc seems to be much slower than previously suggested; and (c) the localized volcanism in the center of the Lesser Antilles arc is related to asthenospheric flow through the tear induced on the subducted slab by major strike-slip faults.
Abstract: The tectonic-geodynamic characteristics of the North African-Arabian region are complicated by interaction of numerous factors. To study this interaction, we primarily used the satellite gravity data (retracked to the Earth's surface), recognized as a powerful tool for tectonic-geodynamic zonation. The applied polynomial averaging of gravity data indicated the presence of a giant, deep quasi-ring structure in the Eastern Mediterranean, the center of which is located under the Island of Cyprus. Simultaneously, the geometrical center of the revealed structure coincides with the Earth's critical latitude of 35?. A quantitative analysis of the obtained gravitational anomaly made it possible to estimate the depth of the upper edge of the anomalous body as 1650?1700 km. The GPS vector map coinciding with the gravitational trend indicates counterclockwise rotation of this structure. Review of paleomagnetic data on the projection of the discovered structure into the earth's surface also confirms its counterclockwise rotation. The analysis of the geoid anomalies map and seismic tomography data commonly approve presence of this deep anomaly. The structural and geodynamic characteristics of the region and paleobiogeographic data are consistent with the proposed physical-geological model. Comprehensive analysis of petrological, mineralogical, and tectonic data suggests a relationship between the discovered deep structure and near-surface processes. The revealed geological deep structure sheds light on specific anomalous effects in the upper crustal layer, including the high-intensity Cyprus gravity anomaly, counterclockwise rotation of the Mesozoic terrane belt, configuration of the Sinai plate, and the asymmetry of sedimentary basins along the continental faults.
Nature Communications, Vol. 12, 1, doi:10.10.1038/ s41467-021-21657-8 8p. Pdf
Mantle
geophysics - seismic
Abstract: ediments play a key role in subduction. They help control the chemistry of arc volcanoes and the location of seismic hazards. Here, we present a new model describing the fate of subducted sediments that explains magnetotelluric models of subduction zones, which commonly show an enigmatic conductive anomaly at the trenchward side of volcanic arcs. In many subduction zones, sediments will melt trenchward of the source region for arc melts. High-pressure experiments show that these sediment melts will react with the overlying mantle wedge to produce electrically conductive phlogopite pyroxenites. Modelling of the Cascadia and Kyushu subduction zones shows that the products of sediment melting closely reproduce the magnetotelluric observations. Melting of subducted sediments can also explain K-rich volcanic rocks that are produced when the phlogopite pyroxenites melt during slab roll-back events. This process may also help constrain models for subduction zone seismicity. Since melts and phlogopite both have low frictional strength, damaging thrust earthquakes are unlikely to occur in the vicinity of the melting sediments, while increased fluid pressures may promote the occurrence of small magnitude earthquakes and episodic tremor and slip.
Abstract: The Deccan volcanic province (DVP) witnessed a massive outpouring of flood basalts of ?2 million km3 volume, at ?65 Ma, in less than a Myr. The volcanic eruption is concomitant with crustal extension, lithospheric thinning and magma influx beneath the major rift systems namely Cambay, Narmada, and Kutch. In this study, we investigate the anisotropic and isotropic variations within the crust and upper mantle beneath the northwestern DVP by estimating the shear wave velocity (VSV, VSH, and VSoigt) and radial anisotropy (?oigt) models using the Surface Wave Tomography technique. A joint inversion of the regionalized Rayleigh and Love wave group velocities is performed, using the genetic algorithm approach. Our results reveal different intracrustal layers, lid, and a low?velocity zone (LVZ). This LVZ comprises of a uniform asthenospheric low?velocity layer (LVL) of average VSV 4.44 km/s and VSH 4.47 km/s, and another LVL below, of average VSV 4.45 km/s and VSH 4.41 km/s. Furthermore, the LVZ represents a negative anomaly with reference to different global models (AK135, STW105, PREM, and S2.9EA). A negative ?oigt is observed in the LVZ, indicating dominance of vertical flow. This could be related to presence of partials melts, volatile materials and/or a thermal anomaly. We also identified the Moho (?34-40 km) and lithosphere?asthenosphere boundary (?84-123 km). The low VS values, negative ?oigt and a thin lithosphere (?84 km) in the vicinity of Gulf of Cambay affirm the presence of a plume head beneath it, in concurrence with the hypothesis of Indian Plate?Reunion plume interaction.
Physics of the Earth and Planetary Interiors, Vol. 316, 106710 13p. Pdf
Mantle
geophysics - seismic
Abstract: While the lateral limits of tectonic plates are well mapped by seismicity, the bottom boundary of the lithosphere, the uppermost rigid layer of the Earth comprising both crust and shallow mantle, remains elusive. The lithosphere is usually viewed as consisting of cold, rigid, internally undeformed blocks that translate coherently. The base of the lithosphere, designated as the lithosphere-asthenosphere boundary (LAB), may thus be characterised by different physical fields: temperature, viscosity, strain rate and velocity. The LABs as defined by these different fields are investigated here using thermo-mechanical models of plate and upper mantle dynamics, either in a transient subduction or in a steady-state plate-driven set-up. Mantle material is modelled as homogeneous in composition with a viscosity that depends on temperature, pressure and strain rate. In such a system, the thermo-mechanical transition between lithosphere and asthenosphere occurs over a finite depth interval in temperature, strain rate and velocity. We propose that the most useful dynamical LAB is defined as the base of a “constant-velocity” plate (i.e. the material translating at constant horizontal velocity). The bottom part of this plate deforms at strain rates comparable to those in the underlying asthenosphere mantle: the translating block is not fully rigid. Thermal structure exerts a major control on this dynamical LAB, which deepens with increasing plate age. However, the surface plate velocity, the asthenospheric flow geometry and magnitude also influence the depth of the dynamical LAB, as well as the thickness of the deformed region at the base of the constant-velocity plate. The mechanical transitions from lithosphere to asthenosphere adjust when mantle dynamics evolves. The dynamical and thermo-mechanical LABs occur within a thermal lithosphere-asthenophere gradual transition, similar to the one imaged by geophysical proxies. The concept of a constant-velocity plate can be extended to a constant-velocity subducting slab, which also deforms at its borders and drags the surrounding mantle. This framework is relevant to quantify mass transport within the Earth's mantle.
Geochemistry, Geophysics, Geosystems, 10.1029/2020G C009512 17p. Pdf
Europe, Greenland
geophysics - seismic
Abstract: Measurements of seismic anisotropy (the direction-dependent variation in seismic wavespeed) provide useful information about the orientation of deformation in the Earth. We measured seismic anisotropy using shear waves refracted through the outer core and recorded by stations in Greenland. Due to new stations and data, this study includes more measurements of the effects of anisotropy than previously possible. We show that a model with two layers of anisotropy explains dominant patterns in the fast vibration direction of the shear waves as a function of the angle at which they approach each station. We suggest that the shallow layer reflects coherent deformation in the continental lithosphere of Greenland due to its history of plate collisions and that the lower layer reflects deformation in the asthenospheric mantle induced by the motion of the plate above or a second layer of lithospheric anisotropy.
Abstract: The eastern Indian shield consists of Archaean Singhbhum Craton and Proterozoic Chhotanagpur Gneissic Complex sandwiching the Singhbhum Mobile Belt. Since the cratonization of the Singhbhum Craton in Archaean, the growth of the eastern Indian shield took place in time and space through tectono-magmatic processes. The stability of cold and thick lithosphere is fundamental to long-term survival of cratons, whereas the geophysical studies have detected the lithosphere-asthenosphere boundary (LAB) under the eastern Indian shield at depths too shallow to be called stable. We analysed the terrestrial Bouguer gravity anomaly, and satellite-based free-air anomaly, geoid undulation, and elevation data to ascertain the 2D lithospheric density structure across the region. Our density model illustrates that the density inhomogeneity exists in the crust across the three tectonic domains of the eastern Indian shield. The derived crustal model shows an upper and lower crustal density variation from 2740 to 2770 kg/m3, and from 2930 to 2940 kg/m3, respectively, and a reasonably smooth Moho at 37-41 km depth. Towards the north, the Moho undulates from 40 to 43 km under the foreland Ganga basin, whereas in the south, it varies from 38 to 30 km under the Eastern Ghats Mobile Belt and lastly moves to ~20 km in the Bay of Bengal. In the southern part of the Singhbhum Craton, an undissipated lithospheric mantle root is found at a depth of ~150 km. Otherwise, the LAB shallows to ~132 km in the northern Singhbhum Craton and Singhbhum Mobile Belt and then thickens to about 135-140 km depth beneath the Chhotanagpur Gneissic Complex. The foreland Ganga basin toward the extreme north is characterized by a more in-depth LAB lying at a depth of over 200 km. The LAB, in the Bay of Bengal, is at a depth of 112-125 km, except for the Kolkata coast (135 km). Moderate crustal density difference in various crustal domains, as well as an almost smooth crust-mantle boundary at 37-40 km depth, suggests the effect of substantial mafic-ultramafic crustal intrusion and together with the thin (135-140 km) lithosphere reinforces the evidence of thermo-chemical processes that controlled the lithospheric modification in the eastern Indian shield.
Abstract: We investigate the structure of the continental lithosphere (tectonic plate) in regions that have had negligible tectonic activity, such as mountain building, for the past 500 million years. The internal structure of the lithosphere in these regions can be indicative of the ancient processes that first formed continents. Due to challenges in methodology, layering within the upper 150 km of the continental lithosphere is poorly understood. We carefully process earthquake data to avoid problems that previous studies encountered. We observe layering in 50% of the ancient continental regions. Most of this layering can be explained by the presence of minerals that have lower seismic velocities than the surrounding rock because they have been altered by fluids during the formation of the continent. In regions closer to more recent tectonic activity, some layering has stronger seismic velocity decreases, indicating the effects of more recent alteration. We also find that layering is more prevalent in the continental regions that last experienced tectonic activity no later than 1.6 billion years ago. This corresponds with a global transition in the depth to which the subducting lithosphere carries fluids into the mantle, indicating that subduction has a key role in generating layering in the ancient continental lithosphere.
Abstract: From west to east along the Sunda-Banda arc, convergence of the Indo-Australian plate transitions from subduction of oceanic lithosphere to arc-continent collision. This region of eastern Indonesia and Timor-Leste provides an opportunity for unraveling the processes that occur during collision between a continent and a volcanic arc, and it can be viewed as the temporal transition of this process along strike. We collected a range of complementary geological and geophysical data to place constraints on the geometry and history of arc-continent collision. Utilizing ?4 yr of new broadband seismic data, we imaged the structure of the crust through the uppermost mantle. Ambient noise tomography shows velocity anomalies along strike and across the arc that are attributed to the inherited structure of the incoming and colliding Australian plate. The pattern of anomalies at depth resembles the system of salients and embayments that is present offshore western Australia, which formed during rifting of east Gondwana. Previously identified changes in geochemistry of volcanics from Pb isotope anomalies from the inner arc islands correlate with newly identified velocity structures representing the underthrusted and subducted Indo-Australian plate. Reconstruction of uplift from river profiles from the outer arc islands suggests rapid uplift at the ends of the islands of Timor and western Sumba, which coincide with the edges of the volcanic-margin protrusions as inferred from the tomography. These findings suggest that the tectonic evolution of this region is defined by inherited structure of the Gondwana rifted continental margin of the incoming plate. Therefore, the initial template of plate structure controls orogenesis.
Abstract: The central region of the Kaapvaal craton is relatively understudied in terms of its lithospheric mantle architecture, but is commonly believed to be significantly impacted by post-Archean magmatism such as the ca. 2056 Ma Bushveld large igneous event. We investigate a collection of 17 eclogite xenoliths from the Cretaceous Palmietfontein kimberlite at the Western Limb of the Bushveld Complex for their mineralogical compositions (major and trace elements, plus Fe3+ contents), as well as stable oxygen and radiogenic Pb isotopic compositions to gain further insights into the nature and evolution of the central Kaapvaal cratonic mantle lithosphere. New U/Pb age determinations on mantle-derived zircon yield a magma emplacement age of ca. 75 Ma for the Palmietfontein Group-1 kimberlite, which means that the entrained eclogite xenoliths may record a protracted metasomatic history from the Proterozoic through to most of the Phanerozoic eon. Garnet ?18O values of up to 6.9‰ and positive Eu anomalies for the bulk rocks suggest seawater-altered oceanic crustal protoliths for the Palmietfontein eclogite xenolith suite, which is typical for the eclogitic components of the Kaapvaal root and other cratonic mantle sections worldwide. However, several features of the Palmietfontein eclogites are commonly not observed in other mantle-derived eclogite xenolith suites. Firstly, the samples studied yield relatively low equilibration pressures and temperatures between 2.7 and 4.5 GPa and 740-1064 °C, indicating a relatively shallow residence between 90 and 150 km depths. Secondly, euhedral coarse amphibole is present in several eclogite nodules where it is in equilibrium with ‘touching’ garnet, supporting eclogite residence within the amphibole stability field at uppermost lithospheric mantle conditions. Thirdly, primary omphacitic clinopyroxene is often overgrown by diopside, and is significantly enriched in incompatible trace elements. The clinopyroxene is also characterized by elevated 206Pb/204Pb of 17.28-19.20 and 207Pb/204 Pb of 15.51-16.27, and these Pb isotopic compositions overlap with those of Mesozoic Group-2 kimberlites from the Kaapvaal craton. Our results show that eclogites reside at ~85 km depth beneath the central Kaapvaal craton as part of a layer that corresponds to an approximately 50 km thick seismically-detected mid-lithospheric discontinuity. Mid-lithospheric discontinuities have been interpreted as metasomatic fronts formed by focussed crystallization of hydrous mineral phases from enriched volatile-bearing melts, and as such the strongly overprinted amphibole-bearing eclogite xenoliths from Palmietfontein may represent a physical expression of such seismically anomalous metasomatic layer at mid-lithospheric depth. Our Pb isotope data suggest that the focussed metasomatism can be attributed to volatile-rich melts reminiscent of potassic Group-2 kimberlites, which have been invoked in MARID-style metasomatic overprinting of the lower lithospheric mantle beneath the western Kaapvaal craton. However, the relatively low fO2 recorded by the Palmietfontein eclogites (minimum FMQ-4.5) suggests that the metasomatism at mid-lithospheric depth was less protracted compared to the more intensive and oxidizing metasomatism typically observed near the base of cratonic mantle roots. While it is possible that Proterozoic magmatic events were responsible for the focussed mid-lithospheric metasomatism of the Kaapvaal mantle, on the basis of the Pb isotope constraints the Palmietfontein eclogites were most likely overprinted during ca. 120 Ma Group-2 kimberlite magmatism.
Abstract: Structure within the Earth is best studied in three dimensions and using several coincident overlays of diverse information with which one can best see where unusual properties match up. Here we use regional surfaces causing discontinuities in seismic waves a few hundred kilometers deep in the Earth, intersected and thus calibrated by rebuilt rock columns using rare rock samples erupted to the surface in two locations. Electrically conductive regions can be mapped using natural (magnetotelluric) currents. East- and west-dipping seismic discontinuity surfaces match surface structures that developed about 1.8 billion years ago marginal to the Superior crustal block. Surfaces dipping to the southeast and northwest match some boundaries between crustal blocks that are over 2.5 billion years old, but many such crustal boundaries trend more east-west. Conductive rocks appear more commonly above these discontinuity surfaces where gas-rich fluids apparently flowed and that the discontinuities somehow filtered these fluids. The mismatch in orientation and dip between the most ancient deep and exposed structures suggests that plate tectonic processes operating today differed earlier than 2.5 billion years ago.
Geochemistry, Geophysics, Geosystems, 10.1029/2021GC009660 11p. Pdf
United States, Kansas
geophysics- seismic
Abstract: Wastewater injection tied to oil and gas development has induced earthquakes across a broad swath of the central US. These earthquakes have almost exclusively occurred in the crystalline rocks, many kilometers below the main disposal zone in this portion of the country, the Arbuckle Group aquifer. However, the hydrologic properties of these deep crystalline rocks are not well known and are needed to accurately model pressure transmission from injection wells to faults in the deep basement and related earthquake hazards. Newly compiled pressure data, from wastewater disposal wells in Kansas, provide an opportunity to constrain these properties. In this study, we construct a detailed, three-dimensional geological model for the Arbuckle and basement, based on data from >400 wells covering south-central Kansas. We use the model to simulate injection and pressures from more than 300 wells. The model indicates that Arbuckle pressures increased by 1.1 MPa in high injection rate areas and overpressures of <0.1 MPa may be the cause of seismicity in the basement. The simulation results also yield the likely range in hydrologic properties for the crystalline basement and suggests large-scale properties of the basement are enhanced by hidden networks of faults and fractures.
Geochemistry, Geophysics, Geosystems, 10.1029/2021GC009710 20p. Pdf
United States, Canada
geophysics- seismic
Abstract: We probe the properties of upper mantle rocks beneath northeastern North America using the observations of seismic waves from distant earthquakes. We examine signals of converted P-S waves that originate from locations of rapid vertical or directional changes in seismic velocities. These abrupt velocity boundaries are thought to originate from rock deformation, variations in composition, temperature, or melt content. The sharp transitions detectable by this method are compositionally more plausible within the cold tectonic plate than within the hot convecting asthenosphere. Previous studies in this region that analyzed the same type of seismic data report boundaries with sharp downward reduction in seismic velocities between the depths of 60 and 100 km. Their widespread distribution and local consistency with seismic velocity models was used as evidence of them marking the transition between the cold tectonic plate and the hot convecting asthenosphere. Here we expand our search to other types of boundaries and find numerous examples at much greater depths (down to ?185 km). These deeper boundaries primarily reflect changes in directional variation of seismic velocities (anisotropy). The distribution of our deep boundaries broadly agrees with lithospheric thickness estimates in global upper mantle models that consider seismic, gravity, and heat flow data.
Communications Earth & Environment, Vol. 2, 10p. Pdf doi:10.1038/s43247-021-00276-7.
Australia
geophysics - seismics
Abstract: The Australian continental crust preserves a rich geological history, but it is unclear to what extent this history is expressed deeper within the mantle. Here an investigation of Quasi-Love waves is performed to detect scattering of seismic surface waves at mantle depths (between 100-200?km) by lateral gradients in seismic anisotropy. Across Australasia 275 new observations of Quasi-Love waves are presented. The inferred scattering source and lateral anisotropic gradients are preferentially located either near the passive continental margins, or near the boundaries of major geological provinces within Australia. Pervasive fossilized lithospheric anisotropy within the continental interior is implied, on a scale that mirrors the crustal geology at the surface, and a strong lithosphere that has preserved this signal over billions of years. Along the continental margins, lateral anisotropic gradients may indicate either the edge of the thick continental lithosphere, or small-scale dynamic processes in the asthenosphere below.
Abstract: The two most abundant minerals in the Earth’s lower mantle are bridgmanite and ferropericlase. The bulk modulus of ferropericlase (Fp) softens as iron d-electrons transition from a high-spin to low-spin state, affecting the seismic compressional velocity but not the shear velocity. Here, we identify a seismological expression of the iron spin crossover in fast regions associated with cold Fp-rich subducted oceanic lithosphere: the relative abundance of fast velocities in P- and S-wave tomography models diverges in the?~1,400-2,000 km depth range. This is consistent with a reduced temperature sensitivity of P-waves throughout the iron spin crossover. A similar signal is also found in seismically slow regions below?~1,800 km, consistent with broadening and deepening of the crossover at higher temperatures. The corresponding inflection in P-wave velocity is not yet observed in 1-D seismic profiles, suggesting that the lower mantle is composed of non-uniformly distributed thermochemical heterogeneities which dampen the global signature of the Fp spin crossover.
Journal of the Geological Society of India, Vol. 97, 10, pp. 1169-1189.
India
geophysics - seismics
Abstract: The article reviews the history and accomplishments of CSIR-NGRI over the past 60 years, related to elucidating the seismic structure of the crust and lithospheric mantle of the Indian shield. Extensive investigations have been carried out in diverse geological and tectonic provinces of India, employing seismic reflection, refraction/wide-angle reflection and passive seismology to decipher (a) the evolution of the Indian plate through geological time, (b) hazard and its mitigation and (c) accumulation and disposition of natural resources. These endeavours entailed the application and development of state-of-the-art methodologies. Synthesis of the results from active and passive seismology reveals that the thickness of the crust varies between 28 and 65 km in the Kachchh and Aravalli regions respectively, consistent with their evolutionary histories. The thickest crust is observed in the western Dharwar craton (WDC) and the shallowest lies in the west coast. The crust in the shield region is mostly thicker, while it is thin beneath the rift zones. Results from coincident reflection and wide-angle seismic reflection studies broadly suggest a three-layered crust with magmatic underplating. Interestingly, the seismic sections traversing the Aravalli fold belt, central Indian suture zone, Dharwar craton and Southern Granulite Terrain (SGT) depict paleo-collision and subduction environments. The diverse character of the Moho, crustal fabrics and structure in different geological provinces indicate that contrasting tectonic environments might have influenced their evolution and support the hypothesis that plate tectonic processes were operative since Neoarchean. The thickness of the lithosphere estimated from receiver functions varies from 80 to 140 km. An undulation in the Lithosphere Asthenosphere Boundary reveals evidence for a flexure on a regional scale, owing to the continental collision of the Indian and Asian plates. However, the lithospheric thickness derived from surface wave dispersion studies is somewhat larger, ranging from 100 to 250 km, with some body wave tomographic studies suggesting it to be ?400 km, in consonance with the concept of Tectosphere. The thickness values derived from both the methods agree at a few locales such as the Eastern Dharwar Craton, SGT, Cambay, Singhbhum and western DVP. However, a broad disagreement prevails in WDC and northern part of the Indian shield where surface wave tomography reveals the thickness of lithosphere to be 140 to 200 km.
Abstract: New results are presented on the features of the deep velocity structure of two of the three main tectonic blocks that make up the Kola region-Murmansk and Belomorskii-by the P receiver function technique. The research is based on data from the broadband seismic stations Teriberka and Kovda. The results are compared with the models obtained by mutual inversion of PRF and SRF using the data from the stations Apatity and Lovozero. It is shown that the crust has a two-layer structure with the border at a depth of 11 km under the Murmansk block and at a depth of 15 km under the Kola and Belomorskii blocks. The crust thickness of the Murmansk, Belomorskii, and Kola blocks are 35, 33, and 40 km, respectively. The presence of the MLD was revealed in all tectonic structures analyzed for the first time, with a top at a depth of about 70 km for the Murmansk and Belomorskii blocks and 90 km for the Kola block and a bottom at 130-140 km for all structures.
Journal of Asian Earth Sciences, in press available 10.1016/j.jseaes.2021.105092 99 p. Pdf
India
geophysics - seismics
Abstract: The Dharwar Craton of southern India is an important stable cratonic province of the world with complex geology and tectonic settings. Extensive studies provide insights of crustal velocity structure for the tectonic and geodynamic evolution of this Archean craton. This region has experienced several tectonically disturbed zones like Chitradurga Shear Zone (CSZ), Bababudan Shear Zone (BSZ) and Closepet Granites (CG). We have developed a comprehensive geologically plausible tectonic model using both P- and S-wave velocity structures to image major structural elements like shear zones and decipher the compositional distinctions of different rock assemblages of Western Dharwar Craton (WDC) and Eastern Dharwar Craton (EDC) part using 3-C wide-angle seismic data acquired along the 200-km long Perur-Chikmagalur deep seismic profile. The tectonic model show large compositional changes of subsurface rocks with anomalous high , , , Poisson’s ratio () and density () forming a major tectonic divide or suture zone called CSZ between EDC and WDC blocks. Significant crustal thinning (37-41 km) is observed due to Moho upwarping towards the Neo-Archean EDC block mainly composed of felsic granites and granodiorites. The WDC block show relatively thick crust (48-50 km) due to mafic underplating and mantle plume activity below CSZ forming Meso-Archean greenschist-facies-gneisses with dominant mafic/ultra-mafic compositions. Hence, crustal velocity, density, heat-flow, geology and geochronology studies support a plume-arc model with evidence of thick magmatic underplating of the lower-crust, complex subduction and development of highly strained shear zones like CSZ as suture juxtaposing EDC and WDC blocks.
Abstract: The mantle transition zone (MTZ) of the Earth lies between 410 and ?1000 km in depth and has a key role in mantle convection processes. In particular, the discontinuity at 660 km and its associated endothermic mineralogical transformation can slow or inhibit the passage of matter between the upper and lower mantle. The MTZ thus acts as a boundary layer within the mantle. The depth variations of radial and azimuthal seismic anisotropies enable the detection of boundary layers within the mantle. However, the 3D imaging is difficult due to the lack of sensitivity of surface waves of fundamental modes, and the poor global coverage of this depth range by body-wave data. We present a new 3D general anisotropy model (both radial and azimuthal anisotropies) of the mantle down to 1200 km in depth using surface-wave overtone datasets. We find that there is little seismic anisotropy in most of the MTZ, except below subduction zones around the Pacific Ocean and, more surprisingly, in a large area beneath eastern Eurasia where the Pacific subducting plate is stagnant. Seismic anisotropy is usually associated with intense deformation processes but also possibly to water transportation or to fine layering. This significant anisotropy in this part of MTZ might reveal a large water ‘reservoir’ associated with hydrous minerals or a strong stratification. It reflects a complex history beneath central Asia, where the Tethys, Izanagi and Pacific plates appear to have strongly interacted during the last 100 My, having subducted in orthogonal directions under the Asian continent, with the Tethys plate descending into the lower mantle, and the Izanagi plate remaining stagnant in the MTZ. The Asian continent is the only region in the world where subducting slabs originating from different plates can interact. This unique slab distribution might explain why some plates descend while others remain in the lower transition zone.
Abstract: The rifted continental margins of Mozambique provide excellent examples of continental passive margins with a significant structural variability associated with magmatism and inheritance. Despite accumulated knowledge, the tectonic structure and nature of the crust beneath the South Mozambique Coastal Plain (SMCP) are still poorly known. This study interprets high-resolution seismic reflection data paired with data from industry-drilled wells and proposes a structural model of the Limpopo transform margin in a magma-rich context. Results indicate that the Limpopo transform margin is characterized by an ocean-continent transition that links the Beira-High and Natal valley margin segments and represents the western limit of the continental crust, separating continental volcano-sedimentary infilled grabens from the oceanic crust domain. These basins result from the emplacement of the Karoo Supergroup during a Permo-Triassic tectonic event, followed by an Early Jurassic tectonic and magmatic event. This latter led to the establishment of steady-state seafloor spreading at ca.156 Ma along the SMCP. A Late Jurassic to Early Cretaceous event corresponds to formation of the Limpopo transform fault zone. Which accommodated the SSE-ward displacement of Antarctica with respect to Africa. We define a new type of margin: the magma-rich transform margin, characterized by the presence of voluminous magmatic extrusion and intrusion coincident with the formation and evolution of the transform margin. The Limpopo transform fault zone consists of several syn-transfer and -transform faults rather than a single transform fault. The intense magmatic activity was associated primarily with mantle dynamics, which controlled the large-scale differential subsidence along the transform margin.
Geophysical Research Letters, Vol. 48, e2021GLo94784
Mantle
geophysics -seismics
Abstract: Earth's inner core (IC) anisotropy-the directional dependence of seismic wave speed in the IC-contains essential information of deep Earth's structure and dynamics. It results from a preferred alignment of iron crystals related to the formation and post-formation dynamics of the IC. Many studies have investigated the IC anisotropy observed for compressional waves. In contrast, possible anisotropy for the inner-core shear waves remains elusive. This study presents a new class of inner-core shear-wave anisotropy observations based on recent advances in earthquake coda-correlation wavefield. We find that the coda-correlation feature I2-J, sensitive to the inner-core shear-wave speed, exhibits variable timing and amplitude for sampling the IC in different directions. Quantitatively, inner-core shear waves travel faster for at least ?5 s in directions oblique to the Earth's rotation axis than directions parallel to the equatorial plane. The simplest and most plausible explanation for our observations is the inner-core shear-wave anisotropy with a strength of ?0.8% or higher. We can rule out at least one of the body-centered-cubic iron models in the IC, although the other models are not distinguishable.
Watremez, L., Leroy, S., d'Acremont, E., Roche, V., Evain, M., Lepretre, A., Verrier, F., Aslanian, D., Dias, N., Afilhado, A., Schnurle, P., Castilla, R., Despinois, F., Moulin, M.
Abstract: A variety of structures results from the interplay of evolving far-field forces, plate kinematics, and magmatic activity during continental break-up. The east Limpopo transform margin, offshore northern Mozambique, formed as Africa and Antarctica separated during the mid-Jurassic period break-up of the Gondwana supercontinent. The nature of the crust onshore has been discussed for decades in an effort to resolve issues with plate kinematic models. Two seismic refraction profiles with coincident multichannel seismic reflection profiles allow us to interpret the seismic velocity structures across the margin, both onshore and offshore. These seismic profiles allow us to (a) delineate the major regional crustal domains; (b) identify widespread indications of magmatic activity; and (c) map crustal structure and geometry of this magma-rich transform margin. Careful examination of the profiles allows us to make the following observations and interpretations: (a) on land, continental crust is overlain by a >10-km thick volcano-sedimentary wedge related to an early rifting stage, (b) offshore, thick oceanic crust formed due to intense magmatic activity, and between the two (c) a 50-60-km wide transform zone where the crustal structures are affected by intense magmatic activity and faulting. The prominent presence of intrusive and extrusive igneous units may be attributed to the combination of a deep-seated melting anomaly and a trans-tensional fault zone running through thinned lithosphere that allowed melt to reach the surface. A comparison of the crustal thinning along other transform margins shows a probable dependence with the thermal and/or tectonic history of the lithosphere.
Earth and planetary Science Letters, Vol. 579, 117343, 11p. Pdf
Australia
geophysics- seismics
Abstract: Plate tectonics, including rifting, subduction, and collision processes, was likely to have been different in the past due to the secular cooling of the Earth. The northeastern part of the West Australian Craton (WAC) has a complex Archean and Paleoproterozoic tectonic history; therefore, it provides an opportunity to study how subduction and collision processes evolved during the emergence of plate tectonics, particularly regarding the assembly of Earth's first supercontinent, Columbia. Because the northeastern boundary of the WAC and the southwestern boundary of the North Australian Craton (NAC) are covered by the Phanerozoic Canning Basin, the regional tectonic evolution has remained enigmatic, including how many tectonic elements were assembled and what may have driven rifting and subsequent collision events. Here, we use new passive-source seismic modeling to identify a seismically distinct segment of the lithosphere, the Percival Lakes Province, which lies east of the Pilbara Craton and is separated by two previously unknown southeast-trending lithosphere scale Paleoproterozoic sutures. We interpret that the northeastern suture, separates the Percival Lakes Province from the NAC, records the amalgamation of the WAC with the NAC. The southwestern suture separates the PLP from the reworked northeastern margin of the Pilbara Craton, including the East Pilbara Terrane and the Rudall Province. A significant upper mantle dipping structure was identified in the southwestern suture, and we interpret it to be a relic of subduction that records a previously unknown Paleoproterozoic collision that pre-dated the amalgamation of the WAC and NAC. By comparing our findings with previously documented dipping features, we show that the Paleoproterozoic collisions are seismically distinguishable from their Phanerozoic counterparts.
Journal of South American Earth Sciences, Vol. 111, 19p. Pdf
South America, Brazil
geophysics - seismics
Abstract: The Orocaima SLIP consists of an association of acid-intermediate volcanic-plutonic rocks. The volcanic rocks were generated in explosive eruptions through low eruptive columns, probably associated with fissural volcanism in the north of the Amazonian Craton, Brazil, between 2.0 and 1.98 Ga. It generated ignimbrites, whose facies (volcanic breccia rich in lithic, lapilli-tuff and lithic lapilli-tuff) show the proximity of the source. The extensive area of ca. 200.000 km2 of ignimbrite, rhyolite and dacite deposits, as well as the age range (2.0-1.98 Ga) and geochemical signatures suggest that the Orocaima volcano-plutonism may correspond to one of the oldest silicic LIPs in the world. The silicic volcanism is essentially subaerial and characterized by high-grade ignimbrites (densely welded) and subordinate lava, the ages of which indicate the longevity of the volcanic event in the Orosirian. They are included in the Surumu Group and comprise rocks with high-K calc-alkaline affinities and were emplaced in a subduction-related setting, similar to the rocks that extend through Venezuela, Guyana and Suriname (Cuchivero-Surumu-Iwokrama-Dalbana metavolcanic belt - CSID). The occurrence of mafic fragments disseminated in volcanic and granitic rocks in the north of Roraima, Brazil and in other segments of the CSID belt suggests the coexistence of acid and basic magmas. Except for one sample (?Nd(t) = -2.3), the Nd isotopic data of analyzed Surumu Group volcanic rocks yielded positive ?Nd(t) values (0.5-4.48; TDM = 2.0-2.47 Ga), suggesting generation from magmas derived from the mantle or from the melting of new juvenile crust. The Orocaima volcanism bears no evidence of involvement of Archean sources in the generation of the rocks. Thus, the Orocaima volcano-plutonism may represent one of the most significant ignimbrite eruption events during the Palaeoproterozoic in the world.-
Journal of South American Earth Sciences, Vol. 111, 10p. Pdf
South America, Brazil
geophysics - seismics
Abstract: The crustal thickness in South America has been mostly determined using seismological techniques. However, because these techniques provide point constraint or profile-specific results, the crustal thickness maps become especially dependent on both the number and spatial distribution of seismological stations. In the Amazonian Craton, the extensive forest cover restricts the number of existing stations, not allowing to elaborate a solely seismological crustal thickness model with homogeneous data coverage. Therefore, to overcome this difficulty, this work proposes a crustal thickness model for the Amazonian Craton developed based on the Parker-Oldenburg method and the Global Geopotential Model called GECO, considering the relationships between wavelengths and depths of the investigation sources. Furthermore, the developed iterative process allowed to determine the average depth of the crust-mantle interface, the density contrast at the interface, and the minimum and maximum frequencies used in the signal filtering process, making the model more robust for defining the used constants. The average crustal thickness of the Amazonian Craton was estimated as 40.25 km, with a standard deviation of the differences of 4.91 km, compared to crustal thickness defined by the seismological data. The estimated model shows great consistency with the data set used while allowing important inferences about craton compartmentation. Also, the geological provinces displayed an N-S connecting trend under the Amazonas, Solimões, and Acre basins, correlating the Guyana Shield with the Central Brazil Shield. Additionally, we observed various tectonic cycles acting on the craton while significantly modifying the structure of the provinces, possibly removing cratonic roots and rejuvenating the crust in older provinces.
Abstract: The seismological properties of Earth’s inner core are key to understanding its composition, dynamics and growth history. Within the inner core, fine-scale heterogeneity has previously been identified from backscattering of high-frequency compressional waves. Here we use historical earthquake and explosion data from the Large Aperture Seismic Array, USA, between 1969 and 1975 to build a 3D map of heterogeneity from the inner-core boundary to 500?km depth and determine the geographical distribution of the scatterers across the 40% of the inner core that is visible to the array. Our model has two regions of strong scattering, one beneath eastern Asia and the other beneath South America, both located where past local surveys have identified scattering. We suggest that these loci of strong, fine-scale heterogeneities may be related to random alignments of small, inner-core crystals due to fast freezing. These areas, which have been identified as having high attenuation and lie beneath colder areas of the core-mantle boundary, potentially provide constraints on the dynamics of the inner core and the motions in the outer core, with downwelling in the mantle and outer core possibly associated with strong scattering and inner-core heterogeneity.