Next Article in Journal
Effect of Hemp Seed Oil Addition on the Rheological Properties of Dough and Bread
Next Article in Special Issue
Estimating Full-Coverage PM2.5 Concentrations Based on Himawari-8 and NAQPMS Data over Sichuan-Chongqing
Previous Article in Journal
Setting Up and Operating Electric City Buses in Harsh Winter Conditions
Previous Article in Special Issue
Investigation of Petroleum Hydrocarbon Fingerprints of Water and Sediment Samples of the Nestos River Estuary in Northern Greece
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On the Geochemistry of Major and Trace Elements Distribution in Sediments and Soils of Zarafshon River Valley, Western Tajikistan

by
Djamshed A. Abdushukurov
1,
Daler Abdusamadzoda
1,2,
Octavian G. Duliu
2,3,4,*,
Inga Zinicovscaia
2,5 and
Pavel S. Nekhoroshkov
2
1
Institute of Water Problems, Hydropower and Ecology of Academy of Science, 14a Ainy Str., Dushanbe 734042, Tajikistan
2
Joint Institute for Nuclear Research, 6, Joliot Curie Str., 141980 Dubna, Russia
3
Department of Structure of Matter, Faculty of Physics, University of Bucharest, Earth and Atmospheric Physics and Astrophysics, 405, Atomistilor Str., 077125 Magurele, Romania
4
Geological Institute of Romania, 1, Caransebes Str., 012271 Bucharest, Romania
5
Department of Nuclear Physics, Horia Hulubei National Institute for R&D in Physics and Nuclear Engineering, 30 Reactorului Str., P.O. Box MG-6, 077125 Magurele, Romania
*
Author to whom correspondence should be addressed.
Appl. Sci. 2022, 12(6), 2763; https://doi.org/10.3390/app12062763
Submission received: 7 February 2022 / Revised: 26 February 2022 / Accepted: 3 March 2022 / Published: 8 March 2022
(This article belongs to the Special Issue Monitoring and Analysis of Environmental Pollution)

Abstract

:
To assess the geochemical features of sedimentary material of Zarafshon river, (Western Tajikistan) catchment basin, the mass fractions of 38 major and trace elements were determined by Instrumental Neutron Activation Analysis (INAA) in 2 × 116 paired samples of sediments and soils collected along the Zarafshon River and its main tributaries from the sources to Tajik—Uzbek border. At each collecting location, the distance between sediments and soils’ sampling was no greater than 10 m allowing the studying of the interrelation between sediments and soils. This evidenced a significant similarity between paired soils and sediments’ samples, including the potentially contaminating elements As, Sb and Hg, whose mass fractions in some places were significantly higher than for the Upper Continental Crust (UCC) and North American Shale Composite (NASC), suggesting a common provenience. At the same time, the distribution of major, as well as of incompatible trace elements, Sc, Zr, REE, Th, and U, in spite of geological diversity of the Zarafshon river catchment basin, suggest a possible felsic origin of investigated material.

1. Introduction

The Zarafshon River originates at the Zarafshon glacier, at an altitude of 2775 m, at the junction of the Turkestan and Zarafshon ranges, both of them belonging to the western sector of the Pamir-Alay system [1]. The initial part of the river, about 300 km long, lies in a narrow, deep valley. On the southern left bank, flowing between the Turkestan and Zarafshon ranges, it receives the Yagnob, Arthuch and Mogiyon rivers, as well as many small tributaries (Figure 1) [1,2,3]. After the city of Panjakent, the river follows a flat region, crosses the Tajikistan-Uzbekistan border, passes Samarkand and Bukhara, to finally disappear in the desert in the vicinity of the city of Karakul, without reaching the Amu Darya River [3,4]. Its catchment basin, with an area of 17,700 km2, covers diverse geological formations spread from the alpine zone, 3200–3500 m above sea level, to less than 150 m elevation westward towards the city of Karakul. In the Zarafshon mountain gorges, below 1500 m there are small areas of semi-savannahs as well as mountain forests. Here, regosols are prevalent [5], while in the western sector of Zarafshon Valley luvisols formed on loess deposits are common. Moreover, in mountain gorges, terraces of alluvial deposits and outflow cones of lateral tributaries serve as the basis for soil components (Figure 1a, Table A1).
The “Zarafshon” name, which in Farsi means “spreader of gold”, reflects the existence of ore fields along the river valley which, together with other outcrops of nonferrous ores makes it a mineral rich province requiring a detailed investigation.
Although the total population of Zarafshon valley does not exceed 400,000 inhabitants, mainly populating small villages, the existence of non-ferrous ore outcrops could represent the main source of down-stream contamination, previously analyzed in [6].
For this reason, the distribution of major and rock-forming elements as well as of incompatible lithophile elements such as Sc, Ti, V, Cr, Zr, lanthanides (LN) Th or more soluble U together with siderophile Co and Ni were chosen for this study. All of them play a significant role in geochemistry and pedology as tracers of the origin of natural material which, disintegrated during weathering and erosion, represents the main constituents of sediments and soils [7,8,9,10]. Besides them, there are other trace elements such as V, Cr, Co, Ni, Zn, As, Sb, or Hg, intensively used in diverse industrial processes, and, the presence of which, at levels significantly exceeding the natural ones, are related to an anthropogenic contamination which could be industrial [11,12] or due to urban as well as rural agglomerations [13]. According to [14], soil is a result of the combined action of weather, geographic relief, organisms and human activity on the parent mineral substrate, all of them continuously interacting. The actual soil is no older than Pleistocene, i.e., 2.59 Ma almost coincident with the onset of Northern hemisphere glaciation.
In turn, sediments represent a fragmented mineral substrate under the action of weathering and erosion and which is transported at long distances by the action of water, wind or ice and finally deposited. With respect to soil, sediments consist of unconsolidated clastic material with a reduced amount of organic matter. When consolidated by natural cementation processes, sediments turn into sedimentary rocks such as siltstone, sandstone, shale, etc.
Both sediments and soils preserve the fingerprint of the parental mineral substrate, which, under the action of external factors, evolve in a different manner. Collected from the same places and investigated in tandem, these sedimentary materials could bring a bonus of information concerning the local geochemistry as well as their interrelationship.
In this regard, in the period of spring-autumn, between 2015 and 2018, the Institute of Water Problems, Hydropower and Ecology of the Academy of Science of Tajikistan organized a field expedition to collect the most representative samples of sediments and soils from the Tajik sector of Zarafshon catching basin. For a better characterization of the soil–sediments interrelation, we collected from the same places distributed along the Zarafshon River and its main tributaries, including both sediments and soils samples (Figure 1a). In the case of presumed contaminated locations, mainly nonferrous mines and processing plants, samples were collected upstream and downstream at distances of no more than 10 to 20 m from them.
This permitted the collecting of 116 sediments, as well as an equal number of soils samples to be analyzed at the Neutron Activation Analysis Section of the Frank Laboratory of Neutron Physics (FLNP) of the Joint Institute for Nuclear Research (JINR), Dubna, Russian Federation. For a better precision and accuracy, we have used the Instrumental Neutron Activation Analysis (INAA), in both its thermal and epithermal variants [15,16]. It is worth mentioning that INAA was chosen due its capacity to determine the mass fractions of up to 40 elements, 38 in the case of this study, without any preliminary processing prone to induce systematic errors (Table 1) [16].

2. Hypothesis and Research Objectives

As presented before, the Zarafshon catchment basin covers a diversity of geological formations of which the origin dates back to Oligocene. This was when, as mentioned in [22], a Jurassic basin filled with continental sediments, under the pressure exerted by the Hindustan continental plate, and generated the Pamir-Alay system. Therefore, it is to be expected that the materials of which the Zarafshon valley is made up to contain a significant proportion of crustal matter. In time, due to weathering, this tends to be transported by the river and its tributaries as suspended load and deposited on river beds and river banks, as the mineral constituents of soils.
Accordingly, the main goals of this study consist of:
(i) to verify to what extent the geochemistry of depositional material, i.e., sediments and soils, can be related to crustal material, better approximated by UCC [17] and NASC [18], PAAS [8] in the case of soils, and SMWR [19], as well as AWSS [20] for sediments;
(ii) to evidence any interrelationship between the geochemistry of sediments and adjacent soils, as it is reflected by the distribution of major rock forming, as well as incompatible elements.
It should be noted that we have paid special attention to those elements whose presence is not related to any industrial activities such as mining, as the presence and distribution of the Presumably Contaminating Elements (PCE) on the Zarafshon catchment basin were previously analyzed in [6].
The results of our study performed under these circumstances will be further presented and discussed.

3. Geological Setting

