Impacts of Large-Scale Land-Use Change on the Uptake of Polycyclic Aromatic Hydrocarbons in the Artificial Three Northern Regions Shelter Forest Across Northern China

Spatiotemporal patterns and deposition of organophosphate esters (OPEs) in air, foliage and litter in a subtropical forest of South China

Recent studies revealed the un-negligible impact of airborne organophosphate esters (OPEs) on phosphorus (P)-limited ecosystems. Subtropical forests, the global prevalence P-limited ecosystems, contain canopy structures that can effectively sequester OPEs from the atmosphere. However, little is known about the behavior and fate of OPEs in subtropical forest ecosystem, and the impact on the P cycling in this ecosystem. OPE concentrations in the understory air (at two heights), foliage, and litterfall were investigated in a subtropical forest in southern China. The median ∑OPE concentrations were 3149 and 2489 pg/m3 in the upper and bottom air, respectively. Foliage exhibited higher ∑OPE concentrations (median = 386 ng/g dry weight (dw)) compared to litter (median = 267 ng/g dw). The air OPE concentrations were ordered by broadleaved forest > mixed forest > coniferous forest, which corresponds to the results of canopy coverage or leaf area index. The spatial variation of OPEs in foliage and litter was likely caused by the leaf surface functional traits. Higher OPE concentrations were found in the wet season for understory air while in the dry season for foliage and litter, which were attributed to the changes in emission sources and meteorological conditions, respectively. The inverse temporal variation suggests the un-equilibrium partitioning of OPEs between leaf and air. The OPE concentrations during the litter-incubation presented similar temporal trends with those in foliage and litter, indicating the strong interaction of OPEs between the litter layer and the near-soil air, and the efficient buffer of litter layer played in the OPEs partitioning between soil and air. The median OPEs-associated P deposition fluxes through litterfall were 270, 186, and 249 μg P/m2·yr in the broadleaved, mixed, and coniferous forests, respectively. Although the fluxes accounted for approximately 0.2% of the total atmospheric P deposition, their significance to this P-limited ecosystem may not be negligible.

Space-specific flux estimation of atmospheric chemicals from point sources to soil

Predicting chemical flux to soil from industrial point sources accurately at a regional scale has been a significant challenge due to high uncertainty in spatial heterogeneity and quantification. To address this challenge, we developed an innovative approach by combining California Air Resources Board Puff (CALPUFF) and mass balance models, leveraging their complementary strengths in quantitative accuracy and spatial precision. Specifically, CALPUFF was used to predict the polycyclic aromatic hydrocarbons (PAHs) flux to soil due to industrial sources. Additionally, the spatial distribution coefficient of PAHs flux (e.g., si for spatial unit i) was calculated by neural network and combined with the mass balance model to obtain the results of total PAHs fluxes, which were then combined with the results predicted by CALPUFF to effectively estimate the contribution of industrial sources to soil PAHs flux. Taking a petrochemical industry region located in Zhejiang province, China as a case study, results showed the input Phenanthrene (Phe) and Benzo(a)pyrene (BaP) fluxes predicted by CALPUFF were generally lower than those by the mass balance model, with slightly different distribution patterns. CALPUFF results, based on 36 industrial sources, partially represent those of the mass balance model, which includes all sources and pathways. It was suggested that industrial sources contributed 49%-89% and 65%-100% of soil Phe and BaP, respectively across the study area. The average Phe flux from point sources by deposition averaged 2.68 mg m-2∙a-1 in 2021, accounting for approximately 60% of the total Phe flux to soil. The average BaP flux from point sources by deposition averaged 0.0755 mg m-2∙a-1, accounting for only 0.1%-3.65% of the total BaP flux to soil. Thereby, our approach fills up a gap between the relevance to point sources and the accuracy of deposition quantification in estimating chemical flux from specific point sources to soil at a regional scale.

Polycyclic aromatic hydrocarbons (PAHs) in air, foliage, and litter in a subtropical forest: Spatioseasonal variations, partitioning, and litter-PAH degradation

