Thermal Upwelling Research Articles

Imaging upper lithospheric structures of the Benue Trough and adjoining basement areas in Nigeria and Cameroon from satellite gravity data

This study investigates the structures of upper lithosphere in the Benue Trough (BT) and its abutting basements in Nigeria and Cameroon through detailed analyses and interpretation of satellite Bouguer gravity anomaly (BA) data. Matched bandpass filtering, edge enhancements (tilt- and theta-derivatives), 3D inversion, and 2D modelling methods were applied to the BA data. The result from matched bandpass filtering reveals four gravity layer sources with their respective average depths, namely near-surface geological bodies (c. 0.6274 km), depth to basement (sediment thickness; c. 5.8921 km), intra-crustal boundary (c. 12.7538 km) and Moho (c. 27.980 km). Tilt- and theta-derivatives of the BA imaged numerous regional tectonic fabrics/structures exhibiting ENE – WSW, NNE –SSW, and NE – SW major trends coupled with N – S, NW – SE and E − W minor trends across the area, with the major regional tectonic trends mostly dominant in the BT and Adamawa Plateau/Cameroon Volcanic Line (AP/CVL) regions. Positive tilt cum positive theta derivatives of BA indicate that high-density (dense) geological bodies occur within the underlying sediments and crustal basements, particularly within the BT. 3D inversion result reveals a conspicuous and elongated (c. > 600 km long, 120 km wide) NE – SW orientated shallow (c. 23–30 ± 3 km) Moho structure beneath the BT and a prominent elongated, large (c. > 550 km long, 115 km wide) NE – SW striking deep (c. 34–44 ± 3 km) Moho structure below the AP/CVL region in the Western Cameroon Domain/Adamawa-Yade Domain (WCD/AYD) of Cameroon. Also, it shows that the Northern Nigerian Basement Complex (NNBC) and Eastern Nigerian Basement Complex (ENBC) regions are characterised by highly variable Moho morphology and depths (c. 26 to 36 ± 3 km). Generated 3D Moho model represents the first detailed Moho model for this studied area. The result from 2-D forward modelling of a NW – SE representative profile, constrained by known outcrop regional geology and available nearby seismic/teleseismic information, indicates the presence of two variable sedimentary sub-basins (c. 2–5 ± 0.4 km deep), a basement ridge and shallow Moho (c. 23–27 km ± 2 km) in BT region with occurrences of shallow intrusive bodies (e.g. granites), dense lower crustal body and deep Moho (c. 30–40 ± 2 km) underneath the basement regions, particularly in the AYD region. It also highlights that the underlying crust in BT is thinnest while thickest in basement regions of NNBC, ENBC, WCD (particularly AP/CVL area) and AYD along the NW – SE profile. The Moho depth values obtained in this study are within ±3.5 km with all those previously published from teleseismic works in the studied area. The imaged major NE – SW and ENE –WSW orientated tectonic fabrics/structures of the Precambrian crustal basement may have strongly influenced the location, emplacement and orientations of the BT, AP/CVL and Younger Granite complexes. Also, the BT developed in response to mechanical stretching (extension) of the lithosphere driven by regional far-field extensional stresses originating from remote plate tectonic/boundary forces coupled with some active-related effects (e.g. thermal upwelling of mantle materials) caused by changes in the dynamics of the sub-lithospheric mantle during the rifting in the Jurassic – Cretaceous.

Ambient Noise Tomography of Northern Borneo Reveals Evidence of Subduction and Post‐Subduction Processes

