Crustal Thickness Research Articles

Insight into intracontinental deformation and crustal reworking of ancient continents implied by crustal structure imaging of the Yangtze Craton

The crustal imprints from multistage tectonic activities in cratons offer valuable insights into continental evolution. Utilizing seismic data from two densely deployed, nearly perpendicular linear arrays, and a newly developed stepwise joint inversion of depth-domain receiver function and surface wave dispersion, we constructed a detailed crustal layering model for the Yangtze Craton. Our analysis revealed elongated double velocity reversal zones as salient features of the crust, which likely record ancient crustal reworking and juvenile crustal growth associated with Neoproterozoic rift-related magmatic processes. The interlayering of low- and high-velocity structures may contribute to the enduring stability of the Yangtze Craton. Additionally, superimposed layers separated by east-dipping interfaces and abrupt changes in crustal thickness in the boundary belts surrounding the Yangtze Craton document the crust's structural responses to intracontinental deformation during continent assembly.

Present-day incipient fault coalescence at a relay zone (Jiloca extensional basin, Spain): Evidence from instrumental seismicity

The relay zones between NW-SE to NNW-SSE striking faults of the Jiloca graben (Iberian Chain) mostly show distributed along-strike fault and fracture patterns. The latter are chiefly controlled by the Late Pliocene-Quaternary regional stress field, and secondarily respond to local controls from inherited structures. Such fracture patterns contrast with the classical models of transverse connecting faults controlled by relay kinematics. North of the Concud fault trace, at the relay zone with the Sierra Palomera fault, an unusually high seismic activity has been noticed since 2014, with magnitudes up to M = 3.5. Upgrading of the National Seismic Network allowed obtaining such new detailed records, while the installation of a new seismometer by the IGN within the study area has improved the reliability of focal depth data since 2017. A high-precision absolute relocation of seismicity from 01/01/2000 to 30/05/2022 has been carried out. The results show that (i) the epicentres are significantly clustered along a nearly N-S trending band, and (ii) the focal depths range from 0 to 14 km, in good agreement with the thickness of the brittle crust. This 3D spatial distribution of seismicity is interpreted as a consequence of activation of either a single fault or a fault zone, nearly vertical and N-S striking. Such structural setting is consistent with the surficial fracture patterns observed at both map and outcrop scale: NNW-SSE and NNE-SSW oriented faults and fractures, orthogonal to the ENE-WSW to ESE-WNW regional σ3 trajectories, together with NW-SE trending ones controlled by inherited contractive faults. The present-day seismic activity suggests that along-strike, incipient fault propagation at the relay zone between the Concud and Sierra Palomera faults is currently operating under the control of the remote stress field.

Seismic Evidence for Velocity Heterogeneity Along ∼40 Ma Old Oceanic Crustal Segment Formed at the Slow‐Spreading Equatorial Mid‐Atlantic Ridge From Full Waveform Inversion of Ocean Bottom Seismometer Data

AbstractIn slow spreading environments, oceanic crust is formed by a combination of magmatic and tectonic processes. Using full waveform inversion applied to active‐source ocean bottom seismometer data, we reveal the presence of a strong lateral variability in the 40–48 Ma old oceanic crust formed at the slow‐spreading Mid‐Atlantic Ridge in the equatorial Atlantic Ocean. Over a 120 km‐long section between the St Paul fracture zone (FZ) and the Romanche transform fault (TF), we observe four distinct 20–30 km long crustal segments. The segment affected by the St Paul FZ consists of three layers, an ∼2 km thick layer with a P‐wave velocity <6 km/s, a 1.5 km thick middle crust with a velocity of 6–6.5 km/s, and an underlying layer where velocity is ∼7 km/s, representing the lower crust. The segment associated with an abyssal hill morphology contains a high velocity of ∼7 km/s at 2–2.5 km below the basement, indicating the presence of primitive gabbro or serpenized peridotite. The segment associated with a low basement morphology seems to have 5.5–6.5 km/s velocity starting near the basement extending down to ∼4 km depth, indicating chemically distinct crust. The segment close to the Romanche TF, a velocity 4.5–5 km/s near the seafloor increasing to 7 km/s at 4 km depth indicates a magmatic origin. The four distinct crustal segments have a good correlation with the overlying seafloor morphology. These observed strong crustal heterogeneities could result from alternate tectonic and magmatic processes along the ridge axis, possibly modulated by thermal and/or chemical variations in the mantle during their formation along the ridge segment.

