Crustal Thickness Map Research Articles

Quantifying Crustal Growth in Arc‐Backarc Systems: Gravity Inversion Modeling of the Lau Basin

AbstractThe formation and evolution of arc‐backarc systems govern crustal production in some of the most volcanically and hydrothermally active environments on Earth. This study presents the first complete three‐dimensional density model of the active arc‐backarc system in the southwest Pacific comprising the Lau Basin and Tofua arc. Seafloor density and crustal thickness maps reveal changes in crustal composition and growth rates throughout the basin and along the volcanic arc. Crustal thickness varies significantly between the different centers of accretion (i.e., assemblages), resulting from seafloor spreading and subsurface melt accumulation below volcanic fields. Volumetric growth rates were calculated for each assemblage, corresponding to their respective contribution to basin expansion. The highest crustal density and growth rates are thought to be related to a mantle‐derived melt source entering the basin from the north around the edge of the subducting Pacific Plate. This study shows that the inverse modeling approach can be applied to global gravity data sets to characterize and quantify the density and thickness of the crust anywhere in the oceans.

Mercury's Crustal Porosity as Constrained by the Planet's Bombardment History

AbstractKnowing the structure of the crust is critical to understanding a planet's geologic evolution. Crustal thickness inversions rely on bulk density estimates, which are primarily affected by porosity. Due to the absence of high‐resolution gravity data, Mercury's crustal porosity has remained unknown. Here, we use a model that was calibrated to the Moon to relate Mercury's impact crater population and long‐wavelength crustal porosity in the cratered terrains. Therein, porosity is created by large impacts and then decreased as the surface ages due to pore compaction by smaller impacts and overburden pressure. Our models fit independent gravity‐derived porosity estimates in the northern regions, where data is well resolved. Porosity in the cratered terrains is found to be 9%–18% with an average and standard deviation of 13% 2%, indicating lunar‐like crustal bulk densities of 2,565 70 kg from which updated crustal thickness maps are constructed.

Regional crustal structure of Indonesia from receiver functions

Characterizing the crustal structure of Indonesia is important to gain a better understanding of its geodynamic evolution and improve seismic hazard assessments in the area. However, a unified crustal model of the entire Indonesian region and its surroundings is lacking. We present new maps of crustal thickness and bulk Vp/Vs ratio in Indonesia and the surrounding area that are obtained using P-wave receiver functions at 36 seismic stations from several permanent regional networks. The measured crustal thickness varies from ∼24 km to ∼38 km. The thickest crust, ∼38 km, is beneath Flores Island, southern Maluku, and neighboring northernmost Australia, whereas the thinnest crust, ∼24 km, is found under eastern Malaysia. Thus, crustal thickness varies by ∼14 km (from ∼24 km to ∼38 km) despite the small changes in elevation at the measurement points. The Vp/Vs ratios are 1.79±0.11, with high values (>1.85) found along the Banda-Sunda arc-trench system. We attribute these high values to: (1) the presence of mafic island arc and oceanic crust and (2) partial melting within this volcanic region, which causes a larger decrease in S-wave velocities compared with P-wave velocities. The comparison of the seismic properties of Indonesian island arc crust, particularly the Vp/Vs ratio, with laboratory measurements and the petrology of the exhumed Talkeetna island arc, Alaska, allows us to infer the crustal composition of Indonesian island arc crust.

A High‐Resolution Shear Velocity Model of the Crust and Uppermost Mantle Beneath Westernmost Mediterranean Including Radial Anisotropy

AbstractUsing seismic data from 1,186 stations deployed across the westernmost Mediterranean, we construct a high‐resolution 3‐D radially anisotropic model from a joint inversion of receiver functions and Rayleigh and Love wave dispersions. The Rayleigh and Love data are extracted from both ambient noise interferograms and earthquake waveforms, and a new three‐station ambient noise interferometry method is used to further improve the data coverage. The obtained crustal thickness map reproduces independently the Moho depth maps compiled from previous data, showing crustal roots beneath the Pyrenean‐Cantabrian range and the Gibraltar Arc and thicker crust beneath the Variscan Iberian Massif than in the area affected by Alpine orogenesis. TheVsvcrustal model outlines the Iberian Massif as a high shear wave velocity block and shows extremely low velocities in the Gulf of Cadiz and Gibraltar Arc. At mantle levels, a sharp boundary between the Aquitanean Basin and the Massif Central is imaged, with lowVsvbeneath the Massif Central probably reflecting a remaining signature of the magma. To the south of Iberia, the geometry of the Alboran slab is captured by the model, while beneath the Atlas Mountains, widespread lowVsvand positive radial anisotropy is observed, favoring the edge‐driven convection model explaining the lithospheric thinning. The most relevant contribution of this work is mapping, for the first time at this scale, the radial anisotropy anomalies at crustal level, with higher values observed in the Alpine Iberia, dominated by present‐day extensional tectonics.