AbstractThe region of northern Borneo in South East Asia sits within a post‐subduction setting formed by the recent termination of two sequential but opposed subduction systems. In this study we use seismic data from a recent temporary array deployment to image the crustal velocity structure beneath northern Borneo using a two‐stage Bayesian trans‐dimensional tomography scheme, in which period dependent phase velocity maps are first generated, and then used to build a 3‐D shear wave model through a series of 1‐D inversions. In the second stage, we also apply an Artificial Neural Network to solve the 1D inverse problem, which results in a smoother 3‐D model compared to the TransD approach without sacrificing data fit. Our shear wave velocity model reveals a complex crustal structure. Under the Crocker Range, a heterogeneous velocity structure likely represents remnants of early Miocene subduction, including underthrust continental crust from subsequent continent‐continent collision. In the east we observe high velocities that are interpreted to be igneous rocks in the crust generated by melting due to mid Miocene Celebes Sea subduction and later decompression melting as well as a low velocity zone that could represent underthrust sediment or duplexes from Celebes Sea subduction. A low velocity zone in the lower crust is present in a region of apparent crustal thinning. Our preferred explanation for this anomaly is remnant thermal upwelling within a failed rift that represents the on‐shore continuation of the extension of the Sulu Sea, most likely caused by rollback of the Celebes Sea slab.

Cenozoic Volcanic Provinces and Shallow Asthenospheric Volumes in the Circum‐Mediterranean: Evidence From Magmatic Geochemistry, Seismic Tomography, and Integrated Geophysical‐Petrological Thermochemical Inversion

AbstractWe compile an extensive catalog comprising geochemistry and ages of Cenozoic volcanic provinces in the Mediterranean region, distinguishing between three groups according to the geochemistry of magmatic rocks: intraplate (IVP), subduction‐related (SRVP), and mixed‐origin volcanic provinces (MVP; intraplate with subduction imprint). In order to relate their spatial distribution to properties of the lithosphere‐asthenosphere system, we determine temperature‐depth profiles by integrated geophysical‐petrological inversion at representative locations, using Rayleigh and Love wave phase velocities, heat flow, rock densities, and accounting for thermochemical conditions. The results confirm the occurrence of thin lithosphere (<100 km) above warm asthenosphere (>1,300°C) in areas of low shear‐wave velocities in the shallow upper mantle. Nine shallow asthenospheric volumes (SAVs) between 70 and 300 km depths are identified, forming a partly interconnected belt across the Circum‐Mediterranean. A remarkable colocation exists between the SAVs and the intraplate and mixed‐origin volcanic provinces (IMVPs). Whereas dense networks of IMVPs have formed above the SAVs, IMVPs are absent in areas of thick mantle lithosphere. Magmatic activity in IMVPs at 60–70 Ma indicates that several SAVs existed already in the Paleogene (Central European, Adriatic, Western Mediterranean, and Moesian SAVs). The formation of SAVs is related either to asthenospheric upwelling caused by slab rollback and back‐arc extension (Aegean‐Anatolian, Moesian, Pannonian, Western Mediterranean SAVs), or to thermal upwelling (Adriatic, Central European, Middle East [MEA], North African, Rhine‐Rhone SAVs), with some of the latter coupled partly with continental rifting (Central European, MEA, North African, Rhine Rhone SAVs).

Topography of the 410 and 660 km discontinuities beneath the Tibetan plateau and adjacent areas

To investigate the three-dimensional spatial distribution of subducted oceanic slab segments and their consequential effect on the thermal and hydrous composition of the mantle transition zone (MTZ) beneath the entire Tibetan Plateau, we conducted a comprehensive analysis of approximately 655,000 P-to-s receiver functions (RFs), obtained from 735 broadband seismic stations. These RFs were utilized to delineate the 410 km (d410) and 660 km (d660) seismic discontinuities, which represent the upper and lower boundaries of the MTZ, respectively. The RFs were grouped within circular bins with a radius of 1 degree and stacked to image the discontinuities. The mean apparent depths derived based on the 1-D IASP91 Earth model for the d410 and d660 across the entire study area are 412.2±8.3 km and 668.9±8.5 km respectively, and the MTZ thickness is 256.5±6.9 km. The observed apparent depths underwent subsequent correction utilizing multiple velocity models. Several areas in central Tibet exhibit a normal d410 and an anomalously deep d660, which can be attributed to the combined effect of the negative thermal anomaly and dehydration associated with subducted slab segments that have penetrated at least to the d660 depth. The anomalous thickening of the MTZ beneath the southeastern Tibetan Plateau surrounding the Tengchong volcanic field can be explained by the dehydration of the subducted Indian Slab. Significant thinning of the MTZ associated with the deepening of the d410 beneath the western Tian Shan may indicate active thermal upwelling originating from the MTZ.