Due to its location between two important mountain ranges—Turkestan and Zarafshon–the Zarafshon catchment’s basin covers a variety of lithologies of which clastic material represents the main component of sediments and bordering soils.
According to [22], the actual Turkestan and Zarafshon ranges evolved during Oligocene from a Jurassic basin filled with continental, shallow-water marine and lagoonal sediments of which crustal origin seems to influence its present geochemistry. This was the main reason for which we have chosen UCC [17] and NASC [18] as the most appropriate descriptor of investigated material.
At present, Djidjikrut, Fondarya, both tributaries of the Yagnob, as well as Mogiyon and its tributaries have the most important contribution to the mineral composition of Zarafshon sedimentary material [23,24].
The Djidjikrut River flowing in the Ayni district passes through clays, marls, limestones, dolomites and salts, significantly influencing the mineral composition of water and neighboring sediments and soils.
At the confluence with the Zarafshon (Yagnob river), is located the Djidjikrut ore field related to Zarafshon-Gissar mercury-antimony belt, and dating to Middle and Upper Paleozoic.
Antimonite, cinnabar and metacinnabarite are the main minerals of the Djidjikrut ore field, while pyrite and marcasite with an admixture of Co, Ni, Cu, As, Sb, Au, Te and Bi, realgar and orpiment, as well as hematite and sphalerite are common minerals too. For this reason, at the confluence of the Djidjikrut and Yagnob rivers, on the northern slope of the Gissar range, the Anzob mining and processing plant is operating [5] with the final release of Hg-Sb concentrate, an important source of downstream contamination (Figure A1).
The Yagnob River flows westward through the narrow mountain gorge along the southern edge of the Zarafshon mountains, being the first tributary of the Fondarya River. The total length of the Yagnob river is 114 km with a basin area of 1660 km2, while the mineral composition of sediments and soils of its upper reaches is still poorly investigated. Downstream, near the confluence with Fondarya, the Fon-Yagnob coal deposit occupies almost the entire left-bank site between the Djidjikrut and Pasrud tributaries [25]. Here, due to exothermic oxidation of the sub-surface coal layers, rich in pyrite and marcasite, a self-ignition process contaminated neighboring locations. Here, the paralava, a product of complete melting and (partial) recrystallization of local rocks generated by the high temperature attained during underground coal combustion, consists of almost melted sediments enclosing coal seams [26,27,28]. It is supposed that the underground coal fire may contribute to diagenesis of some REE minerals such as godovikovite, efremovite, maskanite, and so on [27,28,29,30].
The Fondarya River with a significant catchment basin of 3230 km2 is formed by the confluence of the Yagnob and Iskandardarya, which then crosses the Zarafshon ridge and flows into the Zarafshon river. In addition, the Fondarya is fed by the waters of the Sarvoda river which originates from the Alovdin lakes, as well as the lower and upper Kumarg tributaries, Chore and Pete (Figure 1a).
It should be noted that in the upper flows of the Iskandardarya river operates the mining and processing plant “Tolko-Gold”, which extracts and processes Au, Ag, Hg, and Sb ore as well as fluorite. The mining operation is also carried out at the gold deposit in the Kumargi Bolo gorge. The Chore River gorge is classified as gold-quartz deposits at the Turkestan-Chore zones of Central Tajikistan [31]. This deposit was discovered in 1960 and intensively exploited during 1974–1984 mainly for pyrite and arsenopyrite which represented more than 95% of all ore minerals together with small amounts of antimonites, chalcopyrite, sphalerite, pyrrhotite, and others [32].
The Pete tributary flows between the Pete-Takfon ore field hosting Takfon, Pete, Simich, and other deposits. In the geological structure of this ore field, Paleozoic formations are widely developed, which are represented by Silurian (shales, limestones, dolomites with interlayers of quartzites), Devonian (siliceous and micaceous shales with interlayers of cherts and limestones), Carboniferous (limestones with horizons of clay shales, sandstones) and Cong deposits. To the south of the ore field area, the complex of Paleozoic formations is overlapped by sub-platform sediments of the Mesozoic. The Igneous formations are represented by the Pete intrusions of granodiorites, dikes of granodiorite-porphyry, and lamprophyres. There are also tin and silver, lead, and zinc-containing ores [33].
The Mogiyon river (Figure 1a), with a length of 68 km and a basin area of 1100 km2, is one of the largest Zarafshon tributaries. From the source to the village of Mogiyon in a narrow valley, the river flows in a northern direction, cutting through the Zarafshon ridge with a deep transverse gorge. Further on, having received its largest tributary–the Shing River–it continues to flow in a narrow rapids channel with a width of up to 10 m. Nine kilometers above the city of Panjakent, leaving the gorge and passing through the village of Sujina, it then falls into the Zarafshon River (Figure 1). Between the Shing and Mogiyon interfluve is located the Taror ore field, which belongs to the Western Zarafshon-Gissar metallogenic zone [34].
However, along with gold-sulfide mineralization, silver-tin-base metal mineralization, localized in its upper and middle horizons, is also developed. The silver deposits in the Taror ore field belong to the Mirkhant deposit and its mineralization is predominantly represented by a silver-tin-polymetallic association spatially combined with sulfide, showing the closest correlation between Ag and Sb [27]. It is worth noting that in the Tarror gold ore deposit operates the Joint Venture “Zarafshon” whose main product consists of Ag-Au alloy [35].
In addition, along Zarafshon river and some of its main tributaries there is an active mining and metallurgical industry, based on the extraction and processing of non-ferrous metals. Among them, the Anzob mining and processing plant (AMPP), Taror gold ore plant (TGOP), Upper Kumarg and Kanchoch AGOK and TGOP specialized in the production of antimony-mercury concentrates and, respectively, gold. All of them could be regarded as a potential source of anthropogenic contamination which needs continuous monitoring extended to both the sediments and neighboring soils [35].
In this regard, the investigated Zarafshon valley together with some of its tributaries are the object of controversy between scientists and environmental experts representing governmental and non-governmental organizations. For this reason, it is not entirely correct to assess the environmental state of a particular region without paying a special attention to its geochemical features.

4. Materials and Methods

4.1. Sampling

An equal number of 116 sediments and 116 soils samples were collected in the investigated area (Figure 1a), using an AMS 12 Multi Stage sediments Sampler (AFMS Inc., American Falls, ID, USA) equipped with tubular glass sampler. Unconsolidated columns of sediment, 20 to 30 cm long, were handpicked from the riverbeds during a high-water period. As at altitudes greater than 2000 m, the soil cover is less developed, so that soil columns were no longer than 10 cm. For our study we retained the top soil (horizon A), as we were also interested to evidence recent contamination processes [6].
We collected soil samples always nearby the river banks. After collecting, any traces of vegetable matter such as fragments of grass and grass roots, leaves, wood pieces, as well as fragments of rock were carefully removed.
From each sampling point, three sub-samples of collected material, e.g., sediment or soil, covering an area of about 10 m2 were gathered and carefully mixed to form a composite sample. Any possible cross-contamination was avoided by washing the glass sampler before any sampling with a solution of diluted hydrochloric acid with deionized water and soaking them with disposable cotton serves.
To evidence any correlation between sediments’ and soils’ geochemistry, sediment and soil samples were collected in the same places, with two exceptions, and the distance between the sediments and the soils’ sampling location never exceeded 5–10 m (Table 1) and [36]. The same distances separated mines, adits or tunnels of the downstream collecting places.

4.2. Samples Preparation

Samples consisting of silty grayish sand without organic detritus were placed into the plastic container with cooled refrigerant and kept at 4 °C. In the laboratory, all samples were dried at 70 °C at constant weight and ground. Finally, the dried material was sieved using a 0.425 mm (42 mesh) sieve, homogenized and sent to the Frank Laboratory of Neutron Physics (FLNP) of the Joint Institute for Nuclear Research (JINR) for further INAA.
Here, 10 g of each sample were homogenized again for 15 min in a planetary ball mill (PULVERISETTE 6, Fritsch Laboratory Instruments GmbH, Idar-Oberstein, Germany) at 400 rpm. To increase the precision and accuracy of measurements, six aliquots of about 0.3 g were selected from each sample and independently analyzed.

4.3. INAA Measurements and Quality Control