Influence of Lateral Plume Channel on the Evolution of Rift Arms of the Afar Triple Junction: Constraints From 3‐D Gravity Interpretation

AbstractThe Afar region presents a unique example of a juvenile R‐R‐R triple junction with the arms evolved to various stages of maturity. The lithosphere in this area has been modified heterogeneously by the nearby Afar mantle plume. To gain a comprehensive and consistent understanding of the link between crustal geometry and geodynamics of the region and facilitate comparisons across different tectonic elements, we developed a high‐resolution crustal thickness map of the region that resolves major tectonic features. The calculation involved a well‐constrained 3‐D inversion of upper mantle‐corrected crustal Bouguer gravity anomalies computed with the help of tomography‐derived densities. Our results match earlier localized crustal thickness estimates, illustrating the method's efficacy. A thinner crust (∼20–26 km) is revealed in the Afar compared to adjacent continental regions. The Main Ethiopian Rift has the thickest crust in the central part of the rift (∼38–40 km) and thinner to the north (∼30 km) and south (∼35 km). The crustal thickness along the Gulf of Aden changes from east (∼5 km) to west (∼17 km) in the west, consistent with a westward transition from seafloor spreading to continental rifting. The southern Red Sea is characterized by asymmetric crustal thickness, with a thicker crust (∼17 km) on the eastern flank compared to the western flank (∼11 km), which we attribute to the influence of a sub‐lithospheric channel from the Afar mantle plume. We propose a model describing the channel leading to asymmetric features of the southern Red Sea, such as crustal accretion, marginal topography, and volcanism.

Crustal architecture and metallogeny associated with the Paleo-Tethys evolution in the Eastern Kunlun Orogenic Belt, Northern Tibetan Plateau

The Eastern Kunlun Orogenic Belt (EKOB) in the Northern Tibet Plateau hosts a wide variety of metal deposits related to the Late Paleozoic to Mesozoic magmatism. In this study, we investigate the spatio-temporal distribution of the Late Paleozoic to Mesozoic granitic rocks and associated metal deposits in the EKOB and provide a comprehensive compilation of the geochronological, geochemical and isotopic data on these rocks. We compute regional zircon Hf isotope and crustal thickness maps from the data, based on which a comprehensive model is proposed involving subduction (ca. 270–240 Ma), continental collision (ca. 240–224 Ma), and post-collisional extension (ca. 224–200 Ma) for the Late Paleozoic to Mesozoic Paleo-Tethys evolution in the EKOB.Zircon Hf isotopic and crustal thickness mapping of Late Paleozoic to Mesozoic magmatic rocks was carried out to evaluate their spatio-temporal and genetic links with the regional metallogeny. The polymetallic Fe-skarn and porphyry Cu (Mo) deposits in the EKOB are located above the Moho uplift region, featuring a comparatively thin crust. Granites associated with porphyry Cu (Mo) and polymetallic Fe skarn mineralization are commonly characterized by high εHf(t) and younger TDMc values, whereas granite related to Cu-Mo-Sn skarn deposits exhibit more variable εHf(t) values, TDMc ages, and the crust thickness, which suggest that more crustal materials contributed to the formation of Cu-Mo-Sn skarn deposits than those for porphyry Cu (Mo) and polymetallic Fe skarn mineralization. In contrast, vein-type Au deposits are located primarily where the Moho surface displays a depression, i.e., where the continental crust is relatively thick. The magmatic rocks associated with Au mineralization are characterized by low εHf(t) and high TDMc values, representing reworked ancient crustal components, similar to those associated with porphyry Mo and epithermal Ag-Pb-Zn-(Au) deposits. Our study indicates that the emplacement of magmatic-hydrothermal deposits was controlled by the crustal structure and magma sources.

A Positive Feedback Between Crustal Thickness and Melt Extraction for the Origin of the Martian Dichotomy.