Effects of the types and ages of vegetation restoration on land surface roughness in the Eastern Hobq Desert, Northern China

Vegetation restoration has an important remodelling effect on near-surface characteristics, and consequent changes in land surface roughness (LSR) are key influences on soil wind erosion processes. However, the effects of vegetation restoration types and ages on LSR and the underlying mechanisms are not fully understood. In this study, the sand-fixing vegetation restoration area of the Hobq Desert was examined in comparison to a bare sand control area. The LSR of four artificial vegetation types (Salix psammophila, Caragana korshinskii, Artemisia ordosica, and Populus simonii) with restoration age of 36 years, and Salix psammophila and Caragana korshinskii after different periods of restoration (20, 28, 36 and 45 years) were measured using Structure-from-Motion (SfM) photogrammetry. Near-surface characteristics that may affect LSR were also measured. The results showed that vegetation restoration was associated with a 230–409 % higher LSR compared to the control site (1.74 mm). LSR in the different vegetation restoration areas were ranked, from high to low, as follows: AO (8.85 mm) > CK (7.89 mm) > SP (6.70 mm) > PS (6.61 mm). LSR also increased with time since restoration, with the greatest rate of increase during the first 20 years. The change of LSR is mainly affected by the change of near-surface characteristics, with the direct effects of biological crust thickness (0.331), litter thickness (0.289), soil bulk density (−0.239), and clay content (0.171) being significant. Stem diameter, litter density, biological crust coverage, and soil organic matter affected LSR indirectly, mainly through acting on the above factors. Finally, LSR was effectively estimated based on biological crust thickness, litter thickness, and soil bulk density (R2 = 0.904). The research results will help to further deepen the understanding of the influence mechanism of vegetation restoration on LSR, and provide scientific basis and practical reference for vegetation ecological restoration in similar areas.

Relationships between high-K magmas and Cu–Au porphyry deposits in the Macquarie Arc, Australia

The Macquarie Arc is a long-lived volcanic arc terrane which began in the latest Cambrian and continued until the early Silurian. Early Silurian high-K suites associated with Cu–Au porphyry mineralization are typically considered shoshonitic; however, comparisons with well-studied examples of shoshonites in continental arcs and associated settings suggest few Macquarie Arc rocks would classify as shoshonitic. Notably, the Au-rich Cadia Intrusive Complex, which is marginally shoshonitic, and most high-K to ‘shoshonitic’ Macquarie Arc rocks are essentially high-K calc-alkaline, analogous to ‘primitive’ shoshonites generated within intraoceanic arc settings. Arc magmatism in the Macquarie Arc is consistently primitive and mantle derived prior to an increase in siliciclastic input and crustal thickening during the Silurian Benambran Orogeny, consistent with prior isotopic studies. Magma source regions can thus be interpreted as varying spatially and temporally due to variable enrichment from slab and mantle derived components, reflecting a maturing and potentially migrating arc system.

The geodynamics of plume-influenced mid-ocean ridges: insights from the Foundation Segment of the Pacific-Antarctic Ridge

