<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article  PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd"><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="3.0" xml:lang="en" article-type="research article"><front><journal-meta><journal-id journal-id-type="publisher-id">JEP</journal-id><journal-title-group><journal-title>Journal of Environmental Protection</journal-title></journal-title-group><issn pub-type="epub">2152-2197</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/jep.2024.155031</article-id><article-id pub-id-type="publisher-id">JEP-133512</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Earth&amp;Environmental Sciences</subject></subj-group></article-categories><title-group><article-title>
 
 
  Soil Profile Concentration Distributions of Antimony and Bismuth across Southeastern Missouri (USA)
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Michael</surname><given-names>Aide</given-names></name><xref ref-type="aff" rid="aff1"><sub>1</sub></xref></contrib></contrib-group><aff id="aff1"><label>1</label><addr-line>Department of Agriculture, Southeast Missouri State University, Cape Girardeau, MO, USA</addr-line></aff><pub-date pub-type="epub"><day>24</day><month>05</month><year>2024</year></pub-date><volume>15</volume><issue>05</issue><fpage>537</fpage><lpage>551</lpage><history><date date-type="received"><day>15,</day>	<month>April</month>	<year>2024</year></date><date date-type="rev-recd"><day>27,</day>	<month>May</month>	<year>2024</year>	</date><date date-type="accepted"><day>30,</day>	<month>May</month>	<year>2024</year></date></history><permissions><copyright-statement>&#169; Copyright  2014 by authors and Scientific Research Publishing Inc. </copyright-statement><copyright-year>2014</copyright-year><license><license-p>This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/</license-p></license></permissions><abstract><p>
 
 
  Recent research has suggested that increased industrial and technological utilization of antimony and bismuth necessitates greater research to determine the soil and water chemistry and the environmental risks associated with these elements. The near-total soil profile concentrations of antimony and bismuth were determined for key soil series across southeastern Missouri. The antimony concentrations ranged from 0.65 to 0.08 mg kg&lt;sup&gt;&amp;#8722;&lt;/sup&gt;&lt;sup&gt;1&lt;/sup&gt;, whereas the bismuth soil profile concentrations ranged from 0.92 to 0.03 mg kg&lt;sup&gt;&amp;#8722;&lt;/sup&gt;&lt;sup&gt;1&lt;/sup&gt;. Most pedons showed antimony concentrations ranging from 20 to 30 mg kg&lt;sup&gt;&amp;#8722;&lt;/sup&gt;&lt;sup&gt;1&lt;/sup&gt;, whereas bismuth concentrations were commonly 10 to 20 mg kg&lt;sup&gt;&amp;#8722;&lt;/sup&gt;&lt;sup&gt;1&lt;/sup&gt;. For soils having argillic horizons, antimony and bismuth concentrations were greater for the illuvial horizons than the eluvial horizons, whereas Entisols, Inceptisols, and one Vertisol showed rather uniform antimony and bismuth concentrations, features paralleling the soil texture distribution. Both antimony and bismuth showed significant correlations with iron.
 