A North/South difference in crustal thickness is likely at the origin of the Martian dichotomy in topography. Recent crustal thickness maps were obtained by inversion of topography and gravity data seismically anchored at the InSight station. On average, the Martian crust is 51-71km thick with a southern crust thicker by 18-28km than the northern one. The origin of this crustal dichotomy is still debated although the hypothesis of a large impact is at present very popular. Here, we propose a new mechanism for the formation of this dichotomy that involves a positive feedback between crustal growth and mantle melting. As the crust is enriched in heat-producing elements, the lid of a one-plate planet is hotter and thinner where the crust is thicker, inducing a larger amount of partial melt below the lid and hence a larger rate of melt extraction and crustal growth. We first demonstrate analytically that larger wavelength perturbations, that is, hemispherical perturbations, grow faster because smaller wavelengths are more attenuated by thermal diffusion. We then use a parameterized thermal evolution model with a well-mixed mantle topped by two different lids characterized by their thermal structures and thicknesses to study the growth of the crust in the two hemispheres. Our results demonstrate that this positive feedback can generate a significant crustal dichotomy.

Lithospheric and Slab Configurations From Receiver Function Imaging in Northwestern South America, Colombia

AbstractThe northwestern corner of South America presents a highly complex geodynamical setting including the convergence between South America, Nazca, and Caribbean plates; which also entails heterogeneity and complexity in the lithospheric structure. Here, we present depth‐migrated common‐conversion point (CCP) teleseismic receiver function profiles carried out using seismological data between 2011 and 2020 from the Colombian Geological Service, combined with some stations of observatories and international networks. These profiles provide better constraints on the lateral interpretation of the main lithospheric discontinuities in this region. We provide an updated crustal thickness map displaying some differences from previous studies, especially in the northernmost part of the cordilleran system. Our results also show some intracrustal features such as the regional detachment surface beneath the Eastern Cordillera and the presence of deep melt reservoirs under the Central Cordillera and the Paipa‐Iza complex in the Eastern Cordillera. We also discuss the lithospheric thickness and the potential existence of a highly hydrated and serpentinized mantle wedge under the forearc over the slab segment south of the Caldas tear. Furthermore, the CCP profiles further support the subduction system segmentation under northwestern South America, with at least five slab segments, two of them overlapping between ∼5.5 and ∼8°N, three belonging to the Nazca plate, and two having a Caribbean origin. The Bucaramanga nest appears to be linked to the interaction between the two northmost slab segments, while the Cauca cluster might be related to a slab dehydration process of a highly fractured zone.

The Neoproterozoic basement of the Parnaíba Basin (NE Brazil) from combined geophysical-geological analysis: A missing piece of the western Gondwana puzzle

A regional Pre-Silurian unconformity marks the erosive planar base of the Phanerozoic Parnaíba Basin (PB), and indistinctly truncates key tectonic units for the understanding of the Neoproterozoic evolution of western Gondwana in northeastern Brazil. We present a comparative analysis of geophysical and geological datasets in the PB and propose a new tectonic configuration for its pre-Silurian basement, composed of different terranes amalgamated during the Brasiliano orogeny. Using integrated seismic interpretation and gravity modelling, constrained by an updated crustal thickness map, well data and a compilation of recent geophysical studies, we have identified two main blocks in the center of the basin. They represent pre-Brasiliano inliers, surrounded by Brasiliano mobile belts, and assigned to two major crustal building blocks of western Gondwana: the Amazonian-West Africa (AWB; to the west/north) and the Central African blocks (CAB; to the east/northeast). The Grajaú block belongs to the AWB and is characterized by low gravity anomaly, thicker crust (41–45 km), transparent seismic character of the basement and a high velocity lower crust. The Teresina block belongs to the CAB and is characterized by slightly thinner crust (39–41 km), higher values of gravity and magnetic anomalies and by the presence of mid-crustal reflectors (MCR), observed in seven seismic lines and here interpreted as a remnant of a paleosuture zone between the two blocks. Along a NE-SW 500-km seismic and gravity profile, we interpreted the MCR as composed of crustal-scale thrust faults verging westwards and defining the Barra do Corda belt, which was formed by the closure of the Goiás-Pharusian ocean along the Transbrasiliano-Kandi corridor and subsequent collision of the Grajaú and Teresina blocks. This belt deforms the eastern margin of the Ediacaran Riachão foreland basin (RB), observed in seismic and well data beneath the SW portion of PB and overlying the southern Grajaú block. The western margin of RB is bounded by eastwards verging thrust faults, interpreted as a zone of back-thrusts and thinned crust (∼36 km) in the eastern prolongation of the Araguaia belt beneath PB. To the east, the limit between the Teresina block and the Borborema Province is marked by the NE-SW Transbrasiliano Shear Zone, along which Neoproterozoic mylonites were recovered from wellbores, and a system of narrow pull-apart basins is interpreted in the seismic data, triggered by Early Paleozoic transtensional reactivations of the TBSZ, prior to the deposition of the basal Ordovician-Silurian sequence of PB. This complex collisional tectonic setting of PB basement contradicts the idea of a stable cratonic block beneath it, inciting new formation models that account for the observed crustal and lithospheric heterogeneities.