The intersection of the Foundation Plume with the Pacific-Antarctic Ridge is a key location in global geodynamics where a mantle plume is approached by and interacting with a fast-spreading mid-ocean ridge. Here, we discuss a comprehensive major and trace element and Sr-Nd-Pb isotope dataset of new and existing samples from the young Foundation Seamount Chain (<5 Ma) and adjacent section of the Pacific-Antarctic Ridge. We use the geochemistry of axial, off-axis and intraplate lavas to map the spatial extent of plume dispersal underneath the ridge as well as the internal zonation of the upwelling plume. We show that the unusual length, increased crustal thickness and occurrence of silicic rocks on the axis of the Foundation Segment are the direct result of plume being tapped by the axial melting zone. We demonstrate that the plume is not homogeneous but shows a HIMU-like (high time-integrated 238U/204Pb) OIB (Ocean Island Basalt) component characterized by 206Pb/204Pb of up to 20.5 in its center and a more EM1-like (Enriched Mantle one) OIB component characterized by low U/Pb and 206Pb/204Pb but high Rb/Nb and 87Sr/86Sr towards its edges. Plume entrainment leads to a high magma supply rate that fosters the formation of silicic rocks and triggers the lengthening of the segment over time. However, plume dispersal is not symmetric as the geochemical tracers for the OIB component are extending <100 km northwards but >300 km southwards. We relate this to the current plate tectonic framework in which the obliquity between the migrating ridge and the absolute plate motions induces a sub-axial asthenospheric flow that preferentially channels plume material southwards.

Pleistocene glacial advances and exposure age scatter in the Olympus Range, Antarctica: A study of cosmogenic 36Cl/3He in dolerites and 10Be in sandstones

In three cirques in the western Olympus Range of the McMurdo Dry Valleys, Antarctica, previous advances of cirque glaciers are recorded by a sequence of three drifts in each of the cirques. We dated drift limits and the deposits on modern glaciers in two of these cirques, Dean and Dipboye, via cosmogenic 3He in pyroxene from 41 dolerite boulders, 36Cl in pyroxene from 12 of those dolerites, and 10Be in quartz from 11 sandstone boulders. Exposure age scatter is high on all deposits. The 3He exposure ages across all deposits range from ∼35 to ∼2300 ka and 10Be exposure ages range from ∼7 to ∼435 ka. Coupled 36Cl/3He from dolerites support constant exposure with erosion for nine of the 12 samples, while the other three might have experienced complex exposure-burial histories. Due to the mesa-butte topography and slow bedrock erosion rates, nuclide inheritance is the primary cause of age scatter in dolerites, accounting for >1 Myr of exposure age error. Mean exposure ages from sandstones are 2–7 times younger than those from dolerites for the same deposits, indicating that inheritance is less common in sandstones in this region. Weathering analyses of sandstone boulders show an increase in average siliceous crust thickness and rock strength with deposit age, an example of case hardening. Based on both relative and exposure age dating, drift age increases with distance from the modern glaciers in both Dean and Dipboye cirques, with three advances during the past <700 ka. However, due to high exposure age scatter, it cannot be determined if the three drifts are temporally correlated across the two cirques and therefore the drifts might record different glacial advances in Dean Cirque vs. Dipboye Cirque despite the apparent stratigraphic correlation of the drifts. This study has implications for drift depositional processes of cold-based glaciers and the importance of source-bedrock lithology and geomorphology on nuclide inheritance in Antarctica.

Seismic mapping of the central and southern segments of the Tanlu fault zone using P-wave receiver functions

The central and southern segments of the Tanlu fault zone, located in the collision boundary between the Yangtze Plate and the North China Craton, underwent complex tectonic deformation associated with the Pacific Plate subduction. The crustal thickness and Vp/Vs ratio are important parameters for comprehending tectonic evolution and geodynamic processes. By integrating a newly dense seismic array, our results based on P-wave receiver function analyses reveal the crust thickness varies significantly across the Tanlu fault zone. The Yangtze Plate is characterized by a thinner crust, contrasting a thicker crust imaged in the North China Craton, indicating the Tanlu fault zone is a prominent structural block boundary. The Vp/Vs ratio within the Tanlu fault zone is visibly higher than its surroundings, probably correlating with abundant fluids activity along the fault zone, as a result of the subduction of the Pacific Plate. Additionally, we observe a gradually deepening Moho discontinuity beneath the Sulu orogenic belt with a low-angle dip, presenting direct seismic evidence for supporting the ancient Yangtze Plate underthrusting beneath the Sulu orogenic belt.

