Nothing Special   »   [go: up one dir, main page]

You seem to have javascript disabled. Please note that many of the page functionalities won't work as expected without javascript enabled.
 
 
Sign in to use this feature.

Years

Between: -

Subjects

remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline

Journals

Article Types

Countries / Regions

Search Results (111)

Search Parameters:
Keywords = Songnen Plain

Order results
Result details
Results per page
Select all
Export citation of selected articles as:
25 pages, 7570 KiB  
Article
Evaluating Maize Residue Cover Using Machine Learning and Remote Sensing in the Meadow Soil Region of Northeast China
by Zhengwei Liang, Jia Du, Weilin Yu, Kaizeng Zhuo, Kewen Shao, Weijian Zhang, Cangming Zhang, Jie Qin, Yu Han, Bingrun Sui and Kaishan Song
Remote Sens. 2024, 16(21), 3953; https://doi.org/10.3390/rs16213953 - 23 Oct 2024
Viewed by 608
Abstract
The management of crop residues in farmland is crucial for increasing soil organic matter and reducing soil erosion. Identifying the regional extent of crop residue cover (CRC) is vital for implementing conservation tillage and formulating agricultural subsidy policies. The Google Earth Engine (GEE) [...] Read more.
The management of crop residues in farmland is crucial for increasing soil organic matter and reducing soil erosion. Identifying the regional extent of crop residue cover (CRC) is vital for implementing conservation tillage and formulating agricultural subsidy policies. The Google Earth Engine (GEE) and remote sensing images from 2019 to 2023 were used to obtain spectral characteristics before the maize seedling stage in Northeast China, followed by constructing the CRC estimation models using machine learning algorithms. To avoid the impact of multicollinearity among data, three machine learning algorithms—ridge regression (RR), partial least squares regression (PLSR), and least absolute shrinkage and selection operator (LASSO)—were employed. By comparing the accuracy of these methods, the most accurate model was determined and applied to subsequent CRC estimation. Based on the estimated CRC and Conservation Technology Information Center definitions of tillage practices, the conservation tillage mapping was completed, and the spatiotemporal distribution characteristics were thoroughly analyzed. The following findings were demonstrated: (1) the PLSR-based model outperformed RR (Pearson’s correlation coefficient (r) = 0.8875, R2 = 0.7877, RMSE = 6.99%) and LASSO (r = 0.8903, R2 = 0.7926, RMSE = 6.88%) with higher accuracy (r = 0.9264, R2 = 0.8582, RMSE = 4.93%). (2) Over the five years, the average no-tillage (NT) proportion in the study area was 15.9%, reduced tillage (RT) was 17.8%, and conventional tillage (CT) was 66.3%. In 2020 and 2022, NT rates were significantly higher at 27.5% and 15.5%, while RT were 15.7% and 30.0%, respectively. (3) Compared to the Sanjiang and Liaohe Plains (RT = 1907 km2 and 1336 km2, and NT = 559 km2 and 585 km2, respectively), the Songnen Plain exhibited higher conservation tillage rates (where RT was 3791 km2 and NT was 1265 km2). This provides crucial scientific evidence for the management and planning of conservation tillage, thereby optimizing farmland production planning, enhancing production efficiency, and promoting the development of sustainable agricultural production systems. Full article
Show Figures

Figure 1

Figure 1
<p>(<b>a</b>) Geographic location of Northeast China. (<b>b</b>) Location of the three plains in Northeast China. (<b>c</b>) Distribution of sample points and soil types.</p>
Full article ">Figure 2
<p>(<b>a</b>) The five-point sampling method. (<b>b</b>) Field sampling data. (<b>c</b>) Classification result.</p>
Full article ">Figure 3
<p>(<b>a</b>–<b>e</b>) Time windows in Northeast China from spring 2019 to spring 2023, respectively.</p>
Full article ">Figure 4
<p>Workflow diagram for the mapping of maize tillage practices.</p>
Full article ">Figure 5
<p>The correlations between the MRC and (<b>a</b>) NDTI; (<b>b</b>) STI.</p>
Full article ">Figure 6
<p>The learning curves that present both R<sup>2</sup> and RMSE for the three models: (<b>a</b>) RMSE for RR learning curve; (<b>b</b>) R<sup>2</sup> for RR learning curve; (<b>c</b>) RMSE for LASSO learning curve; (<b>d</b>) R<sup>2</sup> for LASSO learning curve; (<b>e</b>) RMSE for PLSR learning curve; (<b>f</b>) R<sup>2</sup> for PLSR learning curve.</p>
Full article ">Figure 7
<p>The relationships between the measured and predicted MRC: (<b>a</b>) RR; (<b>b</b>) LASSO; (<b>c</b>)PLSR.</p>
Full article ">Figure 8
<p>The spatial distribution of MRC across the study area: (<b>a</b>–<b>e</b>) for spring 2019 to spring 2023. Note: AS = Anshan; BC = Baicheng; BS = Baishan; BX = Benxi; CY = Chaoyang; CF = Chifeng; DL = Dalian; DQ = Daqing; DXAL = Daxinganling; DD = Dandong; FS = Fushun; FX = Fuxin; HEB = Haerbin; HG = Hegang; HH = Heihe; HLBE = Hulunbeier; HLD = Huludao; JX = Jixi; JL = Jilin; JMS = Jiamusi; JZ = Jinzhou; LYL = Liaoyang; LYJ = Liaoyuan; MDJ = Mudanjiang; PJ = Panjin; QTH = Qitaihe; QQHE = Qiqihaer; SYL = Shenyang; SYS = Shuangyashan; SP = Siping; SYJ = Songyuan; SH = Suihua; TLH = Tieling; TH = Tonghua; TLN = Tongliao; XA = Xingan; YB = Yanbian; YC = Yichun; YK = Yingkou; CC = Changchun.</p>
Full article ">Figure 9
<p>(<b>a</b>–<b>e</b>) Spatial distribution of tillage methods in Northeast China from spring 2019 to spring 2023, respectively. (<b>f</b>) Proportions of the different tillage methods used in Northeast China.</p>
Full article ">Figure 10
<p>Spatial distribution of tillage methods in the (<b>a</b>) Songnen Plain, (<b>b</b>) Sanjiang Plain, and (<b>c</b>) Liaohe Plain in spring 2023. (<b>d</b>) Area under the different tillage methods in the three plains.</p>
Full article ">
18 pages, 2901 KiB  
Article
Comparative Study of Back-Propagation Artificial Neural Network Models for Predicting Salinity Parameters Based on Spectroscopy Under Different Surface Conditions of Soda Saline–Alkali Soils
by Yating Jing, Xuelin You, Mingxuan Lu, Zhuopeng Zhang, Xiaozhen Liu and Jianhua Ren
Agronomy 2024, 14(10), 2407; https://doi.org/10.3390/agronomy14102407 - 17 Oct 2024
Viewed by 654
Abstract
Soil salinization typically exerts a highly negative influence on soil productivity, crop yields, and ecosystem balance. As a typical region afflicted by soil salinization, the soda saline–alkali soils in the Songnen Plain of China demonstrate a clear cracking phenomena. Nevertheless, the overall spectral [...] Read more.
Soil salinization typically exerts a highly negative influence on soil productivity, crop yields, and ecosystem balance. As a typical region afflicted by soil salinization, the soda saline–alkali soils in the Songnen Plain of China demonstrate a clear cracking phenomena. Nevertheless, the overall spectral response to the cracked soil surface has scarcely been studied. This study intends to study the impact of salt parameters on the soil cracking process and enhance the spectral measurement method used for cracked salt-affected soil. To accomplish this goal, a controlled desiccation cracking experiment was carried out on saline soil samples. A gray-level co-occurrence matrix (GLCM) was calculated for the contrast (CON) texture feature to measure the extent of cracking in the dried soil samples. Additionally, spectroscopy measurements were conducted under different surface conditions. Principal component analysis (PCA) was subsequently performed to downscale the spectral data for band integration. Subsequently, the prediction accuracy of back-propagation artificial neural network (BP-ANN) models developed from the principal components of spectral reflectance was compared for different salt parameters. The results reveal that salt content is the dominant factor determining the cracking process in salt-affected soils, and that cracked soil samples had the highest model prediction accuracy for different salt parameters rather than uncracked blocks and 2 mm comparison soil samples. Furthermore, BP-ANN prediction models combining spectral response and CON were further developed, which can significantly enhance the prediction accuracy of different salt parameters with R2 values of 0.93, 0.91, and 0.74 and a ratio of prediction deviation (RPD) of 3.68, 3.26, and 1.72 for soil salinity, electrical conductivity (EC), and pH, respectively. These findings provide valuable insights into the cracking mechanism in salt-affected soils, thereby advancing the field of hyperspectral remote sensing for monitoring soil salinization. Furthermore, this study also aids in enhancing the design of spectral measurements for saline–alkali soils and is also helpful for local soil remediation with supporting data. Full article
(This article belongs to the Special Issue Crop Improvement and Cultivation in Saline-Alkali Soils)
Show Figures

Figure 1