Before irradiation at the IBR-2 reactor of the JINR FLNP in Dubna, each aliquot was packed in polyethylene bags for the determination of short-living isotopes (three aliquots) and in aluminum foil (three aliquots) for long-living ones. It is worth mentioning that the mass fractions of short-living isotopes were determined by thermal neutron activation while in the case of long-living isotopes we have used epithermal neutrons. More details concerning irradiation time and neutron flux densities in each irradiation channel can be found in [37].
After irradiation, the gamma spectrum of each sample was recorded using a HPGe detector with a resolution of 1.9 keV for the 60Co 1332 keV line, analyzed by means of a Genie 2000 Canberra software and processed by a proprietary software [37]. This permitted determination of the mass fraction of each of the 38 considered elements together with the associated CSU [21]. In all cases, the CSU was estimated by taking into account the statistical error, influence of measurement geometry, detector efficiency, the content, as well as the error of each element for all Certified Standard Material (CSM) utilized for calibration. To avoid any environmental contamination, all resulting radioactive materials were deposited in a permanent low activity repository.
The quality control was assured in the case of sediments and soils by a simultaneous analysis of Standard Reference Materials (SRM) provided by the National Institute for Standard and Technologies (NIST), i.e., 1633c Coal fly ash, 667 Estuarine sediment, 2710 Montana Soil and 1547 Peach leaves for short half-life time isotopes. In the case of long life-time isotopes, the best results were obtained using–2709 Trace elements in soil, 1632c–Trace elements in coal, 690CC–Calcareous soil, 2709a–San Joaquin soil, 1632c–Trace elements in coal and SRM-AGV2–Andesite.
Further, by using all of the above-mentioned SRM a Group of Standard Sample (GSS) [37,38] was realized which permitted the choice of the most appropriate SRM lines to minimize errors in determining the mass fractions of all investigated elements. It is worth mentioning that GSS is a proprietary software which, besides helping to determine the mass fractions of chosen elements, it allowed the calculation of the extent to which the mass fractions of SRM corresponded to certified ones [37,38]. Following this procedure, the accuracy quantified by means of standard deviation corresponding to each group of three aliquots was lower than 10%. In the case of SRM, the CSU varied between 3 and 15%, higher in the case of REE, but never greater than 20%. It is worth mentioning that for all 38 investigated elements, the detection limit was of 1 mg/kg and lower [37,38].
To illustrate the sensitivity and accuracy of INAA measurements following the CSS software, in Table A2 we have reproduced, as an example, the experimentally determined mass fractions of SRM 2709. Similar determinations were completed in the case of all utilized SRM.

5. Results and Discussion

A complete set of data concerning the mass fractions of all 38 elements in the investigated samples are provided in the Mendeley repository [36], while Table 1 reproduces only the mass fractions’ average values together with corresponding values of the UCC [17], NASC [18], PAAS [8], SMWR [19], as well as AWSS [20].
In interpreting compositional data of river sediments and adjacent soils, the geochemistry of local formations usually plays a determinant role. Additionally, due to a variable mobility, different mineral fractions are mixed and deposited along river beds depending on water speed, grain sizes and specific gravity. In this regard, Figure 2 illustrates the mass fractions of investigated elements in sediments normalized to UCC (Figure 2) [16], and in soils normalized to SMWR (Figure 2) [19]. The UCC [17] was chosen as an appropriate descriptor of the continental crust, while SMWR [19] drew on the lithogenesis theory [39], according to which, the mobilization of suspended and dissolved materials in drainage basins represent the initial stage of pedogenesis.
Here, a significant fraction of chemical elements including highly mobile alkaline and calc-alkaline ones were transported as suspended sediments and deposited along river beds and banks. Therefore, other descriptors such as NASC [17], AWSS [19] or PAAS [8] showed to be almost identical.
According to Table 1 and [36], the elemental composition of sediments and soils appears, excepting chalcophile Zn, As, Sb ad Hg, closer to UCC [17] and SMWR [19]. After a careful analysis of the experimental data presented in [36], mass fractions of potentially contaminating elements As, Sb and Hg showed increased values only for samples collected in the vicinity of mines and adits located on the Djidjikrut, Kanchoch, Chore and Mogiyon tributaries (Figure 1a, Figure 2 and Figure A1).
This finding is well evidenced in Figure 3a, where the As, Sb and Hg standard deviations appear disproportionately high due to the fact that in a very few places, such as mines or mine adits, their mass fractions exceeded the UCC [17] one by two to three orders of magnitude. Excepting these, the average mass fractions of investigated elements were shown to be relatively closer in both sediments and soils, as suggested by ANOVA analysis (Table 2) and by both biplots reproduced in Figure 3. In addition, a Spearman’s correlation coefficient of 0.83 sustains this observation (Figure 3b).
Figure 3a evidenced also that the ratio between the same elements in sediments and soil slightly varies around one, excepting the above mentioned PCE As, Sb and Hg and, to a lesser extent, Na of which the average mass fraction was slightly higher in sediments than in soils, as will be discussed in the next section.

5.1. Major Elements

The average values of the mass fractions of major, rock-forming elements excepting for MgO and CaO showed to be relatively closer to UCC [17], NASC, [18], PAAS [8], or other sedimentary systems such as SMWR [19], or Average World Suspended Sediments AWSS [20] (Table 1).
Silica represents the major constituent of all investigated samples, its mass fraction of 34.64 wt.% to 84 wt.% (62.8 ± 11.9 wt.% average) for sediments and of 31.53 wt.% to 75.12 wt.% (59.4 ± 8.8 wt.% average) for soils. On discriminating diagrams Na2O + K2O vs SiO2 (Figure 4a), SiO2-Al2O3-Na2O + K2O + CaO (Figure 4b) and K2O-Al2O3-Na2O + CaO (Figure 4c) the composition of investigated materials points towards a somewhat felsic to intermediate origin with a small contribution of mafic material as was the case in the Pete intrusion [5,24] (Table 1). Its average mass fraction was closer to the UCC [17] but with only a few exceptions in the sediments collected in the vicinity of the Anzob tunnel, Sarvoda 1, and Sarvoda 2 as well as soils from Chore 1 and 3, Djidjikrut 1, Obi Sara 1, Saritag 2, as well as Sarvoda 1 to Sarvoda 7 (Figure 1a and Figure A1). In all locations, the lower content of silica is due to an increased content of calcium oxide of which mass fractions varied between 26 and 47 wt.% (Table 1).
Aluminum oxide showed in sediments an average mass fraction 9.91 ± 3.2 wt.%, smaller than the soils’ one of 12.92 ± 3.7 wt.%, both of them lower than the UCC [17] and NASC [17] (Table 1). At the same time, the sodium oxide has an average mass fraction in sediments and especially in soils smaller than the corresponding UCC [17] (Table 1), suggesting a low content of plagioclase, and by consequence predominance of a sedimentary substrate (Figure 4a,b). This finding is illustrated by the discriminating diagram SiO2-Al2O3-Na2O + K2O + CaO (Figure 4b) where some points are shifted towards lower SiO2 content, which could be due to an increased content of CaO whose mass fraction in some places reached 28% [36]. It is worth remarking that CaO correlates positively only with MgO, in sediments as well as in soils (Table A3), suggesting the presence of dolomite CaMg(CO3)2. Besides these, the K-A-CN ternary diagram suggests a weak influence of weathering on both sediments and soils, expressed by a reduced fraction of total alkali–calc-alkali oxides Na2O + CaO of which normalized values appear smaller than 0.25 (Figure 4c).
The mass fractions of the other major elements oxides, i.e., TiO2, FeO, and MnO were closer to the UCC [17] values in conformity with the previous suppositions concerning a rather felsic origin of investigated material (Table 1).

5.2. Trace Elements