Regional Geophysics of the Caribbean and Northern South America: Implications for Tectonics

AbstractThe Caribbean plate is an enclosed oceanic basin whose formation and evolution are controversial. In the most commonly accepted model, the Caribbean plate is mainly composed of the Caribbean Large Igneous Province (CLIP) and the buoyant characteristic of this oceanic plateau resisted subduction and allowed an eastward migration to its present position north of South America. In this study, we integrate a broad range of geophysical and geomorphological data to define structural elements and present‐day tectonics of the Caribbean plate and the surrounding region. We present a Bouguer gravity anomaly map and a new crustal thickness map that documents large areas of normal‐thickness oceanic crust within the Venezuela and Colombia basins of the Caribbean plate. Selected cross sections of seismicity and P‐wave anomalies from a seismic tomographic model depict the present‐day geometry of subducting oceanic plates within the Caribbean region. We observe that rather than resisting subduction, as expected for the thick crust of a buoyant large igneous province, the subduction of the Caribbean plate can be traced to a depth of 600 km beneath NW South America. This, together with the crustal thickness map, implies that a significant area of the Caribbean plate, including the subducted portion, is composed of normal‐thickness oceanic crust. As proposed by the Pacific origin model, the Caribbean plate likely migrated eastward from the Pacific Ocean as an oceanic plate mostly with normal‐thickness crust and limited portions of the crust thickened by hot spot volcanism (CLIP).

Crustal Thinning From Orogen to Back‐Arc Basin: The Structure of the Pannonian Basin Region Revealed by P‐to‐S Converted Seismic Waves

AbstractWe present the results of P‐to‐S receiver function analysis to improve the 3D image of the sedimentary layer, the upper crust, and lower crust in the Pannonian Basin area. The Pannonian Basin hosts deep sedimentary depocentres superimposed on a complex basement structure and it is surrounded by mountain belts. We processed waveforms from 221 three‐component broadband seismological stations. As a result of the dense station coverage, we were able to achieve so far unprecedented spatial resolution in determining the velocity structure of the crust. We applied a three‐fold quality control process; the first two being applied to the observed waveforms and the third to the calculated radial receiver functions. This work is the first comprehensive receiver function study of the entire region. To prepare the inversions, we performed station‐wise H‐Vp/Vs grid search, as well as Common Conversion Point migration. Our main focus was then the S‐wave velocity structure of the area, which we determined by the Neighborhood Algorithm inversion method at each station, where data were sub‐divided into back‐azimuthal bundles based on similar Ps delay times. The 1D, nonlinear inversions provided the depth of the discontinuities, shear‐wave velocities and Vp/Vs ratios of each layer per bundle, and we calculated uncertainty values for each of these parameters. We then developed a 3D interpolation method based on natural neighbor interpolation to obtain the 3D crustal structure from the local inversion results. We present the sedimentary thickness map, the first Conrad depth map and an improved, detailed Moho map, as well as the first upper and lower crustal thickness maps obtained from receiver function analysis. The velocity jump across the Conrad discontinuity is estimated at less than 0.2 km/s over most of the investigated area. We also compare the new Moho map from our approach to simple grid search results and prior knowledge from other techniques. Our Moho depth map presents local variations in the investigated area: the crust‐mantle boundary is at 20–26 km beneath the sedimentary basins, while it is situated deeper below the Apuseni Mountains, Transdanubian and North Hungarian Ranges (28–33 km), and it is the deepest beneath the Eastern Alps and the Southern Carpathians (40–45 km). These values reflect well the Neogene evolution of the region, such as crustal thinning of the Pannonian Basin and orogenic thickening in the neighboring mountain belts.