Forest canopies play a vital role in scavenging airborne semi-volatile organic compounds. The present study measured polycyclic aromatic hydrocarbons (PAHs) in the understory air (at two heights), foliage, and litterfall in a subtropical rainforest (the Dinghushan mountain) in southern China. ∑₁₇PAH concentrations in the air ranged from 2.75 to 44.0 ng/m³ (mean = 8.91 ng/m³), showing a spatial variation depending on the forest canopy coverage. Vertical distributions of the understory air concentrations also indicated PAH inputs from the above-canopy air. The concentrations of PAHs in fresh litter (with a mean of 261 ± 163 ng/g dry weight (dw)) were slightly lower than those in the foliage (362 ± 291 ng/g dw). Unlike the stable air PAH concentrations for most of the time of the year, the temporal variations of foliage and litter concentrations were remarkable but generally similar. Higher or comparable leaf/litter-air partition coefficients (K_LA) in fresh litter compared with living K_LA in leaves suggest that the forest litter layer is an efficient storage media for PAHs. Degradation of three-ring PAHs in litter under the field conditions follows first-order kinetics (R² = 0.81), while the degradation is moderate for four-ring PAHs and insignificant for five- and six-ring PAHs. The yearly net cumulative deposition of PAHs through forest litterfall in the whole Dinghushan forest area over the sampling year was about 1.1 kg, 46% of the initial deposition (2.4 kg). This spatial variations study provides the results of in-field degradation of litter PAHs and makes a quantitative assessment of the litter deposition of PAHs, deducing their residence dynamics in the litter layer in a subtropical rainforest.

Atmospheric benzo[a]pyrene in the Yangtze River Delta, China: pollution level and lung cancer risk in 2016 and future predictions

The Yangtze River Delta (YRD) has undergone widespread polycyclic aromatic hydrocarbon (PAH) pollution. In this study, we simulated the spatial distribution of atmospheric benzo[a]pyrene (BaP, the most carcinogenic PAH) in the YRD in 2016 and 2030 under different emission scenarios using a 3-D atmospheric transport model and evaluated the lung cancer risks posed by BaP during the study years. The purpose of this study is to suggest targeted policy recommendations for policy-makers to mitigate BaP pollution through numerical simulation. Our results showed that the average BaP concentration in the YRD was 0.30 ng/m3 in 2016; however, a significant spatial variation was observed, with the highest BaP concentration in Shanghai (0.59 ng/m3). The population-weighted incremental lifetime lung cancer risk (PILCR) was 6.67 × 10-6 in 2016, whereas it ranged from 2.70 × 10-6 to 1.05 × 10-5 in 2030 under the five emission scenarios. A higher future population density in the YRD region could increase lung cancer risk. In all scenarios, Shanghai had the highest number of lung cancer cases (range: 208-476). The results suggest that BaP pollution could be effectively improved through the synergistic effect of reducing activity levels and improving technology. Finally, we provide specific suggested pollution control strategies (e.g., accelerating the use of clean energy in rural areas) for atmospheric BaP in the YRD.

Long-term spatiotemporal variation and lung cancer risk of atmospheric polycyclic aromatic hydrocarbons (PAHs) in the Yangtze River Delta, China

The Yangtze River Delta (YRD), which is the most developed region in China, suffers from atmospheric polycyclic aromatic hydrocarbons (PAH) pollution. However, the long-term spatiotemporal variation of atmospheric PAHs and the lung cancer risk caused by PAH exposure in the YRD remain unclear. Herein, we simulated the daily atmospheric concentration of benzo[a]pyrene (BaP, the most carcinogenic PAH) from 2001 to 2016 using an atmospheric transport model. During this period, the atmospheric BaP concentration showed a general trend of first increasing and then decreasing (average BaP concentration = 0.50 ± 0.12 ng/m³) and was highest in 2005 (0.72 ng/m³). Moreover, the BaP concentration in Jiangsu and Shanghai was 5.17- and 4.98-fold higher than that in Zhejiang. BaP pollution was severe in Jiangsu during the winter. The average area proportion of BaP exceeding the national standard in winter in Jiangsu was 69.09%. The population-weighted incremental lifetime cancer risk from 2001 to 2016 ranged 6.67 × 10^-6-1.50 × 10^-5, and the excess lung cancer cases ranged 1054-2130. Compared with 2005, excess lung cancer cases in the YRD decreased by 49.49% in 2016. Reducing BaP pollution in winter in Jiangsu is crucial for reducing lung cancer risk in the YRD.

Signatures of Indian endosulfan usage in China's environment

Endosulfan, as an organochlorine pesticide (OCPs), was widely used in agriculture. As the largest endosulfan user country in the world and adjacent to China, it is interesting to know to what extent the endosulfan usage of India could affect the environment in China. In this study, we established gridded endosulfan usage, atmospheric emission, and soil residue inventories in 2010 based on its application in different crops in China and India. We employed an atmospheric transport model CanMETOP to simulate atmospheric and soil concentrations, as well as dry and wet deposition flux of α- and β-endosulfan. Results were used to assess the signatures of Indian endosulfan usage in the China's environment. In 2010, endosulfan usage, atmospheric emissions, and highest soil residue in China were 3083.9, 1312.7, and 587.5 tonnes, and 3204.8, 1441.4, and 463.7 tonnes in India, respectively. The spatial distribution of modeled atmospheric and soil concentrations, and dry deposition fluxes of endosulfan were in line with its usage but wet deposition fluxes were mainly identified in Southern China and Sichuan basin with heavy rainfall, especially for α-endosulfan. Endosulfan tended to transport from India to Tibetan plateau, Yunnan-Kweichow Plateau, and some provinces in southern China under the Indian Summer Monsoon regime. Due to its stronger volatility, α-endosulfan posed a more significant impact on China's environment via the atmospheric transport from India compared to β-endosulfan. Although rainfall during Indian Summer Monsoon reduced endosulfan levels in the air during its journey from India to China, it was observed that Indian endosulfan usage in 2010 contributed more than 50% of atmospheric concentrations and 30% of soil concentrations of α-endosulfan in some regions in Tibetan plateau.