More data concerning the origin of investigated sedimentary materials can be inferred analyzing the distribution of trace elements and especially, as mentioned in the Introduction, of those considered incompatible and insoluble such as Sc, Zr, La, Ce, Nd, Sm, Eu, Gd, Tm, Yb, Th as well as U and, to a lower extent, Cr and Ni [6,7,40].
Notably, the felsic rocks are richer in Th and light REE, whereas mafic rocks contain more Sc, Cr, Co, Ni and an increased fraction of heavy REE [41]. In this regard, Sc of average mass fraction is less than 20 mg/kg in felsic rocks but increases to 20–40 mg/kg in mafic rocks, showing it to be one of the most appropriate discriminating elements [42]. In the case of investigated sediments and soils, Sc mass fractions were 9.56 ± 4 and 11.4 ± 4.6 mg/kg, respectively, confirming the previous inference concerning the felsic origin of Zarafshon valley sedimentary material.
The same finding was confirmed by discriminating Th/Co vs. La/Sc bi-plot [43] (Figure 5a) and ternary Sc-La-Th [7,44] diagram (Figure 5b). On both diagrams, the majority of Zarafshon experimental data regarding either sediments or soils are distributed around UCC [17] and NASC [18]. Moreover, on the Th/Co vs. La/Sc bi-plot [43], all Zarafshon data are closely grouped on the felsic sector of discriminating bi-plot reproduced in Figure 5a.
Two relatively closer bi-plots La/Th vs. Hf (Figure 6a) [45] and Th/Sc vs. Zr/Sc (Figure 6b) [8,46,47] provide a new type of information concerning not only the source material, but also the past processes that have been shaping the current position. Indeed, Hf has an increasing tendency to be enriched during the erosion of ancient (meta) sedimentary rocks. The mineral zircon, in turn, showed to be extremely resilient during recycling, so that the higher the Zr mass fraction, the larger the sedimentary material sorting and recycling [48]. In this regard, both bi-plots confirm previous inferences concerning the felsic nature of investigated sediments and soils, as well as reducing the cycle of erosion and recycling, i.e., a relatively young material. Only sediments and soils collected around the Fondarya tributary, the Kumargi Poyon, Yangob and Chore estuaries, showed a mixed felsic/basic source characterized by a higher La/Th ratio simultaneously with a reduced Hf mass fraction. On contrary, the increased mass fraction of Hf together with reduced values of La/Th ratio suggest the presence of old sediment sources in the Mastchohi kuhi, Shahriston and some regions of Panjekent [38].
Rare Earth Elements (REE) represent a group of 17 elements, including 15 LN, Sc and Y, all of them playing a significant role in provenance studies. Traditionally, LN are classified in light LN or light REE (LREE)—La to Sm and heavy REE (HREE)—from Eu to Lu. Since REE are not easily fractionated during sedimentation, they keep the “fingerprints” of the mother rocks making them helpful in establishing their origin [48,49]. Moreover, the Eu, which can exist in two valence states (Eu2+ and Eu3+) in its bivalent state due to a smaller ionic radius, can be substituted by Sr2+ ions in plagioclase feldspars. As this process takes place only in reducing conditions which can be found only in the superior mantle or lower crust, this leads to a distinctly negative anomaly in rocks coming from the upper crust [49,50]. Another distinctive characteristic of LN, useful in evidencing the rock origin, is related to the asymmetry of LREE and HREE in felsic rocks being enriched in LREE with respect to mafic ones, depleted in this category of REE.
The La/Th ratio could be useful in inferring the provenance of sedimentary material, by considering that for the UCC [17] and NASC [18] this ratio is around 2.88 while in the case of the average Mid-Ocean Ridge Basalt (MORB) [51] it is estimated to be between 7.2 and 15.6. In the case of Zarafshon samples, La/Th ratios for sediments and soils were 3.23 ± 1 and 3.03 ± 0.9, respectively (Figure 7a), both of them closer to the UCC ratio of 2.96 [17] and of 2.53 for NASC [18]. Further analysis of the distribution function of mass fraction of eight REE, i.e., La, Ce, Nd, Sm, Eu, Gd, Tb, Tm and Yb, are shown to be closer to the UCC [16] one. This is the same peculiarity that we have noticed in the case of the La/Yb ratio, as well as in the case of Eu negative anomaly (Figure 7b, Table 3).
A possible explanation of this small discrepancy could be related to an increased U content of which average mass fraction were between 3.04 ± 1.38 mg/kg in sediments and 3.8 ± 1.5 mg/kg in soils.
These results suggest the prevalence of the felsic material and sustains the hypotheses that sediments are the main material of the Zarafshon soils.

6. Conclusions

To determine at which extent the geochemistry of Turkestan and Zarafshon ranges is reflected by the adjacent sedimentary material, 116 sediments, and an equal number of soil samples were collected along the Tadjik sector of Zarafshon river and its main tributaries Arthuch, Djidjikrul, Iskadardarya, Mogiyon, Sarvoda, Shing, and Yagnob. All elemental analysis were performed by Instrumental Neutron Activation Analysis given its ability to analyze row material, without any previous treatment such as acid dissolution. This permitted determination of mass fractions of 38 major and trace elements with a sensitivity of 1 mg/kg, in both sediments and soils samples. As for each collecting point, the distance between sediments and soils sampling was no greater than 10 m, it was possible to investigate the interrelationship between sediments and soils geochemistry.
All data were compared and analyzed in correlation with similar data for the Upper Continental Crust (UCC) and North American Shale Composite (NASC). As a general conclusion, the distribution of mass fractions of all studied elements, except As, Sb and Hg showed a remarkable similarity in sediments and soils. In the case of the Presumed Contaminating Elements As, Sb and Hg, their mass fractions in few mines’ adits located around Djidjikrut or Kanchoch tributaries exceeded two to three order of magnitude their natural values. At the same time, mass fractions of both major and trace elements presented an increase variability characterized by coefficients of variation between 15 and 200% which could be attributed to the diversity of geological formation covering the Zarafshon catchment basin.
A detailed analysis of the major element oxides pointed towards a felsic origin of the collected material, with a small contribution of mafic one, as the alkali-silica discriminating diagrams SiO2-Al2O3-Na2O + K2O + CaO suggested. At the same time, the K2O-Al2O3-CaO + Na2O ternary diagram pointed towards a certain degree of weathering processes. Only in some places we have noticed an increased amount of magnesium and calcium oxides, most probable related to local deposits of limestone and dolomite. Relevant information was obtained by analyzing the distribution of low mobility and incompatible Sc, Zr, La, Ce, Nd, Sm, Eu, Gd, Tm, Yb as well as Th and U. Accordingly, Sc presented an average mass fraction of 10 to 11 mg/kg, characteristic of a felsic origin, a hypothesis suggested also by the reciprocal distribution of mass fractions of Sc, Co, La, Hf, or Th, in sediments as well as in soils. The la/Th vs. Hf and Th/Sc vs. Zr/Sc bi-plots suggest the absence of a significant sorting and recycling process of sedimentary material. The REE distribution was shown to be almost identical for sediments and soils and closer to UCC and NASC. The same pattern was noticed in the case of negative Eu anomaly, which was shown to be closer to UCC and NASC.
Given the surface of the investigated area which exceeds 17,700 km2, as well as the diversity of geological formation still insufficiently explored, the presented report should be seen as a baseline set of data for further detailed exploration.

Author Contributions

Conceptualization, D.A.A. and D.A.; methodology, D.A., O.G.D. and I.Z.; software, D.A. and O.G.D.; validation, O.G.D., D.A.A. and I.Z.; formal analysis, D.A. and P.S.N.; investigation, O.G.D. and D.A.; resources, D.A.A. and O.G.D.; data curation, D.A., I.Z. and O.G.D.; writing—original draft preparation, D.A.; writing—review and editing, D.A., O.G.D., I.Z. and D.A.A.; visualization, D.A. and O.G.D.; supervision, D.A.A. and O.G.D.; project administration, D.A.A.; funding acquisition, D.A.A., I.Z. and O.G.D.; All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All experimental data concerning sampling point coordinates, mass fraction of all 38 investigated elements are accessible at: Mendeley Reference Manager as: D.A., D.A.A., O.G.D., I.Z., P.S.N. Mendeley Data, https://doi.org/10.17632/vj24br7gnt.4 (accessed on 30 July 2021).

Acknowledgments