DEEP CRUSTAL STRUCTURE IN NORTHEASTERN EURASIA AND ITS CONTINENTAL MARGINS

The paper reports on the deep geophysical studies performed by the Geological Survey of Russia (VSEGEI) under the international project – Deep Processes and Metallogeny of Northern, Central and Eastern Asia. A model of the deep crustal structure is represented by a set of crustal thickness maps and a 5400-km long geotransect across the major tectonic areas of Northeastern Eurasia. An area of 50000000 km2 is digitally mapped in the uniform projection. The maps show the Moho depths, thicknesses of the main crustal units (i.e. the sedimentary cover and the consolidated crust), anomalous gravity and magnetic fields (in a schematic zoning map of the study area), and types of the crust. The geotransect gives the vertical section of the crust and upper mantle at the passive margin of the Eurasian continent (including submarine uplifts and shelf areas of the Arctic Ocean) and the active eastern continental margin, as well as an area of the Pacific plate.

Integration of gravity, magnetic, and seismic data for subsalt modeling in the Northern Red Sea

Rifts and rifted passive margins are often associated with thick evaporite layers, which challenge seismic reflection imaging in the subsalt domain. This makes understanding the basin evolution and crustal architecture difficult. An integrative, multidisciplinary workflow has been developed using the exploration well, gravity and magnetics data, together with seismic reflection and refraction data sets to build a comprehensive 3D subsurface model of the Egyptian Red Sea. Using a 2D iterative workflow first, we have constructed cross sections using the available well penetrations and seismic refraction data as preliminary constraints. The 2D forward model uses regional gravity and magnetic data to investigate the regional crustal structure. The final models are refined using enhanced gravity and magnetic data and geologic interpretations. This process reduces uncertainties in basement interpretation and magmatic body identification. Euler depth estimates are used to point out the edges of high-susceptibility bodies. We achieved further refinement by initiating a 3D gravity inversion. The resultant 3D gravity model increases precision in crustal geometries and lateral density variations within the crust and the presalt sediments. Along the Egyptian margin, where data inputs are more robust, basement lows are observed and interpreted as basins. Basement lows correspond with thin crust ([Formula: see text]), indicating that the evolution of these basins is closely related to the thinning or necking process. In fact, the Egyptian Northern Red Sea is typified by dramatic crustal thinning or necking that is occurring over very short distances of approximately 30 km, very proximal to the present-day coastline. The integrated 2D and 3D modeling reveals the presence of high-density magnetic bodies that are located along the margin. The location of the present-day Zabargad transform fault zone is very well delineated in the computed crustal thickness maps, suggesting that it is associated with thin crust and shallow mantle.

High‐Resolution Crustal and Uppermost Mantle Structure Beneath Central Mongolia From Rayleigh Waves and Receiver Functions