Enhanced coastal upwelling indices for Moroccan Atlantic Coast, their force–response variability and sensitivity to extreme events between 1993 and 2021

The Eastern Boundary Upwelling Systems (EBUS) are highly productive ecosystems globally, contributing approximately 20% of the world’s fishery resources and bearing significant ecological and socioeconomic importance. These EBUSs are primarily influenced and sustained by seasonal equatorward wind patterns. Along the Atlantic coast of Morocco, there are two distinct and persistently occurring upwelling patterns: a robust one to the south of (26°N) and a weaker one to the north. In this study, we present a coastal upwelling (transport) index based on combined contributions from both Ekman and Geostrophic transports, complemented by a recently refined thermal upwelling index based on sea surface temperature (SST). The seasonality patterns of these indices align with earlier research findings, revealing a sequence of wind-to-chlorophyll-to-SST lag relationship at all latitudes (24°N, 27°N and 32°N) during the period from 1993 to 2021. Subsequently, the two indices and Chlorophyll-a signature and three potential basin-scale driving factors, the North Atlantic Oscillation (NAO), the Azores High and the East Atlantic Oscillation (EAO) are analyzed to explore the interannual variability of the upwelling system and connect them to those of the prevailing wind patterns. Our analysis reveals three prominent periodicity bands: one near-decadal cycle spanning 8–10 years and two interannual cycles of 4–5 years and 2–3 years. It appears that the interannual bands are influenced by variations in the amplitude of the Azores high center pressure and are associated with the EAO, while the near-decadal cycles may be linked to the NAO. Furthermore, our examination of extreme events suggests that high (low) upwelling years tend to coincide with low (high) EAO years.

Constraints from lamprophyre petrogenesis on the timing of Eocene lithospheric thinning and associated rifting of Borneo and Sulawesi

The petrogenetic history of well characterised magmatic rocks may contain invaluable information of paleo-lithospheric thickness and structure that can inform on tectonic evolution. Petrographically distinctive rocks, such as lamprophyres, provide temporal and spatial constraints to build, assess, and verify tectonic models.An alkaline lamprophyre dyke that intrudes the rift-related sediments of the Eocene Tanjung Formation of the Senakin Peninsula in South Kalimantan, Indonesia has been investigated using petrology, whole rock geochemistry and Ar-Ar geochronology.The Senakin lamprophyre dyke is mostly geochemically homogenous but petrographically stratified with vesicle free monchiquite alkaline lamprophyre mineralogy at the base and an amygdale bearing kaersutite phenocryst free upper zone. The major and trace element chemistry is consistent with an intra-plate Ocean Island Basalt like affinity. 40Ar/39Ar dating of hornblende and groundmass aliquots indicates an emplacement age of 37.45 ± 0.10 Ma (2σ), which is consistent with the Late Eocene age of the Tanjung Formation derived from microfossil analysis.The petrology and geochemistry of the Senakin lamprophyre dyke is consistent with magma genesis at > 2GPa beneath a thick lithosphere followed by fractional crystallisation at Moho depths equivalent to 0.5–0.8 GPa. This indicates a much thicker southeastern Kalimantan paleo-lithosphere. Lamprophyre magmatism was initiated by lithospheric extension associated with the opening of the Makassar Strait and also regional thermal upwelling. This extension continued to thin the lithosphere to its current thickness equivalent to 0.5 GPa.

Warm versus cold crust in the Tien Shan orogenic belt revealed by seismic Lg attenuation tomography