Stress heterogeneity in the eastern Tibetan Plateau and implications for the present-day plateau expansion

The eastward expansion of the Tibetan Plateau has resulted in different earthquake types in the eastern Tibetan Plateau, but the mechanism remains unclear. Here, we construct a three-dimensional visco-elastoplastic finite element model considering the topography to investigate the influence of fault geometry and rheological heterogeneity on stress fields. In our best-fitting model, the minimum principal stress is nearly vertical around the southern Huya fault zone, which is adjacent to the Longmen Shan fault zone, due to the significant mid-lower crust lateral rheological heterogeneity, and the thrust stress regime accounts for the reverse fault and thrust-dominated earthquakes. In this scenario, the eastward horizontal motion of the mid-lower crust is obstructed and facilitates thrust faulting, suggesting the limited eastward expansion of the Tibetan Plateau. In contrast, the northern Huya fault zone, one of the terminal branches of the East Kunlun fault, accommodates the continuous eastward extrusion of the East Kunlun fault, where the stress regime under a more homogenized crust favors the strike-slip faulting process, along with the dominant strike-slip earthquakes. Moreover, the best-fitting of stress regime explains the thrust-dominated 2008 Ms. 8.0 Wenchuan and 2013 Ms. 7.0 Lushan earthquakes on the Longmen Shan fault zone. Combining geophysical and geodetic observations and model analyses, we propose that the hybrid deformation mode in the eastern Tibetan Plateau is accommodated by upper crustal shear and thrusting deformation and mid-lower crustal thickening driven by the gravitational potential energy gradient. Our results elucidate the mechanism for differences in strong historical earthquakes and, more importantly, isolate the effect of fault geometry from those of heterogeneous viscosity on crustal deformation and stress heterogeneity in the eastern Tibetan Plateau.

A Pn-wave spectral inversion technique based on trust region reflective algorithm

This paper suggests utilizing the trust region reflective (TRR) algorithm to address the inverse problem related to seismic spectrum analysis of Pn waves. By applying this approach, it becomes feasible to simultaneously deduce the seismic moment, corner frequency, tomographic model of Pn-wave attenuation, and site responses while incorporating appropriate constraints based on prior knowledge to ensure solution accuracy. The efficacy of this method has been validated using synthetic data. Specifically, the recorded Pn waves from four underground explosions in North Korea were employed to evaluate the performance of the proposed technique. The derived seismic moment and corner frequency for these events were found to align with the characteristics of the seismic data, earthquake magnitude calculations, and empirical data. The resulting tomographic model of Pn-wave attenuation revealed a spatial distribution of high and low attenuation, indicating significant lithospheric heterogeneity in northeastern China. Areas with low Q0 and low η, such as the Xialiaohe Basin and southwestern Changbai Mountain, suggest a high-temperature environment in the upper mantle cap layer or the presence of molten substances. Conversely, regions with high Q0 and high η, like the northern Changbai Mountain and the vicinity of the eastern Tanlu Fault Zone, characterized by high Pn-wave velocity and thick crust, indicate minimal modification or destruction in the lithosphere. Furthermore, site responses were determined for 75 seismic stations in northeastern China, with their characteristics preliminarily analyzed and interpreted in conjunction with geological context. The iterative process based on the TRR algorithm for Pn-wave spectral inversion proposed in this study demonstrates robust convergence and enhances the fitting accuracy of the observed spectrum. The inversion outcomes from this spectral technique yielded synthetic Pn-wave spectra that closely matched the observed spectra of the four underground explosions in North Korea across various frequency bands for over 90 % of the stations.

Effects of Different Shade Treatments on the Epidermal Wax Deposition of Hosta Genotypes with Different Glaucousness of Leaf Surface