Figure 1
<p>Study area and distribution of sampling points.</p>
Full article ">Figure 2
<p>Measurements of some soil salinity parameters in this study: (<b>a</b>) soil suspension; (<b>b</b>) pH measurement; (<b>c</b>) EC measurement.</p>
Full article ">Figure 3
<p>The pre-processing process of crack image standardization. (<b>a</b>) Standard photograph of cracked soil sample; (<b>b</b>) calibration plate image; (<b>c</b>) colorful crack image; (<b>d</b>) grayscale crack image; (<b>e</b>) binary crack image; (<b>f</b>) inversion of binary crack image.</p>
Full article ">Figure 4
<p>Measurement process of spectral reflectance of soil samples under different surface conditions. (<b>a</b>) cracked sample as a whole; (<b>b</b>) local non-cracked block area; (<b>c</b>) comparison sample with a particle size of 2 mm.</p>
Full article ">Figure 5
<p>A simple schematic diagram of the BP-ANN.</p>
Full article ">Figure 6
<p>Reflectance curves under different surface conditions: (<b>a</b>) 2 mm comparison soil samples; (<b>b</b>) uncracked blocks; (<b>c</b>) overall cracked soil samples; (<b>d</b>) reflectance of a typical soil sample.</p>
Full article ">Figure 7
<p>Principal component scores under different surface states: (<b>a</b>) 2 mm comparison soil samples; (<b>b</b>) uncracked blocks; (<b>c</b>) overall cracked soil samples.</p>
Full article ">Figure 8
<p>Correlation coefficient curve between spectral reflectance and salt parameters: (<b>a</b>) 2 mm comparison soil sample; (<b>b</b>) uncracked blocks; (<b>c</b>) cracked soil samples.</p>
Full article ">Figure 9
<p>Scatter points between measured and predicted salt parameters. (<b>a1</b>–<b>a4</b>) Total salinity; (<b>b1</b>–<b>b4</b>) EC; (<b>c1</b>–<b>c4</b>) pH; (<b>a1</b>,<b>b1</b>,<b>c1</b>) the 2 mm comparison soil samples; (<b>a2</b>,<b>b2</b>,<b>c2</b>) uncracked blocks; (<b>a3</b>,<b>b3</b>,<b>c3</b>) cracked soil samples; (<b>a4</b>,<b>b4</b>,<b>c4</b>) combined with CON.</p>
Full article ">
15 pages, 2583 KiB  
Article
Phosphogypsum with Rice Cultivation Driven Saline-Alkali Soil Remediation Alters the Microbial Community Structure
by Guanru Lu, Zhonghui Feng, Yang Xu, Fachun Guan, Yangyang Jin, Guohui Zhang, Jiafeng Hu, Tianhe Yu, Mingming Wang, Miao Liu, Haoyu Yang, Weiqiang Li and Zhengwei Liang
Plants 2024, 13(19), 2818; https://doi.org/10.3390/plants13192818 - 8 Oct 2024
Viewed by 748
Abstract
The improvement of saline-alkali land plays a key role in ensuring food security and promoting agricultural development. Saline soils modifies the response of the soil microbial community, but research is still limited. The effects of applying phosphogypsum with rice cultivation (PRC) on soil [...] Read more.
The improvement of saline-alkali land plays a key role in ensuring food security and promoting agricultural development. Saline soils modifies the response of the soil microbial community, but research is still limited. The effects of applying phosphogypsum with rice cultivation (PRC) on soil physicochemical properties and bacterial community in soda saline-alkali paddy fields in Songnen Plain, China were studied. The results showed that the PRC significantly improved the physicochemical properties of soil, significantly reduced the salinity, increased the utilization efficiency of carbon, nitrogen, and phosphorus, and significantly increased the activities of urease and phosphatase. The activities of urease and phosphatase were significantly correlated with the contents of total organic carbon and total carbon. A redundancy analysis showed that pH, AP, ESP, HCO3, and Na+ were dominant factors in determining the bacterial community structure. The results showed that PRC could improve soil quality and enhance the ecosystem functionality of soda saline-alkali paddy fields by increasing nutrient content, stimulating soil enzyme activity, and regulating bacterial community improvement. After many years of PRC, the soda-alkali soil paddy field still develops continuously and healthily, which will provide a new idea for sustainable land use management and agricultural development. Full article
Show Figures

Figure 1

Figure 1
<p>(<b>A</b>) Urease activity analyses of soil bacterial microorganisms. (<b>B</b>) Invertase activity analyses of soil bacterial microorganisms. (<b>C</b>) Phosphatase activity analyses of soil bacterial microorganisms. CK, P4, P5, and P6 indicate unmodified soda saline-alkali paddy fields and modification for 4, 5, and 6 years, respectively. Different letters indicate significant differences between different treatments (ANOVA, LSD test, <span class="html-italic">p</span> &lt; 0.05). (<b>D</b>) Correlations between soil physicochemical characteristics and soil enzyme activity. The color of the correlation coefficient ranges from dark red (positive correlation) to dark blue (negative correlation). The color intensity is proportional to the correlation coefficient. The size of the grid is proportional to the absolute value of the correlation. The thickness of the connection between nodes and each environmental factor indicates the degree of correlation. The thicker the connection is, the stronger the correlation is. The color of the line between the node and the environment factor represents the significance value. SWC, soil water content; pH, potential of hydrogen; EC, electric conductivity; TOC, total organic carbon; TN, total nitrogen; TP, total phosphorus; AP, available phosphorus; NH<sub>4</sub>-N, nitrate; NO<sub>3</sub>-N, ammonium; ENa, exchangeable sodium ion; ESP, soil alkalinity. *, <span class="html-italic">p</span> &lt; 0.05; **, <span class="html-italic">p</span> &lt; 0.01; ***, <span class="html-italic">p</span> &lt; 0.001.</p>
Full article ">Figure 2
<p>(<b>A</b>) Venn diagram representing the number of shared and unique amplicon sequence variations (ASVs) in the bacterial communities of the different experimental sites. Each ellipse represents a treatment, the overlapping area between ellipses represents shared ASVs between groups, and the number of each block represents the number of ASVs contained within the block. Microbial internal transcribed spacer amplicon sequences were clustered using the “dada2” algorithm to obtain non-monad ASVs. (<b>B</b>) Linear discriminant analysis (LDA) score. Enriched taxa with an LDA score &gt; 4.2 are shown in the histogram. The ordinate is the classification unit with significant differences between groups, and the horizontal coordinate is a bar chart to visually display the LDA analysis logarithm scores of each classification unit. The taxa are sorted by score value size to describe their specificity within the sample grouping. Longer lengths indicate more significant differences in the taxon, and the color of the bar plot indicates the sample group corresponding to the taxon with the highest abundance. (<b>C</b>) Relative abundances of bacteria from the experimental sites at the phylum level. (<b>D</b>) Relative abundances of bacteria from the experimental sites of the first 20 most abundant class. CK, P4, P5, and P6 indicate unmodified soda saline-alkali paddy fields and modification for 4, 5, and 6 years, respectively.</p>
Full article ">Figure 3
<p>(<b>A</b>) Simpson index analyses of soil bacterial microorganisms. (<b>B</b>) Shannon index analyses of soil bacterial microorganisms. (<b>C</b>) Faith_pd index analyses of soil bacterial microorganisms. (<b>D</b>) Pielou_e index analyses of soil bacterial microorganisms. CK, P4, P5, and P6 indicate unmodified soda saline-alkali paddy fields and modification for 4, 5, and 6 years, respectively. The ends of the whiskers represent the minima and maxima, the bottom and top of the box are the first and third quartiles, respectively, and the line inside the box is the median. * indicates significant differences at the 0.05 level (ANOVA, LSD test, <span class="html-italic">p</span> &lt; 0.05).</p>
Full article ">Figure 4
<p>(<b>A</b>) Nonmetric multidimensional scaling plots of bacterial community structures based on Bray–Curtis distance. Each ellipse represents a treatment. Each ellipse represents a treatment. (<b>B</b>) Redundancy analysis of soil physicochemical characteristics and the microbial community structure of bacteria. This figure presents the scores of samples and significantly varied soil physicochemical characteristics on the first two axes. (<b>C</b>) Pearson correlation analysis between the relative abundance of bacterial at the phylum level and environmental factors (soil physicochemical characteristics and soil ions content). The color of the correlation coefficient ranges from dark red (positive correlation) to dark blue (negative correlation). The color intensity is proportional to the correlation coefficient. SWC, soil water content; pH, potential of hydrogen; EC, electric conductivity; TOC, total organic carbon; TN, total nitrogen; TP, total phosphorus; AP, available phosphorus; NH<sub>4</sub>-N, nitrate; NO<sub>3</sub>-N, ammonium; ENa, exchangeable sodium ion; ESP, soil alkalinity. *, <span class="html-italic">p</span> &lt; 0.05; **, <span class="html-italic">p</span> &lt; 0.01.</p>
Full article ">Figure 5
<p>This model illustrates the mechanism of the improvement of phosphogypsum with rice cultivation on soil quality and bacterial community in paddy field. PRC significantly reduced soil salinity (pH, EC) and increased soil fertility (TOC, TN, TP, NH<sub>4</sub>-N, Urease, Phosphatase). In addition, bacterial diversity and community composition also changed significantly, Faith_pd index increased significantly, and the relative abundance of Proteobacteria, Bacteroidetes, Gammaproteobacteria, and Bacteroidia increased. The relative abundance of Actinobacteria, Acidobacteria, Chloroflexi, Alphaproteobacteria, and Anaerolineae decreased. The reduction in salt pressure and the increase in nutrients provide more suitable conditions for bacterial reproduction and habitat and promote the activity of beneficial bacteria. This helps to improve the ability of bacterial decomposition, strengthen material circulation, and improve soil quality. The red arrow up and the blue arrow down represent the increase and decrease in the corresponding parameter. TOC, total organic carbon; TN, total nitrogen; TP, total phosphorus; NH<sub>4</sub>-N, nitrate.</p>
Full article ">
17 pages, 2512 KiB  
Article
The Spatiotemporal Evolution of Extreme Climate Indices in the Songnen Plain and Its Impact on Maize Yield
by Bowen Tang, Fanxiang Meng, Fangli Dong, Hengfei Zhang and Bo Meng
Agronomy 2024, 14(9), 2128; https://doi.org/10.3390/agronomy14092128 - 19 Sep 2024
Viewed by 614
Abstract
Global climate change is intensifying and extreme weather events are occurring frequently, with far-reaching impacts on agricultural production. The Songnen Plain, as an important maize production region in China, faces challenges posed by climate change. This study aims to explore the effects of [...] Read more.
Global climate change is intensifying and extreme weather events are occurring frequently, with far-reaching impacts on agricultural production. The Songnen Plain, as an important maize production region in China, faces challenges posed by climate change. This study aims to explore the effects of climate extremes on maize yield and provide a scientific basis for the adaptation of agriculture to climate change in this region. The study focuses on the spatial and temporal evolution characteristics of climate extremes during the maize reproductive period from 1988 to 2020 in the Songnen Plain and their impacts on maize yield. Daily temperature and precipitation data from 11 meteorological stations were selected and combined with maize yield information to assess the spatial and temporal trends of extreme climate indices using statistical methods such as the moving average and Mann–Kendall (M-K) mutation test. Pearson correlation analysis and a random forest algorithm were also used to quantify the degree of influence of extreme climate on maize yield. The results showed that (1) the extreme heat and humidity indices (TN90p, TX90p, CWD, R95p, R10, and SDII) tended to increase, while the cold indices (TN10p, TX10p) and the drought indices (CDD) showed a decreasing trend, suggesting that the climate of the Songnen Plain region tends to be warmer and more humid. (2) The cold indices in the extreme temperature indices showed a spatial pattern of being higher in the north and lower in the south and lower in the west and higher in the east, while the warm indices were the opposite, and the extreme precipitation indices showed a spatial pattern of being higher in the east and lower in the west. (3). Both maize yield and trend yield showed a significant upward trend. Maize meteorological yield showed a fluctuating downward trend within the range of −1.64~0.79 t/hm2. During the 33 years, there were three climatic abundance years, two climatic failure years, and the rest of the years were normal years. (4) The cold index TN10p and warm indices TN90p and CWD were significantly correlated with maize yield, in which TN90p had the highest degree of positive correlation with yield, and in the comprehensive analysis, the importance of extreme climatic events on maize yield was in the order of TN90p, TN10p, and CWD. This study demonstrates the impact of extreme climate indices on maize yield in the Songnen Plain, providing a scientific basis for local agricultural management and decision-making, helping to formulate response strategies to mitigate the negative effects of extreme climate, ensure food security, and promote sustainable agricultural development. Full article
Show Figures