</p></abstract><kwd-group><kwd>Trace Elements</kwd><kwd> Antimony</kwd><kwd> Bismuth</kwd><kwd> Soils</kwd><kwd> Group 15</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>Group 15 of the Periodic <xref ref-type="table" rid="table">Table </xref>consists of the elements nitrogen (N), phosphorus (P), arsenic (As), antimony (Sb) and bismuth (Bi). In this manuscript, the emphasis will focus on antimony and bismuth. Previously, Aide et al. [<xref ref-type="bibr" rid="scirp.133512-ref1">1</xref>] published soil profile distributions involving arsenic across southeastern Missouri.</p><sec id="s1_1"><title>1.1. Introduction to Antimony</title><p>The ground state electronic configuration for Sb is [Kr] 4d<sup>10</sup>5s<sup>2</sup>p<sup>3</sup> and the covalent radius is 0.141 nm [<xref ref-type="bibr" rid="scirp.133512-ref2">2</xref>] . Filella et al. [<xref ref-type="bibr" rid="scirp.133512-ref3">3</xref>] noted that the primary Sb valence states are Sb<sup>+5</sup>, Sb<sup>+3</sup>, and Sb<sup>−</sup><sup>3</sup>, with Sb<sup>5+</sup> as the more commonly occurring oxidation state. Stability constants and other selected thermodynamic data are limited. Antimony compounds are considered priority pollutants in the United States and European Union and elevated antimony soil and sediment concentrations are primarily related to anthropogenic sources or high arsenic bearing sulfide ores. Stable anionic thiocomplexes include SbS<sub>2</sub><sup>2</sup><sup>−</sup> and SbS<sub>4</sub><sup>3</sup><sup>−</sup> [<xref ref-type="bibr" rid="scirp.133512-ref3">3</xref>] . Antimony has a wide range of uses including the manufacture of semiconductors, diodes, flameproof retardants, lead hardeners, batteries, small arms, tracer bullets, automobile brake linings, and pigments [<xref ref-type="bibr" rid="scirp.133512-ref3">3</xref>] .</p><p>In a detailed survey of the literature, Kabata-Pendias [<xref ref-type="bibr" rid="scirp.133512-ref4">4</xref>] reported the crustal average Sb concentrations are approximately 0.2 mg kg<sup>−</sup><sup>1</sup>, with argillaceous sediments showing slightly larger average concentrations of up to 4 mg kg<sup>−</sup><sup>1</sup>. Hou et al. [<xref ref-type="bibr" rid="scirp.133512-ref5">5</xref>] reported that the average Sb concentration ranges in Japanese soils were 0.83 &#177; 0.32 mg kg<sup>−</sup><sup>1</sup>. Stibnite (Sb<sub>2</sub>S<sub>3</sub>) is the most important antimony ore, followed by valentinite (Sb<sub>2</sub>O<sub>3</sub>). Other sulfide minerals include pyrargyrite, zinkenite, jamesonite, and boulangerite.</p><p>Filella et al. [<xref ref-type="bibr" rid="scirp.133512-ref3">3</xref>] reported an ordering for antimony species toxicity, which is antimonites (Sb<sup>3+</sup>) &gt; antimonates (Sb<sup>5+</sup>) &gt; organoantimonials. Kolesnikov et al. [<xref ref-type="bibr" rid="scirp.133512-ref6">6</xref>] evaluated the ecotoxicity of 23 metals, metalloids, and nonmetals in a Haplic Chernozem and proposed three hazard classes with Sb in class II (intermediate hazard class). Antimony was determined to be bound to relatively immobile Fe and Al oxihydroxides and to a lesser degree as antimony-organic substances [<xref ref-type="bibr" rid="scirp.133512-ref7">7</xref>] .</p><p>Antimony (III) and (V) ions undergo hydrolysis readily in aqueous solutions. In oxic water and soil conditions, the Sb(V) species are the major species [<xref ref-type="bibr" rid="scirp.133512-ref8">8</xref>] . Antimonic acid has been represented as H[Sb(OH)<sub>6</sub>], Sb(OH)<sub>5</sub>, or HSbO<sub>3</sub> and frequently forms polymers as the pH increases. Antimonous acid [Sb(III)] is frequently represented as SbO<sup>+</sup>, Sb(OH)<sub>2</sub><sup>+</sup> in acidic media and as Sb(OH)<sub>4</sub><sup>−</sup> or hydrated SbO<sub>2</sub><sup>−</sup> in basic media. Antimony(III) chloride will form successive chlorocomplexes [<xref ref-type="bibr" rid="scirp.133512-ref8">8</xref>] . Under reducing conditions and in the presence of sulfur, the mineral stibnite (Sb<sub>2</sub>S<sub>3</sub>) will crystalize in acid soils, whereas in alkaline pH levels the mineral SbS<sub>2</sub> replaces stibnite.</p><p>Sb(III) prefers sulfur as a ligand; however, Sb<sup>3+</sup> will form stable complexes with ligands such as citric, lactic, mandelic, and tartaric acids Antimony(V) forms complexes with polyhydric alcohols, polyhydric phenols, and citric, malic, and lactic acids. Antimonic acid (Sb<sub>2</sub>O<sub>5</sub>.nH<sub>2</sub>O or H<sub>3</sub>O<sub>4</sub>Sb) forms complexes with low molecular weight acids that have oxygen-containing functional groups, including humic compounds. Pilarski et al. [<xref ref-type="bibr" rid="scirp.133512-ref9">9</xref>] investigated the adsorption of Sb(III) and Sb(V) with humic acid, demonstrating adsorption was effectively described using Langmuir isotherms. Bagherifam et al. [<xref ref-type="bibr" rid="scirp.133512-ref10">10</xref>] investigated high organic matter soils incubated with antimony and documented that humic substances were responsible for 63% of the extractable Sb. In the humic acid fraction, antimony was largely associated with the low molecular weight fulvic acids and correlated with total organic carbon and nitrogen contents [<xref ref-type="bibr" rid="scirp.133512-ref9">9</xref>] .</p><p>Thanabalasingam et al. [<xref ref-type="bibr" rid="scirp.133512-ref11">11</xref>] investigated antimony adsorption onto hydrous oxides of Mn, Fe, and Al. Capacity values decreased along the sequence MnOOH &gt; Al(OH)<sub>3</sub> &gt; FeOOH and adsorption on each substrate decreased gradually at pH values greater than pH 6. Factors influencing Sb sorption included substrate surface charge, chemical form of Sb and surface interactions. In oxic soil conditions antimony is present as Sb(OH)<sub>6</sub> and in anoxic soil conditions Sb(OH).</p><p>Tang et al. [<xref ref-type="bibr" rid="scirp.133512-ref12">12</xref>] and Vidya et al. [<xref ref-type="bibr" rid="scirp.133512-ref13">13</xref>] noted that mine waste leaching, the weathering of sulfide ores and shooting ranges using Sb adulterated lead-based ammunition are major Sb pathways for impacting soils and aquatic environments. Tang et al. [<xref ref-type="bibr" rid="scirp.133512-ref12">12</xref>] also documented that Sb toxicity in plants reduces i) root and shoot growth, ii) seed germination, and iii) yield potential. Furthermore, antimony induces chlorosis, reduces photosynthetic efficiency, membrane stability, and nutrient uptake, and increases reactive oxygen species. Vidya et al. [<xref ref-type="bibr" rid="scirp.133512-ref13">13</xref>] reported that the majority of plant uptake of Sb is confined to root tissues; however, some of the metalloid is translocated to the shoot. Inhibition of photosynthesis, modified root and leaf patterns, activation of antioxidant systems, and plant membrane disruption are some of the deleterious effects of antimony on plant growth and development. Bowen [<xref ref-type="bibr" rid="scirp.133512-ref14">14</xref>] documented that terrestrial plant uptake of antimony may result in a plant tissue concentration of approximately 60 &#181;g kg<sup>−</sup><sup>1</sup>.</p></sec><sec id="s1_2"><title>1.2. Introduction to Bismuth</title><p>The ground state electronic configuration Bi is [Xe] 4f<sup>14</sup>5d<sup>10</sup>6s<sup>2</sup>p<sup>3</sup> [<xref ref-type="bibr" rid="scirp.133512-ref2">2</xref>] . The covalent radius for Bi is 0.152 nm [<xref ref-type="bibr" rid="scirp.133512-ref2">2</xref>] . In a detailed survey of literature, Kabata-Pendias [<xref ref-type="bibr" rid="scirp.133512-ref4">4</xref>] noted the crustal average Bi concentrations are approximately 0.2 mg kg<sup>−</sup><sup>1</sup>, with argillaceous sediments showing slightly larger concentrations up to 4 mg kg<sup>−</sup><sup>1</sup>. In the United States, Govindaraju [<xref ref-type="bibr" rid="scirp.133512-ref15">15</xref>] reported that soil Bi averages range from 0.03 to 0.69 mg kg<sup>−</sup><sup>1</sup>. Hou et al. [<xref ref-type="bibr" rid="scirp.133512-ref5">5</xref>] reported average Bi concentration ranges in Japanese soils were 0.32 &#177; 0.12 mg kg<sup>−</sup><sup>1</sup>. Hou et al. [<xref ref-type="bibr" rid="scirp.133512-ref5">5</xref>] reported that average Bi distributions among the more dominant chemical fractions were noncrystalline (26%), peroxide extractable organic (26%) and metal organic (19%), residual (17%) and crystalline iron (12%). In Japan, Manaka [<xref ref-type="bibr" rid="scirp.133512-ref16">16</xref>] documented that Sb, and Bi were positively correlated with the amorphous Fe<sub>2</sub>O<sub>3</sub>.</p><p>Murata [<xref ref-type="bibr" rid="scirp.133512-ref17">17</xref>] investigated Bi solubility as influenced by pH and the presence of EDTA, citric acid, tartaric acid, L-cysteine, soil humic acids, and dissolved organic matter. Solution pH and the presence of citric acid, tartaric acid, L-cysteine, and soil humic acids influenced bismuth solubility. Kleja et al. [<xref ref-type="bibr" rid="scirp.133512-ref18">18</xref>] investigated the binding of Bi3+ to organic soil materials and noted that Bi3+ formed organic soil complexes, forming a dimeric Bi3+ complex.</p><p>The objectives of this investigation are: i) to estimate the soil abundances of antimony and bismuth across southeastern Missouri, and ii) to determine their potential soil profile distribution because of eluviation-illuviation and iron oxide abundances.</p></sec></sec><sec id="s2"><title>2. Materials and Methods</title><sec id="s2_1"><title>2.1. Study Area</title><p>The study area in Missouri is located between the Mississippi River and the St. Francois River. The northern section consists of thin to thick loess mantles overlying primarily Precambrian igneous and Ordovician carbonate rocks. The southern portion is in the Mississippi River embayment and consists of floodplains and terraces that have coarse to fine sediments.</p><p>The climate is continental humid. The average daily January temperatures are 2 to 4˚C (35 to 39˚F), whereas the average summer temperatures are 25 to 26˚C (77 to 79˚F). The rainfall is reasonably well distributed, with the total annual precipitation averaging 1.20 m. The remnants of tropical storms from the Gulf of Mexico provide periodic intense rainfall events [<xref ref-type="bibr" rid="scirp.133512-ref19">19</xref>] [<xref ref-type="bibr" rid="scirp.133512-ref20">20</xref>] [<xref ref-type="bibr" rid="scirp.133512-ref21">21</xref>] .</p></sec><sec id="s2_2"><title>2.2. Methods</title><p>Soils were selected from the following soil orders: i) Mollisols, ii) Alfisols, iii) Ultisols, iv) Entisols, v) Inceptisols, and vi) Vertisols. In total, 27 soil series were selected, many with multiple pedons. The soils used in this investigation were routinely characterized: i) to verify that the pedon was a member of the soil series, and ii) to provide routine soil chemical characterization. Laboratory analysis was performed only on the fine earth fraction, that is material finer than 2 mm. Standard routine methods included pH in water, exchangeable cations, total neutralizable acidity, and organic matter content by loss on ignition. These methods were performed by the soil testing laboratory at the University Missouri-Columbia Fisher Delta Center (Portageville, MO). Soil taxonomic classifications were from the United States Department of Agriculture official soil series descriptions [<xref ref-type="bibr" rid="scirp.133512-ref22">22</xref>] .</p><p>An aqua regia digestion was employed to obtain a near total estimation of elemental abundance associated with all but the most recalcitrant soil chemical environments. Homogenized samples (0.75 g) were equilibrated with 0.01 liter of aqua-regia (3 mole nitric acid: 1 mole hydrochloric acid) in a 35˚C incubator for 24 hours. Samples were shaken, centrifuged, and filtered (0.45 &#181;m), with a known aliquot volume analyzed using inductively coupled plasma emission – mass spectrometry. Selected samples were duplicated and known reference materials were employed to guarantee analytical accuracy. A water extraction was performed to recover only the most labile or potentially labile fractions. A hot water extraction involved equilibrating 0.5 g samples in 0.02 L distilled-deionized water at 80˚C for one hour followed by 0.45 &#181;m filtering and elemental determination using inductively coupled plasma emission – mass spectrometry. For the water extraction, selected samples were duplicated, and reference materials were employed to guarantee analytical precession. Simple statistics included mean, standard deviation (STD), coefficient of variation, and linear regression analysis were each performed using Excel.</p></sec></sec><sec id="s3"><title>3. Soil Series Characterization</title><p>Twenty-seven soil series were characterized, involving ten series from the Alfisol order, five series from the Ultisol order, five series from the Entisol order, three series from the Inceptisol order, three series from the Mollisol order and one series from the Vertisol order. Most of the soil series were deep to very deep, whereas a few soil series were shallow to moderately deep. Soil profiles ranged from excessively well-drained to poorly-drained. The soil classifications for the 27-soil series are listed in Appendix 1.</p></sec><sec id="s4"><title>4. Antimony and Bismuth Soil Profile Concentrations</title><p>The mean soil profile antimony concentrations range from 0.65 mg kg<sup>−</sup><sup>1</sup> for the Foley pedon to 0.08 mg kg<sup>−</sup><sup>1</sup> for the Killarney pedon, which are within the typical soil concentration ranges reported by the cited literature. The Irondale (0.53 mg kg<sup>−</sup><sup>1</sup>) and Taumsauk (0.5 mg kg<sup>−</sup><sup>1</sup>) pedon are adjacent to each other and are composed of mass-wasting loess mixed with rhyolite residuum on steep sideslopes, whereas the Knobtop (0.08 mg kg<sup>−</sup><sup>1</sup>) resides on summit positions and has a comparatively thick loess mantle discretely overlying rhyolite (<xref ref-type="table" rid="table">Table </xref>1). Thus, the stark differences are likely attributed to parent material inheritance.</p><p>The mean soil profile bismuth concentrations range from 0.92 mg kg<sup>−</sup><sup>1</sup> for the Wakeland pedon and 0.53 mg kg<sup>−</sup><sup>1</sup> for the Haymond pedon to 0.03 mg kg<sup>−</sup><sup>1</sup> for the Malden pedon and 0.05 mg kg<sup>−</sup><sup>1</sup> for the Clana pedon. Interesting the Wakeland and Haymond are adjacent pedons in a silty-textured floodplain, whereas the adjacent Malden and Clana pedons reside on coarse-textured terrace positions. The Sb and Bi mean and standard deviations for all soil series are in <xref ref-type="table" rid="table">Table </xref>1.</p><p>The mean and standard deviations for Sb resulting from a hot water extraction was intended to estimate the quantity of antimony that is biologically available. Bi concentrations from the hot water extraction were typically below detection limit (0.8 &#181;g kg<sup>−</sup><sup>1</sup>). The Sb hot water extraction values range from 0.27 &#181;g kg<sup>−</sup><sup>1</sup> for the Lilbourn soil series to 5.94 &#181;m kg<sup>−</sup><sup>1</sup> for the Calhoun soil series (<xref ref-type="table" rid="table">Table </xref>2). The mean soil profile aqua regia digestion concentration for the Calhoun soil series was 0.24 mg kg<sup>−</sup><sup>1</sup>, whereas the hot water extraction was 0.0059 mg kg<sup>−</sup><sup>1</sup>, or 2.5% of the total antimony concentration.</p><p>The coefficient of variation for antimony and bismuth may be easily calculated (standard deviation * 100 / mean) for each soil series. The mean and standard deviation of the Sb and Bi coefficients of variation were determined for those soil series having argillic horizons and lacking argillic horizons. Argillic horizons are defined as mineral soil horizons that are characterized by the illuvial accumulation of layer-lattice silicate clays, which is soil horizons having accumulated translocated clay from superimposed soil horizons. The mean and standard deviation for the coefficients of variation for antimony and bismuth for soil series having and lacking argillic horizons are displayed in <xref ref-type="table" rid="table">Table </xref>3. Application of t-test for mean separation demonstrates that antimony means are significantly different between soils having and lacking argillic horizons (α = 0.015), whereas bismuth means are not significantly different between soils having and lacking argillic horizons (α = 0.13).</p><p>The implication is that the greater antimony coefficient of variations for soils having argillic horizons is attributed to the greater antimony concentrations of the clay-enriched illuvial soil horizons. The greater clay contents, with their associated Fe-oxyhydroxides, provide a more intense antimony adsorption capacity. In this manuscript, Entisols, Inceptisols and the Vertisol samples have rather uniform soil profile textures and correspondingly rather uniform adsorption capacities, thus antimony soil profile concentration distributions have smaller variances.</p><table-wrap id="table1" ><label><xref ref-type="table" rid="table">Table </xref>1</label><caption><title> Mean and standard deviation for antimony and bismuth in soils</title></caption><table><tbody><thead><tr><th align="center" valign="middle" ></th><th align="center" valign="middle"  colspan="2"  >Antimony (mg / kg)</th><th align="center" valign="middle"  colspan="2"  >Bismuth (mg / kg)</th></tr></thead><tr><td align="center" valign="middle" >Soil</td><td align="center" valign="middle" >Mean</td><td align="center" valign="middle" >Standard dev</td><td align="center" valign="middle" >Mean</td><td align="center" valign="middle" >Standard dev</td></tr><tr><td align="center" valign="middle" >Alred</td><td align="center" valign="middle" >0.47</td><td align="center" valign="middle" >0.18</td><td align="center" valign="middle" >0.15</td><td align="center" valign="middle" >0.08</td></tr><tr><td align="center" valign="middle" >Amagon</td><td align="center" valign="middle" >0.26</td><td align="center" valign="middle" >0.13</td><td align="center" valign="middle" >0.09</td><td align="center" valign="middle" >0.04</td></tr><tr><td align="center" valign="middle" >Broseley</td><td align="center" valign="middle" >0.20</td><td align="center" valign="middle" >0.05</td><td align="center" valign="middle" >0.07</td><td align="center" valign="middle" >0.03</td></tr><tr><td align="center" valign="middle" >Calhoun</td><td align="center" valign="middle" >0.24</td><td align="center" valign="middle" >0.08</td><td align="center" valign="middle" >0.07</td><td align="center" valign="middle" >0.04</td></tr><tr><td align="center" valign="middle" >Clana</td><td align="center" valign="middle" >0.28</td><td align="center" valign="middle" >0.06</td><td align="center" valign="middle" >0.05</td><td align="center" valign="middle" >0.02</td></tr><tr><td align="center" valign="middle" >Commerce</td><td align="center" valign="middle" >0.33</td><td align="center" valign="middle" >0.06</td><td align="center" valign="middle" >0.17</td><td align="center" valign="middle" >0.05</td></tr><tr><td align="center" valign="middle" >Dubbs</td><td align="center" valign="middle" >0.30</td><td align="center" valign="middle" >0.09</td><td align="center" valign="middle" >0.13</td><td align="center" valign="middle" >0.03</td></tr><tr><td align="center" valign="middle" >Foley</td><td align="center" valign="middle" >0.65</td><td align="center" valign="middle" >0.22</td><td align="center" valign="middle" >0.22</td><td align="center" valign="middle" >0.05</td></tr><tr><td align="center" valign="middle" >Frenchmill</td><td align="center" valign="middle" >0.10</td><td align="center" valign="middle" >0.02</td><td align="center" valign="middle" >0.17</td><td align="center" valign="middle" >0.03</td></tr><tr><td align="center" valign="middle" >Haymond</td><td align="center" valign="middle" >0.19</td><td align="center" valign="middle" >0.03</td><td align="center" valign="middle" >0.53</td><td align="center" valign="middle" >0.30</td></tr><tr><td align="center" valign="middle" >Hildebrecht</td><td align="center" valign="middle" >0.21</td><td align="center" valign="middle" >0.05</td><td align="center" valign="middle" >0.18</td><td align="center" valign="middle" >0.07</td></tr><tr><td align="center" valign="middle" >Irondale</td><td align="center" valign="middle" >0.53</td><td align="center" valign="middle" >0.15</td><td align="center" valign="middle" >0.20</td><td align="center" valign="middle" >0.02</td></tr><tr><td align="center" valign="middle" >Kaintuck</td><td align="center" valign="middle" >0.24</td><td align="center" valign="middle" >0.02</td><td align="center" valign="middle" >0.13</td><td align="center" valign="middle" >0.01</td></tr><tr><td align="center" valign="middle" >Killarney</td><td align="center" valign="middle" >0.09</td><td align="center" valign="middle" >0.02</td><td align="center" valign="middle" >0.11</td><td align="center" valign="middle" >0.01</td></tr><tr><td align="center" valign="middle" >Knobtop</td><td align="center" valign="middle" >0.08</td><td align="center" valign="middle" >0.03</td><td align="center" valign="middle" >0.16</td><td align="center" valign="middle" >0.04</td></tr><tr><td align="center" valign="middle" >Lilbourn</td><td align="center" valign="middle" >0.26</td><td align="center" valign="middle" >0.09</td><td align="center" valign="middle" >0.10</td><td align="center" valign="middle" >0.09</td></tr><tr><td align="center" valign="middle" >Malden</td><td align="center" valign="middle" >0.25</td><td align="center" valign="middle" >0.05</td><td align="center" valign="middle" >0.03</td><td align="center" valign="middle" >0.01</td></tr><tr><td align="center" valign="middle" >Menfro</td><td align="center" valign="middle" >0.38</td><td align="center" valign="middle" >0.10</td><td align="center" valign="middle" >0.20</td><td align="center" valign="middle" >0.05</td></tr><tr><td align="center" valign="middle" >Overcup</td><td align="center" valign="middle" >0.44</td><td align="center" valign="middle" >0.26</td><td align="center" valign="middle" >0.20</td><td align="center" valign="middle" >0.03</td></tr><tr><td align="center" valign="middle" >Portageville</td><td align="center" valign="middle" >0.38</td><td align="center" valign="middle" >0.06</td><td align="center" valign="middle" >0.35</td><td align="center" valign="middle" >0.01</td></tr><tr><td align="center" valign="middle" >Reelfoot</td><td align="center" valign="middle" >0.39</td><td align="center" valign="middle" >0.07</td><td align="center" valign="middle" >0.17</td><td align="center" valign="middle" >0.03</td></tr><tr><td align="center" valign="middle" >Rueter</td><td align="center" valign="middle" >0.32</td><td align="center" valign="middle" >0.08</td><td align="center" valign="middle" >0.10</td><td align="center" valign="middle" >0.04</td></tr><tr><td align="center" valign="middle" >Sharkey</td><td align="center" valign="middle" >0.39</td><td align="center" valign="middle" >0.05</td><td align="center" valign="middle" >0.34</td><td align="center" valign="middle" >0.04</td></tr><tr><td align="center" valign="middle" >Taumsauk</td><td align="center" valign="middle" >0.50</td><td align="center" valign="middle" >0.17</td><td align="center" valign="middle" >0.19</td><td align="center" valign="middle" >0.04</td></tr><tr><td align="center" valign="middle" >Tiptonville</td><td align="center" valign="middle" >0.41</td><td align="center" valign="middle" >0.08</td><td align="center" valign="middle" >0.16</td><td align="center" valign="middle" >0.04</td></tr><tr><td align="center" valign="middle" >Wakeland</td><td align="center" valign="middle" >0.16</td><td align="center" valign="middle" >0.03</td><td align="center" valign="middle" >0.92</td><td align="center" valign="middle" >0.51</td></tr><tr><td align="center" valign="middle" >Wilbur</td><td align="center" valign="middle" >0.39</td><td align="center" valign="middle" >0.09</td><td align="center" valign="middle" >0.13</td><td align="center" valign="middle" >0.04</td></tr></tbody></table></table-wrap><table-wrap id="table2" ><label><xref ref-type="table" rid="table">Table </xref>2</label><caption><title> Mean and standard deviation for antimony water extract</title></caption><table><tbody><thead><tr><th align="center" valign="middle" ></th><th align="center" valign="middle"  colspan="2"  >Antimony (&#181;g / kg)</th></tr></thead><tr><td align="center" valign="middle" >Soil</td><td align="center" valign="middle" >Mean</td><td align="center" valign="middle" >Standard Deviation</td></tr><tr><td align="center" valign="middle" >Amagon</td><td align="center" valign="middle" >3.43</td><td align="center" valign="middle" >0.96</td></tr><tr><td align="center" valign="middle" >Calhoun</td><td align="center" valign="middle" >5.94</td><td align="center" valign="middle" >6.22</td></tr><tr><td align="center" valign="middle" >Irondale</td><td align="center" valign="middle" >4.56</td><td align="center" valign="middle" >3.66</td></tr><tr><td align="center" valign="middle" >Kaintuck</td><td align="center" valign="middle" >1.98</td><td align="center" valign="middle" >0.36</td></tr><tr><td align="center" valign="middle" >Lilbourn</td><td align="center" valign="middle" >0.27</td><td align="center" valign="middle" >0.10</td></tr><tr><td align="center" valign="middle" >Wilbur</td><td align="center" valign="middle" >4.27</td><td align="center" valign="middle" >1.31</td></tr></tbody></table></table-wrap><table-wrap id="table3" ><label><xref ref-type="table" rid="table">Table </xref>3</label><caption><title> The mean and standard deviation for the coefficients of variation for Sb and Bi for soil series having and lacking argillic horizons</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Statistic</th><th align="center" valign="middle"  colspan="2"  >Soils Having Argillic Horizons</th><th align="center" valign="middle"  colspan="2"  >Soils Lacking Argillic Horizons</th></tr></thead><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" >Sb</td><td align="center" valign="middle" >Bi</td><td align="center" valign="middle" >Sb</td><td align="center" valign="middle" >Bi</td></tr><tr><td align="center" valign="middle" >Mean</td><td align="center" valign="middle" >30</td><td align="center" valign="middle" >27</td><td align="center" valign="middle" >19</td><td align="center" valign="middle" >39</td></tr><tr><td align="center" valign="middle" >Standard Dev.</td><td align="center" valign="middle" >11</td><td align="center" valign="middle" >15</td><td align="center" valign="middle" >7</td><td align="center" valign="middle" >25</td></tr></tbody></table></table-wrap><p>Standard Dev is standard deviation.</p></sec><sec id="s5"><title>5. Relationship of Antimony and Bismuth with Iron</title><p>The antimony and bismuth concentrations for the Alred and Rueter pedons were pooled given their profile similarities and their adjacent geographic locations. The antimony and bismuth concentrations show positive linear relationships with respect to iron (<xref ref-type="fig" rid="fig1">Figure 1</xref> and <xref ref-type="fig" rid="fig2">Figure 2</xref>). Similarly, the Menfro pedons show positive antimony and bismuth linear relationships with respect to iron (<xref ref-type="fig" rid="fig3">Figure 3</xref> and <xref ref-type="fig" rid="fig4">Figure 4</xref>). Two features are prominent: i) the iron concentration variance is predicated on greater Fe concentrations associated with the phyllosilicate enriched argillic horizons, and ii) antimony and bismuth concentrations are correlated with iron.</p><p>The soil profile distribution of antimony and bismuth in the Menfro pedons (<xref ref-type="fig" rid="fig5">Figure 5</xref> and <xref ref-type="fig" rid="fig6">Figure 6</xref>) illustrate the importance of the argillic horizon’s antimony and bismuth adsorption capacity. Similarly, the antimony soil profile distribution demonstrates greater concentrations for the illuvial horizons (<xref ref-type="fig" rid="fig7">Figure 7</xref>); however, the portageville pedon demonstrates rather uniform antimony and bismuth concentrations (<xref ref-type="fig" rid="fig8">Figure 8</xref> and <xref ref-type="fig" rid="fig9">Figure 9</xref>).</p></sec><sec id="s6"><title>6. Discussion</title><sec id="s6_1"><title>6.1. Assessment of Antimony and Bismuth Soil Concentrations across Southeastern Missouri</title><p>The antimony concentrations ranged from very low concentrations to 0.65 mg kg<sup>−</sup><sup>1</sup>, whereas the bismuth soil profile concentrations ranged from very low concentrations to 0.92 mg kg<sup>−</sup><sup>1</sup>. For the Alfisols and Ultisols, the antimony soil profile distributions show greater concentrations in the illuvial horizons than the eluvial horizons, whereas for the antimony soil profile distributions for the majority of the Entisols, Inceptisols, Mollisols and the solitary Vertisol show rather uniform concentrations. Zhao et al. [<xref ref-type="bibr" rid="scirp.133512-ref23">23</xref>] noted that adsorption was an important phenomenon influencing antimony accumulation, with clay and iron oxyhydroxides important soil substrates.</p><p>The magnitude of the antimony and bismuth concentrations are within the concentration ranges documented for pristine, non-impacted (geogenic) soils [<xref ref-type="bibr" rid="scirp.133512-ref4">4</xref>] [<xref ref-type="bibr" rid="scirp.133512-ref23">23</xref>] [<xref ref-type="bibr" rid="scirp.133512-ref24">24</xref>] [<xref ref-type="bibr" rid="scirp.133512-ref25">25</xref>] [<xref ref-type="bibr" rid="scirp.133512-ref26">26</xref>] . In their review, Bolan et al. [<xref ref-type="bibr" rid="scirp.133512-ref24">24</xref>] did report that background and average antimony concentration values were 0.3 to 8.6 mg kg<sup>−</sup><sup>1</sup>. Thus, the sampled soils across southeastern Missouri are not inferred to be impacted by mining or other anthropogenic activities.</p></sec><sec id="s6_2"><title>6.2. Potential Environmental Risks Associated with Antimony and Bismuth across Southeastern Missouri Are Presently Extremely Limited</title><p>Tang et al. [<xref ref-type="bibr" rid="scirp.133512-ref26">26</xref>] recently described sources of soil antimony pollution and its environmental impact, and proposed remediation techniques, noting use of ferrous sulfate, phosphate amendments, clay minerals and biochar for adsorption, selected bacterial communities, and phytoremediation. Antimony pollution sources include mining, pharmaceutical manufacturing, vehicle emissions, plastic waste leaching, shooting ranges, and others [<xref ref-type="bibr" rid="scirp.133512-ref23">23</xref>] [<xref ref-type="bibr" rid="scirp.133512-ref25">25</xref>] [<xref ref-type="bibr" rid="scirp.133512-ref26">26</xref>] . The maximum concentration of antimony for drinking water is 6 ppb by the United States Environmental Protection Agency and 20 ppb by World Health Organization [<xref ref-type="bibr" rid="scirp.133512-ref26">26</xref>] . Vidya et al. [<xref ref-type="bibr" rid="scirp.133512-ref13">13</xref>] reviewed the influence of antimony on plant growth and development, including i) stunted growth, ii) reduced photosynthesis and accumulation biomass, iii) generation of reactive oxygen species, and lipid peroxidation. However, the literature addressing antimony and bismuth as environmental risks is a relatively recent activity.</p><p>Bolan et al. [<xref ref-type="bibr" rid="scirp.133512-ref24">24</xref>] noted that no definitive review of the biogeochemistry of antimony has completely described antimony mobilization, bioavailability, toxicity, and threats to environmental and human health. Bolin et al. [<xref ref-type="bibr" rid="scirp.133512-ref24">24</xref>] significantly reviewed recent literature of the biogeochemical processes influencing soil antimony, noting that clay minerals and the oxyhydroxides of Mn, Al and Fe were implemental in regulating antimony soil chemistry.</p></sec></sec><sec id="s7"><title>7. Conclusion</title><p>The soil profile concentrations of antimony and bismuth were determined for key soil series across southeastern Missouri. The antimony concentrations ranged from 0.65 to 0.08 mg kg<sup>−</sup><sup>1</sup>, whereas the bismuth soil profile concentrations ranged from 0.92 to 0.03 mg kg<sup>−</sup><sup>1</sup>. For soils having argillic horizons, antimony and bismuth concentrations were greater for the illuvial horizons than the eluvial horizons. Both elements show significant correlations with iron.</p></sec><sec id="s8"><title>Conflicts of Interest</title><p>The author declares no conflicts of interest regarding the publication of this paper.</p></sec><sec id="s9"><title>Cite this paper</title><p>Aide, M. (2024) Soil Profile Concentration Distributions of Antimony and Bismuth across Southeastern Missouri (USA). Journal of Environmental Protection, 15, 537-551. https://doi.org/10.4236/jep.2024.155031</p></sec><sec id="s10"><title>Appendix 1. Soil Series Characterization</title></sec></body><back><ref-list><title>References</title><ref id="scirp.133512-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">Aide, M.T., Beighley, D. and Dunn, D. (2013) Soil Profile Arsenic Concentration Distributions in Missouri Soils Having Cambic and Argillic Soil Horizons. &lt;i&gt;Soil S&lt;/i&gt;&lt;i&gt;e&lt;/i&gt;&lt;i&gt;diment Contamination&lt;/i&gt;, 23, 313-327. &lt;br&gt;https://doi.org/10.1080/15320383.2014.831028</mixed-citation></ref><ref id="scirp.133512-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">Lee, J.D. (1991) Concise Inorganic Chemistry. Chapman Hall, London.</mixed-citation></ref><ref id="scirp.133512-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">Filella, M., Belzile, N. and Chen, Y.W. (2002) Antimony in the Environment: A Review Focused on Natural Waters: I. Occurrence. &lt;i&gt;Earth&lt;/i&gt;-&lt;i&gt;Science Reviews&lt;/i&gt;, 57, 125-176. &lt;br&gt;https://doi.org/10.1016/S0012-8252(01)00070-8</mixed-citation></ref><ref id="scirp.133512-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">Kabata-Pendias, A. (2011) Trace Elements in Soils and Plants. CRC Press, Boca Raton. &lt;br&gt;https://doi.org/10.1201/b10158</mixed-citation></ref><ref id="scirp.133512-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">Hou, H., Takamatsu, T., Koshikawa, M.K. and Hosomi, M. (2005) Concentrations of Ag, In, Sn, Sb and Bi, and Their Chemical Fractionation in Typical Soils in Japan. &lt;i&gt;European Journal of Soil Science&lt;/i&gt;, 57, 214-227. &lt;br&gt;https://doi.org/10.1111/j.1365-2389.2005.00731.x</mixed-citation></ref><ref id="scirp.133512-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">Kolesnikov, S., Minnikova, T., Kazeev, K., Akimenko, Y. and Evstegneeva, N. (2022) Assessment of the Ecotoxicity of Pollution by Potentially Toxic Elements by Biological Indicators of Haplic Chernozem of Southern Russia (Rostov Region). &lt;i&gt;Water&lt;/i&gt;,&lt;i&gt; Air&lt;/i&gt;,&lt;i&gt; &amp; Soil Pollution&lt;/i&gt;, 233, Article No. 18. &lt;br&gt;https://doi.org/10.1007/s11270-021-05496-3</mixed-citation></ref><ref id="scirp.133512-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">Lintschinger, J., Michalke, B., Schulte-Hostede, S. and Schramel, P. (1997) Studies on Speciation of Antimony in Soil Contaminated by Industrial Activity. &lt;i&gt;Intern&lt;/i&gt;&lt;i&gt;a&lt;/i&gt;&lt;i&gt;tional Journal of Environmental Analytical Chemistry&lt;/i&gt;, 72, 11-25. &lt;br&gt;https://doi.org/10.1080/03067319808032641</mixed-citation></ref><ref id="scirp.133512-ref8"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">Filella, M., Belzile, N. and Chen, Y.W. (2002) Antimony in the Environment: A Review Focused on Natural Waters: II. Relevant Solution Chemistry. &lt;i&gt;Earth&lt;/i&gt;-&lt;i&gt;Science Reviews&lt;/i&gt;, 59, 265-285. &lt;br&gt;https://doi.org/10.1016/S0012-8252(02)00089-2</mixed-citation></ref><ref id="scirp.133512-ref9"><label>9</label><mixed-citation publication-type="other" xlink:type="simple">Pilarski, J., Waller, P. and Pickering, W. (1995) Sorption of Antimony Species by Humic Acid. &lt;i&gt;Water&lt;/i&gt;,&lt;i&gt; Air&lt;/i&gt;,&lt;i&gt; &amp; Soil Pollution&lt;/i&gt;, 84, 51-59. &lt;br&gt;https://doi.org/10.1007/BF00479588</mixed-citation></ref><ref id="scirp.133512-ref10"><label>10</label><mixed-citation publication-type="other" xlink:type="simple">Bagherifam, S., Brown, T.C., Bagherifam, S. and Baglieri, A. (2023) Sequential Extraction of Labile and Recalcitrant Fractions of Soil Organic Matter: A Case Study Focusing on Antimony (Sb) in Humic Acids, Fulvic Acids and Humin Fractions of Long-Term Aged Contaminated Soils. &lt;i&gt;Environmental Pollution&lt;/i&gt;, 327, Article ID: 121610. &lt;br&gt;https://doi.org/10.1016/j.envpol.2023.121610</mixed-citation></ref><ref id="scirp.133512-ref11"><label>11</label><mixed-citation publication-type="other" xlink:type="simple">Thanabalasingam, P. and Pickering, W.F. (1990) Specific Sorption of Antimony (III) by the Hydrous Oxides of Mn, Fe, and Al. &lt;i&gt;Water&lt;/i&gt;,&lt;i&gt; Air&lt;/i&gt;,&lt;i&gt; &amp; Soil Pollution&lt;/i&gt;, 49, 175-185. &lt;br&gt;https://doi.org/10.1007/BF00279519</mixed-citation></ref><ref id="scirp.133512-ref12"><label>12</label><mixed-citation publication-type="other" xlink:type="simple">Tang, H., Meng, G., Xiang, J., Mahmood, A., Xiang, G., Sana, U., Liu, Y. and Huang, G. (2022) Toxic Effects of Antimony in Plants: Reasons and Remediation Possibilities&amp;#8212;A Review and Future Prospects. &lt;i&gt;Frontiers in Plant Science&lt;/i&gt;, 13, Article ID: 1011945. &lt;br&gt;https://doi.org/10.3389/fpls.2022.1011945</mixed-citation></ref><ref id="scirp.133512-ref13"><label>13</label><mixed-citation publication-type="other" xlink:type="simple">Vidya, S.-N., Shetty, R., Vaculikova, M. and Vaculik, M. (2022) Antimony Toxicity in Soils and Plants and Mechanisms of Its Alleviation. &lt;i&gt;Environmental and Exper&lt;/i&gt;&lt;i&gt;i&lt;/i&gt;&lt;i&gt;mental Botany&lt;/i&gt;, 202, Article ID: 104996. &lt;br&gt;https://doi.org/10.1016/j.envexpbot.2022.104996</mixed-citation></ref><ref id="scirp.133512-ref14"><label>14</label><mixed-citation publication-type="other" xlink:type="simple">Bowen, H.J.M. (1979) Environmental Chemistry of Elements. Academic Press, Cambridge.</mixed-citation></ref><ref id="scirp.133512-ref15"><label>15</label><mixed-citation publication-type="other" xlink:type="simple">Govindaraju, K. (1994) 1994 Compilation of Working Values and Sample Description for 383 Geostandards. &lt;i&gt;Geostandards Newsletter&lt;/i&gt;, 18, 1-158. &lt;br&gt;https://doi.org/10.1046/j.1365-2494.1998.53202081.x-i1</mixed-citation></ref><ref id="scirp.133512-ref16"><label>16</label><mixed-citation publication-type="other" xlink:type="simple">Manaka, M. (2006) Amount of Amorphous Materials in Relationship to Arsenic, Antimony, and Bismuth Concentrations in a Brown Forest Soil. &lt;i&gt;Geoderma&lt;/i&gt;, 136, 75-86. &lt;br&gt;https://doi.org/10.1016/j.geoderma.2006.02.002</mixed-citation></ref><ref id="scirp.133512-ref17"><label>17</label><mixed-citation publication-type="other" xlink:type="simple">Murata, T. (2010) Bismuth Solubility through Binding by Various Organic Compounds and Naturally Occurring Soil Organic Matter. &lt;i&gt;Journal of Environmental Science and Health&lt;/i&gt;,&lt;i&gt; Part A&lt;/i&gt;, 45, 746-753. &lt;br&gt;https://doi.org/10.1080/10934521003651465</mixed-citation></ref><ref id="scirp.133512-ref18"><label>18</label><mixed-citation publication-type="other" xlink:type="simple">Kleja, D.B., Gustafsson, J.P., Kessler, V. and Persson, I. (2022) Bismuth(III) Forms Exceptionally Strong Complexes with Natural Organic Matter. &lt;i&gt;Environmental Science &lt;/i&gt;&lt;i&gt;&amp; Technology&lt;/i&gt;, 56, 3076-3084. &lt;br&gt;https://doi.org/10.1021/acs.est.1c06982</mixed-citation></ref><ref id="scirp.133512-ref19"><label>19</label><mixed-citation publication-type="other" xlink:type="simple">Festervand, D.F. (1981) Soil Survey of Cape Girardeau, Scott and Mississippi Counties, Missouri. Produced in Cooperation with the United States Department of Agriculture, United States Forest Service, and the University Missouri-Columbia. Printed Washington DC.</mixed-citation></ref><ref id="scirp.133512-ref20"><label>20</label><mixed-citation publication-type="other" xlink:type="simple">Brown, B.C. and Childress, J.D. (1985) Soil Survey of Ste. Genevieve County, Missouri. Produced in Cooperation with the United States Department of Agriculture, United States Forest Service, and the University Missouri-Columbia. Printed Washington DC.</mixed-citation></ref><ref id="scirp.133512-ref21"><label>21</label><mixed-citation publication-type="other" xlink:type="simple">Festervand, D.F. (1986) Soil Survey of Perry County, Missouri. Produced in Cooperation with the United States Department of Agriculture, United States Forest Service, and the University Missouri-Columbia. Printed Washington DC.</mixed-citation></ref><ref id="scirp.133512-ref22"><label>22</label><mixed-citation publication-type="other" xlink:type="simple">USDA Official Soil Series Descriptions. &lt;br&gt;https://www.nrcs.usda.gov/resources/data-and-reports/official-soil-series-descriptions-osd </mixed-citation></ref><ref id="scirp.133512-ref23"><label>23</label><mixed-citation publication-type="other" xlink:type="simple">Zhao, S., Shi, T., Terada, A. and Riya, S. (2022) Evaluation of Pollution Level, Spatial Distribution, and Ecological Effects of Antimony in Soils of Mining Areas: A Review. &lt;i&gt;International Journal of Environmental Research and Public Health&lt;/i&gt;, 20, Article No. 242. &lt;br&gt;https://doi.org/10.3390/ijerph20010242</mixed-citation></ref><ref id="scirp.133512-ref24"><label>24</label><mixed-citation publication-type="other" xlink:type="simple">Bolan, N., Kumar, M., Singh, E., Kumar, A., Singh, L., Kumar, S., Santhirasekaram, K., Hoang, A., El-Naggar, A., Vithanage, M., Sarkar, B., Wijeysekara, H., Diyabalanage, S., Sooriyakumar, P., Vinu, A., Wang, H., Kirkham, M., Shaheen, S., Rinklebe, J. and Siddique, K. (2022) Antimony Contamination and Its Risk Management in Complex Environmental Settings: A Review. &lt;i&gt;Environment International&lt;/i&gt;, 158, Article ID: 106908. &lt;br&gt;https://doi.org/10.1016/j.envint.2021.106908</mixed-citation></ref><ref id="scirp.133512-ref25"><label>25</label><mixed-citation publication-type="other" xlink:type="simple">Tschan, M., Robinson, B.H. and Schulin, R. (2009) Antimony in the Soil-Plant System&amp;#8212;A Review. &lt;i&gt;Environmental Chemistry&lt;/i&gt;, 6, 106-115. &lt;br&gt;https://doi.org/10.1071/EN08111</mixed-citation></ref><ref id="scirp.133512-ref26"><label>26</label><mixed-citation publication-type="other" xlink:type="simple">Tang, H., Hassan, M.U., Nawaz, M., Yang, W., Liu, Y. and Yang, B. (2023) A Review on Sources of Soil Antimony Pollution and Recent Progress on Remediation of Antimony Polluted Soils. &lt;i&gt;Ecotoxicology and Environmental Safety&lt;/i&gt;, 266, Article ID: 115583. &lt;br&gt;https://doi.org/10.1016/j.ecoenv.2023.115583</mixed-citation></ref></ref-list></back></article>