Epidermal wax is strategically situated at the interface between plants and air; therefore, it plays a key role in plants’ interactions with their surroundings. It is also unstable and susceptible to light intensity. Hosta plants are shade-loving herbs with admirable flowers and leaves. Hosta ‘Halcyon’ and Hosta ensata F. Maek. are two species of Hosta with a glaucous and a glossy appearance, respectively. Light intensity can affect the composition of epicuticular wax on the leaf surface, which influences the leaf color phenotype and ornamental value. In this paper, the crystal micromorphology, content, and components of epicuticular wax on the leaves of two species of Hosta under different light conditions (10%-, 30%-, 50%-, 70%-, and 100%-intensity sunlight, relative light intensity (RLI)) have been studied using pot experiments. The results indicate that the epicuticular wax crystals of H. ‘Halcyon’ and H. ensata are tubular and platelet-like, respectively. The wax crystals of H. ‘Halcyon’ melted and formed a thick crust under 100% RLI, and those of H. ensata melted and formed a thick crust under 70% and 100% RLI conditions. The primary ingredients of the epicuticular wax of the two species of Hosta contained primary alcohols, alkanes, fatty acids, and esters; β-diketones were only detected in H. ‘Halcyon’. The quantity of epicuticular wax of H. ‘Halcyon’ reduced at first and then increased with an RLI increase, achieving its lowest value at 50% RLI, but that of H. ensata declined little by little. The amounts of C28 primary alcohols, C31 alkanes, and C18 fatty acids were significantly higher than those of other carbon atoms in the two genotypes of Hosta. The C31β-diketones content decreased with the increase in light intensity, which caused the white frost phenotype to gradually weaken in H. ‘Halcyon’.

Driven magmatism and crustal thinning of coastal southern China in response to subduction

Abstract. The late Mesozoic igneous rocks along the coastal South China Block (SCB) exhibit complex parental sources involving a depleted mantle, subducted sediment-derived melt, and melted crust. This period aligns with the magmatic flare-up and lull in the SCB, debating with the compression or extension in coastal region. Our study employs numerical models to investigate the dynamics of the ascent of underplating magma along the Changle–Nan'ao Belt (CNB), simulating its intrusion and cooling processes while disregarding the formational background. The rheological structure of the lithospheric mantle significantly influences magma pathways, dictating the distribution of magmatism. This work reveals that the ascent of magma in the presence of faults is considerably faster than in the absence of faults, and contemporaneous magmatic melts could produce different cooling and diagenetic processes. Additionally, the influence of pre-existing magma accelerated the emplacement of underplating magma. The magma beneath the fault ascended rapidly, reaching the lower crust within 20 million years, with a cooling rate of approximately ∼ 35 °C Myr−1. Conversely, the thickened magma took 40–50 million years to ascend to the lower crust, cooling at a rate of ∼ 10 °C Myr−1. In contrast, magma without thickening and fault would take a considerably longer time to reach the lower crust. The ascent of magma formed a mush-like head, contributing to magmatic circulation beneath the crust and decreasing crustal thickness. Multiphase magmatism increases the geothermal gradient, reducing lithospheric viscosity and promoting underplating magma ascent, leading to magmatic flare-ups and lulls. Our findings suggest that the Cretaceous magmatism at different times in the coastal SCB may be associated with the effects of lithospheric faults under similar subduction conditions. Boundary compression forces delay magma ascent, while rising magma induces a significant circulation, decreasing the crustal thickness of the coastal SCB. This study provides new insights into the complex interplay of magmatic processes during subduction, emphasizing the role of lithospheric structure in shaping the temporal and spatial evolution of coastal magmatism.

Low-velocity middle-and-lower crustal materials blocked by the red river fault in the SE margin of the Tibetan plateau