Figure 1

Figure 1
<p>Study area and distribution of stations.</p>
Full article ">Figure 2
<p>Time variation trend for extreme climate indices during the maize growth period. (<b>a</b>) TN10p; (<b>b</b>) TN90p; (<b>c</b>) TX10p; (<b>d</b>) TX90p; (<b>e</b>) CDD; (<b>f</b>) CWD; (<b>g</b>) R95p; (<b>h</b>) R10; (<b>i</b>) SDII.</p>
Full article ">Figure 2 Cont.
<p>Time variation trend for extreme climate indices during the maize growth period. (<b>a</b>) TN10p; (<b>b</b>) TN90p; (<b>c</b>) TX10p; (<b>d</b>) TX90p; (<b>e</b>) CDD; (<b>f</b>) CWD; (<b>g</b>) R95p; (<b>h</b>) R10; (<b>i</b>) SDII.</p>
Full article ">Figure 3
<p>Trend and mean spatial distribution of extreme climate indices changes during maize growth period from 1988 to 2020. (<b>a</b>) TN10p; (<b>b</b>) TN90p; (<b>c</b>) TX10p; (<b>d</b>) TX90p; (<b>e</b>) CDD; (<b>f</b>) CWD; (<b>g</b>) R95p; (<b>h</b>) R10; (<b>i</b>) SDII.</p>
Full article ">Figure 4
<p>Changes in maize production in the Songnen Plain from 1988 to 2020. (<b>a</b>) yield; (<b>b</b>) trend yield; (<b>c</b>) meteorological yield; (<b>d</b>) relative yield.</p>
Full article ">Figure 5
<p>Correlation coefficients of extreme climate indices during maize growth period in Songnen Plain.</p>
Full article ">Figure 6
<p>Polynomial fitting curve of extreme climate indices and maize yield in Songnen Plain from 1988 to 2020. (<b>a</b>) TN10p; (<b>b</b>) TN90p; (<b>c</b>) TX10p; (<b>d</b>) TX90p; (<b>e</b>) CDD; (<b>f</b>) CWD; (<b>g</b>) R95p; (<b>h</b>) R10; (<b>i</b>) SDII.</p>
Full article ">Figure 7
<p>Results of random forest algorithm weight assessment of extreme climate indices during maize fertility in Songnen Plain.</p>
Full article ">
20 pages, 10669 KiB  
Article
Spatial and Temporal Variations in Soil Organic Matter and Their Influencing Factors in the Songnen and Sanjiang Plains of China (1984–2021)
by Hongju Zhao, Chong Luo, Depiao Kong, Yunfei Yu, Deqiang Zang and Fang Wang
Land 2024, 13(9), 1447; https://doi.org/10.3390/land13091447 - 6 Sep 2024
Viewed by 570
Abstract
Soil organic matter (SOM) is essential for assessing land quality and enhancing soil fertility. Understanding SOM spatial and temporal changes is crucial for sustainable soil management. This study investigates the spatial and temporal variations and influencing factors of SOM content in the Songnen [...] Read more.
Soil organic matter (SOM) is essential for assessing land quality and enhancing soil fertility. Understanding SOM spatial and temporal changes is crucial for sustainable soil management. This study investigates the spatial and temporal variations and influencing factors of SOM content in the Songnen Plain (SNP) and Sanjiang Plain (SJP) of Heilongjiang Province, China, based on high-precision SOC content data (RMSE = 4.84 g/kg−1, R2 = 0.75, RPIQ = 2.43) from 1984 to 2021, with geostatistical analyses and geodetector models. This study aims to quantitatively reveal and compare the long-term spatial and temporal characteristics of SOM changes and their influencing factors across these two plains. The results show that SOM content in both plains has decreased over the past 37 years. In the SNP, the average SOM decreased from 48.61 g/kg to 45.6 g/kg, representing a reduction of 3.01 g/kg, or a 6.10% decrease; SOM decreased spatially from northeast to southwest, covering 63.1% of the area. In the SJP, the average SOM declined from 48.41 g/kg to 44.31 g/kg, a decrease of 4.1 g/kg, or an 8.50% decrease; no pronounced spatial pattern was observed, but the declining area comprises 67.49%. Changing SOM hotspots are concentrated in southern SNP and central and northwestern SJP, showing clear heterogeneity across counties. Geodetector model analysis indicates annual mean temperature as the primary driver of SOM variations in SNP; while elevation is the main driver in SJP, the combined explanatory power of multiple factors surpasses individual ones. There is a positive correlation between SOM and temperature in SNP, and policy protection positively influences SOM in both plains. These findings provide insights into the differential protection of SOM in SNP and SJP. Full article
Show Figures

Figure 1