SUMMERY Due to the far-field effect of the India–Eurasia collision, the Tien Shan orogenic belt has been undergoing reactivation and modification. Two end-member models of the geodynamic mechanisms are (1) surface uplift due to crustal shortening caused by lithospheric compression and (2) mountain formation resulting from thermal upwelling of asthenospheric mantle materials generated by lithospheric subduction. However, the topography along the Tien Shan orogenic belt changes significantly, and the deep structure and dynamic process are quite different beneath the Tien Shan orogenic belt from both geological and geophysical observations. Therefore, the reactivation and modification of the Tien Shan orogenic belt are likely influenced by both geodynamic mechanisms, which also generate various thermal anomalies in the crust. Seismic Lg-wave attenuation is very sensitive to crustal material composition and status and can point to the presence of partial melting within the crust resulting from mantle upwelling. In this study, we develop a high-resolution Lg-wave attenuation model between 0.05 and 10.0 Hz in Northwest China and use lateral attenuation variations to infer thermal structures in the crust. The central Tien Shan is characterized by prominent low-QLg anomalies, whereas relatively high-QLg distributions are imaged beneath the eastern and western Tien Shan. The surface uplift and crustal deformation are mostly related to the upwelling of hot mantle materials in the central Tien Shan and are likely induced by lithospheric compression in the eastern and western Tien Shan. However, low-Q anomalies are observed in the junction between the Pamir Plateau and western Tien Shan, indicating that the uplift in the south of the western Tien Shan is related to thermal subduction-induced upwelling and intrusion into the crust due to the collision between the Indian and Eurasian plates. The Kazakh Shield, characterized by pronounced high-QLg values, is likely a cold and hard terrane, and hence blocks the far-field effect of the India–Eurasia collision.

Topographic variations of mantle transition zone discontinuities in the Java subduction zone and adjacent regions from receiver function analyses

An unprecedented number of broadband seismic data from 165 stations were utilized to constrain the topographic variations in seismic discontinuities at the top and bottom of the mantle transition zone (MTZ) in the Java subduction zone. After moveout corrections under a nonplane wavefront assumption and time-to-depth conversions based on the IASP91 Earth model, the depth undulations of two discontinuities that border the MTZ (d410 and d660) were obtained by stacking P receiver functions that were grouped into circular bins with a radius of 1°. The apparent depths exhibit systematically spatial variations that range from 385 to 427 km for the d410 and 608 to 700 km for the d660, which result in a MTZ thickness of 217–309 km. The mean values are 406 ± 12, 665 ± 19, and 258 ± 18 km, suggesting thicker than normal MTZ by an average of 8 km. The major factor that accounts for the systematic depth variations across the study area originates from the subducted slab that perturbs the entirety of the MTZ, which induces a negative thermal anomaly, high velocity zone, and abnormal water presence concentrated in the lower MTZ of western Java. In Lombok and Sumbawa Islands, after the thermal and velocity heterogeneities from the subducted slab were considered, the remnant d660 uplift is found due to a significant thermal anomaly with a magnitude of up to +200 K near the d660. In eastern Java and adjacent regions, advective thermal upwelling associated with toroidal mantle flow that traversed a slab window in central Java induces a low velocity anomaly with a magnitude that is greater than the high velocity anomaly caused by the subducted slab, which apparently depressed the d410 and d660 by the same amount.

Continental tapering of South America caused asymmetric non-uniform opening of the South Atlantic Ocean

South Atlantic opening has been typically modelled as being related to symmetric and static thermal upwelling and seafloor spreading that drive divergent continental drift of South America and Africa. Comparative analyses, however, show that South Atlantic opening is asymmetric and non-uniform. For neither asymmetric nor non-uniform opening are the underlying mechanisms clear. Here I use geological and geophysical data to inform analytical modelling, revealing that westward drifting and southward tapering of the South American continent have controlled the asymmetry and the non-uniformity in South Atlantic opening. I interpret that the asymmetric non-uniform seafloor spreading caused the ridge and hotspots to migrate, leaving behind non-linear seamount trails that are indicative of the speed of hotspot migration rather than direction of plate movement. The findings point towards a chain reaction from continental drifting, through seafloor spreading to ridge-hotspot interaction, which is instrumental in understanding the geodynamics for global plate tectonics.