The ongoing collision between the Indian and Eurasian plates propels the eastward movement of the Tibetan plateau (TP), leading to substantial crustal deformation around the southern Sichuan-Yunnan block (SYB). Using ambient noise data from multiple temporary seismic arrays and permanent stations, we construct a high-resolution regional crustal azimuthally anisotropic Vs model in the SYB. Our new model reveals two significant low-velocity anomalies with strong azimuthal anisotropy near the block boundary faults in the middle-and-lower crust. The extensive low-velocity anomalies around the middle-south segment of the Xiaojiang Fault (XJF) possibly result from partial melting due to spontaneous deformation caused by crustal thickening and increased felsic components, as well as the superimposition of shear heating faults and local upwelling asthenosphere. The N‒S trending low-velocity anomaly at the northwest end of the Red River Fault (RRF) may be associated with weak material migration from the TP, potentially serving as a conduit for mantle upwelling. The azimuthal anisotropy along the block boundary faults exhibits spatial variations linked to segmented distortion resulting from southeastward crustal movement and various geological activities. A key finding is that the crustal channelized low-velocity along the XJF is clearly blocked by the RRF, instead of going through. Notably, the azimuthal anisotropy in the E‒W direction, observed above the Moho and at depths deeper than 30 km in the intersection end, implies the potential intrusion of localized mantle materials into the lower crust. Therefore, lithospheric deformation is significantly affected by block boundary faults and the properties of the crust and mantle.

Subaerial crust emergence hindered by phase-driven lower crust densification on early Earth.

Earth owes much of its dynamic surface to its bimodal hypsometry, manifested by high-riding continents and low-riding ocean basins. The thickness of the crust in the lithosphere exerts the dominant control on the long-wavelength elevations of continents. However, there is a limit to how high elevations can rise by crustal thickening. With continuous crustal thickening, the mafic lower crust eventually undergoes a densifying phase transition, arresting further elevation gain-an effect clearly observed in modern orogenic belts. On early Earth, lower crust densification should also limit how high a thickening crust can rise, regardless of the thickening mechanisms. We suggest that lower crust densification combined with a thicker oceanic crust in the Archean may have limited the whole-Earth topographic relief to 3 to 5 kilometers at most-half that of the present day. Unless the oceans were far less voluminous, limited relief would inevitably lead to a water world on early Earth.

Neoproterozoic High-K calc-alkaline granites forming the large Agudos Grandes and Socorro batholiths, SE Brazil: Crystallization conditions and tectonic implications