AbstractIn this study, I present a high‐resolution 3‐D Vs model of the crust and uppermost mantle beneath central Mongolia, based on seismic data recorded at the Central Mongolia Seismic Experiment array and some other networks. The model is constructed by a Bayesian Monte Carlo inversion that jointly interprets Rayleigh wave phase speeds and receiver functions. The Rayleigh wave dispersion measurements are derived from both ambient noise and teleseismic earthquake data. To improve ambient noise data coverage, a recently developed three‐station interferometry (Zhang et al., 2020, https://doi.org/10.1093/gji/ggaa046) technique is applied to produce high‐quality three‐station interferograms that augment the ambient noise data base. There are several prominent geological and tectonic features captured by the model, including (1) the Khövsgöl rift system is imaged as a high‐speed anomaly in the upper crust, whose southern edge follows the Bulnay fault. (2) The estimated crustal thickness map identifies relatively thick crust beneath the Khövsgöl rift system, the Hangay Dome, and the Gobi‐Altai Mountains. (3) A crustal high‐velocity anomaly is imaged beneath the Taryatu‐Chulutu volcanic field, which is interpreted as a solidified intrusive magma body. (4) Dispersed low‐velocity anomalies in the mantle identify locations with possible partial melting. The melting zones could be related to small‐scale asthenospheric upwelling and convective removal of lithosphere, which produces uplift of the Hangay Dome.

The Thermal Evolution of Mercury Over the Past ∼4.2 Ga as Revealed by Relaxation States of Mantle Plugs Beneath Impact Basins

AbstractWe investigated the relaxation states of 11 large impact basins on Mercury based on an updated crustal thickness map, finding that the pre‐Tolstojan basins have comparable instead of varied relaxation states, suggesting that Moho temperature (TempMoho ) did not decrease substantially from ∼4.2 to 3.8 Ga. At the same time, mantle uplift beneath the Caloris basin is the least degraded, therefore implying a sharp decrease of TempMoho ∼ 3.8 Ga. These findings contrast with our thermal evolution models that predict a fast decrease of TempMoho between ∼4.2 and 3.8 Ga. Therefore, the discrepancies in the cooling rate suggest that the relatively elevated bombardment history between ∼4.2 and 3.8 Ga might have input additional energy to Mercury and substantially decreased the cooling rate.

Crustal thickness in the Black Sea and surrounding region, estimated from the gravity data

In this research, the crustal thickness (Moho depth) variations in the Black Sea and surrounding region were calculated from the regional gravity anomaly data using an empirical equation. The data were obtained from the open source of Bureau Gravimétrique International (BGI) in digitized form. The gravity anomaly values are changing from −80 mGal in the onshore Black Sea region in Turkey to 260 mGal in the deeper parts of the sea. Maximum gravity anomaly indicates a consistent pattern in the offshore Black Sea and anomaly pattern follows the shape of the shoreline. The pattern is consistent with the deeper parts in the marine area and the eastern anomaly extends into the Russian territory between NW of Georgia and SE of the Crimean Peninsula. The crustal thickness map indicates almost the same trends. In this map, the minimum crustal thickness is about 10 km around the mid of the Shatsky Ridge, close to the coastal regions of SE Russia and the second thinnest (12 km) place is located around the Western Black Sea Basin. The maximum thicknesses are represented with the 34–36 km contours in the onshore Eastern Black Sea region and east of Anatolia. On the other hand, the maximum crustal thickness to the north can reach up to 28 km around the Sea of Azov and to the east. The thickness variations are generally related with the main tectonic trends in this region. The Analytic Signal (AS) and maxima points of the horizontal gradients (maxspots) maps constructed from the gravity anomalies are also prepared in this research to see if these anomaly trends follow the tectonic lines. The regions of maximum and minimum crustal thicknesses are well-consistent with the AS anomalies and boundaries of the bifocal thin crustal sections coinciding with the Eastern and Western Black Sea basinal parts are surrounded with the maxspots. The northern thick crustal region around the southern margin of the Indolo-Kuban Basin is also surrounded by the maxspots. It is possible to claim that the maxspots map from the horizontal gradient of Bouguer anomalies reflects the boundaries of the main tectonic units in the Black Sea Region. Particularly, southern boundary of the Scythian Platform, southern boundary of the Shatsky Ridge and Great Caucasus thrust are distinctive and easy to follow in the maxspots map.

Lithospheric structures of and tectonic implications for the central–east Tibetan plateau inferred from joint tomography of receiver functions and surface waves

SUMMARYWe present an updated joint tomographic method to simultaneously invert receiver function waveforms and surface wave dispersions for a 3-D S-wave velocity (Vs) model. By applying this method to observations from ∼900 seismic stations and with a priori Moho constraints from previous studies, we construct a 3-D lithospheric S-wave velocity model and crustal-thickness map for the central–east Tibetan plateau. Data misfit/fitting shows that the inverted model can fit the receiver functions and surface wave dispersions reasonably well, and checkerboard tests show the model can retrieve major structural information. The results highlight several features. Within the plateau crustal thickness is >60 km and outwith the plateau it is ∼40 km. Obvious Moho offsets and lateral variations of crustal velocities exist beneath the eastern (Longmen Shan Fault), northern (central–east Kunlun Fault) and northeastern (east Kunlun Fault) boundaries of the plateau, but with decreasing intensity. Segmented high upper-mantle velocities have varied occurrences and depth extents from south/southwest to north/northeast in the plateau. A Z-shaped upper-mantle low-velocity channel, which was taken as Tibetan lithospheric mantle, reflecting deformable material lies along the northern and eastern periphery of the Tibetan plateau, seemingly separating two large high-velocity mantle areas that, respectively, correspond to the Indian and Asian lithospheres. Other small high-velocity mantle segments overlain by the Z-shaped channel are possibly remnants of cold microplates/slabs associated with subductions/collisions prior to the Indian–Eurasian collision during the accretion of the Tibetan region. By integrating the Vs structures with known tectonic information, we derive that the Indian slab generally underlies the plateau south of the Bangong–Nujiang suture in central Tibet and the Jinsha River suture in eastern Tibet and west of the Lanchangjiang suture in southeastern Tibet. The eastern, northern, northeastern and southeastern boundaries of the Tibetan plateau have undergone deformation with decreasing intensity. The weakly resisting northeast and southeast margins, bounded by a wider softer channel of uppermost mantle material, are two potential regions for plateau expansion in the future.

Bulk Structure of the Crust and Upper Mantle beneath Alaska from an Approximate Rayleigh-Wave Dispersion Formula

AbstractWe introduce a method for estimating crustal thickness and bulk crustal and upper-mantle shear-wave velocities directly from high-quality measurements of fundamental-mode Rayleigh-wave dispersion in the period range from 10 to 40 s. The method is based on an approximate Rayleigh-wave dispersion formula and provides fast results with minimal model parameterization. We apply the method to Rayleigh-wave phase maps in Alaska to reveal first-order structure in a region that had not been systematically and densely instrumented prior to the Transportable Array (TA). To demonstrate the consistency of the results, we also apply the same method to existing Rayleigh-wave phase maps derived from TA data in the conterminous United States, where crustal and upper mantle structures are better known. We contrast features observed in maps of crustal thickness and bulk shear-wave velocity between the Cascadia and Alaska-Aleutian subduction zones to highlight differences in the two regions. Our results show that, contrary to conventional wisdom, first-order information on the location of major depth discontinuities (e.g., the Moho) can be extracted in a fast, straightforward manner from measurements of Rayleigh-wave dispersion alone.

Crustal Thickness Variation Across the Sea of Marmara Region, NW Turkey: A Reflection of Modern and Ancient Tectonic Processes

AbstractThe Marmara region in Turkey is an important geological setting, both from a tectonic and a seismic hazard/risk perspective. We present a new map of crustal thickness variation across this complex region to better understand the interplay of past and present tectonic processes that have formed present‐day structure. Maps of crustal thickness are created using Ps converted phases and receiver function (RF) analysis of earthquakes recorded at all publicly available seismic stations and stations in the national monitoring network (run by AFAD Disaster and Emergency Management Authority Turkey). RFs are converted from time to depth using a local 3‐D full‐waveform tomographic model and are combined in multiphase common conversion point stacks. Direct P to S converted arrivals and associated multiples are mapped to produce continuous maps of the Moho discontinuity. Results show Moho depths ranging from 26–41 km with a regional trend of westward thinning reflecting the effects of the extensional regime in western Anatolia and the neighboring Aegean Sea. The thinnest crust is observed beneath the western end of the Sea of Marmara, attributed to transtensional basin opening. A distinct region of increased crustal thickness bounded by the West Black Sea Fault in the west, and the northern strand of the North Anatolian Fault in the south, defines the ancient crustal terrane of the Istanbul Zone. Isostatic arguments indicate that the thickened crust and lower elevation in the Istanbul Zone require it to be underlain by thicker lithosphere, a conclusion that is consistent with its hypothesized origin near the Odessa shelf.

PHILCRUST3.0: New constraints in crustal growth rate computations for the Philippine arc

A revised crustal thickness map for the Philippine island arc, PHILCRUST3.0, has been generated using new gravity data from the EGM2008. The gravity-derived crustal thickness is compared with crustal thickness estimates obtained using the Sr/Y ratios. PHILCRUST3.0 shows that crustal thickness in the Philippines varies from 13 to 33 km. The Philippine arc appears to be at its thinnest in areas underlain by ophiolites and young magmatic arcs. The resulting crustal thickness values are used to revisit the growth rates of the Philippine arc from magmatic and amagmatic contributions. Arc magmatism contributes more significantly to the growth of Philippine crust with rates from specific volcanoes varying from 5 to 35 km3/km/m.y. Despite being relatively young, magmatism associated with the volcanoes of the East Philippine Arc displayed a higher growth rate (32 km3/km/m.y.) compared to volcanoes belonging to older magmatic arcs in the archipelago. This is attributed to the more hydrated oceanic lithosphere subducting beneath the Philippine Trench. PHILCRUST3.0 provides the most complete crustal thickness map of the Philippines to date. This will enable the computation of crustal growth rate in areas within the Philippines where point gravity or geochemical data is not available.