The research was performed within the framework of the Cooperation Agreement between the Institute of Water Problems, Hydropower and Ecology of Academy of Sciences of Tajikistan and Sector of Neutron Activation Analysis and Applied Research of Frank Laboratory of Neutron Physics of JINR. OGD wishes to acknowledge his contribution was provided within the Cooperation Protocol no. 4920-4-20/22 between the University of Bucharest and the Joint Institute of Nuclear Research, Dubna, Russian Federation. The authors are grateful to the staff of the SNAAPI of the Frank Laboratory of Neutron Physics of JINR for carrying out the neutron activation analysis of samples. We thank two anonymous reviewers for their comments and useful suggestions.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Figure A1. The map of Zarafshon catchment basin showing locations of the main sources of anthropogenic contamination [6].
Figure A1. The map of Zarafshon catchment basin showing locations of the main sources of anthropogenic contamination [6].
Applsci 12 02763 g0a1
Table A1. Sampling points coordinates, altitude as well as the local type of soil.
Table A1. Sampling points coordinates, altitude as well as the local type of soil.
Sampling PointLongitudeLatitudeAltitude (m)Soil
Djidjikrut 139.10914268.6883892726Regosol
Anzob tunnel adits39.10885368.6873142703Regosol
Djidjikrut 239.12658368.6829442574Regosol
Djidjikrut 339.13794468.6636392444Regosol
Djidjikrut 439.14680668.6517782398Regosol
Djidjikrut 539.16452868.6418612318Regosol
Djidjikrut 6 (Anzob ore-mining combine)39.17633368.6274172252Regosol
Djidjikrut 7 (Anzob ore-mining combine)39.19952868.6332781742Regosol
Yagnob 139.21733369.0211112351Regosol
Yagnob 239.21680668.9193332215Regosol
Yagnob 339.17275668.8744142151Regosol
Yagnob 4 (before Anzob confluence}39.16208368.8452222090Regosol
Estuary of Anzob river to the Yagnob39.15030468.8500062144Regosol
Yagnob 5 (after Anzob confluence)39.18966768.7587222070Regosol
Yagnob 639.19297268.6934721888Regosol
Yagnob 739.19997268.6401391739Regosol
Yagnob 839.19427868.5982781706Regosol
Yagnob 939.18747268.5434441657Regosol
Saritag 139.02852868.2712782477Regosol
Saritag 239.05641768.3464442194Regosol
Konchoch39.04216768.3327502361Luvisol
Kanchoch mine adit39.03952868.3378892457Luvisol
Iskandardarya 139.10136168.4035832015Regosol
Iskandardarya 239.13402868.4535281820Regosol
Iskandardarya 3 (estuary of river)39.18655668.5349441658Regosol
Fondarya 139.19638968.5360281644Regosol
Alovdin Lake 139.23764268.2542782780Regosol
Alovdin Lake 239.23747268.2610832780Regosol
Sarvoda 139.25019468.2688892620Regosol
Sarvoda 239.26294468.3109172488Cambisol
Sarvoda 339.24627868.3557502358Cambisol
Sarvoda 439.25211168.4050832166Gleysol
Sarvoda 539.22661168.4410281958Gleisol
Sarvoda 639.22319468.4846111808Gleisol
Sarvoda 7—estuary of river39.21988968.5276391660Regosol
Fondarya 239.22516768.5302501626Regosol
Pete 1 (right tributary to Fondarya)39.24863968.5617781895Gleisol
Pete 2 (right tributary to Fondarya)39.24936168.5498891710Gleisol
Fondarya 339.25880668.5440001589Regosol
Kumargi Bolo 139.27686168.5730562029Regosol
Kumargi Bolo 239.27544468.5475831690Regosol
Chore 1—upper of river39.31363968.4796942148Regosol
Chore 2—adit39.30894468.5061941863Regosol
Chore 3—adit39.30863968.5140561753Regosol
Chore 4—adit39.30822268.5286941577Regosol
Chore 5—estuary of river39.30608368.5329171580Regosol
Fondarya 439.30552868.5338611580Cambisol
Fondarya 539.34475068.5543331466Cambisol
Kumargi poyon 139.32594468.6283893629Gleysol
Kumargi poyon 239.36316768.6026941782Gleysol
Kumargi poyon 339.36396768.5610831434Gleysol
Fondarya 639.38052868.5495561389Cambisol
Zarafshon 139.45266770.4261002732Regosol
Zarafshon 2 (Dihisor)39.42619470.3075832612Regosol
Zarafshon 3 (Dihisor)39.44447270.1897502517Regosol
Zarafshon 4 (Paldorak)39.45758370.0761942423Regosol
Zarafshon 5 (Paldorak)39.42641769.9634442330Regosol
Zarafshon 6 (Sabak)39.44763969.8411942272Regosol
Zarafshon 7 (Madrushkat)39.45316769.7238612267Regosol
Rovut– tributary to Zarafshon39.44919469.6763062242Regosol
Zarafshon 8 (Madrushkat)39.44035869.6013242158Regosol
Madrushkat– tributary to Zarafshon39.44694469.5634722123Regosol
Tributary 2—to Zarafshon (Isizi poyon)39.43988969.5423612075Regosol
Zarafshon 9 Isizi poyon)39.43122269.4836391999Regosol
Tributary 3—to Zarafshon (Isizi poyon)39.41575069.4186941985Regosol
Zarafshon 10 (Isizi poyon)39.41330669.3768331927Regosol
Zarafshon 11 (Pastigov)39.41383369.2691941919Regosol
Kallakhona–tributary to Zarafshon39.41819469.2648611913Regosol
Oburdon (Pastigov)39.42561169.1672501875Regosol
Zarafshon 12 (Pastigov)39.41705669.1568331820Regosol
Zarafshon 13 (Shamtuch)39.40122269.0605001742Regosol
Zarafshon 14 (Veshab)39.40283368.9540281660Regosol
Zarafshon 15 (Shavadki poyon)39.39147268.8533061619Regosol
Zarafshon 16 (Rarz)39.38063968.7733611508Regosol
Zarafshon 17 (Fatmev)39.38586168.6738611462Regosol
Zarafshon 18 (Sangiston)39.38725068.5702501405Regosol
Shahriston 1, tributary right site of tunnel39.52080668.5533892728Regosol
Shahriston 2, tributary left site of tunnel39.51997268.5548892744Regosol
Shahriston 3–estuary39.44058368.5413891682Luvisol
Zarafshon 1939.41945568.5093111366Regosol
Zarafshon 2039.43919468.4415831306Regosol
Zarafshon 2139.43919468.4415831284Regosol
Zarafshon 2239.44080668.2732501226Regosol
Zarafshon 2339.44347268.0877501206Regosol
Artuch 1 (Alpine camp)39.27488968.1391672142Luvisol
Artuch 239.33527868.0877501573Luvisol
Artuch 339.39780668.0432501267Regosol
Zarafshon 2439.47166767.9918331113Gleysol
Zarafshon 2539.46647267.8871391072Gleysol
Obi Sara 139.37927867.8800281775Regosol
Obi Sara 239.41588967.8656391315Regosol
Obi Sara 339.46838967.8728331079Regosol
Zarafshon 2639.49730667.7850561035Gleysol
Lake 1 (Khazorchashma)39.11325067.8552502393Regosol
Lake 2 (Marguzor)39.14097267.8595282106Regosol
Lake 339.16608367.8389721888Regosol
Lake 439.19727867.8193331811Regosol
Lake 539.20356467.8106501718Regosol
Lake 639.21005667.8050281701Regosol
Lake 739.22030667.8008061610Regosol
Shing 139.24080667.7980281565Regosol
Shing 239.27366767.8027781382Regosol
Shing 3 industrial zone39.30691767.7837781282Regosol
Mogiyon 139.11261167.7336392383Gleysol
Mogiyon 239.14775067.7370002144Gleysol
Mogiyon 339.21519467.7119441722Gleysol
Mogiyon 439.27372267.6886391474Gleysol
Mogiyon 539.31558367.7701671250Luviysol
Mogiyon 6 after junction with Shing river39.33958367.7703061222Regosol
Mogiyon 739.37877867.7594721180Regosol
Mogiyon 839.41944467.7660831099Regosol
Mogiyon 939.45655667.7358891041Regosol
Mogiyon 1039.49508367.711861994Regosol
Zarafshon 2739.50288967.678861981Regosol
Zarafshon 2839.49672867.552247934Gleysol
Zarafshon 2939.52950067.433860883Gleysol
Table A2. The certified as well as experimentally determined mass fractions of the NIST 2709 CRM elements.
Table A2. The certified as well as experimentally determined mass fractions of the NIST 2709 CRM elements.
ElementCertified (w%)Measured (w%)ElementCertified (mg/kg)Measured (mg/kg)
Al7.377.11 ± 0.16Sb1.551.6 ± 0.06
Ca1.911.99 ± 0.09Ba979 951 ± 29
Fe3.363.45 ± 0.07Cr130121 ± 9.2
Mg1.461.38 ± 0.03Co12.813.6 ± 0.3
K2.112.15 ± 0.06Mn529546 ± 19
Si30.330.95 ± 0.41Sr239230 ± 8
Na1.221.26 ± 0.03V110121 ± 11
Ti0.340.3 ± 0.01Zr195201 ± 50
Table A3. The matrix of Spearman’s correlation coefficients of major elements (as oxides) in soils and sediments. All correlations significant at p < 0.05 are represented in red (positive correlation) or blue (negative correlation) color.
Table A3. The matrix of Spearman’s correlation coefficients of major elements (as oxides) in soils and sediments. All correlations significant at p < 0.05 are represented in red (positive correlation) or blue (negative correlation) color.
Sediments
SiO2TiO2Al2O3FeOMnO2MgOCaOK2O
TiO2−0.251
Al2O3−0.3020.840
FeO−0.3400.7580.819
MnO2−0.0180.6210.6390.510
MgO−0.8040.3810.3570.3820.013
CaO−0.752−0.195−0.184−0.153−0.1910.577
K2O−0.1550.6120.6270.7280.3920.268−0.260
Na2O0.2040.4190.4100.2760.313−0.066−0.4350.262
Soils
SiO2TiO2Al2O3FeOMnO2MgOCaOK2O
TiO2−0.081
Al2O3−0.1190.645
FeO−0.2910.4450.348
MnO2−0.0090.4420.4560.300
MgO−0.6710.062−0.0350.1530.022
CaO−0.668−0.436−0.527−0.123−0.3560.463
K2O−0.2130.5430.6260.4590.2550.108−0.276
Na2O−0.0810.5360.4030.4220.378−0.011−0.2680.431