I-type, high-K calc-alkaline (HKCA) granites are the main components of two large Neoproterozoic batholiths that occur in different tectonic domains in SE Brazil (Agudos Grandes (emplaced in the central Ribeira Fold Belt) and Socorro (emplaced in the Socorro-Guaxupé Nappe). We investigated representative samples from both batholiths using detailed petrography by optical and electronic microscopy and chemical analysis of main minerals (amphibole, biotite and plagioclase) by EPMA aiming to obtain reliable estimates of magmatic temperatures and pressures. Granites from both batholiths show broadly similar whole-rock compositions and crystallized from oxidized magmas at ∼ NNO+1 to +2. However, they were emplaced at different depths (Agudos Grandes: 3–5 kbar; Socorro: 5–6 kbar; Al-in-hornblende geobarometry) and crystallized at slightly, but consistently different temperatures (respectively, lower and greater than 800 °C). The lower temperatures of the Agudos Grandes granites may account for the presence of titanite as a main accessory mineral, in contrast with its absence, or occurrence as thin rims in ilmenite in the Socorro granites. The pressure estimates are consistent with the metamorphic conditions of the country rocks and suggest that magma emplacement was influenced by a combination of density control and evolving tectonic stress regimes, with the dominantly more mafic rocks of Socorro emplaced at greater depths. On the other hand, higher Sr/Y ratios and REE fractionation shown by the Agudos Grandes whole-rocks may result from magma equilibration in a thicker crust (>60 km, as contrasted with ∼50 km for Socorro) at the time of magma generation. However, the possible presence of cumulate components (as suggested, for instance, by high apatite and hornblende saturation temperatures) requires that these contrasts are tested in the future by proxies that are independent of whole-rock composition (e.g., trace-element chemistry of minerals).

Magma storage conditions of Lascar andesites, central volcanic zone, Chile

Abstract. Lascar volcano, located in northern Chile, is among the most active volcanoes of the Andean Central Volcanic Zone (CVZ). Its activity culminated in the last major explosive eruption in April 1993. Lascar andesites which erupted in April 1993 have a phase assemblage composed of plagioclase, clinopyroxene, orthopyroxene, Fe–Ti oxides, and rhyolitic glass. To better constrain storage conditions and mechanisms of magmatic differentiation for andesitic magmas in a thick continental crust, crystallization experiments were performed in internally heated pressure vessels at 300 and 500 MPa, in the temperature (T) range of 900–1050 °C, at various water activities (aH2O) and oxygen fugacities (logfO2 between QFM+1.5 and QFM+3.3 at aH2O =1; QFM is quartz–fayalite–magnetite oxygen buffer). The comparison of experimental products with natural phase assemblages, phase compositions, and whole-rock compositions was used to estimate magma storage conditions and to reconstruct the magma plumbing system. We estimate that Lascar two-pyroxene andesitic magmas were stored at 975±25 °C, 300±50 MPa, and logfO2 of QFM+1.5±0.5, under H2O-undersaturated conditions with 2.5 wt % to 4.5 wt % H2O in the melt. The geochemical characteristics of the entire suite of Lascar volcanics indicates that a fractionating magmatic system located at a depth of 10–13 km is periodically replenished with less evolved magma. Some eruptive stages were dominated by volcanic products resulting most probably from the mixing of a mafic andesitic magma with a felsic component, whereas compositional variations in other eruptive stages are better explained by crystal fractionation processes. The relative importance of these two mechanisms (mixing vs. crystal fractionation) may be related to the amount and frequency of magma recharge in the reservoir.

Metamorphic and chronological records of early Permian continental collision in Beni Bousera granulites (Internal Rif belt): implications for the tectonic evolution of northern Morocco during Pangea construction

This work presents petrological data obtained from shielded mineral inclusions within large garnets from the Beni Bousera metamorphic unit (internal Rif belt, northern Morocco) combined with in situ U‒Th‒Pb dating of monazite inclusions. In the considered Beni Bousera metapelites, the occurrence of mineral inclusions of kyanite + rutile + plagioclase + K-feldspar + quartz ± graphite trapped in garnet cores is diagnostic of high-pressure granulite-facies metamorphism. Geothermobarometry combined with thermodynamic modelling yields minimal P ‒ T conditions of approximately 1.37–1.5 GPa and 780–795°C, consistent with the development of palaeogradients during crustal thickening in collision belts. A concordant U/Pb age of 281 ± 3 Ma is obtained from associated monazite inclusions in the same garnet cores. These new petrochronological data attest to the early Permian continental collision in northern Morocco. The obtained dataset is constrained with petrological and geochronological parameters available on other neighbouring segments of the Variscan belt. Within the framework of recently unified full plate reconstruction models, the presented tectonic scenario suggests that all these orogenic segments were built during the late Carboniferous–early Permian crustal thickening at the northwestern Gondwana margin in response to the convergence of Laurentia and Gondwana during the late stages of Pangaea construction.

Crustal structure of the Bushveld complex, South Africa from 1D shear wave velocity models: Evidence for complex-wide crustal modification

Thirty-nine 1D shear wave velocity profiles, obtained by jointly inverting receiver functions and Rayleigh wave group velocities, are used to investigate the crustal structure of the Bushveld Complex in northern South Africa. Data from teleseismic earthquakes recorded on broadband seismic stations between 1997 and 1999 and 2015–2020 were used to compute P-wave receiver functions. Rayleigh wave group velocities between 5 and 30 s period were obtained from an ambient noise tomography and combined with group velocities between 30 and 60 s period from a published continental-scale surface wave tomography model. Moho depths of 45–47.5 km are found under the center of the complex compared to 40 km thick crust, on average, surrounding the complex, indicating ∼5–7 of crustal thickening. The bottom ∼10 km or more of the lower crust across much of the Bushveld Complex has a Vs ≥ 4.0 km/s, consistent with a mafic composition, whereas in most areas around the margins of the complex the thickness of the mafic lower crust is much less than 10 km. In the upper crust higher velocity structure (Vs > 3.6 km/s) above 15 km depth underlain by lower velocity structure is seen in many locations, suggesting the presence of mafic/ultramafic layering. These results favor the continuous-sheet model for the structure of the Bushveld Complex because the ensemble of 1D models is characterized by three diagnostic features consistent with that model: (1) thicker crust under the center of the complex than away from the complex; (2) a greater thickness of high-velocity (i.e., mafic) layering in the lower crust under the complex compared to away from the complex; (3) high-velocity (i.e., mafic/ultramafic) layering in the upper crust beneath much of the complex. The lack of upper crustal mafic/ultramafic layering beneath some parts of the complex is consistent with the post-emplacement tectonic and magmatic history of the complex.

Rapid Impact Crater Relaxation Caused by an Insulating Methane Clathrate Crust on Titan

Titan’s impact craters are hundreds of meters shallower than expected, compared to similar-sized craters on Ganymede. Only 90 crater candidates have been identified, the majority of which have low certainty of an impact origin. Many processes have been suggested to shallow, modify, and remove Titan’s craters, including fluvial erosion by liquid from rainfall, aeolian sand infill, and topographic relaxation induced by insulating sand infill. Here we propose an additional mechanism: topographic relaxation due to an insulating methane clathrate crustal layer in Titan’s upper ice shell. We use finite element modeling to test whether a clathrate crust 5, 10, 15, or 20 km thick could warm the ice shell and relax craters to their currently observed depths or remove them completely. We model the viscoelastic evolution of crater diameters 120, 100, 85, and 40 km, with two initial depths based on depth−diameter trends of Ganymede’s craters. We find that all clathrate crustal thicknesses result in rapid topographic relaxation, despite Titan’s cold surface temperature. The 5 km thick clathrate crust can reproduce nearly all of the observed shallow depths, many in under 1000 yr. A 10 km thick crust can reproduce the observed depths of the larger craters over geologic timescales. If relaxation is the primary cause of the shallow craters, then the clathrate thickness is likely 5–10 km thick. Topographic relaxation alone cannot remove craters; crater rims and flexural moats remain. To completely remove craters and reproduce the observed biased crater distribution, multiple modification processes must act together.

Controlling factors of Riedel shear spacing in the simple shear mode of strike-slip fault: Insights from sandbox models

Strike-slip faults generally develop Riedel shears (R-shears), which exhibit parallel and evenly-spaced distribution characteristics. However, the factors controlling the R-shear spacing in strike-slip faults are still unclear. The influence of material properties such as internal friction angle and cohesion, basal friction, and the thickness of brittle layers (T) on the R-shear spacing (S) are investigated using analogue models in this paper. Research findings indicate that the internal friction angle of the material and the thickness of the brittle layer have a significant impact on the R-shear spacing, with the thickness of the brittle layer directly determining the R-shear spacing as evidenced by their linear correlation. In comparison, cohesion and basal friction have insignificant effects on R-shear spacing. Based on this, experiments were carried out using various thicknesses of brittle layers (specified materials) to investigate the impact of the brittle layer thickness on the R-shear spacing, and Particle Image Velocimetry (PIV) is used to analyze the distribution pattern of R-shear development at each stage. The results indicate that fractures occur in regions where the vorticity field alternates between positive and negative values, and as the evolution progresses, the maximum strain gradually converges towards the center of the deformation zone, leading to a reduction in the activity of the R-shear, while the spacing of the R-shear remains unaltered. The normalized (S/T) results indicate that the experimental value of 1.32 aligns with natural laws and is very close to the normalized value of the natural faults, which is 1.24. It can be inferred that the thickness of the seismogenic crust within the range of the Altyn Tagh Fault is 40.9–43.5 km.