Yuze Zhoua,
Ming Dou*ab,
Yan Zhangbc,
Kaizi Ningd and
Yuxuan Lib
aSchool of Water Conservancy Engineering, Zhengzhou University, Zhengzhou, China
bSchool of Ecology and Environment, Zhengzhou University, Zhengzhou, China. E-mail: dou_ming@163.com
cFarmland Irrigation Research Institute, Chinese Academy of Agricultural Sciences, Xinxiang, China
dFaculty of Science, Monash University, Melbourne, Australia
First published on 15th July 2024
Microplastic (MPs) pollution has become a global issue, with particular concern regarding MPs in soil. To determine the characteristics of MPs in agricultural production areas and their impact on soil physicochemical properties, soil samples were collected from different land use types in the North China Plain. Layered sampling was conducted and the soil physicochemical properties were determined. A novel image recognition method based on fluorescence staining was proposed for the batch analysis of MPs in the study area. Together with the results of the soil physicochemical properties, the impact of MPs on soil physicochemical properties was analyzed and evaluated. The results showed that the soil MPs abundance in this agricultural area was moderate to low compared to other agricultural areas, with a larger proportion of particle-type and fragment-type MPs smaller than 10 μm. The soil MPs were predominantly composed of polyvinyl chloride (PVC) and polypropylene (PP). MPs abundance was higher in farmland and forest land than in vegetable fields. The impact of MPs on soil physicochemical properties was mainly manifested in the changes in soil structure due to the different MPs characteristics. Apart from abundance, the type of MPs was found to be the main factor affecting soil bulk density, with particle size and shape influencing the soil aggregate structure. MPs may effect the pH values of sandy and loamy soils, primarily by altering the soil porosity and water holding capacity, but also by increasing the area and duration of contact between the soil medium and external water sources. This study revealed the MPs characteristics in agricultural areas as well as the pathways by which they can impact soil physicochemical properties.
Environmental significanceMicroplastics (MPs) pollution has become a global issue, with particular concern regarding MPs in soil. But most studies have focused on the abundance and distribution characteristics of MPs in soil, with little discussion on the impact of MPs on soil physicochemical properties. And the traditional MPs identification analysis methods suffer from issues of low efficiency and high costs. Therefore, a novel image recognition method based on fluorescence staining was proposed for the batch analysis of MPs in the study area. According to the results of the soil physicochemical properties, the impact of MPs on soil physicochemical properties was analyzed and evaluated. |
Soil is generally considered to be a reservoir for MPs.13–15 In Europe, it is estimated that 1270–2130 tons of MPs enter the soil each year through various pathways. For example, in soil samples from the French soil quality monitoring network,16 rural vegetable fields in southeastern Spain,17 floodplain soils in Switzerland,18 and composted vegetable fields in southern Netherlands,19 the MPs abundances were reported to be 80 ind kg−1, 2130 ± 950 ind kg−1, 593 ind kg−1, and 67.34–1109 ind kg−1, respectively, with polyethylene (PE) and polystyrene (PS) being the main components.
In China, soils containing MPs have been reported in farmland areas in the hilly and gully region of Yan'an,20 farmland soils in Chengdu,21 farmland in the Hetao Irrigation District,22 farmland in the karst plateau of Guizhou,23 farmland in the rocky desertification area of Guizhou,24 the Yellow River Delta,25 farmland in Shouguang City,26 soil in vegetable bases in Xianyang City,27 farmland in suburban Tianjin,28 farmland on the Qinghai-Tibet Plateau,29 farmland in the Hangzhou Bay and southeast coastal plains,30 the Danjiangkou Reservoir area of the South-to-North Water Diversion Project,31 as well as the cities of Harbin32 and Shihezi.33 MPs abundances can reach up to 28100 items per kg, with the main components including polyethylene terephthalate (PET), polypropylene (PP), PE, and PS. Most previous studies have focused on the abundance and distribution characteristics of MPs in soil, with little discussion on the impact of MPs on soil physicochemical properties. The methods of extracting and observing MPs in these studies are also varied.
The distribution and environmental impacts of MPs have been more extensively studied in the atmosphere and aquatic environments than in soils because the extraction and sampling of MPs from air and water are relatively simple and the resulting particles are easy to identify.34–38 The extraction and identification of MPs are quite difficult due to the influences of inorganic substances (such as roots and minerals) and organic matter.39–41 At present, density extraction is a common method for extracting microplastics. After modification and improvement, the recovery rate of MPs by density extraction is generally above 70%.42 Commonly used extractants include NaCl, ZnCl2, NaBr and NaI.43–46 NaCl is the most commonly used, but the density of NaCl extracts is only 1.2 g cm−3, it can only extract MPs with lower density such as PP and PE. The density of ZnCl2, NaBr and NaI extracts are above 1.55 g cm−3, which are higher than most common types of MPs. But the cost of using NaI as an extract is high. ZnCl2 is corrosive, and in order to maximize the density of the ZnCl2 extract, additional acid solution is usually required, which will seriously hinder the extraction process.43,47
Existing methods for MPs identification and counting primarily rely on visual identification, Raman spectroscopy and FT-IR.48–52
Raman spectroscopy is commonly coupled with microscopy techniques and can be used to identify MPs with sizes larger than 1 μm. However, the detection process is time-consuming and can be easily influenced by the spontaneous fluorescence of soil organic matter.53,54 FT-IR is only suitable for detecting microplastics larger than 20 μm and is prone to interference from organic matter.55,56 Fourier Transform Infrared Microspectroscopy combines the advantages of FT-IR and microscopy, enabling the batch identification of MPs within a certain range. However, due to its high cost, its popularity is low.57
As one of China's three major plains, the North China Plain is an important hub for population aggregation and grain production. Consequently, considerable attention has been given to its ecological environment and food security. MPs, as a novel pollutant in soil, are widely distributed in the farmland soils of the North China Plain. According to existing research, the presence of MPs have already had an obvious impact on the soil ecology and food security of the North China Plain.58 Therefore, it is necessary to conduct field sampling surveys to clarify the extent and process of changes in the physicochemical properties of soil in the agricultural areas of the North China Plain caused by MPs. This will provide a basis for the subsequent management of soil microplastic pollution. Given the poor identification accuracy and high detection costs associated with existing MPs detection methods, a fluorescence staining-based image recognition method was proposed in this study. This method utilizes images observed under a fluorescence microscope, integrates image processing software to merge and extract the boundaries of MPs fluorescence images in different bands, and then uses geographic information system (GIS) software to batch analyze the quantity and characteristics of MPs in processed images. Simultaneously, the physicochemical properties of soils in the study area were tested, the characteristics of MPs were studied, and the impacts of MPs on the soil environment were analyzed. The results provide a reference for the management of MPs pollution in agricultural soils.
The main land use types in the study area are farmland and woodland. The crops in the farmland are mainly wheat and corn. According to the field investigation, no plastic film is used. The average annual volume of groundwater irrigation amounts to approximately 0.37 × 104 m3 per ha, while the average annual fertilizer application (including both base fertilizer and topdressing) stands at around 750 kg per ha. The crops in the forest are mainly economic trees such as birch and fruit trees. Except for the black plastic film used at the bottom of the fruit trees, the other crops basically do not use plastic film. The permanent population in study area is about 8300 people. There are three large-scale building ruins in the study area, as well as an east-west viaduct under construction and a supporting temporary accommodation area. The project has lasted for two years. These are all potential sources of microplastic pollution.
The study area experiences moderate precipitation, but with distinct interannual variability. The multi-year average precipitation is 639.4 mm (1956–2019), with most (53%) concentrated from July to September. The least precipitation occurs from December to February of the following year, constituting less than 5% of the annual precipitation. The multi-year average evaporation is 1822.9 mm (1956–2019), which is primarily concentrated in the months of April, May, and June.
Sample | East longitude | North latitude | Type of soil | Land use |
---|---|---|---|---|
1# | 113.8328 | 34.8446 | Sand | Vegetable plot |
2# | 113.8498 | 34.8394 | Loam | Vegetable plot |
3# | 113.8600 | 34.8568 | Silt | Farmland |
4# | 113.8929 | 34.8705 | Loam | Woodland |
5# | 113.8854 | 34.8594 | Loam | Farmland |
6# | 113.9074 | 34.8519 | Loam | Woodland |
7# | 113.8824 | 34.8453 | Sand | Farmland |
8# | 113.8767 | 34.8343 | Silt | Woodland |
9# | 113.8997 | 34.8312 | Sand | Vegetable plot |
10# | 113.9082 | 34.8277 | Silt | Woodland |
Additionally, some soil samples were placed in aluminum boxes and brought back to the laboratory for soil pH determination. After air-drying, the samples were ground and passed through a sieve to obtain soil samples smaller than 2 mm. A 10 g soil sample was weighed and placed into a 50 ml centrifuge tube, then 25 ml of deionized water was added and stirred for 1 min. After 30 min of settling, the pH of the supernatant was measured to determine the soil pH by pH instrument (PHSJ-4F, INESA INSTRUMENT).
The soil bulk density at each sampling site was determined using a drying method. Undisturbed soil samples from the surface layer at each site were collected in the upper, middle, and lower layers using a ring cutter (Φ50.46 × 50 mm), and then stored in sealed aluminium boxes.
A 25 g undisturbed soil sample was used to determine the distribution of soil aggregates using the wet sieving method. The soil was placed on a set of nylon sieves with mesh sizes of 2, 0.25, and 0.053 mm. After adding distilled water to submerge the soil to a depth of 3 cm, the samples were soaked for 5 min. Then, the sieves were slowly moved up and down 50 times in 3 cm increments. The samples collected on the upper part of each sieve were washed into beakers, separating the soil aggregates into macroaggregates (>2 mm, 2–0.25 mm) and microaggregates (0.25–0.05 and <0.05 mm). After drying, the aggregates were weighed to determine their mass by analytical balance (JD400-3, NANBEI INSTRUMENT).
MPs were extracted from soil samples using a saturated sodium bromide solution density separation method. The collected soil samples were ground and sieved through a 2 mm mesh sieve. A 20 g sieved soil samples was then added to 200 ml of saturated sodium bromide solution, stirred vigorously on a magnetic stirrer (ZNDL-6S, HONGHUAYIQI) for 1 h, and left to settle for 24 h to extract the supernatant. Hydrogen peroxide (H2O2) and deionized water were added to dilute the sodium bromide and hydrogen peroxide concentrations to 8.8% and 10%, respectively. Subsequently, all samples were digested in a constant temperature water bath at 70 °C for 24 h. According to statistics, the fluorescence intensity of chitin and cellulose stained in 10% (w/v) H2O2 solution decreased by an average of 53.7 and 52.5%, respectively, compared with the samples digested without adding H2O2.71 A Nile red stock solution standard solution (1 mg ml−1) was added to control the Nile red concentration in the supernatant to 10 μg ml−1, and the samples were further heated for an additional 0.5 h. The MPs in the supernatant were filtered and extracted using a 0.45 μm nylon filter (0.45 μm, Onion) membrane in a vacuum filtration device (SCJ-10, SUPO), and the filter membrane was placed in a Petri dish and allowed to dry for a subsequent analysis.69 To prevent external contamination, the separation, extraction, and observation of MPs were conducted in a clean and sealed environment.
The flow of MPs observations analysis is shown in Fig. 2. The filter membranes obtained through density separation were placed under a fluorescence microscope, and photographs of the membranes were taken under different fluorescence wavelengths. Using ImageJ software, the fluorescence images from each wavelength were stitched and fused together. The brightness threshold was used to extract the boundaries of each MPs in the images and eliminate noise and natural light reflections outside the selected MPs boundaries. Subsequently, ArcGIS 10.2 software was used to reclassify the processed images by converting raster files to feature files to obtain the boundaries of each MPs. A statistical analysis of the number, area, size, and shape of the MPs at each sampling site was then conducted.
The quantity of MPs was obtained by counting the number of polygon features in the shapefile of each sampling site. The abundance value was used to represent the MPs concentrations in the soil at each sampling site.
(1) |
The size and shape of MPs were identified by calculating the aspect ratio and the degree of fitting to minimum circumscribed circle in the shapefile.
(2) |
(3) |
Recognizing the shape information of images through shape parameters including aspect ratio is widely used in the field of medical tumor recognition and other image recognition research,74,75 but it is less used in MPs recognition. This study attempts to use aspect ratio and fitting degree to identify the shape of MPs. According to the area of each MPs, the minimum aspect ratio of the outer rectangle, and the minimum circumscribed circle area data calculated by ArcGIS 10.2. First, MPs shapes were filtered based on aspect ratios. MPs with P > 2 were considered to be fiber-type MPs. D is calculated for each MPs whose aspect ratio is less than 2. The smaller D is, the greater the deviation from the standard circle is. Based on this criterion, MPs were classified into fragments (P < 2, D < 0.8), pellets (P < 2, D > 0.8), and fibers (P > 2.5) (Fig. 3).
Fig. 3 Observation results of MPs with different shapes (a1, b1, c1 and d1) – natural light; (a2, b2, c2 and d2) – fluorescence excitation. |
The IBM SPSS statistics software was used for data processing. A single-factor analysis of variance (ANOVA) was used to conduct a significance analysis of the differences in MPs abundance among various sampling sites and different soil depths. Graphs were generated using Excel 2019 and Origin 2022 software.
Shape parameters were analysed for 15 samples collected from sampling sites 1 to 5 during the process of extracting MPs. Thresholds were determined for each parameter and some results of threshold analysis are presented in Table 2. From Table 2, it can be observed that pellets have smaller aspect ratios and higher degrees of fitting above 0.8 compared to other shapes such as fibers which have aspect ratios greater than 2 but lower degrees of fitting below 0.5. Fragments essentially appear as irregular polygons without specific shape characteristics. Based on these thresholds, MPs were classified into fragments (P < 2, D < 0.8), pellets (P < 2, D > 0.8), and fibers (P > 2).
To validate the identified thresholds, 3 fluorescence images obtained from sampling site 6 were used (Fig. 4). It is evident that the selected thresholds effectively identify different shapes of MPs and enable simultaneous identification and analysis within various ranges including particle size, quantity, and other features.
The MPs identification method from soil used in this study is more efficient than the traditional visual counting method and less expensive than FT-IR and Raman analysis methods.
The MPs abundance was highest at sampling site 10, while sampling site 1 had the lowest abundance. In forested areas, the MPs abundance ranged from 1650 to 19350 ind kg−1, with an average abundance of 10350 particles per kg. Significant differences in MPs abundance were observed among the sampling sites (p < 0.05). MPs abundance was highest at sampling site 10, while it was lowest at sampling site 4. In farmland areas, the MPs abundance ranged from 2300 to 24300 ind kg−1, with an average abundance of 8822 ind kg−1. MPs abundance was highest at sampling site 3 and lowest at sampling site 5. In vegetable fields, the MPs abundance ranged from 2500 to 14550 ind kg−1, with an average abundance of 5583 ind kg−1. MPs abundance was highest at sampling site 9 and lowest at sampling site 1.
In terms of sampling depth, the MPs abundance in the study area decreased with increasing depth. The highest MPs abundance was observed in the 0–10 cm soil samples, with an average abundance of 12165 ind kg−1. Sampling site 3 had the highest MPs abundance in surface soil, while sampling site 1 had the lowest MPs abundance in surface soil. In the 10–20 cm soil samples, the MPs abundance was lower compared to the surface soil, with an average abundance of 7550 ind kg−1. Sampling site 10 had the highest MPs abundance, while sampling site 5 had the lowest MPs abundance. The lowest MPs abundance was observed in the 20–30 cm soil samples, with an average abundance of 5670 ind kg−1. Sampling site 10 had the highest MPs abundance, while sampling site 7 had the lowest MPs abundance.
The MPs particle sizes in the study area soil could be largely classified into categories of <10, 10–50, 50–100, 100–200, and >200 μm (Fig. 6a). The <10 μm MPs particles had the highest abundance, ranging from 950 to 13000 ind kg−1, with an average abundance of 5040 ind kg−1. The highest abundance was found in the topsoil of sampling site 3, which had an abundance of 13000 ind kg−1, while the lowest abundance was detected in the bottom layer samples of sampling sites 2 and 4, which both had an abundance of 950 ind kg−1. As the MPs particle sizes increased, their abundance decreased accordingly. MPs particles >200 μm were the least widely distributed size class in the study area soil, ranging from 0 to 300 ind kg−1, with an average abundance of 97 ind kg−1. The highest abundance was found in the middle layer soil of sampling site 10, with an abundance of 300 ind kg−1, while the lowest abundance was detected in the middle layer samples of sampling sites 1, 2, 3, and 4, where no >200 μm MPs particles were detected.
The MPs in the soils of the study area were mainly fragments, particles, and fibers (Fig. 6b). Particle-type MPs had the highest abundance in the study area soil, followed by fragments, while fibers accounted for the lowest proportion. The abundance of particle-type MPs ranged from 650 to 16300 ind kg−1, with an average abundance of 3915 ind kg−1. The highest abundance of particle-type MPs was found in the bottom layer soil samples of sampling site 10, while the lowest abundance was detected in the topsoil of sampling site 1. Fragment-type MPs had an abundance ranging from 600 to 9600 ind kg−1, with an average abundance of 3692 ind kg−1. The highest abundance of fragment-type MPs was found in the bottom layer soil samples of sampling site 10, while the lowest abundance was detected in the topsoil of sampling site 1. Fiber-type MPs had an abundance ranging from 50 to 3750 ind kg−1, with an average abundance of 855 ind kg−1. The highest abundance of fiber-type MPs was found in the topsoil of sampling site 3, while the lowest abundance was detected in the bottom layer soil samples of sampling site 1.
The MPs in the soils of the study area mainly consisted of PVC, PET, PE, polyamide (PA), and PP (Fig. 6c). MPs consisting of PVC had the highest abundance, followed by PE and PA, while MPs consisting of PP had the lowest abundance. The abundance of PVC MPs ranged from 850 to 15450 ind kg−1, with an average abundance of 5782 ind kg−1. The highest abundance of PVC MPs was found in the topsoil samples of sampling site 8, while the lowest abundance was detected in the topsoil of sampling site 3. The abundance of PE MPs ranged from 150 to 16200 ind kg−1, with an average abundance of 1262 ind kg−1. The highest abundance of PE MPs was found in the topsoil samples of sampling site 3, while the lowest abundance was detected in the bottom layer soil of sampling site 7. The abundance of PA MPs ranged from 0 to 4650 ind kg−1, with an average abundance of 825 ind kg−1. The highest abundance of PA MPs was found in the bottom layer soil samples of sampling site 10, while the lowest abundance was detected in the bottom layer soil of sampling site 7. The abundance of PET MPs ranged from 0 to 3900 ind kg−1, with an average abundance of 320 ind kg−1. The highest abundance of PET MPs was found in the topsoil samples of sampling site 3, while the lowest abundance was detected in the soil of sampling sites 2 and 4. The abundance of PP MPs ranged from 0 to 2100 ind kg−1, with an average abundance of 273 ind kg−1.
Fig. 7 The relationship of MPs abundance to soil physicochemical properties ((a) pH, salinity content, (b) bulk density, moisture content). |
From the Fig. 7b, it can be observed that the salinity content in the soil of the study area ranged from 0 to 0.1, with an average of 0.04. Site 3 had the highest salinity content, while site 2 had the lowest, with the salinity content at the other sites generally ranging from 0.01 to 0.06, indicating relatively low salinity levels. There was considerable spatial variation in the salinity distribution. The salinity of the surface soil ranged from 0 to 0.1, with site 2 having the lowest salinity and site 3 having the highest. In mid-level soil, the salinity content ranged from 0 to 0.1, with sites 2 and 9 having the lowest salinity and site 3 the highest. In bottom-level soil, the salinity content ranged from 0 to 0.1, with site 9 having the lowest salinity and sites 3 and 5 the highest. The vertical distribution of the salinity content at each site indicated relatively minor differences, with no apparent fluctuations from the surface to the bottom layers.
The soil moisture content in the study area ranged from 6.2 to 100%, with an average of 44.4%, indicating significant spatial variability (p < 0.001). The moisture content of surface soil ranged from 9.8 to 100%, with an average of 42.1%. Site 2 had the lowest soil moisture content, while site 3 had the highest. The moisture content of the middle soil layer ranged from 7 to 100%, with an average of 45.1%. Site 9 had the lowest moisture content, while site 3 had the highest. The moisture content of the bottom soil layer ranged from 6.2 to 100%, with an average of 45.9%. Site 9 had the lowest moisture content, while site 3 has the highest. There was obvious overall variation in the longitudinal distribution of the moisture content at each site, with sites 5 and 7 having the largest longitudinal variations. Site 3 was close to a fishpond and had recently been subjected to heavy water irrigation. It had a moisture content of 100% across all three layers.
The distribution of soil water-stable aggregates at each sampling site is shown in Fig. 8. The particle size distribution of soil water-stable aggregates at the various sampling sites revealed that aggregates smaller than 0.25 mm, including microaggregates and fine silt and clay aggregates, constituted the majority (66%) of particles, while larger aggregates (>0.25 mm) accounted for the other 34%. In the surface soil aggregates at each site, microaggregates were dominant in sites 1, 2, 3, 4, and 9, with larger aggregates being less prevalent. The remaining sites were primarily characterized by larger aggregates. In the middle soil layer, microaggregates were predominant in sites 1, 2, 3, 4, 7, 9, and 10, while larger aggregates were less common. In the bottom soil layer, sites 5, 6, and 8 were characterized by a higher proportion of larger aggregates, while the rest of the samples were dominated by microaggregates.
The distribution of soil microaggregates (<0.25 mm) at each sampling site is shown in Fig. 8. The particle size distribution of soil microaggregates mainly consisted of microaggregates and fine, clayey microaggregates (<0.25 mm), accounting for 66% of the total, while larger microaggregates (>0.25 mm) accounted for the remaining 34%. In the surface soil, microaggregates dominated in sites 1, 2, 3, 4, and 9, with a smaller proportion of larger aggregates. The other sites were primarily composed of larger aggregates. In the middle soil layer, sites 1, 2, 3, 4, 7, 9, and 10 had a higher proportion of microaggregates, with smaller proportions of larger aggregates. In the bottom soil layer, sites 5, 6, and 8 had a higher proportion of larger aggregates, while the remaining sites were predominantly composed of microaggregates.
The composition of soil macroaggregates (>0.25 mm) at each site was primarily dominated by particles in the size range of 2–0.25 mm, followed by particles larger than 2 mm. In sites 9 and 10, the proportion of particles larger than 2 mm in the soil macroaggregate structure was higher than that of particles in the 2–0.25 mm range. In the surface soil, the microaggregates in sites 3 and 8 were mainly composed of macroaggregates in the 2–0.25 mm range, while those in sites 4, 7, and 10 were primarily composed of macroaggregates larger than 2 mm, with a uniform distribution of both particle sizes in other sites. In the middle soil layer, the microaggregates in sites 3, 5, 8, 9, and 10 were mainly composed of macroaggregates in the 2–0.25 mm range, while those in sites 6 and 7 were predominantly composed of macroaggregates larger than 2 mm, with a uniform distribution of both particle sizes in the other sites. In the bottom soil layer, the microaggregates in sites 1, 4, 5, 8, and 9 were primarily composed of macroaggregates in the 2–0.25 mm range, with a uniform distribution of both particle sizes in the other sites.
Land use | Location | Depth (cm) | Abundance (ind kg−1) | References |
---|---|---|---|---|
Vegetable plot | Wuhan | 0–5 | 320–12560 | 59 |
Shanghai | 0–6 | 62.5 ± 12.97–78 ± 12.91 | 60 | |
Wuhan | 0–5 | 43000–62000 | 61 | |
South Korea | 0–5 | 1880 ± 1563 | 62 | |
Yunnan | 0–10 | 8180–18100 | 63 | |
Cropland | Shanxi | 0–10 | 1430–3410 | 13 |
Guangdong | 0–20 | 9450–9520 | 11 | |
Yunnan | 0–30 | 900–40800 | 64 | |
Chile | 0–25 | 600–10400 | 65 | |
Forest | Wuhan | 0–5 | 96000–690000 | 61 |
Nanyang | 0–40 | 654–15161 | 66 | |
South Korea | 0–5 | 371–3443 | 67 | |
Shanxi | 0–10 | 1360–4960 | 68 | |
Yunnan | 0–10 | 7100–42960 | 63 |
Compared with the other sampling sites (Fig. 5) sites 1 and 2 had a lower MPs abundance. Both of these sites were vegetable fields, situated far from roads, indicating a minimal influence from human activity. Additionally, the sampling occurred in winter (December), suggesting the prolonged abandonment of these plots due to seasonal effects, with no traces of plastic mulching. In contrast to the other sites, the source of MPs pollution at these sites was likely atmospheric deposition.
Sites 3 and 10 had the highest MP abundance and were located in farmland and forest land, respectively. These sites were adjacent to roads, indicating a larger influence from human activity and an abundance of MPs sources compared to sites 1 and 2. Site 3 consisted mainly of farmland, largely corn, while site 2 in the forest was primarily populated with poplar trees. Both sites had low planting densities and tall mature crops, exposing a large area of the soil surface and intercepting exogenous MPs more effectively. Wheat was cultivated in sites 5 and 7 with a high planting density and short plants, resulting in a lower MPs abundance compared to site 3 despite the similar land use.
Sites 4 and 8 were distant from roads. Fruit trees were grown at both sites with ground coverings of plastic film and straw, resulting in a poor external MPs supply and lower MPs abundance compared to site 2, which was also forest land.
In the study area, vegetable fields are characterized by small plot sizes, remote locations, and a strong seasonality in planting, with minimal use of plastic film. Consequently, the main source of MPs supply is primarily atmospheric deposition, resulting in a lower MPs abundance compared to other land use types. In forest land, the taller plants facilitate the interception of exogenous MPs, while the lower planting densities imply larger exposed surface areas, leading to a higher MPs abundance. Similarly, MPs abundance in farmland is influenced by the crop types. Compared to wheat, the taller corn plants with a lower planting density are more conducive to the accumulation of exogenous MPs. Additionally, based on field observations, there was no widespread use of plastic film coverings, indicating that exogenous MPs were the main source of soil MPs in the study area.
According to the test results of 30 soil samples in the study area, the proportion of granular MPs is relatively high, which may be due to the presence of several large building ruins in the study area, as well as the viaduct project that has been under construction for two years and the supporting accommodation facilities. This has produced a large number of foam particles, plastic microbeads of toiletries, and mechanical parts in the study area, all of which are potential sources of granular MPs.
The main factor influencing bulk density was soil type. There were three types of soil in the study area: sandy soil (samples 1, 7, 9), loamy soil (samples 2, 4, 5, 6), and silty sand soil (samples 3, 8, 10). Among the different soil types, the higher the abundance of MPs, the lower the bulk density. This was more pronounced when PVC and PP, which have relatively low densities, constituted the highest proportion of MPs. This indicates that the presence of MPs can cause changes in soil bulk density. Compared to soil, MPs have a lower density. Therefore, as the abundance of MPs in soil increased, the mass of soil per unit volume decreased accordingly. The influence of MPs on soil bulk density was therefore mainly reflected in both the abundance and types of MPs present in the soil.
Compared to the other two soil types, the upper layer of the silty sand soil had obvious higher MPs concentrations than the middle and lower layers. This indicates that silty sand soil was better able to retain MPs than the other soil types, and its higher soil density implied a lower soil porosity. Additionally, in silty sand soil, the distribution of MPs revealed a larger proportion of MPs with sizes below 50 μm, while those above 200 μm constituted a smaller proportion. This confirmed that it was difficult for larger MPs to enter soils characterized by smaller pores and lower frequencies of anthropogenic disturbance (such as tillage, sowing, and fertilization). The proportion of fibers in sandy soil was higher than in the other two soil types. Yan X. and O'Connor D. suggested that the high soil porosity and increased soil fissures resulting from alternating wet and dry conditions in sandy soils facilitated the downward movement of MPs from the soil surface.77,78 Considering that the sandy soil sampling sites in the study area were located within vegetable fields, where human activity is relatively low, the higher abundance of fiber-type MPs in this region was likely due to differences in soil physical properties.
The abundance of fiber-type MPs was positively correlated with soil bulk density, with higher levels observed in locations with a greater soil bulk density. This may be because fiber-type MPs typically have longer longitudinal lengths and smaller transverse diameters, facilitating their downward movement in soils with small pore sizes.79 In contrast, fragment and film-type MPs do not possess these characteristics, making it difficult for them to migrate downward in soils with high bulk density.
The correlation between soil moisture content and MPs abundance was not clear, which may be due to the greater influence of environmental factors than MPs on soil moisture content. However, in terms of the vertical distribution of moisture, locations with a higher MPs abundance in the middle and lower layers exhibited a smaller decrease in soil moisture content compared to the upper layers. This phenomenon was more pronounced in sandy and loamy soils, suggesting an impact of MPs on soil water retention. Luo et al. reported that the vertical distribution of the soil moisture content in the surface soil (0–30 cm depth) generally displayed a decreasing trend,80 while Chen et al. found that the distribution of the soil moisture content before and after irrigation (0–30 cm depth) did not change obviously, and generally followed the trend of an increasing moisture content with depth.81 However, in this study, the soil moisture content at sites 4, 7, and 9 displayed a decreasing trend with increasing soil depth. These three sites consisted of sandy and loamy soils, with the MPs contents in the upper layer ranging from 11575 to 14475 ind kg−1. There was an obvious decrease in MPs abundance with increasing depth. This might be because the MPs affected the soil structure, leading to increased water retention in the upper soil layers. Due to the reduced water supply, the moisture content subsequently decreased in the middle and lower soil layers.
Sandy soil has a loose structure and is usually dominated by microaggregates. In this study, it was found that the proportion of aggregates in sandy soil increased with the increase in MPs abundance. Silty soil, due to its small particle size and weak cohesive forces, was less affected by MPs abundance. Several locations with a high proportion of soil aggregates were found to contain a large amount of fiber-type MPs and large fragment-type MPs. Small-sized MPs can also adhere to soil particles, as evidenced by the increase in the proportion of microaggregates in the size range of 0.25–0.05 mm with the increase in the fraction of small MPs. The proportion of fragment-type MPs was larger than the proportion of particle-type MPs, resulting in a noticeable increase in the proportion of microaggregates in the size range of 0.25–0.05 mm. This may be because fragment-type and fiber-type MPs have a larger specific surface area and contain more surface charge sites, which can effectively stabilize the soil structure.82
The study area is 35 km2. According to the collected data, there are 6 soil points in the study area, with a monitoring density of about 0.17 ind km−2; there are 8 groundwater monitoring wells, with a monitoring well density of about 0.23 ind km−2. The density of monitoring wells is relatively large, basically covering the entire study area. At the same time, according to the test results of soil samples in the study area, the soil in the study area has been polluted by MPs to varying degrees. Therefore, the arithmetic mean of pH calculated at each monitoring point in the study area using the collected data can basically represent the overall level of soil and groundwater pH in the study area. For the North China Plain, the frequency of mulch film use and irrigation water volume of cultivated land before 2000 were both low, and the degree of development and utilization of groundwater and land resources was low. The changes in soil, precipitation and groundwater pH during this period can basically represent the changes in pH under the natural state of the North China Plain. Therefore, this study selected the average pH calculation value before 2000 as the background value of pH, and used it as a basis to evaluate the impact of MPs brought by human activities after 2000 on the soil in the study area.
According to calculations, the background pH value of soil in the study area is 8.5, the background pH value of precipitation is 6.31, and the background pH value of groundwater is 7.52. The data revealed that from 2000 to 2020, the soil pH in the study area changed from 8.5 to 8.55, the precipitation pH changed from 6.31 to 7.46, and the groundwater pH changed from 7.52 to 7.38 (Fig. 9). There were distinct differences in the trend around 2015. From 2000 to 2015, soil and precipitation pH displayed a decreasing trend, while groundwater pH displayed an increasing trend. From 2015 to 2020, the soil, precipitation, and groundwater pH all displayed a decreasing trend.
The pH of shallow groundwater in the study area has remained relatively stable since 1980, while the pH of precipitation has increased by approximately 1.15 over the same period, exhibiting a trend similar to that of soil pH. The main irrigation water sources for farmland in the study area are atmospheric precipitation and groundwater. Therefore, the increase in soil pH may be related to the increase in precipitation pH. Using the average soil pH value in 2015 as a reference, the range of the pH variation at each sampling site was calculated. In conjunction with the distribution of MPs, it was observed that the two sites with large variations (sites 7 and 9) were located in farmland and vegetable fields with sandy soils, respectively. These sites also had relatively high MPs concentrations of 7450 and 8608 ind kg−1, respectively. This suggests that in areas with a loose soil structure, high MPs concentrations may have an influence on the quality of the water supply source, subsequently affecting the soil pH. This could be attributed to the increase in the concentration of small MPs particles, mainly fragment-type, which tend to loosen the soil structure in sandy soil-dominated farmland. This assists with precipitation infiltration, increases the soil water retention capacity by creating larger pores, enhances the contact between precipitation and soil, and consequently leads to a greater impact of external precipitation on soil chemical properties.
The range of soil pH variation in forested areas was relatively low compared to the variation in 2015. Despite the relatively high abundance of MPs in forested areas within the research area, the predominant soil type was sandy loam, characterized by a compact soil structure, limiting the impact of MPs on the soil structure.
The correlation between soil salinity and MPs abundance in the research area was weak. Li et al. suggested that the primary factors influencing soil salinity are meteorological conditions (such as precipitation and evaporation) and human activities (such as vegetation coverage), which could lead to changes in the soil moisture content, thereby affecting salinity levels.83 This relationship needs to clarified through the continuous monitoring of soil sample indicators such as soil moisture content. Because this study only conducted soil sampling on a single occasion, it was not possible to further investigate the relationship between MPs abundance and soil salinity.
This study also analyzed the changes in soil physical and chemical properties caused by MPs. MPs could impact the soil structure, bulk density, moisture content, and aggregate distribution. In addition to abundance, the type and size of MPs were key factors affecting soil bulk density, while the size and shape of MPs influenced the soil aggregate structure. Additionally, MPs altered the soil water retention, thereby affecting the longitudinal distribution of soil moisture. In contrast, MPs had less noticeable effects on soil pH and salinity than the soil physical properties. MPs had a more noticeable impact on the pH of sandy and loamy soils, primarily through altering soil porosity and water retention. The main change was an increase in the area and duration of contact between soil media and external water sources, thereby altering soil pH. The study did not conclusively determine the influence of MPs on soil salinity. Future studies should involve the continuous sampling and testing of the MPs abundance and soil physicochemical properties in the study area to further investigate the relationship between MPs and soil chemical properties.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d4em00242c |
This journal is © The Royal Society of Chemistry 2024 |