Receiver Function Mapping of the Mantle Transition Zone Beneath the Tian Shan Orogenic Belt

AbstractThis study provides new constraints on the mantle transition zone (MTZ) structure under the Tian Shan orogenic belt in the context of the double‐sided subduction of the Junggar lithosphere and the Tarim lithosphere. The 410‐ and 660‐km discontinuities bordering the MTZ under the Tian Shan and its adjacent areas are mapped by stacking 68,361 receiver functions from 4,122 events recorded by 100 broadband seismic stations. Regional and large‐scale 3‐D velocity models are used to explain the effect of the laterally heterogeneous velocity anomalies. We identify a thickened MTZ of about 9.3 km under the eastern Tian Shan and the Darbut belt. This thickening indicates a lower temperature, which is correlated with the broken‐off subducted lithosphere or the delamination of the eastern Tian Shan lithosphere. We document the thinning of the MTZ caused by the depressed 410‐km discontinuity and the uplifted 660‐km discontinuity across the Tian Shan orogenic belt, Junggar Basin and Tarim Basin. This thinning corresponds to a + 100 K thermal anomaly in the MTZ, and we suggest that it may have been produced by thermal upwelling originating from the lower mantle. The slightly thinned MTZ beneath the Altai Mountains was caused by the uplift of these two discontinuities, which is related to the high‐velocity anomalies in the upper mantle and MTZ. In contrast, the thinned MTZ below the easternmost segment of the Tian Shan orogenic belt, formed by the depression of 6.5 km for the 410‐km discontinuity may indicate the presence of small‐scale mantle upwelling.

Fingerprinting secondary mantle plumes

Many vertical seismic velocity anomalies observed below different parts of the Eurasian plate are rooted in the transition zone between the upper and lower mantle (410–660 km), forming so-called secondary plumes. These anomalies are interpreted as the result of thermal effects of large-scale thermal upwelling (primary plume) in the lower mantle or deep dehydration of fluid-rich subducting oceanic plates. We present the results of thermo-mechanical numerical modelling to investigate the dynamics of such small-scale thermal and chemical (hydrous) anomalies rising from the lower part of the Earth's upper mantle. Our objective is to determine the conditions that allow thermo-chemical secondary plumes of moderate size (initial radius of 50 km) to penetrate the continental lithosphere, as often detected in seismo-tomographic studies. To this end, we examine the effect of the following parameters: (1) the compositional deficit of the plume density due to the presence of water and hydrous silicate melts, (2) the width of the weak zone in the overlying lithosphere formed because of plume-induced magmatic weakening and/or previous tectonic events, and (3) a tectonic regime varied from neutral to extensional. In our models, secondary plumes of purely thermal origin do not penetrate the overlying plate, but flatten at its base, forming “mushroom”-shaped structures at the level of the lithosphere-asthenosphere boundary. On the contrary, plumes with enhanced density contrast due to a chemical (hydrous) component are shown to be able to pass upwards through the lithospheric mantle to shallow depths near the Moho when (1) the compositional density contrast is ≥ 100 kg m−3 and (2) the width of the lithospheric weakness zone above the plume is ≥ 100 km. An extensional tectonic regime facilitates plume penetration into the lithosphere but is not mandatory. Our findings can explain observations that have long remained enigmatic, such as the “arrow”-shaped zone of low seismic velocities below the Tengchong volcano in south-western China and the columnar (“finger”-shaped) anomaly within the lithospheric mantle discovered more than two decades ago beneath the Eifel volcanic fields in north-western Germany. It appears that a chemical component is a characteristic feature not only of conventional hydrous plumes located over presently downgoing oceanic slabs, but also of upper mantle plumes in other tectonic settings.

Heat Flow and Thermal Source of the Xi’an Depression, Weihe Basin, Central China

The Xi’an Depression of the Weihe Basin, located in the transition zone where the North China, Qinling and Yangtze plates collide with each other, is an important area of geothermal energy utilization in China. Studies of heat flow and thermal sources are of great significance to the exploration and development of geothermal resources in this area. In this paper, six temperature logs boreholes, and 14 thermal conductivity samples have been used to study the geothermal gradient and terrestrial heat flow in the area. The results show that the geothermal gradients of Xi’an Depression range from 20.8 C/km to 49.1 C/km, with an average of 31.7 ± 5.0 C/km. The calculated heat flow is 59.4–88.6 mW/m2, and the average value is 71.0 ± 6.3 mW/m2, which indicates a high thermal background in the area. The high anomalous zones are near the Lintong-Chang’an Fault zone in the southeast, the Weihe Fault in the north, and the Fenghe Fault in the central Xi’an Depression. These deep and large faults not only control the formation of the Xi’an Depression but also provide an important channel for the circulation of groundwater and affect the characteristics of the shallow geothermal distribution. The temperature of the Moho in the Xi’an Depression ranges from 600 to 700°C, and the thermal lithosphere thickness is about 90–100 km. The characteristics of lithospheric thermal structure in Xi’an Depression indicate that the higher thermal background in the study area is related to lithospheric extension and thinning and asthenosphere thermal material upwelling.

Mesoscale patterns in barnacle reproduction are mediated by upwelling-driven thermal variability

Environmental variables are known to regulate the reproductive output of marine intertidal organisms, but typically these variables are studied as averages and interpreted at a macroscale level. Along 200 km of coast in NW Iberia, great variability in the reproductive activity of the stalked barnacle Pollicipes pollicipes was found among 7 different locations. We found the highest number of broods reported to date in this species and suggest a more realistic method of predicting reproductive success that takes into account sea surface temperature (SST). At these same locations, we studied 13 yr of SST, thermal upwelling index (UI), ground skin temperature and chlorophyll a (chl a) satellite data, using spectral analyses to partition their temporal variability over ecologically relevant time scales. SST played the most relevant role as an environmental driver, explaining 48% of the variability in the proportion of breeding individuals (BI), but the SST-BI goodness-of-fit decreased sharply northwards. Variance-partitioning analyses indicated that cycles between 20 and 100 d in SST and UI were more important southwards, which is consistent with a latitudinal gradient in upwelling intensity and frequency along this coast. Thus, we found better biophysical coupling towards the south, where shorter fluctuation time scales in SST match the gonadal development period (around 1 mo). This pattern may explain the spatial variability in the strength of association between key environmental variables and the reproductive cycle of coastal species along their distribution range.

Phosphate Scavenging During Lava‐Seawater Interaction Offshore of Kīlauea Volcano, Hawaii

AbstractThe 2018, subaerial eruption of Kīlauea volcano, Hawaii, resulted in a 5‐km‐long stretch of coastline that actively drained lava into the ocean. Nutrients were added to the surrounding ocean through the dissolution of basaltic rock and thermal upwelling of deep water, thereby fueling a large phytoplankton bloom. Lava‐impacted, surface seawater had high suspended particle loads, and concentrations of chlorophyll, silicic acid, phosphate (Pi), nitrate, and iron that were elevated up to 12, 36, 5, 960, and 1,400 times, respectively, above the background oligotrophic levels. Widespread precipitation of iron oxyhydroxides (Feox) led to extensive scavenging of the dissolved Pi pool, similar to what occurs along mid‐ocean ridge hydrothermal systems. This scavenging transformed a “fertilization” event into a Pi sink near the coast of the ocean entry; however, nutrient data from outside the bloom suggest that Pi could also desorb from the Feox as it is dispersed into the open ocean. From lava quench experiments, we estimate that the hydration state of the Feox precipitate (H2O/Fe) was 5.2–5.7, and that the equilibrium partition coefficient of Pi into Feox (solid/liquid) was 106. In addition, 33Pi radiotracer incubations were used to differentiate between biotic and abiotic uptake of Pi at Kīlauea's ocean entry. These findings are important for understanding modern‐day volcanic fertilization events, modeling nutrient dynamics during major events in Earth history (such as oxygenation of the atmosphere and the formation of large igneous provinces), and predicting the marine response to greater continental weathering in a warming climate.

Effects of Thermal Water Upwelling on Microclimate Change in the High Geo-Temperature Roadway

Deep-circling thermal water upwelling and trickling to high geo-temperature roadway obviously alter the microclimate in mines, which brings difficulty to the prediction of airflow temperature and humidity. This is the basis of air-conditioning cooling load calculation. The heat and mass transfer between trickling water and airflow is rather complicated. Moreover, humid air exhibits the accumulation effect of heat and humidity in the long-distance flow process. In this paper, an apparatus was designed and developed to explore the influence of thermal water trickling on the airflow thermal parameters of a section of roadway (1L–39L, in which 1L–9L is the trickling section). The results show the following (1) With the rise of trickling water temperature, the total enthalpy difference of dry air in the roadway increases within a small range and that of humid air goes up nonlinearly. Besides, the increase of trickling water flow rate has an insignificant effect on the sensible heat of the airflow, while it plays a notable role in increasing the latent heat of the airflow. (2) High trickling water temperature results in a higher growth rate of humidity ratio at 19L than those at 29L and 39L in the early stage of thermal water trickling. Meanwhile, sensible heat exchange, which becomes strong after thermal water trickles for over 30 min, complicates the enthalpy difference variation rates of wet air at the three measuring points. (3) The three measuring points in the 19L–39L section all display a process of enthalpy growth with time. In the case of point 39L, the enthalpy difference of humid air surges sharply when the trickling water temperature is 80 ° C or the flow rate is 200 ml/min. The research results boast some reference value for thermal water management and microclimate change forecasting after the airflow passes through a trickling roadway.

Early Silurian granitic rocks and associated enclaves as evidence of rapid cooling in a cognate magma system: the case of the Xuehuading–Panshanchong pluton, South China Block

AbstractThe study of enclaves in granitic plutons provides fundamental information on the petrogenesis of their host rocks. Here we combine U–Pb zircon ages, petrography, geochemistry and Nd–Hf isotope composition to investigate the origin of dioritic–granodioritic enclaves and their host granodiorites and biotite granites in the Xuehuading–Panshanchong area, which is a pivotal site to study the Palaeozoic intracontinental orogenic processes of the South China Block. Obtained ages indicate that the host rocks were formed in early Silurian time (c. 432 Ma). The enclaves are fine grained, but with mineral assemblages similar to their hosts and contain amphibole, biotite and plagioclase. All rocks have fractionated rare earth element patterns ((La/Yb)N = 2.86–8.16), except for one biotite granite that has a concave rare earth element pattern ((La/Yb)N = 1.50). Most rocks are depleted in Ta–Nb–Ti, and have negative Eu anomalies and ϵNd(t) (–8.86 to –5.75) and zircon ϵHf(t) (–13.30 to –4.11, except for one, –39.08). We interpret that the enclaves were formed at the borders of magma-ascending conduits, where the mafic mineral crystallization was enhanced by rapid cooling. Conversely, the biotite granites were produced by fractional crystallization from a related granodiorite magma. The sample with a concave rare earth element pattern may have been influenced by hydrothermal fluid–melt interaction. Geochemical modelling suggests that the granodiorites were likely generated by disequilibrium melting of heterogeneous amphibolites in the middle–lower crust. Considering the geological data for the Palaeozoic magmatic rocks in the South China Block, we propose that the Xuehuading–Panshanchong magmatism was likely triggered by piecemeal removal of the thickened lithospheric root and subsequent thermal upwelling of mantle, without a mantle-derived magma contribution to the granites.

Topography of the 410 and 660 km Discontinuities Beneath the Cenozoic Okavango Rift Zone and Adjacent Precambrian Provinces