Figure 1
<p>Location of the study area.</p>
Full article ">Figure 2
<p>Distribution and ranking of factors in the study area, where “<b>1</b>” represents the areas that belong to the farm management, “<b>2</b>” represents the areas not under farm management; “<b>3</b>” represents the dry field areas, and “<b>4</b>” represents the paddy field areas.</p>
Full article ">Figure 3
<p>SOM content obtained by processing SOC content based on the ML-CNN model for 1984–2021, where (<b>a</b>–<b>d</b>) show the SOM content of SNP in 1984–1990, the SOM content of SNP in 2016–2021, the SOM content of SJP in 1984–1990, and the SOM content of SJP in 2016–2021, respectively.</p>
Full article ">Figure 4
<p>Spatial distribution pattern of SOM. (<b>a</b>) The spatial proportion of SOM with different contents (1984–1990); (<b>b</b>) The spatial proportion of SOM with different contents (2016–2021).</p>
Full article ">Figure 5
<p>(<b>a</b>) SNP and SJP SOM changes (1984–2021); (<b>b</b>) the Cv index of SOM in SNP and SJP (1984–2021); (<b>c</b>) the variation in annual SOM in SNP and SJP (1984–2021); and (<b>d</b>) the regional change proportion.</p>
Full article ">Figure 6
<p>Spatial distribution of cold hotspots of SOM change in SNP and SJP counties.</p>
Full article ">Figure 7
<p>Typical areas of SOM changes in SNP and SJP counties, where “1, 2, 3” are Longjiang, Zhaozhou, and Wuchang in SNP; “4, 5, 6” are Fuyuan, Youyi, and Hulin in SJP, respectively.</p>
Full article ">Figure 8
<p>Single factor detection results of drivers for SNP and SJP.</p>
Full article ">Figure 9
<p>Interaction detection results of drivers for SNP and SJP.</p>
Full article ">Figure 10
<p>Differences in dominant factors of spatial and temporal changes in SOM between SNP and SJP.</p>
Full article ">Figure 11
<p>Temporal changes in SOM in the pilot areas carried out in Heilongjiang Province.</p>
Full article ">
16 pages, 5236 KiB  
Article
Effects of Organic Fertilizer and Biochar on Carbon Release and Microbial Communities in Saline–Alkaline Soil
by Pengfei Zhang, Ziwei Jiang, Xiaodong Wu, Nannan Zhang, Jiawei Zhang, Siyuan Zou, Jifu Wang and Shuying Zang
Agronomy 2024, 14(9), 1967; https://doi.org/10.3390/agronomy14091967 - 31 Aug 2024
Viewed by 1097
Abstract
Climate change and aridification have increased the risk of salinization and organic carbon loss in dryland soils. Enrichment using biochar and organic fertilizers has the potential to reduce salt toxicity and soil carbon loss. However, the effects of biochar and organic fertilizers on [...] Read more.
Climate change and aridification have increased the risk of salinization and organic carbon loss in dryland soils. Enrichment using biochar and organic fertilizers has the potential to reduce salt toxicity and soil carbon loss. However, the effects of biochar and organic fertilizers on CO2 and CH4 emissions from saline soils in dryland areas, as well as their microbial mechanisms, remain unelucidated. To clarify these issues, we performed a 5-month incubation experiment on typical soda-type saline soil from the western part of the Songnen Plain using five treatments: control treatment (CK), 5% urea (U), straw + 5% urea (SU), straw + 5% urea + microbial agent (SUH), and straw + 5% urea + biochar (SUB). Compared with the SU treatment, the SUH and SUB treatments reduced cumulative CO2 emissions by 14.85% and 35.19%, respectively. The addition of a microbiological agent to the SU treatment reduced the cumulative CH4 emissions by 19.55%, whereas the addition of biochar to the SU treatment increased the cumulative CH4 emissions by 4.12%. These additions also increased the relative abundances of Proteobacteria, Planctomycetes, and Ascomycota. Overall, the addition of biochar and organic fertilizer promoted CO2 emissions and CH4 uptake. This was mainly attributed to an improved soil gas diffusion rate due to the addition of organic materials and enhanced microbial stress due to soil salinity and alkalinity from the release of alkaline substances under closed-culture conditions. Our findings have positive implications for enhancing carbon storage in saline soils in arid regions. Full article
(This article belongs to the Section Agricultural Biosystem and Biological Engineering)
Show Figures

Figure 1

Figure 1
<p>Cumulative CO<sub>2</sub> emissions (<b>a</b>) and cumulative CH<sub>4</sub> emissions (<b>b</b>) in different treatment groups (n = 3). CK, control treatment; U, urea; SU, straw + urea; SUH, straw + urea + microbial agent; SUB, straw + urea +biochar.</p>
Full article ">Figure 2
<p>Relative abundances of bacteria (<b>a</b>) and fungi (<b>c</b>) at the phylum level (&gt;1%) for all treatment groups on day 150. Nonmetric multidimensional scaling (NMDS) analysis of bacteria (<b>b</b>) and fungi (<b>d</b>). Each graph is grouped and connected based on the samples from each treatment group. CK, control treatment; U, urea; SU, straw + urea; SUH, straw + urea + microbial agent; SUB, straw + urea + biochar.</p>
Full article ">Figure 3
<p>Normalized heatmap analysis of predicted abundances of carbon degradation and CH₄ oxidation functional enzymes derived from soil bacterial (<b>a</b>) and fungal (<b>b</b>) sequencing data following cultivation experiments. CK, control treatment; U, urea; SU, straw + urea; SUH, straw + urea + microbial agent; SUB, straw + urea + biochar. Enzymes include α-amylase (EC 3.2.1.1), glucoamylase (EC 3.2.1.3), α-glucosidase (EC 3.2.1.20), isoamylase (EC 3.2.1.68), glycogen phosphorylase (EC 2.4.1.1), pullulanase (EC 3.2.1.41), cyclodextrin glycosyltransferase (EC 2.4.1.19), exocellobiohydrolase (EC 3.2.1.91), β-glucosidase (BG; EC 3.2.1.21), cellulase (EC 3.2.1.4), xylanase (EC 3.2.1.8), β-mannosidase (EC 3.2.1.25), α-L-arabinosidase (EC 3.2.1.55), β-xylosidase (EC 3.2.1.37), hemicellulase (EC 3.1.1.73), chitinase (EC 3.2.1.14), chitobiase (EC 3.2.1.132), α-N-acetylglucosaminidase (EC 3.2.1.50), particulate methane monooxygenase (pMMO; EC 1.14.18.3), laccase (LA; EC 1.10.3.2), and α-D-glucuronidase (EC 3.2.1.20).</p>
Full article ">Figure 4
<p>Heatmap of Spearman’s correlation analysis between CO<sub>2</sub> emissions (<b>a</b>) and CH₄ emissions (<b>b</b>) with microbial diversity indices and soil physicochemical properties at the phylum-level SAC, soil additives with different characteristics; CCO<sub>2</sub>, cumulative CO<sub>2</sub> emissions; CCH<sub>4</sub>, cumulative CH<sub>4</sub> emissions; CN<sub>2</sub>O, cumulative N<sub>2</sub>O emissions; BOS, bacterial observed species index; FChao1, fungal Chao1 index; FOS, fungal observed species index. * <span class="html-italic">p</span> &lt; 0.05, ** <span class="html-italic">p</span> &lt; 0.01, and *** <span class="html-italic">p</span> &lt; 0.001.</p>
Full article ">Figure 5
<p>Structural equation modeling (SEM) based on the effects of SAC, soil physicochemical properties, and fungal alpha diversity index on cumulative CO<sub>2</sub> emissions (<b>a</b>) and cumulative CH<sub>4</sub> emissions (<b>b</b>) in saline–alkaline soil samples. SEM-based standardized total effect on cumulative CO<sub>2</sub> emissions (<b>c</b>) and cumulative CH<sub>4</sub> emissions (<b>d</b>). Blue and red lines indicate significant positive and negative correlations, respectively (<span class="html-italic">p</span> &lt; 0.05), and dashed lines indicate a potential nonsignificant path. Numbers on the arrows indicate standardized path coefficients (* <span class="html-italic">p</span> &lt; 0.05, ** <span class="html-italic">p</span> &lt; 0.01, and *** <span class="html-italic">p</span> &lt; 0.001). Black double arrows indicate the covariance between the exogenous variables. R<sup>2</sup> denotes the total variance of the dependent variables explained by the model. SAC, soil additives with different characteristics; FC, fungal Chao1 index; FOS, observed fungal species index; FP, fungal Pielou evenness index.</p>
Full article ">
20 pages, 17564 KiB  
Article
Spatiotemporal Dynamics and Evolution of Grain Cropping Patterns in Northeast China: Insights from Remote Sensing and Spatial Overlay Analysis
by Guoming Du, Le Han, Longcheng Yao and Bonoua Faye
Agriculture 2024, 14(9), 1443; https://doi.org/10.3390/agriculture14091443 - 24 Aug 2024
Cited by 1 | Viewed by 950
Abstract
Understanding the spatiotemporal patterns and driving mechanisms of cropping patterns’ evolution tailored to local conditions is crucial for the effective allocation of black soil in northeast China and the advancement of agricultural development. This study utilized the Google Earth Engine platform to extract [...] Read more.
Understanding the spatiotemporal patterns and driving mechanisms of cropping patterns’ evolution tailored to local conditions is crucial for the effective allocation of black soil in northeast China and the advancement of agricultural development. This study utilized the Google Earth Engine platform to extract the spatial distribution data of major grain crops in northeast China for the year 2022. Using crop classification data from 2000 to 2022, the spatial overlay analysis method identified cropping pattern types based on spatial and temporal changes. The primary cropping patterns identified were continuous maize cropping, maize–soybean rotation, mixed cropping, and continuous soybean cropping. Simultaneously, this research constructed three distinct crop periods: Period I (2000–2002), Period II (2010–2012), and Period III (2020–2022). Over three periods, these patterns covered 94.73%, 88.76%, and 86.39% of the area, respectively. The evolution of the dominant cropping pattern from Period I to Period II involved the transition from continuous soybean cropping to continuous maize cropping, while from Period II to Period III, the main shift was from continuous maize cropping to maize–soybean mixed cropping. From a spatial perspective, since Period I, maize has increasingly replaced soybean as the dominant crop, with continuous maize cropping expanding northward and continuous soybean cropping contracting. The maize–soybean rotation area also migrated northward, particularly in the core area of the Songnen Plain, evolving mostly into continuous maize cropping. Maize cropping areas exhibited significant regional characteristics, being densely distributed in the Sanjiang Plain and Liaohe Plain, and along major tributaries in northeast China. Consequently, the interplay of the natural environment, economic policies, and agricultural technologies drove these changes. The findings offer valuable insights for optimizing cropping patterns and developing rotation systems in northeast China. Full article
Show Figures

Figure 1

Figure 1
<p>An overview of the study area.</p>
Full article ">Figure 2
<p>The technological path of data acquisition and processing.</p>
Full article ">Figure 3
<p>Research framework for the spatiotemporal dynamics and evolution of grain cropping patterns in northeast China.</p>
Full article ">Figure 4
<p>Mechanisms of recognizing cropping patterns. MCC, continuous maize cropping; RCC, continuous rice cropping; SCC, continuous soybean cropping; MSR, maize–soybean rotation; MST, maize–soybean 2-year rotation; MMS, maize–maize–soybean 3-year rotation; SSM, soybean–soybean–maize 3-year rotation; OP, other cropping patterns.</p>
Full article ">Figure 5
<p>Evolution of the structure of grain cropping in northeast China, 2000–2002, 2010–2012, and 2020–2022.</p>
Full article ">Figure 6
<p>Spatial layout of the cropping structure for 2000–2002, 2010–2012, and 2020–2022.</p>
Full article ">Figure 7
<p>Distribution map of cropping patterns in northeast China from Period I to Period III.</p>
Full article ">Figure 8
<p>Schematic diagram of cropping patterns’ transitions during Periods I to III.</p>
Full article ">Figure 9
<p>Spatial trajectories of changes in cropping patterns between Periods I and III.</p>
Full article ">Figure 10
<p>The evolutionary drivers of cropping patterns in northeast China.</p>
Full article ">
19 pages, 9669 KiB  
Article
Research on Ground-Based Remote-Sensing Inversion Method for Soil Salinity Information Based on Crack Characteristics and Spectral Response
by Xiaozhen Liu, Zhuopeng Zhang, Mingxuan Lu, Yifan Wang and Jianhua Ren
Agronomy 2024, 14(8), 1837; https://doi.org/10.3390/agronomy14081837 - 20 Aug 2024
Viewed by 567
Abstract
The precise quantification of soil salinity and the spatial distribution are paramount for proficiently managing and remediating salinized soils. This study aims to explore a pioneering methodology for forecasting soil salinity by combining the spectroscopy of soda saline–alkali soil with crack characteristics, thereby [...] Read more.
The precise quantification of soil salinity and the spatial distribution are paramount for proficiently managing and remediating salinized soils. This study aims to explore a pioneering methodology for forecasting soil salinity by combining the spectroscopy of soda saline–alkali soil with crack characteristics, thereby facilitating the ground-based remote-sensing inversion of soil salinity. To attain this objective, a surface cracking experiment was meticulously conducted under controlled indoor conditions for 57 soda saline–alkali soil samples from the Songnen Plain of China. The quantitative parameters for crack characterization, encompassing the length and area of desiccation cracks, together with the contrast texture feature were methodically derived. Furthermore, spectral reflectance of the cracked soil surface was measured. A structural equation model (SEM) was then employed for the estimation of soil salt parameters, including electrical conductivity (EC1:5), Na+, pH, HCO3, CO32−, and the total salinity. The investigation unveiled notable associations between different salt parameters and crack attributes, alongside spectral reflectance measurements (r = 0.52–0.95), yet both clay content and mineralogy had little effect on the cracking process due to its low activity index. In addition, the presence of desiccation cracks accentuated the overall spectral contrast of salt-affected soil samples. The application of SEMs facilitated the concurrent prediction of multiple soil salt parameters alongside the regression analysis for individual salt parameters. Nonetheless, this study confers the advantage of the swift synchronous observation of multiple salt parameters whilst furnishing lucid interpretation and pragmatic utility. This study helps us to explore the mechanism of soil salinity on the surface cracking of soda saline–alkali soil in the Songnen Plain of China, and it also provides an effective solution for quickly and accurately predicting soil salt content using crack characteristics, which also provides a new perspective for the hyperspectral measurement of saline–alkali soils. Full article
Show Figures

Figure 1

Figure 1
<p>Distribution of research area and sampling points.</p>
Full article ">Figure 2
<p>The standardized measurements of photography and spectroscopy of cracked soil samples. (<b>a</b>). Measurement process; (<b>b</b>). Crack pattern; (<b>c</b>). Image of calibration plate; (<b>d</b>). Area of spectral irradiation.</p>
Full article ">Figure 3
<p>The preprocessing process of crack images. (<b>a</b>). Cropped crack image; (<b>b</b>). Grayscale crack image; (<b>c</b>). Binary crack image; (<b>d</b>). Skeletonized crack image.</p>
Full article ">Figure 4
<p>SEM schematic diagram.</p>
Full article ">Figure 5
<p>Spectral characteristics of 11 soil samples with different soil salinity.</p>
Full article ">Figure 6
<p>Correlograms between reflectance and main soil properties of soil samples.</p>
Full article ">Figure 7
<p>The results of fitting measured and predicted values of a single dependent variable by SEM.</p>
Full article ">Figure 8
<p>The results of fitting measured and predicted values of a multiple dependent variables by SEM.</p>
Full article ">
19 pages, 7143 KiB  
Article
Potential Reduction of Spatiotemporal Patterns of Water and Wind Erosion with Conservation Tillage in Northeast China
by Fahui Jiang, Xinhua Peng, Qinglin Li, Yongqi Qian and Zhongbin Zhang
Land 2024, 13(8), 1219; https://doi.org/10.3390/land13081219 - 6 Aug 2024
Viewed by 923
Abstract
Conservational tillage (NT) is widely recognized globally for its efficacy in mitigating soil loss due to wind and water erosion. However, a systematic large-scale estimate of NT’s impact on soil loss reduction in Northeast, China’s primary granary, remains absent. This study aimed to [...] Read more.
Conservational tillage (NT) is widely recognized globally for its efficacy in mitigating soil loss due to wind and water erosion. However, a systematic large-scale estimate of NT’s impact on soil loss reduction in Northeast, China’s primary granary, remains absent. This study aimed to investigate the spatial and temporal variability of soil erosion under NT compared to conventional tillage (CT) in the black soil region and to analyze the underlying mechanisms driving these erosions. The Revised Universal Soil Loss Equation (RUSLE) and the Revised Wind Erosion Equation (RWEQ) models were employed, incorporating previously published plot/watershed data to estimate the potential reduction of water and wind erosion by NT in this region. Results indicated that under CT practices, water- and wind-induced soil losses were widely distributed in the arable land of Northeast China, with intensities of 2603 t km−2 a−1 and 34 t km−2 a−1, respectively. Furthermore, the erosive processes of water and wind erosion were significantly reduced by 56.4% and 91.8%, respectively, under NT practices compared to CT. The highest efficiency in soil conservation using NT was observed in the mountainous regions such as the Changbai Mountains and Greater Khingan Mountains, where water erosion was primarily driven by cropland slopes and wind erosion was driven by the wind speed. Conversely, the largest areas of severe erosion were observed in the Songnen Plain, primarily due to the significant proportion of arable land in this region. In the plain regions, water-induced soil loss was primarily influenced by precipitation, with light and higher levels of erosion occurring more frequently on long gentle slopes (0–3°) than on higher slope areas (3–5°). In the temporal dimension, soil loss induced by water and wind erosion ceased during the winter under both tillage systems due to snow cover and water freezing in the soil combined with the extremely cold climate. Substantial reductions were observed under NT from spring to autumn compared to CT. Ultimately, the temporal and spatial variations of soil loss under CT and NT practices were established from 2010 to 2018 and then projected onto a cropland map of Northeast China. Based on this analysis, NT is recommended as most suitable practice in the southern regions of Northeast China for maintaining soil health and crop yield production, while its suitability decreases in the northern and eastern regions. Full article
(This article belongs to the Topic Slope Erosion Monitoring and Anti-erosion)
Show Figures

Figure 1

Figure 1
<p>The distribution of dryland and the five ecological regions (I–V) in Northeast China.</p>
Full article ">Figure 2
<p>Spatial distribution of factors in the RUSLE model under different tillage practices in 2018, including rainfall erosivity, R (<b>A</b>); soil erodibility, K (<b>B</b>); slope length and steepness, LS (<b>C</b>); crop covering, C (<b>D</b>); soil frozen factor, F (<b>E</b>); and protected effect of conservation tillage, P (<b>F</b>).</p>
Full article ">Figure 3
<p>Spatial distribution of factors in the RWEQ model under different tillage practices in 2018, including wind erosivity, <span class="html-italic">W<sub>f</sub></span> (<b>A</b>); soil wetness, SW (<b>B</b>); snow cover depth, SD (<b>C</b>); soil erodibility factor, K′ (<b>D</b>); soil curst factor, SCF (<b>E</b>); soil roughness, RN (<b>F</b>); soil frozen factor, F′ (<b>G</b>); and straw protection of conservation tillage, P′ (<b>H</b>).</p>
Full article ">Figure 4
<p>Spatial distribution of the water (<b>A</b>–<b>C</b>) and wind (<b>D</b>–<b>F</b>) erosion under different tillage practices in Northeast China’s dryland in 2018. The five ecological regions are the Greater Khingan Mountains (I), the Songnen Plain (II), the Liao River Plain (III), the Changbai Mountains (IV), and the Sanjiang Plain (V).</p>
Full article ">Figure 5
<p>Areas of water (<b>A</b>) and wind (<b>B</b>) erosion under different tillage practices exhibited several grades in Northeast China in 2018. The six erosion grades include tolerable (0–200 t km<sup>−2</sup> a<sup>−1</sup>), slight (200–2500 t km<sup>−2</sup> a<sup>−1</sup>), moderate (2500–5000 t km<sup>−2</sup> a<sup>−1</sup>), severe (5000–8000 t km<sup>−2</sup> a<sup>−1</sup>), very severe (8000–15, 000 t km<sup>−2</sup> a<sup>−1</sup>), and destructive erosion (&gt;15, 000 t km<sup>−2</sup> a<sup>−1</sup>). The five ecological regions are the Greater Khingan Mountains (I), the Songnen Plain (II), the Liao River Plain (III), the Changbai Mountains (IV), and the Sanjiang Plain (V).</p>
Full article ">Figure 6
<p>Annually change in soil loss via water (<b>A</b>) and wind (<b>B</b>) erosion under different tillage practices in Northeast China’s dryland.</p>
Full article ">Figure 7
<p>Monthly change in soil loss via water (<b>A</b>) and wind (<b>B</b>) erosion under different tillage practices in Northeast China’s dryland.</p>
Full article ">Figure 8
<p>The effect of various factors on soil water (<b>A</b>) and wind (<b>B</b>) erosion under different tillage practices in Northeast China’s dryland in 2018.</p>
Full article ">Figure 9
<p>Relationship between soil water erosion and slope under different tillage practices in Northeast China’s dryland in 2018. The slope of the arable land ranged 0–18° in the areas, with 0–5° representing plain areas and 5–18° representing mountainous regions.</p>
Full article ">
26 pages, 15638 KiB  
Article
The Spatiotemporal Correlation between Human Activity Intensity and the Evolution of Ecosystem Service Value in the Songnen Plain, China
by Xinxin Guo, Yang Yang, Yi Zhang, Mohsen Kalantari, Jiali Sun, Weize Sun, Guofeng Guan and Guoming Du
Land 2024, 13(8), 1158; https://doi.org/10.3390/land13081158 - 28 Jul 2024
Viewed by 936
Abstract
For the main grain-producing areas worldwide that balance multi-tasks of grain production, ecological protection, and economic development, quantitatively revealing the correlation between human activity intensity (HAI) and ecosystem service value (ESV) is conducive to formulating adapted ecological protection policies and promoting the coordinated [...] Read more.
For the main grain-producing areas worldwide that balance multi-tasks of grain production, ecological protection, and economic development, quantitatively revealing the correlation between human activity intensity (HAI) and ecosystem service value (ESV) is conducive to formulating adapted ecological protection policies and promoting the coordinated development of the regional economy, society, and ecosystem. In this paper, we took the Songnen Plain of China as an example, employed a modified Equivalent Factor Method (integrating socio-economic data, the normalized difference vegetation index (NDVI), and land use data), and the HAI Assessment Model (based on the data of land use, night-time light, and population spatial distribution) to measure the ESV and HAI in the Songnen Plain of China for the years 1990, 1995, 2000, 2005, 2010, 2015, and 2020. We further applied the standard deviational ellipse method, the coupled coordination degree model, and the bivariate spatial autocorrelation models to reveal the spatiotemporal dynamics and correlation characteristics of ESV and HAI. The results showed the following: (1) Temporally, the ESV declined from 950.96 billion yuan in 1990 to 836.31 billion yuan in 2015, and then increased to 864.60 billion yuan in 2020, with the total loss attributed to the significant decline in the ESV of the natural ecosystem. Spatially, the ESV in the western and northeastern regions was relatively high, with a significant increase in the northeast. (2) HAI showed an upward trend from 1990 to 2020, while the high HAI levels gradually shrank after reaching the peak in 2000. Low HAI levels were mainly distributed in the northeast and southwest, aligning with the ecological space, while high HAI levels were distributed in the middle and southeast. (3) The interaction between ESV and HAI was marked by a negative correlation, transitioning from conflict to coordination. The spatial pattern of HAI and ESV showed H (HAI)-L (ESV) and L-H clustering, with H-H and L-L scattered distributions. This study contributes to providing a framework, methods, and suggestions for the sustainable planning and utilization of land and ecological protection in order to offer scientific references for the Songnen Plain, other major grain-producing areas, and related studies. Full article
(This article belongs to the Section Land, Biodiversity, and Human Wellbeing)
Show Figures

Figure 1

Figure 1
<p>Location of the Songnen Plain.</p>
Full article ">Figure 2
<p>The technical framework of this article.</p>
Full article ">Figure 3
<p>The proportion of the Songnen Plain’s different ecosystem service values in 1990–2020.</p>
Full article ">Figure 4
<p>Value changes of ecosystem service functions during 1990–2020.</p>
Full article ">Figure 5
<p>The spatial distribution of ESV and its changes in the Songnen Plain from 1990 to 2020. (Notes: <b>A</b>–<b>G</b> denotes the ESV of study area. <b>H</b> represents the changing trend of ESV).</p>
Full article ">Figure 6
<p>Standard deviation ellipse of ESV and its center of gravity change in the Songnen plain (1990–2020).</p>
Full article ">Figure 7
<p>The proportion of human activity intensity.</p>
Full article ">Figure 8
<p>Spatial distribution of human activity intensity in the Songnen Plain from 1990 to 2020.</p>
Full article ">Figure 9
<p>The coupling degree and coupling coordination degree between ESV and HAI from 1990 to 2020 in the Songnen Plain.</p>
Full article ">Figure 10
<p>Songnen plain ESV and HAI Moran’s I.</p>
Full article ">Figure 11
<p>LISA clustering map of HAI and ESV (<b>a</b>) and relative proportions of different clustering types to the grid (<b>b</b>).</p>
Full article ">Figure 12
<p>Relationship between human activity and ecosystem service (modified and integrated through [<a href="#B13-land-13-01158" class="html-bibr">13</a>,<a href="#B75-land-13-01158" class="html-bibr">75</a>]).</p>
Full article ">
20 pages, 5896 KiB  
Article
A Comparative Study of Different Dimensionality Reduction Algorithms for Hyperspectral Prediction of Salt Information in Saline–Alkali Soils of Songnen Plain, China
by Kai Li, Haoyun Zhou, Jianhua Ren, Xiaozhen Liu and Zhuopeng Zhang
Agriculture 2024, 14(7), 1200; https://doi.org/10.3390/agriculture14071200 - 21 Jul 2024
Viewed by 975
Abstract
Hyperspectral technology is widely recognized as an effective method for monitoring soil salinity. However, the traditional sieved samples often cannot reflect the true condition of the soil surface. In particular, there is a lack of research on the spectral response of cracked salt-affected [...] Read more.
Hyperspectral technology is widely recognized as an effective method for monitoring soil salinity. However, the traditional sieved samples often cannot reflect the true condition of the soil surface. In particular, there is a lack of research on the spectral response of cracked salt-affected soils despite the common occurrence of cohesive saline soil shrinkage and cracking during water evaporation. To address this research, a laboratory was designed to simulate the desiccation cracking progress of 57 soda saline–alkali soil samples with different salinity levels in the Songnen Plain of China. After completion of the drying process, spectroscopic analysis was conducted on the surface of all the cracked soil samples. Moreover, this study aimed to evaluate the predictive ability of multiple linear regression models (MLR) for four main salt parameters. The hyperspectral reflectance data was analyzed using three different band screening methods, namely random forest (RF), principal component analysis (PCA), and Pearson correlation analysis (R). The findings revealed a significant correlation between desiccation cracking and soil salinity, suggesting that salinity is the primary factor influencing surface cracking of saline–alkali soil in the Songnen Plain. The results of the modeling analysis also indicated that, regardless of the spectral dimensionality reduction method employed, salinity exhibited the highest prediction accuracy for soil salinity, followed by electrical conductivity (EC) and sodium (Na+), while the pH model exhibited the weakest predictive performance. In addition, the usage of RF for band selection has the best effect compared with PCA and Pearson methods, which allows salt information of soda saline–alkali soils in Songnen Plain to be predicted precisely. Full article
(This article belongs to the Special Issue Saline–Alkali Land Ecology and Soil Management)
Show Figures

Figure 1

Figure 1
<p>Distribution of research area and sampling points.</p>
Full article ">Figure 2
<p>Preprocessing of a typical crack image. (<b>a</b>) clipped colorful crack image, (<b>b</b>) grayscale crack image, (<b>c</b>) binary crack image, (<b>d</b>) skeletonized image.</p>
Full article ">Figure 3
<p>The schematic diagram for measuring the spectral reflectance of cracked soil surface.</p>
Full article ">Figure 4
<p>Spectral measurement areas for all cracked soil samples.</p>
Full article ">Figure 5
<p>Cross-correlation heat map among different salt parameters.</p>
Full article ">Figure 6
<p>Reflectance curves of all soil samples.</p>
Full article ">Figure 7
<p>The results of 10-fold cross-validation.</p>
Full article ">Figure 8
<p>The top 20 bands were selected based on the RF method. (<b>a</b>) salinity, (<b>b</b>) EC, (<b>c</b>) Na<sup>+</sup>, (<b>d</b>) pH.</p>
Full article ">Figure 9
<p>The score of principal component.</p>
Full article ">Figure 10
<p>Correlation curves between spectral reflectance and four main salt parameters.</p>
Full article ">Figure 11
<p>Fitting results between the measured and estimated salt parameters.</p>
Full article ">
21 pages, 21293 KiB  
Article
Analyzing Spatial Distribution and Influencing Factors of Soil Organic Matter in Cultivated Land of Northeast China: Implications for Black Soil Protection
by Depiao Kong, Nanchen Chu, Chong Luo and Huanjun Liu
Land 2024, 13(7), 1028; https://doi.org/10.3390/land13071028 - 9 Jul 2024
Cited by 2 | Viewed by 948
Abstract
Soil organic matter (SOM) in cultivated land is vital for land quality and food security. This study examines SOM distribution and influencing factors in northeastern China, providing insights for sustainable agriculture. Utilizing 10 m resolution SOM data, the analysis covers regions including the [...] Read more.
Soil organic matter (SOM) in cultivated land is vital for land quality and food security. This study examines SOM distribution and influencing factors in northeastern China, providing insights for sustainable agriculture. Utilizing 10 m resolution SOM data, the analysis covers regions including the Greater and Lesser Khingan Mountains, Liaohe Plain, Sanjiang Plain, Songnen Plain, the northwest semi-arid region, and the low hilly areas of Paektu Mountain. The Geodetector method is employed to assess various influencing factors. The key findings are as follows: (1) The average SOM content in Northeast China (37.70 g/kg) surpasses the national average, is highest in the Greater and Lesser Khingan Mountains (49.32 g/kg), and lowest in the northwest semi-arid region (26.15 g/kg). (2) SOM content is maximized in regions with high altitudes, steep slopes, low temperatures, and moderate precipitation. (3) The annual average temperature is the primary factor influencing SOM distribution, with a combination of temperature and administrative divisions providing better explanatory power. (4) SOM trends vary across protected areas, with slope being critical in semi-humid plains, elevation in arid regions, and no dominant factors identified in the Sanjiang Plain. These findings underscore the need for tailored black soil protection policies to effectively leverage local resources and preserve ecosystem integrity. Full article
Show Figures

Figure 1

Figure 1
<p>Overview of the study area.</p>
Full article ">Figure 2
<p>Raster plot of each influencing factor.</p>
Full article ">Figure 3
<p>Framework of this study.</p>
Full article ">Figure 4
<p>Spatial distribution statistics of SOM content of cultivated land in Northeast China.</p>
Full article ">Figure 5
<p>Spatial distribution statistics of average SOM content of cultivated land on DEM in Northeast China.</p>
Full article ">Figure 6
<p>Spatial distribution statistics of average SOM content in cultivated land over annual average precipitation in Northeast China.</p>
Full article ">Figure 7
<p>Statistics on the spatial distribution of the average SOM content of cultivated land in Northeast China in terms of annual average temperature.</p>
Full article ">Figure 8
<p>Spatial distribution statistics of average SOM content of cultivated land on slope in Northeast China.</p>
Full article ">Figure 9
<p>Strength analysis of the effects of each influencing factor individually (<b>left</b>) and interactively (<b>right</b>) on the spatial distribution of SOM content in cultivated land in Northeast China.</p>
Full article ">Figure A1
<p>Spatial distribution statistics of average SOM content of cultivated land on DEM in six regions.</p>
Full article ">Figure A2
<p>Spatial distribution statistics of average SOM content in cultivated land over annual average precipitation in six regions.</p>
Full article ">Figure A3
<p>Statistics on the spatial distribution of the average SOM content of cultivated land in six regions in terms of annual average temperature.</p>
Full article ">Figure A4
<p>Spatial distribution statistics of average SOM content of cultivated land on slope in six regions.</p>
Full article ">Figure A5
<p>Influencing factors individually (<b>top</b>) and interactively (<b>bottom</b>) on the spatial distribution of SOM content in cultivated land in six regions.</p>
Full article ">Figure A6
<p>Trends in SOM content of cultivated land along annual average precipitation in Northeast China.</p>
Full article ">Figure A7
<p>Trend of SOM content of cultivated land along the slope in SNP.</p>
Full article ">Figure A8
<p>Trend of SOM content of cultivated land along the slope in LHP.</p>
Full article ">Figure A9
<p>Trends of SOM content along elevation in cultivated land in NSA.</p>
Full article ">Figure A10
<p>Trends of SOM content along elevation in cultivated land in GLKM.</p>
Full article ">Figure A11
<p>Trends of SOM content along elevation in cultivated land in LHPM.</p>
Full article ">
19 pages, 12073 KiB  
Article
Analyzing Spatio-Temporal Dynamics of Grassland Resilience and Influencing Factors in the West Songnen Plain, China, for Eco-Restoration
by Gefei Wang, Zhenyu Shi, Huiqing Wen, Yansu Bo, Haoming Li and Xiaoyan Li
Plants 2024, 13(13), 1860; https://doi.org/10.3390/plants13131860 - 5 Jul 2024
Viewed by 635
Abstract
Grassland plays an indispensable role in the stability and development of terrestrial ecosystems. Quantitatively assessing grassland resilience is of great significance for conducting research on grassland ecosystems. However, the quantitative measurement of resilience is difficult, and research on the spatio-temporal variation of grassland [...] Read more.
Grassland plays an indispensable role in the stability and development of terrestrial ecosystems. Quantitatively assessing grassland resilience is of great significance for conducting research on grassland ecosystems. However, the quantitative measurement of resilience is difficult, and research on the spatio-temporal variation of grassland resilience remains incomplete. Utilizing the Global Land Surface Satellite (GLASS) leaf area index (LAI) product derived from MODIS remote sensing data, along with land cover and meteorological data, this paper constructed the grassland resilience index (GRI) in the west Songnen Plain, China, a typical region with salt and alkali soils. This paper analyzed the spatio-temporal changes of the GRI and explored the contribution of climate factors, human activities, and geographical factors to the GRI. The results revealed that from 2000 to 2021, the GRI in the study area ranged from 0.1 to 0.22, with a multi-year average of 0.14. The average GRI exhibited a pattern of high-value aggregations in the north and low-value distributions in the south. Trend analysis indicated that areas with an improved GRI accounted for 59.09% of the total grassland area, but there were still some areas with serious degradation. From 2000 to 2015, the latitude and mean annual temperature (MAT) were principal factors to control the distribution of the GRI. In 2020, the mean annual precipitation (MAP) and MAT played important roles in the distribution of the GRI. From 2000 to 2021, the influence of human activities was consistently less significant compared to geographical location and climate variables. Full article
Show Figures

Figure 1

Figure 1
<p>The location of the study area (<b>a</b>) and the distribution of grasslands (<b>b</b>).</p>
Full article ">Figure 2
<p>Simplified process of ecosystem anomaly [<a href="#B24-plants-13-01860" class="html-bibr">24</a>].</p>
Full article ">Figure 3
<p>Technical roadmap.</p>
Full article ">Figure 4
<p>Variation of the mean annual GRI in the west Songnen Plain from 2000 to 2021.</p>
Full article ">Figure 5
<p>Spatial distribution of the mean GRI in the west Songnen Plain from 2000 to 2021.</p>
Full article ">Figure 6
<p>Spatial distribution (<b>a</b>) and grade (<b>b</b>) of the slope for the GRI from 2000 to 2021.</p>
Full article ">Figure 7
<p>Slope of the MAT (<b>a</b>) and MAP (<b>b</b>); Spatial distribution of the correlation between the MAT (<b>c</b>), MAP (<b>d</b>), and GRI.</p>
Full article ">Figure 8
<p>The RCRs of the influencing factors to the spatial distribution of the GRI from 2000 to 2021.</p>
Full article ">Figure 9
<p>RCRs of factors to the spatial distribution of the GRI in five prefecture-level cities from 2000 to 2020.</p>
Full article ">Figure 10
<p>Clustered box plot of RCRs for factors in five prefecture-level cities from 2000 to 2020.</p>
Full article ">
17 pages, 7180 KiB  
Article
Alfalfa with Forage Crop Rotation Alleviates Continuous Alfalfa Obstacles through Regulating Soil Enzymes and Bacterial Community Structures
by Yanxia Xu, Zhuxiu Liu, Zhongbao Shen, Zhao Yang, Xuepeng Fu, Xiaolong Wang, Shasha Li, Hua Chai, Ruoding Wang, Xiaobing Liu and Junjie Liu
Agronomy 2024, 14(7), 1349; https://doi.org/10.3390/agronomy14071349 - 21 Jun 2024
Viewed by 971
Abstract
Alfalfa is a perennial herbaceous forage legume that is significantly and adversely affected by monocropping. Crop rotation is the most effective measure to overcome continuous cropping obstacles. However, the mechanisms of how bacterial communities are affected and the potential links between these effects [...] Read more.
Alfalfa is a perennial herbaceous forage legume that is significantly and adversely affected by monocropping. Crop rotation is the most effective measure to overcome continuous cropping obstacles. However, the mechanisms of how bacterial communities are affected and the potential links between these effects and cropping systems remain poorly understood. Based on a long-term field experiments with continuous alfalfa and forage crops with alfalfa rotation in the black soil region of the western Songnen Plain in Northeast China, the alterations in soil bacterial community structure using high-throughput sequencing of the 16S rRNA gene and soil chemical properties and enzyme activities were analyzed. The alfalfa–forage oats–silage maize–alfalfa and alfalfa–silage maize–forage oats–alfalfa system significantly increase the levels of total phosphorus and available phosphorus, and promote the activities of acid phosphatase, β-glucosidase, leucine aminopeptidase, and N-acetyl-β-glucosaminidase in comparison to continuous alfalfa. While alfalfa crop rotation did not affect the α-diversity of soil bacteria, it significantly altered the bacterial community composition and structure. Some key taxa are significantly enriched in the crop rotation system soils, including Bacillus, Sphingobium, Paenibacillus, Hydrogenispora, Rubrobacter, Haliangium, and Rubellimicrobium. Additionally, crop rotation with alfalfa increased the stability and complexity of the soil bacterial co-occurrence network. Based on our findings, we recommend promoting the alfalfa–forage oats–silage maize–alfalfa and alfalfa–silage maize–forage oats–alfalfa rotation systems as ideal practices for overcoming the challenges associated with continuous cropping of alfalfa. These systems not only enhance soil nutrient content and enzyme activities but also foster a beneficial microbial community, ultimately improving soil functionality and crop performance. Full article
(This article belongs to the Section Grassland and Pasture Science)
Show Figures

Figure 1

Figure 1
<p>Effects of continuous cropping/crop rotation of alfalfa on soil bacterial diversity. Bacterial Shannon diversity in bulk soil (<b>a</b>) and rhizosphere soil (<b>b</b>). Principal coordinate analysis (PCoA) based on Bray−Curtis distance for differences in the composition of bacterial communities in bulk soil (<b>c</b>) and rhizosphere soil (<b>d</b>). * <span class="html-italic">p</span> &lt; 0.05; *** <span class="html-italic">p</span> &lt; 0.001; NS: not significant.</p>
Full article ">Figure 2
<p>The relative abundance of soil bacterial communities at phylum level. Changes in the relative abundance of bacterial phylum in bulk soil (<b>a</b>) and rhizosphere soil (<b>b</b>) under different alfalfa continuous cropping/crop rotation treatments. Significance levels are as follows: * <span class="html-italic">p</span> &lt; 0.05; ** <span class="html-italic">p</span> &lt; 0.01; *** <span class="html-italic">p</span> &lt; 0.001.</p>
Full article ">Figure 3
<p>The relative abundance of soil bacterial communities at genus level. Changes in the relative abundance of dominant genera in bulk soil (<b>a</b>) and rhizosphere soil (<b>b</b>) under different alfalfa continuous cropping/crop rotation treatments. Significance levels are as follows: * <span class="html-italic">p</span> &lt; 0.05; ** <span class="html-italic">p</span> &lt; 0.01.</p>
Full article ">Figure 4
<p>Certain OTUs are enriched and depressed under alfalfa continuous cropping/crop rotation. Volcano plot was used to describe OTU differences in AOCA (<b>a</b>), ACCA (<b>b</b>), AOOA (<b>c</b>), and ACOA (<b>d</b>) in bulk soil was expressed as BAAAA, BAOCA, BACCA, BAOOA, and BACOA, respectively. The OTU differences in AOCA (<b>f</b>), ACCA (<b>g</b>), AOOA (<b>h</b>), and ACOA (<b>i</b>) in rhizosphere soil was expressed as RAAAA, RAOCA, RACCA, RAOOA and RACOA, respectively. Each point represents an individual OTU, and the position along the y axis represents the abundance fold change compared with the abundance fold change in AAAA. Number of significantly enriched differential OTUs between continuous cropping and rotation treatments in bulk soil (<b>e</b>) and rhizosphere soil (<b>j</b>).</p>
Full article ">Figure 5
<p>Effects of continuous cropping/crop rotation of alfalfa on soil bacterial co-occurrence network. The bacterial co-occurrence network of AAAA (<b>a</b>), AOCA (<b>b</b>), ACCA (<b>c</b>), AOOA (<b>d</b>) and ACOA (<b>e</b>) in bulk soil was expressed as BAAAA, BAOCA, BACCA, BAOOA, and BACOA, respectively. The bacterial co-occurrence network of AAAA (<b>f</b>), AOCA (<b>g</b>), ACCA (<b>h</b>), AOOA (<b>i</b>), and ACOA (<b>j</b>) in rhizosphere soil was expressed as RAAAA, RAOCA, RACCA, RAOOA, and RACOA, respectively. The nodes are colored according to bacterial phylum. Node size indicates the degree of connection. Edge color represents positive (red) and negative (green) correlations. The key taxa are annotated on the network.</p>
Full article ">Figure 6
<p>Redundancy analysis (RDA) of bacterial OTUs and soil parameters. The relationship between bacterial OTUs and soil parameters in bulk soil (<b>a</b>) and rhizosphere soil (<b>b</b>). Significance levels are as follows: * <span class="html-italic">p</span> &lt; 0.05; ** <span class="html-italic">p</span> &lt; 0.01; *** <span class="html-italic">p</span> &lt; 0.001.</p>
Full article ">Figure 7
<p>Heatmap showing the strength of correlations between soil parameters and soil bacterial communities at phylum level in bulk soil (<b>a</b>) and rhizosphere soil (<b>b</b>). Significance levels are as follows: * <span class="html-italic">p</span> &lt; 0.05; ** <span class="html-italic">p</span> &lt; 0.011. Color legend represents the values of correlation coefficients (r). Negative values indicate negative correlations, while positive values indicate positive correlations.</p>
Full article ">
18 pages, 13047 KiB  
Article
The Contribution of Saline-Alkali Land to the Terrestrial Carbon Stock Balance: The Case of an Important Agriculture and Ecological Region in Northeast China
by Lei Chang, Tianhang Ju, Huijia Liu and Yuefen Li
Land 2024, 13(7), 900; https://doi.org/10.3390/land13070900 - 21 Jun 2024
Viewed by 858
Abstract
Saline-alkali land is an important component of terrestrial ecosystems and may serve as a carbon sink but its net contribution to the overall terrestrial carbon sink is unknown. Using methods recommended by the IPCC, this study evaluates the impacts of interconverting saline-alkali and [...] Read more.
Saline-alkali land is an important component of terrestrial ecosystems and may serve as a carbon sink but its net contribution to the overall terrestrial carbon sink is unknown. Using methods recommended by the IPCC, this study evaluates the impacts of interconverting saline-alkali and non-saline-alkali land on terrestrial carbon stocks by measuring two major carbon pools (soil organic carbon and vegetation carbon) in the saline-alkali land of China’s Songnen Plain. Distinct phases in the evolution of the region’s terrestrial carbon stock were delineated, factors contributing to transitions between phases were identified, and the effects of changes in the saline-alkali land carbon stock on the overall terrestrial carbon sink were estimated. Between 2005 and 2020, the region’s saline-alkali land carbon stock initially increased, then declined, and finally increased again. However, the overall terrestrial carbon stock decreased by 0.5 Tg (1 Tg = 1012 g), indicating that the increase in the saline-alkali land carbon stock was due primarily to expansion of the saline-alkali land area. The conversion of non-saline-alkali land to saline-alkali land was a carbon-emitting process; consequently, in areas undergoing saline-alkali land change, the lower carbon density bound was equal to the carbon density of unconverted saline-alkali land and the upper bound was equal to the carbon density of unconverted non-saline-alkali land. In general, changes in the carbon stock of saline-alkali land correlated negatively with changes in the overall terrestrial carbon stock. The conversion of saline-alkali land into grassland and cropland through biochar improvement and the planting of saline-tolerant crops (Leymus chinensis, salt-tolerant rice) has a positive effect on promoting the enhancement of terrestrial carbon stocks. Full article
(This article belongs to the Special Issue Impact of Climate Change on Land and Water Systems)
Show Figures

Figure 1

Figure 1
<p>The study area and its location within China.</p>
Full article ">Figure 2
<p>Carbon stocks, carbon density, and proportion of carbon stocks in saline-alkali land from 2005 to 2020.</p>
Full article ">Figure 3
<p>Spatial variation and extent of saline-alkali land, 2005–2020.</p>
Full article ">Figure 4
<p>Transfer relationship between saline-alkali and non-saline-alkali land.</p>
Full article ">Figure 5
<p>Changes in carbon stocks and carbon density due to conversion of saline-alkali and non-saline-alkali land. In the figure, (<b>a</b>,<b>c</b>,<b>e</b>) represent the changes in carbon stocks after the conversion of saline-alkali and non-saline-alkali land; (<b>b</b>,<b>d</b>,<b>f</b>) represent the comparison of carbon density before and after the conversion of saline-alkali and non-saline-alkali land. 1, 2, 3, 4, 5, 6, and T denote cropland, woodland, grassland, wetland, construction land, underutilized land, and all land, respectively.</p>
Full article ">Figure 6
<p>Increases or decreases in carbon stocks due to conversion of saline-alkali to non-saline-alkali land, 2005–2020. In the figure, green labels indicate an increase in carbon stocks and red labels indicate a decrease in carbon stocks.</p>
Full article ">Figure 7
<p>Comparison of carbon density before and after conversion of non-saline-alkali to saline-alkali land. The labels 1, 2, 3, 4, 5, 6, and T denote cropland, woodland, grassland, wetland, construction land, underutilized land, and total land, respectively.</p>
Full article ">
Back to TopTop