Per- and polyfluoroalkyl substances (PFAS) in the Three-North Shelter Forest in northern China: First survey on the effects of forests on the behavior of PFAS

Per- and polyfluoroalkyl substances (PFAS) are a family of anthropogenic chemicals, that have attracted increasing attention since the early 2000 s. Although forests have been shown to act as a filter and important sink for nonpolar persistent organic pollutants (POPs), relevant reports on PFAS are lacking. Air, soil, and leaf samples were collected inside and outside the forest from two regions of the Three-North Shelter Forest in northern China between 2017 and 2018. Twenty-seven PFAS were analyzed to study the effect of forest on the transport and fate of PFAS. The average ratios of PFAS in the air outside to inside the forest (Q_air) ranged from 2.83 ± 0.78-10.6 ± 3.1. A significant positive correlation was found between Q_air and the n-octanol-air partition coefficient of individual PFAS (p = 0.041). Higher Q_air values for most ionic PFAS were found in broad-leaved forests than in coniferous forests. Soil samples outside the forests showed higher PFAS levels than those inside. The measured concentrations of 8:2 fluorotelomer alcohol, a volatile neutral PFAS, in leaf samples were two orders of magnitude higher than those estimated using the equilibrium leaf-air partition of nonpolar POPs, indicating that it may not fit the case of PFAS with surface activity.

Organic pollutants in protected plain areas: The occurrence of PAHs, musks, UV-filters, flame retardants and hydrocarbons in woodland soils

Protected woodlands are rare and small portions of the plain territory of northern Italy, where agricultural, industrial and urban activities strongly dominate the landscape. Such natural areas are frequently set on river floodplains and are therefore potentially conditioned by the contamination brought by the surface waters. We investigated the occurrence of multiple categories of organic pollutants, including Polycyclic Aromatic Hydrocarbons (PAHs), Musk fragrances, UV-filters, organophosphorus and novel brominated Flame Retardants (FRs) and Total Petroleum Hydrocarbons (TPH) in woodland soils of eight different protected areas. The samples collected in the floodplains of the Po, Adige and Fratta rivers resulted more contaminated, with levels of PAHs up to 633 ng g^-1. Moreover, these samples for the first time revealed the presence of personal care products, primarily 2-ethylhexyl-4-methoxycinnamate (EHMC) and tonalide (AHTN), in soils of protected woodlands, reaching respectively 3.4 ng g^-1 and 5.0 ng g^-1, together with the occurrence of both organophosphorus and brominated FRs, with total concentrations up to 15 ng g^-1. Higher concentrations of hydrocarbons, with TPH in the range 5-65 μg g^-1, were instead reflecting the inputs of long chain n-alkanes from epicuticular waxes more than petrogenic contamination.

Role of low-latitude forests in modulating forest filter effect on a continental scale: Long-term simulation on PCB-153 in Chinese forests

Forests are important compartments influencing the environmental fate of persistent organic pollutants (POPs). To illustrate the effect of forests on the regional cycle of POPs, a level IV fugacity fate and transport model coupled with a detailed dynamic-forest module was applied to simulate the long-term variations of PCB-153 in China, where forest coverage accounts for approximately one fifth of land area. In the scenarios with forests, atmospheric outflow from China was 69% of that in the scenario without forests due to the enhanced storage in soil, degradation, and leaching. Previous studies regarded high-latitude areas, such as the polar region and boreal forests, as environments capable of reducing mobility of PCB-153, and they act as sinks of POPs. This modeling result suggests that tropical and subtropical forests may also play a similar role despite high temperatures favoring volatilization. Unlike boreal forest, the low-latitude forests may reduce the overall lifetime of PCB-153 in China due to enhanced degradation in warmer and moist soils of the tropical and subtropical area. Given that approximately half of the global forests are located in tropical and subtropical regions, they can be important environments influencing the global geochemical cycle and distribution of POPs, hence deserving more scientific attention by modeling and empirical studies.

Mapping 30 m Fractional Forest Cover over China’s Three-North Region from Landsat-8 Data Using Ensemble Machine Learning Methods

The accurate monitoring of forest cover and its changes are essential for environmental change research, but current satellite products for forest coverage carry many uncertainties. This study used 30-m Landsat-8 data, and aggregated 1-m GaoFen-2 (GF-2) satellite images to construct the training samples and used multiple machine learning algorithms (MLAs) to estimate the fractional forest cover (FFC) in China’s Three North Region (TNR). In this study, multiple MLAs were merged to construct stacked generalization (SG) models based on the idea of SG, and the performances of the MLAs in the FFC estimation were evaluated. The results of the 10-fold cross-validation showed that all non-linear algorithms had a good performance, with an R2 value of greater than 0.8 and a root-mean square error (RMSE) of less than 0.05. In the bagging ensemble, the random forest (RF) (R2 = 0.993, RMSE = 0.020) model performed the best and in the boosting ensemble, the light gradient boosted machine (LGBM) (R2 = 0.992, RMSE = 0.022) performed the best. Although the evaluation index of the RF is slightly better than that of the LGBM, the independent validation results show that the two models have similar performances. The model evaluation results of the independent datasets showed that, in the SG model, the performance of the SG(LGBM) (R2 = 0.991, RMSE = 0.034) was better than that of the single or non-ensemble model. Comparing the FFC estimates of our model with those of existing datasets showed that our model exhibited more forest spatial distribution details and higher accuracy in complex landscapes. Overall, in this study, the method of using high-resolution remote sensing (RS) images to extract samples for FFC estimation is feasible. Our results demonstrate the potential of the ensemble MLAs to map the FFC. The research results also show that among many MALs, the RF algorithm is the most suitable algorithm for estimating FFC, which provides a reference for future research.

Apoplastic and symplastic uptake of phenanthrene in wheat roots

The contamination of agricultural crops by polycyclic aromatic hydrocarbons (PAHs) has drawn considerable attention due to their carcinogenicity, mutagenicity, and toxicity. However, the uptake process of PAHs in plant roots has not been clearly understood. In this work, we first study the radial uptake of phenanthrene in hydroculture wheat roots by vacuum-infiltration-centrifugation method. The concentration-dependent kinetics of apoplastic and symplastic uptake at phenanthrene concentrations of 0-6.72 μM for 4 h can be described with the Langmuir and Michaelis-Menten equations, respectively; whereas, their time-dependent kinetics at 5.60 μM phenanthrene for 36 h follow the Elovich equation. The apoplastic and symplastic uptake increases with temperature of 15-35 °C. The apparent Arrhenius activation energies for apoplastic and symplastic uptake are 77.5 and 9.39 KJ mol^-1, respectively. The symplastic uptake accounts for over 55% of total phenanthrene uptake, suggesting that symplast is the dominant pathway for wheat root phenanthrene uptake. Larger volume of symplast in roots and lower activation energy lead to the greater contribution of symplast to total uptake of phenanthrene. Our results provide not only novel insights into the mechanisms on the uptake of PAHs by plant roots, but also the help to optimize strategies for crop safety and phytoremediation of PAH-contaminated soil/water.

Atmospheric removal of PM2.5 by man-made Three Northern Regions Shelter Forest in Northern China estimated using satellite retrieved PM2.5 concentration

Atmospheric removal of PM_2.5 by the Three Northern Regions Shelter Forest (TNRSF) - the so called Green Great Wall (GGW) in northern China through dry deposition process was estimated using a bulk big-leaf model and a vegetation collection model. Decadal trend of PM_2.5 dry deposition flux from 1999 to 2010 was calculated from modeled dry deposition velocity and air concentration retrieved from the satellite remote sensing. Dry deposition velocities of PM_2.5 calculated using the two deposition models increased in many places of the TNRSF over the last decade due to increasing vegetation coverage of the TNRSF. Both increasing deposition velocity due to forest expansion and PM_2.5 atmospheric level contributed to the increasing deposition flux of PM_2.5. The highest atmospheric deposition flux of PM_2.5 was found in the Central-north region covering Beijing-Tianjin-Hebei area, followed by the Northwestern and the Northeastern regions of the TNRSF. While greater collection of PM_2.5 by vegetation was identified in the Northeastern region of the TNRSF due to higher forest coverage over this region, the most significant incline of the PM_2.5 atmospheric removal due to vegetation collection was discerned in the Central-north region because of the most rapid increase in the vegetation coverage in this region. A total mass of 2.85×10⁷t PM_2.5 was estimated to be removed from the atmosphere through dry deposition process over the TNRSF from 1999 to 2010. The two deposition models simulated similar magnitude and spatial patterns of PM_2.5 dry deposition fluxes. Our results suggest that the TNRSF plays a moderate role in PM_2.5 uptake, but enhances PM_2.5 atmospheric removal by 30% in 2010 than in 1980.