References

  1. Prokhorov, A.M. Zarafshon (river in Central Asia). In Great Soviet Encyclopedia, 3rd ed.; Sovetskaya Encyclopaedia: Moscow, Russia, 1968–1978; Volume 30. (In Russian) [Google Scholar]
  2. Panjakent. Map of Panjakert. 1989. Available online: https://mapj42.narod.ru/map1/index016.html (accessed on 1 September 2021).
  3. Urgut. Map of Urgut. 1989. Available online: https://mapj42.narod.ru/map1/index015.html (accessed on 1 September 2021).
  4. Sadikov, A.S.; Mamarasulov, S.M.; Poslavskiy, V.V. Irrigation in Uzbekistan: Current State and Prospects of Irrigation Development in the Amu Darya River Basin; Fan: Tashkent, Uzbekistan, 1979; pp. 57–63. (In Russian) [Google Scholar]
  5. Geoportal. Mineral resources of the Republic of Tajikistan, Central Asia Geoportal. Available online: https://www.geoportal-tj.org/index.php/geology/deposits (accessed on 1 September 2021).
  6. Abdusamadzoda, D.; Abdushukurov, D.A.; Duliu, O.G.; Zinicovscaia, I. Assessment of the toxic metals pollution of soil and sediment in Zarafshon Valley, Northwest Tajikistan (Part II). Toxics 2020, 8, 113. [Google Scholar] [CrossRef]
  7. Bhatia, M.R.; Crook, K.A.W. Trace element characteristics of greywackes and tectonic setting discrimination of sedimentary basins. Contrib. Mineral. Petrol. 1986, 92, 181–193. [Google Scholar]
  8. Taylor, S.R.; McLennan, S.M. The Continental Crust: Its Composition and Evolution: An Examination of the Geochemical Record Preserved in Sedimentary Rocks; Blackwell Scientific Publications: Hoboken, NJ, USA, 1991; p. 312. ISBN 978-0632011483. [Google Scholar]
  9. McLennan, S.M.; Hemming, S.; McDaniel, D.K.; Hanson, G.N. Geochemical approaches to sedimentation, provenance and tectonics. In Processes Controlling the Composition of Clastic Sediments; Johnsson, M.J., Basu, A., Eds.; Geological Society of America: Washington, DC, USA, 1993; Volume 284, p. 285. [Google Scholar] [CrossRef] [Green Version]
  10. Armstrong-Altrin, J.S. Provenance of sands from Cazones, Acapulco, and Bahia Kino beaches, Mexico. Rev. Mex. Cienc. Geol. 2009, 26, 764–782. Available online: http://satori.geociencias.unam.mx/index.php/rmcg/article/view/346 (accessed on 30 July 2021).
  11. Bussan, D.; Harris, A.; Douvris, C. Monitoring of selected trace elements in sediments of heavily industrialized areas in Calcasieu Parish, Louisiana, United States by inductively coupled plasma-optical emission spectroscopy (ICP-OES). Microchem. J. 2019, 144, 51–55. [Google Scholar] [CrossRef]
  12. Belozertseva, I.A.; Milić, M.; Tošić, S.; Saljnikov, E. Environmental Pollution in the Vicinity of an Aluminium Smelter in Siberia. In Advances in Understanding Soil Degradation. Innovations in Landscape Research; Saljnikov, E., Mueller, L., Lavrishchev, A., Eulenstein, F., Eds.; Springer: Cham, Switzerland, 2022; pp. 379–402. [Google Scholar] [CrossRef]
  13. Badawy, W.M.; Sarhan, Y.; Duliu, O.G.; Kim, J.; El Samman, H.; Hussein, A.A.; Frontasyeva, M.; Shcheglov, A. Monitoring of air pollutants using plants and co-located soil—Egypt: Characteristics, pollution, and toxicity impact. Environ. Sci. Pollut. Res. 2021, 29, 21049–21066. [Google Scholar] [CrossRef] [PubMed]
  14. Ponge, J.-F. The soil as an ecosystem. Biol. Fertil. Soils 2015, 51, 645–648. [Google Scholar] [CrossRef]
  15. Zeisler, R.; Vajda, N.; Kennedy, G.; Lamaze, G.; Molnar, G.L. Activation analysis. In Handbook of Nuclear Chemistry, 2nd ed.; Vértes, A., Nagy, S., Klencsár, Z., Eds.; Springer: Berlin/Heidelberg, Germany, 2011; pp. 1555–1620. [Google Scholar] [CrossRef]
  16. Frontasyeva, M.V. Neutron activation analysis in the life sciences. Phys. Part. Nucl. 2011, 42, 332–378. [Google Scholar] [CrossRef]
  17. Rudnick, R.L.; Gao, S. Composition of the continental crust. In Treatise on Geochemistry; Holland, H.D., Turekian, K.K., Eds.; Elsevier-Pergamon: Amsterdam, The Netherlands, 2004; Volume 3, pp. 1–64. ISBN 978-0-08-098300-4. [Google Scholar]
  18. Gromet, L.P.; Haskin, L.Y.; Korotev, R.L.; Dymek, R.F. The North American shale composite: Its compilation, major and trace element characteristics. Geochim. Cosmochim. Acta 1984, 48, 2469–2482. [Google Scholar] [CrossRef]
  19. Savenko, V.S. Chemical composition of sediments load carried by rivers. Geochem. Int. 2007, 45, 816–824. [Google Scholar]
  20. Viers, J.; Dupre, B.; Gaillardet, G. Chemical composition of suspended sediments in World Rivers: New insights from a new database. Sci. Total Environ. 2009, 407, 853–868. [Google Scholar] [CrossRef]
  21. Farrance, I.; Frenkel, R. Uncertainty of Measurement: A Review of the Rules for Calculating Uncertainty Components through Functional Relationships. Clin. Biochem Rev. 2012, 33, 49–75. [Google Scholar] [PubMed]
  22. Leonov, M.G.; Rybin, A.K.; Batalev, V.Y.; Matyukov, V.E.; Shchelochkov, G.G. Tectonic Structure and Evolution of the Hissar–Alay Mountain Domain and the Pamirs. Geotectonics 2017, 51, 566–583. [Google Scholar] [CrossRef]
  23. Vinogradov, P.D. Information collection of All-Union Scientific Research Geological Institute. USSR 1961, 46, 192–207. (In Russian) [Google Scholar]
  24. Mamadvafoev, M.M.; Bakhtdavlatov, R.D.; Mirzoev, B.; Bakhronov, N. Farob ore field (Central Tajikistan): Geology and assessment of perspectives of industrial gold mineralisation. Proc. Acad. Sci. Rep. Tajikistan 2015, 58, 845–852. Available online: https://cyberleninka.ru/article/n/farobskoe-rudnoe-pole-tsentralnyy-tadzhikistan-geologiya-i-otsenka-perspektiv-promyshlennoy-zolotonosnosti (accessed on 30 July 2021). (In Russian).
  25. Vendrov, S.L.; Kalinin, G.P. Surface-water resources of the USSR: Their utilization and study. Sov. Geogr. 1960, 1, 35–49. [Google Scholar] [CrossRef]
  26. Sharygin, V.V.; Sokol, E.V.; Belakovskii, D.I. Fayalite-sekaninaite paralava from the Ravat coal fire (Central Tajikistan). Geol. Geophys. 2009, 50, 703–721. (In Russian) [Google Scholar] [CrossRef]
  27. Novikov, V.P.; Suprychev, V.V. The conditions of modern mineral formation in the process of underground combustion of coal at the Fan-Yagnob deposit. Mineral. Tajikistan 1986, 7, 91–104. [Google Scholar]
  28. Belakovski, D. Die Mineralien der brennenden Kohlefloeze von Ravat in Tadshikistan. Lapis 1990, 15, 21–26. (In German) [Google Scholar]
  29. Mirakov, M.A.; Fayziev, A.R.; Pautov, L.A. Native selenium in the products of underground fire of Fan-Yagnnob coal deposit (Central Tajikistan). Proc. Acad. Sci. Rep. Tajikistan 2017, 60, 456–460. Available online: https://cyberleninka.ru/article/n/samorodnyy-selen-v-produktah-podzemnogo-pozhara-fan-yagnobskogo-ugolnogo-mestorozhdeniya-tsentralnyy-tadzhikistan (accessed on 30 July 2021). (In Russian).
  30. Nasdala, L.; Pekov, I.V. Ravatite, C14H10, a new organic mineral species from Ravat, Tadzhikistan. Eur. J. Miner. 1993, 5, 699–705. [Google Scholar] [CrossRef]
  31. Madjidov, B.S. Organizational and Economic Substantiation of the Development of the Gold Mining Industry in the Region during the Formation of the Precious Metals Market: The Example of the Republic of Tajikistan. Ph.D. Thesis, Sergo Ordjonikidze Russian State University for Geological Prospecting, Moscow, Russia, 1985; p. 198. Available online: https://www.dslib.net/economika-xoziajstva/organizacionno-jekonomicheskoe-obosnovanie-razvitija-zolotodobyvajuwej-otrasli.html (accessed on 30 July 2021). (In Russian).
  32. Novojhilov, Y.I.; Gavrilov, A.M. Gold-Sulphide Deposits in Carbon-Terrigenous Layers; Publish house of Central Scientific Research Geological Prospecting Institute of Nonferrous and Precious Metals: Moscow, Russia, 1999; p. 174. (In Russian) [Google Scholar]
  33. Fayziev, F.A.; Fayziev, A.R.; Kayumarsi, M. On the silver mineralization of Central Tajikistan. Bull. Ural. State Min. Univ. 2017, 4, 18–22. (In Russian) [Google Scholar] [CrossRef] [Green Version]
  34. Blokhina, N.A. Mineralogy, Geochemistry and Formation Conditions of Gold Sulfide Deposits in the Formation of Magnesian Skarns (Central Tajikistan); Danish Press: Dushanbe, Tajikistan, 1984; p. 277. (In Russian) [Google Scholar]
  35. Solikhov, M.M.; Aksenov, A.V.; Vasiliev, A.A.; Karimov, M.A.; Rahmanov, O.B. Search for optimal processing technology for Taror deposit refractory gold ores (Republic of Tajikistan). Proc. Irkutsk. State Tech. Univ. 2018, 22, 182–192. (In Russian) [Google Scholar] [CrossRef]
  36. Abdusamadzoda, D.; Abdushukurov, D.A.; Duliu, O.G.; Zinicovscaia, I.; Nekhoroshkov, P.S. Major and Trace Elements in Sediments and Soils of the Zerafshon Area, North-West Tajikistan. Available online: https://data.mendeley.com/drafts/vj24br7gnt (accessed on 30 July 2021).
  37. Pavlov, S.S.; Dmitriev, A.Y.; Frontasyeva, M.V. Automation system for neutron activation analysis at the reactor IBR-2, Frank Laboratory of Neutron Physics, Joint Institute for Nuclear Research, Dubna, Russia. J. Radioanal. Nucl. Chem. 2016, 309, 27–38. [Google Scholar] [CrossRef] [Green Version]
  38. Zinicovscaia, I.; Duliu, O.G.; Culicov, O.A.; Frontasyeva, M.V.; Rodica, S. Major and trace elements distribution in Moldavian soils. Rom. Rep. Phys. 2018, 70, 701. [Google Scholar]
  39. Strakhov, N.M. Principles of Lithogenesis, 3rd ed.; Tomkeieff, S.I., Hemingway, J.E., Eds.; Springer: Boston, MA, USA, 1995; p. 245. [Google Scholar] [CrossRef]
  40. Wang, W.; Chen, F.; Hu, R.; Chu, Y.; Yang, Y.Z. Provenance and tectonic setting of Neoproterozoic sedimentary sequences in the South China Block: Evidence from detrital zircon ages and Hf-eNd isotopes. Int. J. Earth Sci. 2012, 101, 1723–1744. [Google Scholar] [CrossRef]
  41. Marshak, S. Essentials of Geology, 6th ed.; W.W. Norton & Company: New York, NY, USA, 2019; ISBN 978-0393667523. [Google Scholar]
  42. Norman, J.C.; Haskin, L.A. The geochemistry of Sc: A comparison to the rare earths and Fe. Geochim. Cosmochim. Acta 1968, 32, 93–108. [Google Scholar] [CrossRef]
  43. Armstrong, A.J.S.; Nagarajan, R.; Madhavaraju, J.; Nagarajan, R.; Madhavaraju, J.; Rosalez, H.L.; Lee, Y.I.; Balaram, V.; Cruz, M.A.; Avila, G.R. Geochemistry of the Jurassic and Upper Cretaceous shales from the Molango Region, Hidalgo, eastern Mexico: Implications for source-area weathering, provenance, and tectonic setting. Comp. Ren. Geosci. 2013, 45, 185–202. [Google Scholar] [CrossRef] [Green Version]
  44. Cullers, R.L. Implications of elemental concentrations for provenance, redox conditions, and metamorphic studies of shales and limestones near Pueblo, CO, USA. Chem. Geol. 2002, 191, 305–327. [Google Scholar] [CrossRef]
  45. Floyd, P.A.; Leveridge, B.E. Tectonic environment of the Devonian Gramscatho basin, south Cornwall: Frame-work mode and geochemical evidence from turbiditic sandstones. J. Geol. Soc. 1987, 144, 531. [Google Scholar] [CrossRef]
  46. Lee, Y.I. Geochemistry of shales of the Upper Cretaceous Hayang Group, SE Korea: Implications for provenance and source weathering at an active continental margin. Sedim. Geol. 2009, 215, 1–12. [Google Scholar] [CrossRef]
  47. Nagarajan, R.; Armstrong-Altrin, J.S.; Kessler, F.L.; Hidalgo-Moral, E.L.; Dodge-Wan, D.; Nur, I.T. Provenance and tectonic setting of Miocene siliciclastic sediments, Sibuti Formation, northwestern Borneo. Arab. J. Geosci. 2015, 8, 8549–8565. [Google Scholar] [CrossRef]
  48. McLennan, S.M.; Hemming, S.; McDaniel, D.K.; Hanson, G.N. Geochemical Approaches to Sedimentation, Provenance, and Tectonics; Special paper; Geological Society of America: Boulder, CO, USA, 1993; Volume 284, p. 20. [Google Scholar]
  49. McLennan, S.M. Rare Earth Elements in sedimentary rocks: Influence of provenance and sedimentary processes. In Geochemistry and Mineralogy of Rare Earth Elements; Lipin, B.R., McKay, G.A., Eds.; Mineralogical Society of America, de Gruyter: Berlin, Germany, 2018; Volume 21, pp. 169–200. [Google Scholar] [CrossRef]
  50. McLennan, S.M. Lanthanide Rare Earths. In Encyclopedia of Geochemistry. Encyclopedia of Earth Sciences Series; White, W., Ed.; Springer: Cham, Switzerland, 2018. [Google Scholar] [CrossRef]
  51. Allison, G.; Dalton, C.A.; Langmuir, C.H.; Su, Y. The mean composition of ocean ridge basalts. Geoch. Geoph. Geosys. 2013, 14, 489–518. [Google Scholar] [CrossRef] [Green Version]
  52. Arevako, R.; McDonough, W.F. Chemical variations and regional diversity observed in MORB. Chem. Geol. 2010, 271, 70–85. [Google Scholar] [CrossRef]
  53. McDonough, W.F.; Sun, S.-s. The composition of the Earth. Chem. Geol. 1995, 120, 223–253. [Google Scholar] [CrossRef]