AbstractBy taking advantage of the recent availability of a broadband seismic data set from Networks NR and BX covering the entire country of Botswana, we conduct a systematic receiver function investigation of the topography of the 410 and 660 km discontinuities beneath the incipient Okavango rift zone (ORZ) in northern Botswana and its adjacent Archean‐Proterozoic tectonic provinces in southern Africa. Similar to a previous mantle transition zone (MTZ) discontinuity study using data from a 1‐D profile traversing the ORZ, a normal MTZ thickness is observed in most parts of the study area. This is inconsistent with the existence of widespread positive thermal anomalies in the MTZ and further implies that active thermal upwelling from the lower mantle plays an insignificant role in the initiation of continental rifting. The results also suggest that cold temperature presumably associated with thick cratonic keels has indiscernible influence on the thermal structure of the MTZ. The expanded data set reveals several isolated areas of slight (~10 km or smaller) MTZ thinning. The largest of such areas has a NE‐SW elongated shape and is mostly caused by relative deepening of the 410 km discontinuity rather than shallowing of the 660 km discontinuity. These characteristics are different from those expected for a typical mantle plume. We speculate that the thinner‐than‐normal MTZ may be induced by minor thermal upwelling associated with late Mesozoic‐early Cenozoic lithospheric delamination, a recently proposed mechanism that might be responsible for the high elevation of southern Africa.

On the Change in the Ocean Surface Temperature in the Benguela Upwelling Region. Part I: Season Cycle

Based on the daily satellite imagery data on the ocean surface temperature for 1982–2017, near-surface wind for 1988–2017 and sea-level anomalies for 1993–2017, the seasonal variability of the hydrothermodynamic characteristics of the upper water layer in the vicinity of the Benguela upwelling is investigated. It is shown that the thermal upwelling index averaged over the entire area with lower temperature values does not provide a correct idea of the seasonal course of the water lifting rate in the upwelling zone due to the significant horizontal advection of waters of upwelling origin. The seasonal variation in the vertical velocity of wind origin in the Benguela upwelling zone is characterized by the presence of two extremums from October to March, which is manifested in the dominance of the semiannual harmonic over the annual one. At the same time, the thermal upwelling index in the zone of distribution of upwelling waters is subject to seasonal variability with an annual period.

On the changes in the sea surface temperature in the benguela upwelling region. Part 2: the long-term tendencies

The paper examines the issue on the long-term trends in the sea surface temperature (SST) in the Benguela upwelling zone and their causes using the daily SST satellite data for 1985–2017’s and the daily near-surface wind for 1988–2017”s. It is shown that in the Benguela upwelling region, there is a significant intensification of driving winds in the last 20 yrs. This is accompanied by a decrease of the thermal upwelling index (taking into account the sign of the index or an increase of its absolute values) in the southern part of the Benguela upwelling, but practically does not influence this indicator in its northern part. The likely reason for this difference is the change in the wind field structure, as a result of which there are opposite trends in the magnitude of the vorticity of the tangential wind stress in different parts of the Benguela upwelling. In the southern part of the Benguela upwelling, both the Ekman’s upwelling and the vertical velocities due to the vorticity of the driving wind intensify, while in the northern part the corresponding trends have the opposite signs. This leads to a partial compensation of these two effects in the northern part of the Benguela upwelling. The reason for the change in the wind field structure is the displacement of the center of the Subtropical High to the south-east and the concomitant reversal of the near-surface wind vector in the coastal zone.

On the changes in the sea surface temperature in the Benguela upwelling region. Part 1: season cycle

Based on daily ocean surface temperature values for 1982–2017, near-surface wind for 1988–2017. and sea level anomalies for 1993–2017, obtained from satellite data, the seasonal variability of the hydrodynamic characteristics of the upper water layer in the vicinity of the Benguela upwelling is investigated. It is shown that the thermal upwelling index averaged over the entire area with lower temperature values does not give a correct idea on the seasonal course of the water lifting rate in the upwelling zone due to the significant horizontal advection of waters of upwelling origin. The seasonal variations of the vertical velocity of wind origin in the Benguela upwelling zone are characterized by the presence of two extremes from October to March, which is manifested in the predominance of the semiannual harmonic. At the same time, the thermal upwelling index in the zone of distribution of upwelling waters is subject to seasonal variability with an annual period.