Figure 1. The map of Zarafshon Valley (a); showing its position within Tajikistan territory (b).
Figure 1. The map of Zarafshon Valley (a); showing its position within Tajikistan territory (b).
Applsci 12 02763 g001
Figure 2. Box-plots illustrating the distribution of investigated elements in Zarafshon sediments (lower); and soils (upper). All mass fractions of soils normalized to UCC [17] and sediments to SMWR [19]. The values corresponding to As and Hg are divided by 10 and by 100 in the case of Sb.
Figure 2. Box-plots illustrating the distribution of investigated elements in Zarafshon sediments (lower); and soils (upper). All mass fractions of soils normalized to UCC [17] and sediments to SMWR [19]. The values corresponding to As and Hg are divided by 10 and by 100 in the case of Sb.
Applsci 12 02763 g002
Figure 3. The sediment to soil ratio of mass fractions of investigated elements (red line) as well as the corresponding standard deviation (gray background) (a); and a biplot illustrating the reciprocal distribution of investigated elements in sediments and soils excepting the Zn, As, Sb and Hg, i.e., elements whose mass fractions in some places exceeds the UCC [17] ones by two to three orders of magnitude (b) (Table 1) [36]. Spearman’s correlation coefficient of 0.837 at p < 0.01.
Figure 3. The sediment to soil ratio of mass fractions of investigated elements (red line) as well as the corresponding standard deviation (gray background) (a); and a biplot illustrating the reciprocal distribution of investigated elements in sediments and soils excepting the Zn, As, Sb and Hg, i.e., elements whose mass fractions in some places exceeds the UCC [17] ones by two to three orders of magnitude (b) (Table 1) [36]. Spearman’s correlation coefficient of 0.837 at p < 0.01.
Applsci 12 02763 g003
Figure 4. Discriminating total alkali metal oxides vs. SiO2 biplot (a); ternary SiO2-Al2O3-Na2O + K2O + CaO (b); and K-A-CN (c) diagrams illustrating the preponderance of felsic material in Zarafshon sediments and soils as well as a small fraction of weathered material.
Figure 4. Discriminating total alkali metal oxides vs. SiO2 biplot (a); ternary SiO2-Al2O3-Na2O + K2O + CaO (b); and K-A-CN (c) diagrams illustrating the preponderance of felsic material in Zarafshon sediments and soils as well as a small fraction of weathered material.
Applsci 12 02763 g004
Figure 5. Discriminating bi-plots La/Th vs Hf (a) and Th/Sc vs. Zr/Sc; (b) confirming the felsic origin of Zarafshon sediments and soils as well as the absence of a significant recycling of sedimentary material.
Figure 5. Discriminating bi-plots La/Th vs Hf (a) and Th/Sc vs. Zr/Sc; (b) confirming the felsic origin of Zarafshon sediments and soils as well as the absence of a significant recycling of sedimentary material.
Applsci 12 02763 g005
Figure 6. Discriminating bi-plot Th/Co vs. La/Sc (a); and ternary Sc-La-Th diagram (b) confirming the felsic origin of investigated sediments and soils.
Figure 6. Discriminating bi-plot Th/Co vs. La/Sc (a); and ternary Sc-La-Th diagram (b) confirming the felsic origin of investigated sediments and soils.
Applsci 12 02763 g006
Figure 7. The La vs. Th bi-plot (a) proving that both Zarafshon sediments and soils are significantly closer to UCC [17] and NASC [18] than to MORB [51,52] as well as the eight REE average values of Zarafshon sediments and soils normalized to chondrite [53] together with the corresponding UCC data (b) [16].
Figure 7. The La vs. Th bi-plot (a) proving that both Zarafshon sediments and soils are significantly closer to UCC [17] and NASC [18] than to MORB [51,52] as well as the eight REE average values of Zarafshon sediments and soils normalized to chondrite [53] together with the corresponding UCC data (b) [16].
Applsci 12 02763 g007
Table 1. The mass fraction of investigated elements in sediments and soils as well as the main literature references: UCC [17]; North American Shale Composite NASC [18]; Post Archaean Australian Shale (PAAS) [8]; Suspended Matter of the World’s Rivers (SMWR) [19]; as well as Average World Suspended Sediments (AWSS) [20]. Oxides’ mass fraction expressed in wt.%, in mg/kg for other elements, excepting Hg of which the mass fraction is expressed in μg/kg. CSU represents the Combined Standard Uncertainty [21].
Table 1. The mass fraction of investigated elements in sediments and soils as well as the main literature references: UCC [17]; North American Shale Composite NASC [18]; Post Archaean Australian Shale (PAAS) [8]; Suspended Matter of the World’s Rivers (SMWR) [19]; as well as Average World Suspended Sediments (AWSS) [20]. Oxides’ mass fraction expressed in wt.%, in mg/kg for other elements, excepting Hg of which the mass fraction is expressed in μg/kg. CSU represents the Combined Standard Uncertainty [21].
ElementSedimentCSUSoilsCSUUCCNASCPAASSMWRAWSS
SiO262.7512.6559.399.766.6264.4862.854.7654.34
TiO20.620.220.831.060.640.7210.650.73
Al2O39.93.5212.924.0715.416.92.216.3116.48
FeO3.792.55.782.015.045.718.96.477.47
MnO20.070.070.090.030.10.060.110.150.22
MgO6.323.425.82.722.482.85---2.392.09
CaO10.5910.719.469.063.593.561.33.643.62
K2O2.40.992.830.792.83.993.72.592.04
Na2O1.40.790.90.413.271.151.21.110.96
Sc9.641151415161418
V10156254150097------120129
Cr714090579212411085130
Co11.56.414.57.117.3------1922.5
Ni463265314758555074
Zn102571195467-------130208
As63254846004.828.4---1436.3
Br2.64.31.91.61.60.9---921.5
Rb9241103398412616077.079
Sr14363177100320142200150.0187
Zr220133210101193200210150.0160
Sb59560382700.4------1.42.2
Cs4.55707.53.64.95.2155.26.3
Ba724450784500624636650500522
La28143212313138.23237
Ce53266125636779.66873.6
Nd27153217272733.92932.2
Sm52.56.33.44.75.65.65.86.1
Eu0.90.61.10.811.21.11.41.3
Gd3.93.32.92.44---4.75.65.3
Tb0.70.310.40.70.90.80.80.8
Tm0.30.30.50.50.3---0.40.40.4
Yb1.91.22.21.523.12.82.52.1
Hf5.535.82.65.36.3---4.44
Ta0.80.410.50.91.1214.60.91.3
W2.55.72.82.71.92.1---1.42
Hg2.514.71.14.50.1------0.1---
Th8.94.310.93.510.512.33.21012.1
U31.53.91.72.72.7---2.43.3
Table 2. The results of two-sample ANOVA tests illustrating the probability the major and trace elements being closer in soils and sediments. The results are based on the average mass fractions values reproduced in Table 1.
Table 2. The results of two-sample ANOVA tests illustrating the probability the major and trace elements being closer in soils and sediments. The results are based on the average mass fractions values reproduced in Table 1.
ANOVA Two-Sample Test
ElementsT-Test
Equal Mean
F-Test
Equal Variance
Mann–Whitney
Equal Median
Anderson–Darling
Equal Distribution
Major elements0.990.871.001.00
Trace elements0.770.590.730.97
Major elements: SiO2, TiO2, Al2O3, FeO, MnO2, MgO, CaO, K2O, Na2O
Trace elements: Sc, V, Cr, Co, Ni, Zn, As, Br, Rb, Sr, Zr, Sb, Cs, Ba, La, Ce, Nd, Sm, Eu, Gd, Tb, Tm, Yb, Hf, Ta, W, Hg, Th, U
Table 3. The La/Sm, Gg/Yb, Eu* as well as the Th/U values for sediments and soils. For comparison, the same ratios corresponding to UCC [17] and NASC [18] are reproduced too.
Table 3. The La/Sm, Gg/Yb, Eu* as well as the Th/U values for sediments and soils. For comparison, the same ratios corresponding to UCC [17] and NASC [18] are reproduced too.
La/SmGd/YbEu*Th/U
Sediment4.91 ± 1.81.92 ± 1.30.69 ± 0.273.23 ± 1
Soils5.81 ± 2.41.49 ± 1.1 0.8 ± 0.323.03 ± 0.9
UCC6.61.90.722.96
NASC5.56------2.53
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Abdushukurov, D.A.; Abdusamadzoda, D.; Duliu, O.G.; Zinicovscaia, I.; Nekhoroshkov, P.S. On the Geochemistry of Major and Trace Elements Distribution in Sediments and Soils of Zarafshon River Valley, Western Tajikistan. Appl. Sci. 2022, 12, 2763. https://doi.org/10.3390/app12062763

AMA Style

Abdushukurov DA, Abdusamadzoda D, Duliu OG, Zinicovscaia I, Nekhoroshkov PS. On the Geochemistry of Major and Trace Elements Distribution in Sediments and Soils of Zarafshon River Valley, Western Tajikistan. Applied Sciences. 2022; 12(6):2763. https://doi.org/10.3390/app12062763

Chicago/Turabian Style

Abdushukurov, Djamshed A., Daler Abdusamadzoda, Octavian G. Duliu, Inga Zinicovscaia, and Pavel S. Nekhoroshkov. 2022. "On the Geochemistry of Major and Trace Elements Distribution in Sediments and Soils of Zarafshon River Valley, Western Tajikistan" Applied Sciences 12, no. 6: 2763. https://doi.org/10.3390/app12062763

APA Style

Abdushukurov, D. A., Abdusamadzoda, D., Duliu, O. G., Zinicovscaia, I., & Nekhoroshkov, P. S. (2022). On the Geochemistry of Major and Trace Elements Distribution in Sediments and Soils of Zarafshon River Valley, Western Tajikistan. Applied Sciences, 12(6), 2763. https://doi.org/10.3390/app12062763

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop