[go: up one dir, main page]

 
 
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

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

Article Types

Countries / Regions

remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline

Search Results (1,393)

Search Parameters:
Keywords = multi-source region

Order results
Result details
Results per page
Select all
Export citation of selected articles as:
16 pages, 3735 KiB  
Article
Material Activity in Debris Flow Watersheds Pre- and Post-Strong Earthquake: A Case Study from the Wenchuan Earthquake Epicenter
by Yu Yang, Ming Chen, Yinghua Cai, Chenxiao Tang, Wenli Huang and Chenhao Xia
Water 2024, 16(16), 2284; https://doi.org/10.3390/w16162284 - 13 Aug 2024
Abstract
The 2008 Wenchuan earthquake released vast quantities of loose material, significantly influencing post-earthquake material dynamics, particularly through recurrent debris flow disasters that posed long-term threats to the earthquake-affected area. To explore the transport and involvement of loose materials in debris flow events within [...] Read more.
The 2008 Wenchuan earthquake released vast quantities of loose material, significantly influencing post-earthquake material dynamics, particularly through recurrent debris flow disasters that posed long-term threats to the earthquake-affected area. To explore the transport and involvement of loose materials in debris flow events within earthquake-affected basins, this study focuses on a representative area near the Wenchuan epicenter, creating a multi-temporal database of active landslides and channel materials pre- and post-earthquake, quantitatively assessing material transport and source replenishment in debris flow basins, and categorizing debris flows based on channel material activity, post-earthquake historical activity, and sustainability of activity. This study revealed that pre-earthquake material activity was concentrated in the watershed’s upper regions, while post-earthquake materials were progressively transported from the central to the lower regions, with many small co-seismic landslides ceasing activity. The supply area ratio from active landslides capable of recharging debris flows, i.e., those connected to channels, consistently remained at approximately 72%, with the peak area of channel material activity comprising approximately 2.5% of the total watershed area. Channel material activity areas serve as valuable indicators for hazard assessment in regions lacking historical debris flow data, with the watershed area predominantly determining the sustainability of post-earthquake debris flow activity. Full article
Show Figures

Figure 1

Figure 1
<p>A focus on our study area, which is primarily situated between two major faults: the Maoxian-Wenchuan fault (pink) and the Yingxiu-Beichuan fault (red), with yellow stars marking the location of the Wenchuan earthquake epicenter. Mapped debris flows (DF) are depicted in light green with grey boundary lines, numbered sequentially. The main trunk of the MinJiang River and its major tributaries are highlighted in blue and all mapped sub-catchments flow into this river. Key locations such as Yinxing and Yingxiu are indicated with green dots. Elevation, represented with a gradient color scale ranging from 606 m above sea level (light brown) to 4560 m above sea level (dark brown), is sourced from ALOS-PALSAR RTC with a resolution of 30 m.</p>
Full article ">Figure 2
<p>Example of multi-temporal interpretation. (<b>a</b>) Co-seismic landslides triggered by the Wenchuan earthquake; (<b>b</b>,<b>d</b>,<b>e</b>) post-earthquake landslides; (<b>c</b>,<b>f</b>) active channel materials.</p>
Full article ">Figure 3
<p>Schematic calculation of d<sub>g</sub>. D represents the horizontal projection length from the centroid of each interpreted material to the tangent at the gully outlet, while L denotes the maximum horizontal projection length from the upstream watershed divide to the tangent at the gully outlet.</p>
Full article ">Figure 4
<p>A demonstration of landslide and channel material activity within the watershed. (<b>a</b>) Google satellite imagery of the demonstration area captured in July 2008, December 2014, and October 2019. (<b>b</b>) Photograph documenting the 2015 field investigation of the landslide identified by the purple arrow. This photo was taken by the author. (<b>c</b>) Photograph documenting the 2015 field investigation of the landslide identified by the yellow arrow. This photo was taken by the author.</p>
Full article ">Figure 5
<p>Classification of interpreted active landslides based on their connectivity to the channel, illustrated with an example from Xiaojia Gully.</p>
Full article ">Figure 6
<p>(<b>a</b>) Interpretation across different time periods shows active landslides, including new and remobilized landslides, mapped in blue and red, respectively. Active channel materials are mapped in orange. (<b>b</b>) Variations in the area of active landslides relative to rainfall are shown. Red dots represent active landslide areas, while blue bar charts depict annual rainfall during the flood season. (<b>c</b>) Variations in the area of active channel materials relative to rainfall are shown. Yellow dots represent active channel material areas, while blue bar charts depict annual rainfall during the flood season.</p>
Full article ">Figure 7
<p>Dynamics of active material relative to area across different periods.</p>
Full article ">Figure 8
<p>Active landslides connected to channels over various time periods, with interpretation showing active connected landslides mapped in purple. The maps illustrate the evolution of the area of active connected landslides from 2005 to 2020, highlighting significant increases post-2008 due to the Wenchuan earthquake.</p>
Full article ">Figure 9
<p>(<b>a</b>) The relationship between active landslides and actively connected landslides. (<b>b</b>) The relationship between actively connected landslides and active channel materials.</p>
Full article ">Figure 10
<p>Classification outcomes of debris flows. (<b>a</b>) Utilizing the active channel material area from <a href="#water-16-02284-t003" class="html-table">Table 3</a> and the natural breaks method, the activity levels of 31 debris flow watersheds were categorized into four groups: extremely active (1.01–3 km<sup>2</sup>), highly active (0.33–1.01 km<sup>2</sup>), moderately active (0.12–0.33 km<sup>2</sup>), and lowly active (0–0.12 km<sup>2</sup>). (<b>b</b>) Utilizing the number of post-earthquake events from <a href="#water-16-02284-t003" class="html-table">Table 3</a> and the natural breaks method, post-earthquake debris flow event counts for 31 watersheds were categorized into four activity levels: extremely active (&gt;6 events), highly active (3–6 events), moderately active (1–3 events), and lowly active (1 event). (<b>c</b>) Utilizing the occurrence of events post-2018 from <a href="#water-16-02284-t003" class="html-table">Table 3</a>, post-earthquake debris flow event sustainability levels were categorized into two groups: relatively strong (Yes) and relatively weak (No).</p>
Full article ">
21 pages, 7056 KiB  
Article
Ecological Security Pattern Construction and Multi-Scenario Risk Early Warning (2020–2035) in the Guangdong–Hong Kong–Macao Greater Bay Area, China
by Junjie Ma, Zhixiong Mei, Xinyu Wang, Sichen Li and Jiangsen Liang
Land 2024, 13(8), 1267; https://doi.org/10.3390/land13081267 - 12 Aug 2024
Viewed by 202
Abstract
The effectiveness of ecological security patterns (ESPs) in maintaining regional ecological stability and promoting sustainable development is widely recognized. However, limited research has focused on the early warning of risks inherent in ESPs. In this study, the Guangdong–Hong Kong–Macao Greater Bay Area (GHKMGBA) [...] Read more.
The effectiveness of ecological security patterns (ESPs) in maintaining regional ecological stability and promoting sustainable development is widely recognized. However, limited research has focused on the early warning of risks inherent in ESPs. In this study, the Guangdong–Hong Kong–Macao Greater Bay Area (GHKMGBA) is taken as the study area, and ecological security risk zones are delineated by combining the landscape ecological risk index and habitat quality, and a multi-level ESP is constructed based on the circuit theory. The PLUS model was employed to simulate future built-up land expansion under different scenarios, which were then extracted and overlaid with the multi-level ESP to enable the multi-scenario early warning of ESP risks. The results showed the following: The ESP in the central plains and coastal areas of the GHKMGBA exhibits a high level of ecological security risk, whereas the peripheral forested areas face less threat, which is crucial for regional ecological stability. The ESP, comprising ecological sources, corridors, and pinch points, is crucial for maintaining regional ecological flow stability, with tertiary corridors under significant stress and risk in all scenarios, requiring focused restoration and enhancement efforts. There are significant differences in risk early warning severity within the ESP across various development scenarios. Under the ecological protection scenario, the ESP will have the best early warning situation, effectively protecting ecological land and reducing ecological damage, providing a valuable reference for regional development policies. However, it must not overlook economic development and still needs to further seek a balance between economic growth and ecological protection. Full article
Show Figures

Figure 1

Figure 1
<p>Location of the study area.</p>
Full article ">Figure 2
<p>Spatial distribution of ecological security risk zones in the study area in 2020.</p>
Full article ">Figure 3
<p>(<b>a</b>–<b>c</b>) Spatial distribution of ecological sources (Note: (<b>a</b>) Ecosystem services (<b>b</b>) Ecological sensitivity (<b>c</b>) Ecological sources).</p>
Full article ">Figure 4
<p>Resistance surface.</p>
Full article ">Figure 5
<p>(<b>a</b>,<b>b</b>) Multi-level ecological security patterns in the study area in 2020 (Note: (<b>a</b>) Ecological security patterns (<b>b</b>) Multi-level ecological security patterns).</p>
Full article ">Figure 6
<p>(<b>a</b>–<b>f</b>) Early warning identification of risks to ecological security patterns in 2035 under different scenarios (Note: (<b>a</b>–<b>c</b>) Multi-scenario land use simulation (<b>d</b>–<b>f</b>) Ecological security pattern multi-scenario risk early warning).</p>
Full article ">Figure 7
<p>Spatial and quantitative distribution of vigilance of ecological security patterns in 2035 under different scenarios.</p>
Full article ">Figure 8
<p>Validation of ecological source.</p>
Full article ">
25 pages, 7111 KiB  
Article
Spatial–Temporal Changes in the Distribution of Populus euphratica Oliv. Forests in the Tarim Basin and Analysis of Influencing Factors from 1990 to 2020
by Xuefei Guo, Lijun Zhu, Zhikun Yang, Chaobin Yang and Zhijun Li
Forests 2024, 15(8), 1384; https://doi.org/10.3390/f15081384 - 7 Aug 2024
Viewed by 305
Abstract
Understanding the spatiotemporal evolution patterns of Populus euphratica Oliv. (P. euphratica) forests in the Tarim Basin (TB) and their influencing factors is crucial for regional ecological security and high-quality development. However, there is currently a lack of large-area, long-term systematic monitoring. [...] Read more.
Understanding the spatiotemporal evolution patterns of Populus euphratica Oliv. (P. euphratica) forests in the Tarim Basin (TB) and their influencing factors is crucial for regional ecological security and high-quality development. However, there is currently a lack of large-area, long-term systematic monitoring. This study utilized multi-source medium and high-resolution remote sensing images from the Landsat series and Sentinel-2, applying a Random Forest classification model to obtain distribution data of P. euphratica forests and shrublands in 14 areas of the TB from 1990 to 2020. We analyzed the effects of river distance, water transfer, and farmland on their distribution. Results indicated that both P. euphratica forests and shrublands decreased during the first 20 years and increased during the last 10 years. Within 1.5 km of river water transfer zones, P. euphratica forests more frequently converted to shrublands, while both forests and shrublands showed recovery in low-frequency water transfer areas. Farmland encroachment was most significant beyond 3 km from rivers. To effectively protect P. euphratica forests, we recommend intermittent low-frequency water transfers within 3 km of rivers and stricter management of agricultural expansion beyond 3 km. These measures will help maintain a balanced ecosystem and promote the long-term sustainability of P. euphratica forests. Full article
(This article belongs to the Section Forest Inventory, Modeling and Remote Sensing)
Show Figures

Figure 1

Figure 1
<p>Spatial distribution map of land use types in the TB (Approved map No. GS (2019)1822).</p>
Full article ">Figure 2
<p>The number of images from different sensors during 1990–2020.</p>
Full article ">Figure 3
<p>Schematic representation of the overall methodological workflow.</p>
Full article ">Figure 4
<p>Extraction and classification of <span class="html-italic">P. euphratica</span> forest boundary (1 is Landsat image (red is vegetation); 2 is Sentinel-2 image (red is vegetation); 3 is DJI P4 multi-spectral UAV aerial image (green is vegetation); 4 is the result of random forest classification; 5 for field photos). (A. Up-Tarim River; B. Mid-Tarim River; C. Down-Tarim River; E. Down-Yarkand River; F. Mid-Hotan River; H. Kongque River; I. Kashgar River; J. Cherchen River).</p>
Full article ">Figure 5
<p>Schematic diagram of <span class="html-italic">P. euphratica</span> forest and shrubland area change in the TB (<b>a</b>). 1990; (<b>b</b>). 1995; (<b>c</b>). 2000; (<b>d</b>). 2005; (<b>e</b>). 2010; (<b>f</b>). 2015, (<b>g</b>) 2020 (<b>h</b>) 1990–2020.</p>
Full article ">Figure 6
<p>Distribution of <span class="html-italic">P. euphratica</span> forest (<b>a</b>) and shrublands (<b>b</b>) in different rivers of the TB at different distances from rivers. (A. Up-Tarim River; B. Mid-Tarim River; C. Down-Tarim River; D. Up-Yarkand River; E. Down-Yarkand River; F. Mid-Hotan River; G. Down-Hotan River; H. Kongque River; I. Kashgar River; J. Cherchen River; K. Rivers in the northern Kunlun Mountains; L. Keriya River; N. Sangzhu River).</p>
Full article ">Figure 7
<p>(<b>a</b>) Analysis of the intensity of conversion from <span class="html-italic">P. euphratica</span> forest; (<b>b</b>) analysis of the intensity of conversion from shrublands to <span class="html-italic">P. euphratica</span> forest; (<b>c</b>) distribution of <span class="html-italic">P. euphratica</span> forest and shrublands conversion at different river distances; (<b>d</b>) distribution of <span class="html-italic">P. euphratica</span> forest and shrublands conversion at different rivers in the TB from 1990 to 2020.</p>
Full article ">Figure 8
<p>(<b>a</b>) Mutual conversion of <span class="html-italic">P. euphratica</span> forest, shrublands, and farmland Sankey map; (<b>b</b>) mutual conversion of <span class="html-italic">P. euphratica</span> forest, shrublands, and farmland spatial distribution map; (<b>c</b>) mutual conversion of <span class="html-italic">P. euphratica</span> forest, shrublands, and farmland in different rivers; (<b>d</b>) farmland-to-shrubland conversion in different rivers; (<b>e</b>) mutual conversion of <span class="html-italic">P. euphratica</span> forest, shrublands, and farmland in different distance from rivers; (<b>f</b>) mutual conversion of <span class="html-italic">P. euphratica</span> forest, shrublands, and farmland in different distance from rivers; (<b>f</b>) conversion of farmland to shrubland at different distances from rivers in the TB from 1990 to 2020.</p>
Full article ">Figure 9
<p>(<b>a</b>) Spatial distribution of <span class="html-italic">P. euphratica</span> forest and shrublands inter-conversion under different water transfer frequencies; (<b>b</b>) different frequencies of water transfers on different rivers; (<b>c</b>) <span class="html-italic">P. euphratica</span> and shrublands inter-conversion under different water delivery frequencies of different rivers; (<b>d</b>) <span class="html-italic">P. euphratica</span> and shrublands inter-conversion under different water delivery frequencies at different distances from the rivers in the TB from 1990 to 2020.</p>
Full article ">Figure 10
<p>(<b>a</b>) Spatial distribution map of farmland change and <span class="html-italic">P. euphratica</span> forest and shrubland transformation in the TB from 1990 to 2020; (<b>b</b>) Pearson correlation heatmap of farmland, <span class="html-italic">P. euphratica</span> forest, and shrubland transformation.”→”Indicates the direction of transformation.</p>
Full article ">
25 pages, 6036 KiB  
Article
Research on Improving the Accuracy of SIF Data in Estimating Gross Primary Productivity in Arid Regions
by Wei Liu, Yu Wang, Ali Mamtimin, Yongqiang Liu, Jiacheng Gao, Meiqi Song, Ailiyaer Aihaiti, Cong Wen, Fan Yang, Wen Huo, Chenglong Zhou, Jian Peng and Hajigul Sayit
Land 2024, 13(8), 1222; https://doi.org/10.3390/land13081222 - 7 Aug 2024
Viewed by 259
Abstract
Coupling solar-induced chlorophyll fluorescence (SIF) with gross primary productivity (GPP) for ecological function integration research presents numerous uncertainties, especially in ecologically fragile and climate-sensitive arid regions. Therefore, evaluating the suitability of SIF data for estimating GPP and the feasibility of improving its accuracy [...] Read more.
Coupling solar-induced chlorophyll fluorescence (SIF) with gross primary productivity (GPP) for ecological function integration research presents numerous uncertainties, especially in ecologically fragile and climate-sensitive arid regions. Therefore, evaluating the suitability of SIF data for estimating GPP and the feasibility of improving its accuracy in the northern region of Xinjiang is of profound significance for revealing the spatial distribution patterns of GPP and the strong coupling relationship between GPP and SIF in arid regions, achieving the goal of “carbon neutrality” in arid regions. This study is based on multisource SIF satellite data and GPP observation data from sites in three typical ecosystems (cultivated and farmlands, pasture grasslands, and desert vegetation). Two precision improvement methods (canopy and linear) are used to couple multiple indicators to determine the suitability of multisource SIF data for GPP estimation and the operability of accuracy improvement methods in arid regions reveal the spatial characteristics of SIF (GPP). The results indicate the following. (1) The interannual variation of GPP shows an inverted “U” shape, with peaks values in June and July. The cultivated and farmland areas have the highest peak value among the sites (0.35 gC/m2/month). (2) The overall suitability ranking of multisource SIF satellite products for GPP estimation in arid regions is RTSIF > CSIF > SIF_OCO2_005 > GOSIF. RTSIF shows better suitability in the pasture grassland and cultivated and farmland areas (R2 values of 0.85 and 0.84, respectively). (3) The canopy method is suitable for areas with a high leaf area proportion (R2 improvement range: 0.05–0.06), while the linear method is applicable across different surface types (R2 improvement range: 0.01–0.13). However, the improvement effect of the linear method is relatively weaker in areas with high vegetation cover. (4) Combining land use data, the overall improvement of SIF (GPP) is approximately 0.11%, and the peak values of its are mainly distributed in the northern and southern slopes of the Tianshan Mountains, while the low values are primarily found in the Gurbantunggut Desert. The annual mean value of SIF (GPP) is about 0.13 mW/m2/nm/sr. This paper elucidates the applicability of SIF for GPP estimation and the feasibility of improving its accuracy, laying the theoretical foundation for the spatiotemporal coupling study of GPP and SIF in an arid region, and providing practical evidence for achieving carbon neutrality goals. Full article
(This article belongs to the Special Issue Land-Based Greenhouse Gas Mitigation for Carbon Neutrality)
Show Figures

Figure 1

Figure 1
<p>(<b>a</b>) Specific locations of the Tianshan Mountains, Ulan Usu Station, Ulastai Station, and Kelameili Station in Xinjiang. (<b>b</b>) Schematic representation of the elevations of the study area. (<b>c</b>) Schematic representation of the land use types at the study area.</p>
Full article ">Figure 2
<p>Interannual variation of monthly average GPP at each site in 2020 (excluding nighttime values).</p>
Full article ">Figure 3
<p>(<b>a</b>) The linear regression fitting of 2020 GPP data from three site with corresponding site data of CSIF satellite products. (<b>b</b>) The linear regression fitting of 2020 GPP data from three site with corresponding site data of RTSIF satellite products. (<b>c</b>) The linear regression fitting of 2020 GPP data from three site with corresponding site data of SIF-OCO-005 satellite products. (<b>d</b>) The linear regression fitting of 2020 GPP data from three site with corresponding site data of GOSIF satellite products.</p>
Full article ">Figure 4
<p>The spatial distribution characteristics of annual mean values of multisource SIF satellite products.</p>
Full article ">Figure 5
<p>Responsiveness of multisource SIF satellite products to major influencing factors of GPP (** indicates significance at the 0.5 level). (<b>a</b>) Kelameili Station, desert vegetation area. (<b>b</b>) Ulastai Station, pasture and grassland area. (<b>c</b>) Ulastai Station, pasture and grassland area. (<b>c</b>) Ulan Usu Staion, cultivate land and farmland area.</p>
Full article ">Figure 6
<p>Differences in GPP/SIF values under different weather conditions. (<b>a</b>) Kelameili Station, desert vegetation area. (<b>b</b>) Ulastai Station, pasture and grassland area. (<b>c</b>) Ulastai Station, pasture and grassland area. (<b>c</b>) Ulan Usu Staion, cultivate land and farmland area.</p>
Full article ">Figure 7
<p>Linear fitting graph of 2020 GPP data and RTSIF corresponding station data for each station after improving based on the canopy method. (<b>a</b>) Ulastai Station, pasture and grassland area. (<b>b</b>) Ulan Usu Staion, cultivate land and farmland area.</p>
Full article ">Figure 8
<p>Linear fitting diagram between the 2020 GPP data of each station and the corresponding RTSIF station data after improving based on the linear method. (<b>a</b>) Kelameili Station, desert vegetation area. (<b>b</b>) Ulastai Station, pasture and grassland area. (<b>c</b>) Ulastai Station, pasture and grassland area. (<b>c</b>) Ulan Usu Staion, cultivate land and farmland area.</p>
Full article ">Figure 9
<p>The R<sup>2</sup> fitting values for various sites based on two accuracy improvement methods: canopy and linear.</p>
Full article ">Figure 10
<p>Changes in spatial characteristics of quarterly average values before and after the improvement of SIF satellite product data. (<b>a1</b>–<b>d1</b>) The spatial variation characteristics of the mean values of each season before improvement, (<b>a1</b>) for spring, and so on. (<b>a2</b>–<b>d2</b>) The spatial variation characteristics of the mean values of each season after improvement, (<b>a2</b>) for spring, and so on.</p>
Full article ">
18 pages, 8753 KiB  
Article
Aeolian Environment Regionalization in Xinjiang and Suggestions for Sand Prevention in Typical Areas
by Jie Zhou, Hongjing Ren, Beibei Han, Yazhou Zhao and Haifeng Wang
Land 2024, 13(8), 1215; https://doi.org/10.3390/land13081215 - 6 Aug 2024
Viewed by 342
Abstract
The Xinjiang region is prone to frequent and complex wind and sand disasters, which present a significant challenge to the sustainable development of local areas. This research uses multi-source data to analyze the spatial distribution of the aeolian environment in Xinjiang, establishes a [...] Read more.
The Xinjiang region is prone to frequent and complex wind and sand disasters, which present a significant challenge to the sustainable development of local areas. This research uses multi-source data to analyze the spatial distribution of the aeolian environment in Xinjiang, establishes a four-level zoning scheme, and proposes recommendations for ecological management and engineering and control. Results indicate that (1) Xinjiang’s aeolian environment and its types exhibit spatial heterogeneity. The aeolian environment types display a high concentration in the eastern region and a low concentration in the western region. Furthermore, the aeolian environment types are concentrated in the basin region. Moreover, the aeolian environment types exhibit a meridional distribution pattern. (2) A four-level zoning system for aeolian environments in Xinjiang was developed, comprising two first-level zones, seven s-level subzones, 22 third-level wind zones, and 31 fourth-level subdivisions. (3) A structural model for a highway sand control system is proposed for aeolian environment types of subdivisions, including fixing-based, combined blocking and fixing, wind-blocking and sand-transferring, and combined blocking and fixing–transferring. The aeolian environment regionalization program proposed in this study can be a scientific reference for relevant departments in formulating and implementing sand prevention and control planning. Full article
Show Figures

Figure 1

Figure 1
<p>Location map of the study area. (<b>a</b>) Geographical location of Xinjiang, China; (<b>b</b>) Spatial distribution in sandy areas of Xinjiang; (<b>c</b>) Average wind speed, annual dominant wind direction, and location of main wind zones in sandy areas of Xinjiang.</p>
Full article ">Figure 2
<p>Data processing flowchart.</p>
Full article ">Figure 3
<p>Correlation coefficient and <span class="html-italic">VIF</span> value of each indicator. (The letters X1 to X8 represent each indicator: X1 represents soil organic carbon content, X2 represents soil water content, X3 represents soil sand content, X4 represents land-use type, X5 represents NDVI, X6 represents precipitation, X7 represents potential evapotranspiration, and X8 represents RDP. Same below.)</p>
Full article ">Figure 4
<p>Characteristics of aeolian environment in Xinjiang: (<b>a</b>) spatial distribution of wind-dynamic conditions in Xinjiang; (<b>b</b>) spatial distribution of the sand material base in Xinjiang; (<b>c</b>) spatial distribution of aeolian environment types in Xinjiang; and (<b>d</b>) characteristics of the change in aeolian environment type values with longitude.</p>
Full article ">Figure 5
<p>Xinjiang aeolian environmental regionalization.</p>
Full article ">Figure 6
<p>Suggestions for sand prevention in typical aeolian environment zones.</p>
Full article ">Figure 7
<p>Structural model of road sand protection system under different types of aeolian environments.</p>
Full article ">
24 pages, 282 KiB  
Article
Rural Land Circulation and Peasant Household Income Growth—Empirical Research Based on Structural Decomposition
by Wenwu Zhang, Shunji Zhao, Jinguo Wang, Xinyao Xia and Hongkui Jin
Sustainability 2024, 16(16), 6717; https://doi.org/10.3390/su16166717 - 6 Aug 2024
Viewed by 405
Abstract
How rural land transfer affects the growth of non-agricultural income and the changes in its sources are important research topics. This study uses the micro-data from the China Family Panel Studies (CFPS) spanning from 2014 to 2020 and empirically analyzes the impact of [...] Read more.
How rural land transfer affects the growth of non-agricultural income and the changes in its sources are important research topics. This study uses the micro-data from the China Family Panel Studies (CFPS) spanning from 2014 to 2020 and empirically analyzes the impact of rural land transfer on the growth of non-agricultural income, based on a multi-dimensional decomposition of rural household income structure. This study found that (1) land transfer has a significant promoting effect on the growth of non-agricultural income. Transferring out land is conducive to increasing wage income and transfer income, while transferring in land compensates for the decrease in operating income by achieving a higher operating income, ultimately leading to an increase in total income. (2) The effect of land transfer on the growth of non-agricultural income is higher in the Eastern region than in the Central and Western regions. The higher the education level of family members, the greater the income-increasing effect of land transfer on farmers. (3) Mechanism analysis shows that land transfer increases farmers’ opportunities for migrant work and improves farmers’ operational efficiency, which are the main channels for the growth in non-agricultural income. This study demonstrates that land circulation will promote farmers’ income growth and prosperity through rental income, share cooperation and dividends, labor transfer and wage income, industrial chain extension and value-added income, and policy support and subsidies. Full article
18 pages, 7338 KiB  
Article
Droop Frequency Limit Control and Its Parameter Optimization in VSC-HVDC Interconnected Power Grids
by Han Jiang, Yichen Zhou, Yi Gao and Shilin Gao
Energies 2024, 17(15), 3851; https://doi.org/10.3390/en17153851 - 5 Aug 2024
Viewed by 321
Abstract
With the gradual emergence of trends such as the asynchronous interconnection of power grids and the increasing penetration of renewable energy, the issues of ultra-low-frequency oscillations and low-frequency stability in power grids have become more prominent, posing serious challenges to the safety and [...] Read more.
With the gradual emergence of trends such as the asynchronous interconnection of power grids and the increasing penetration of renewable energy, the issues of ultra-low-frequency oscillations and low-frequency stability in power grids have become more prominent, posing serious challenges to the safety and stability of systems. The voltage-source converter-based HVDC (VSC-HVDC) interconnection is an effective solution to the frequency stability problems faced by regional power grids. VSC-HVDC can participate in system frequency stability control through a frequency limit controller (FLC). This paper first analyses the basic principles of how VSC-HVDC participates in system frequency stability control. Then, in response to the frequency stability control requirements of the sending and receiving power systems, a droop FLC strategy is designed. Furthermore, a multi-objective optimization method for the parameters of the droop FLC is proposed. Finally, a large-scale electromagnetic transient simulation model of the VSC-HVDC interconnected power system is constructed to verify the effectiveness of the proposed droop FLC method. Full article
Show Figures

Figure 1

Figure 1
<p>Double-ended VSC-HVDC System Structure.</p>
Full article ">Figure 2
<p>FLC Control Block Diagram.</p>
Full article ">Figure 3
<p>Slope characteristics of droop FLC of MMC-HVDC.</p>
Full article ">Figure 4
<p>Slope characteristics of the FLC controller with the limiter.</p>
Full article ">Figure 5
<p>Slope characteristics of the FLC controller after active and voltage limiting.</p>
Full article ">Figure 6
<p>Double-layer optimization framework for FLC of VSC-HVDC.</p>
Full article ">Figure 7
<p>Flowchart for solving the multi-objective optimization model of the initial layer of DC FLC.</p>
Full article ">Figure 8
<p>Flowchart for solving the multi-objective optimization model of the second layer of DC FLC.</p>
Full article ">Figure 9
<p>The topology of system S1.</p>
Full article ">Figure 10
<p>The topology of system S2.</p>
Full article ">Figure 11
<p>Frequency response of system S1.</p>
Full article ">Figure 12
<p>Frequency variation at bus 61.</p>
Full article ">Figure 13
<p>Topology of the interconnected system.</p>
Full article ">Figure 14
<p>The power curve of system S1.</p>
Full article ">Figure 15
<p>The power curve of system S2.</p>
Full article ">Figure 16
<p>Voltage profile of the system in steady state case.</p>
Full article ">Figure 17
<p>The frequency response curve of system S1 after parameter optimization.</p>
Full article ">Figure 18
<p>The active power curve of system S1.</p>
Full article ">Figure 19
<p>The active power curve of system S2.</p>
Full article ">Figure 20
<p>DC voltage profile of the system for low frequency problems.</p>
Full article ">Figure 21
<p>The frequency response curve of system S1 before parameter optimization.</p>
Full article ">Figure 22
<p>The frequency response curve of system S2 after parameter optimization.</p>
Full article ">Figure 23
<p>The frequency response curve of system S1 without parameter optimization.</p>
Full article ">Figure A1
<p>Voltage Curves of LCC2 sending and receiving ends.</p>
Full article ">Figure A2
<p>Active power curve for generator 3.</p>
Full article ">Figure A3
<p>Frequency response of system S1 for a power deficit of 600 MW.</p>
Full article ">Figure A4
<p>Voltage Curves of LCC5 sending and receiving ends.</p>
Full article ">Figure A5
<p>Voltage curves for substations 5 and 6.</p>
Full article ">
21 pages, 9889 KiB  
Article
Research on Multi-Source Data Fusion and Satellite Selection Algorithm Optimization in Tightly Coupled GNSS/INS Navigation Systems
by Xuyang Yu, Zhiming Guo and Liaoni Wu
Remote Sens. 2024, 16(15), 2804; https://doi.org/10.3390/rs16152804 - 31 Jul 2024
Viewed by 359
Abstract
With the increase in the number of Global Navigation Satellite System (GNSS) satellites and their operating frequencies, richer observation data are provided for the tightly coupled Global Navigation Satellite System/Inertial Navigation System (GNSS/INS). In this paper, we propose an efficient and robust combined [...] Read more.
With the increase in the number of Global Navigation Satellite System (GNSS) satellites and their operating frequencies, richer observation data are provided for the tightly coupled Global Navigation Satellite System/Inertial Navigation System (GNSS/INS). In this paper, we propose an efficient and robust combined navigation scheme to address the key issues of system accuracy, robustness, and computational efficiency. The tightly combined system fuses multi-source data such as the pseudo-range, the pseudo-range rate, and dual-antenna observations from the GNSS and the horizontal attitude angle from the vertical gyro (VG) in order to realize robust navigation in a sparse satellite observation environment. In addition, to cope with the high computational load faced by the system when the satellite observation conditions are good, we propose a weighted quasi-optimal satellite selection algorithm that reduces the computational burden of the navigation system by screening the observable satellites while ensuring the accuracy of the observation data. Finally, we comprehensively evaluate the proposed system through simulation experiments. The results show that, compared with the loosely coupled navigation system, our system has a significant improvement in state estimation accuracy and still provides reliable attitude estimation in regions with poor satellite observation conditions. In addition, in comparison experiments with the optimal satellite selection algorithm, our proposed satellite selection algorithm demonstrates greater advantages in terms of computational efficiency and engineering practicability. Full article
Show Figures

Figure 1

Figure 1
<p>A schematic diagram of the ECEF navigation coordinate frame and body coordinate frame. The X-, Y-, and Z-axes and text within the same framework are represented in the same color.</p>
Full article ">Figure 2
<p>The WQOSA execution flowchart, where N represents the current total number of visible satellites, and K represents the intended number to be selected.</p>
Full article ">Figure 3
<p>Pictures of installed experimental equipment.</p>
Full article ">Figure 4
<p>(<b>a</b>) A satellite map of the vehicle’s travel trajectory and (<b>b</b>) the variation curves of observable satellites and HDOP during the experiment.</p>
Full article ">Figure 5
<p>Performance comparison between LC and TC schemes. (<b>a</b>) Comparison chart of attitude errors in horizontal and vertical directions. (<b>b</b>) Comparison chart of velocity errors in eastward, northward, and vertical directions.</p>
Full article ">Figure 6
<p>Performance comparison between LC and TC schemes. Comparison chart of positional errors in the three spatial dimensions of longitude, latitude, and altitude.</p>
Full article ">Figure 7
<p>Performance comparison between with VG and without VG. (<b>a</b>) Comparison chart of attitude errors in horizontal and vertical directions. (<b>b</b>) Comparison chart of velocity errors in eastward, northward, and vertical directions. (<b>c</b>) Comparison chart of positional errors in the three spatial dimensions of longitude, latitude, and altitude.</p>
Full article ">Figure 8
<p>Performance comparison chart of attitude errors in horizontal directions with and without VG.</p>
Full article ">Figure 9
<p>(<b>a</b>) The vehicle’s travel trajectory and (<b>b</b>) the variation curves of observable satellites and the GDOP during the experiment.</p>
Full article ">Figure 10
<p>(<b>a</b>) The GDOP variation curve under various numbers of selected satellites, where the GDOP decreases as the number of selected satellites increases. (<b>b</b>) The curve showing the variation in the average best GDOP values for each number of selected satellites.</p>
Full article ">Figure 11
<p>(<b>a</b>) The GDOP curves for the OSA and WQOSA satellite selection, along with the GDOP difference curve. (<b>b</b>) The variation curves of the satellite selection time consumption for the OSA and WQOSA using the experimental data.</p>
Full article ">Figure 12
<p>Every satellite selected in the WQOSA and QOSA has a standard deviation of the pseudo-range error (psr-std). (<b>a</b>) By taking the average of all psr-std values, we plot the mean values of psr-std over the experimental time period in the chart. (<b>b</b>) Additionally, we also plot the mean curve of the navigation error standard deviation (nav-err) for the corresponding time intervals.</p>
Full article ">
30 pages, 12527 KiB  
Article
Strategic Siting of Direct Air Capture Facilities in the United States
by Jason Boerst, Ivonne Pena Cabra, Smriti Sharma, Connie Zaremsky and Arun K. S. Iyengar
Energies 2024, 17(15), 3755; https://doi.org/10.3390/en17153755 - 30 Jul 2024
Viewed by 601
Abstract
Direct air capture (DAC) systems that capture carbon dioxide (CO2) directly from the atmosphere are garnering considerable attention for their potential role as negative emission technologies in achieving net-zero CO2 emission goals. Common DAC technologies are based either on liquid–solvent [...] Read more.
Direct air capture (DAC) systems that capture carbon dioxide (CO2) directly from the atmosphere are garnering considerable attention for their potential role as negative emission technologies in achieving net-zero CO2 emission goals. Common DAC technologies are based either on liquid–solvent (L-DAC) or solid–sorbent (S-DAC) to capture CO2. A comprehensive multi-factor comparative economic analysis of the deployment of L-DAC and S-DAC facilities across various United States (U.S.) cities is presented in this paper. The analysis considers the influence of various factors on the favorability of DAC deployment, including local climatic conditions such as temperature, humidity, and CO2 concentrations; the availability of energy sources to power the DAC system; and costs for the transport and storage of the captured CO2 along with the consideration of the regional market and policy drivers. The deployment analysis in over 70 continental U.S. cities shows that L-DAC and S-DAC complement each other spatially, as their performance and operational costs vary in different climates. L-DAC is more suited to the hot, humid Southeast, while S-DAC is preferrable in the colder, drier Rocky Mountain region. Strategic deployment based on regional conditions and economics is essential for promoting the commercial adoptability of DAC, which is a critical technology to meet the CO2 reduction targets. Full article
(This article belongs to the Section B: Energy and Environment)
Show Figures

Figure 1

Figure 1
<p>Plant configuration for Case 1: NGCC-based L-DAC.</p>
Full article ">Figure 2
<p>Plant configurations for S-DAC Case 2a and Case 2b.</p>
Full article ">Figure 3
<p>Capture rate represented as a function of temperature (°C) and relative humidity (%) for L-DAC Case 1. Data source: An et al. [<a href="#B8-energies-17-03755" class="html-bibr">8</a>].</p>
Full article ">Figure 4
<p>Energy consumption represented as a function of the system’s capture rate for L-DAC Case 1. Data source: An et al. [<a href="#B8-energies-17-03755" class="html-bibr">8</a>].</p>
Full article ">Figure 5
<p>Productivity as a function of temperature and relative humidity for S-DAC Case 2. Data source: Sendi et al. [<a href="#B9-energies-17-03755" class="html-bibr">9</a>].</p>
Full article ">Figure 6
<p>Electricity requirement as a function of productivity under RH = 5–100% (labels correspond to RH value) for S-DAC Case 2. Data source: Sendi et al. [<a href="#B9-energies-17-03755" class="html-bibr">9</a>].</p>
Full article ">Figure 7
<p>Annual temperature and humidity by city using 2022 monthly averages.</p>
Full article ">Figure 8
<p>Decarbonization commitment (enforced or voluntary) at the state level in a future year. Data source: Anderson et al. [<a href="#B32-energies-17-03755" class="html-bibr">32</a>].</p>
Full article ">Figure 9
<p>State-level adoption of wind and solar technologies presented as the market share of wind and solar in the generation market.</p>
Full article ">Figure 10
<p>Comparative annual estimated gross DAC CO<sub>2</sub> capture rates (tCO<sub>2</sub>/yr) for the three cases considered: (<b>a</b>) Case 1: L-DAC; (<b>b</b>) Case 2a: S-DAC—NGCC and Case 2b: S-DAC—Grid-Powered (Case 2a and Case 2b have the same gross CO<sub>2</sub> capture rates (t-CO<sub>2</sub>/yr) and are represented in one single map. Their capacities are scaled up from the NETL DAC case studies baselines to match Case 1’s design capacity for comparison [<a href="#B12-energies-17-03755" class="html-bibr">12</a>]).</p>
Full article ">Figure 11
<p>Comparative annual estimated net DAC CO<sub>2</sub> removal rates (tCO<sub>2</sub>/yr) for the three cases considered: (<b>a</b>) Case 1: L-DAC; (<b>b</b>) Case 2a: S-DAC—NGCC; (<b>c</b>) Case 2b: S-DAC—Grid-Powered.</p>
Full article ">Figure 11 Cont.
<p>Comparative annual estimated net DAC CO<sub>2</sub> removal rates (tCO<sub>2</sub>/yr) for the three cases considered: (<b>a</b>) Case 1: L-DAC; (<b>b</b>) Case 2a: S-DAC—NGCC; (<b>c</b>) Case 2b: S-DAC—Grid-Powered.</p>
Full article ">Figure 12
<p>First year break-even cost of storage for saline storage formations for a nominal flow of 1.4 Mtpa of CO<sub>2</sub>. Currency is USD.</p>
Full article ">Figure 13
<p>Cost of carbon T&amp;S (USD/tCO<sub>2</sub>) for an assumed MFR of 1.4 Mtpa.</p>
Full article ">Figure 14
<p>Comparative estimated cost of DAC CO<sub>2</sub> removal (USD/tCO<sub>2</sub>) for the three cases considered. (<b>a</b>) Case 1: L-DAC; (<b>b</b>) Case 2a: S-DAC—NGCC; (<b>c</b>) Case 2b: S-DAC—Grid-Powered.</p>
Full article ">Figure 14 Cont.
<p>Comparative estimated cost of DAC CO<sub>2</sub> removal (USD/tCO<sub>2</sub>) for the three cases considered. (<b>a</b>) Case 1: L-DAC; (<b>b</b>) Case 2a: S-DAC—NGCC; (<b>c</b>) Case 2b: S-DAC—Grid-Powered.</p>
Full article ">Figure 15
<p>Impact of concentration on COCR.</p>
Full article ">Figure 16
<p>Policy and market score for 50 states and the District of Columbia used to develop the DAC favorability index.</p>
Full article ">Figure 17
<p>Composite favorability index for (<b>a</b>) L-DAC Case 1 and (<b>b</b>) S-DAC Cases 2a and 2b.</p>
Full article ">Figure 18
<p>Comparative assessment of overall favorability by technology type (L-DAC and S-DAC) with a NGCC power source. This comparative assessment includes Case 1 and Case 2a, as both use NGCC as the power source. Case 2b is excluded due to the difficulty in directly comparing NGCC-related emissions with grid-related emissions.</p>
Full article ">Figure A1
<p>Map of 319 saline reservoirs in FE/NETL CO<sub>2</sub> Saline Storage Cost Model (2017) geologic database. Source: NETL [<a href="#B18-energies-17-03755" class="html-bibr">18</a>,<a href="#B19-energies-17-03755" class="html-bibr">19</a>].</p>
Full article ">Figure A2
<p>Hypothetical T&amp;S infrastructure for potential DAC facilities across the U.S. (71 cities). Map created using [<a href="#B45-energies-17-03755" class="html-bibr">45</a>].</p>
Full article ">
30 pages, 12891 KiB  
Article
Evaluation of GPM IMERG Early, Late, and Final Run in Representing Extreme Rainfall Indices in Southwestern Iran
by Mohammad Sadegh Keikhosravi-Kiany and Robert C. Balling
Remote Sens. 2024, 16(15), 2779; https://doi.org/10.3390/rs16152779 - 30 Jul 2024
Viewed by 275
Abstract
The growing concerns about floods have highlighted the need for accurate and detailed precipitation data as extreme precipitation occurrences can lead to catastrophic floods, resulting in significant economic losses and casualties. Integrated Multi-satellitE Retrievals for the Global Precipitation Measurement (GPM IMERG) is a [...] Read more.
The growing concerns about floods have highlighted the need for accurate and detailed precipitation data as extreme precipitation occurrences can lead to catastrophic floods, resulting in significant economic losses and casualties. Integrated Multi-satellitE Retrievals for the Global Precipitation Measurement (GPM IMERG) is a commonly used high-resolution gridded precipitation dataset and is recognized as trustworthy alternative sources of precipitation data. The aim of this study is to comprehensively evaluate the performance of GPM IMERG Early (IMERG-E), Late (IMERG-L), and Final Run (IMERG-F) in precipitation estimation and their capability in detecting extreme rainfall indices over southwestern Iran during 2001–2020. The Asfezari gridded precipitation data, which are developed using a dense of ground-based observation, were utilized as the reference dataset. The findings indicate that IMERG-F performs reasonably well in capturing many extreme precipitation events (defined by various indices). All three products showed a better performance in capturing fixed and non-threshold precipitation indices across the study region. The findings also revealed that both IMERG-E and IMERG-L have problems in rainfall estimation over elevated areas showing values of overestimations. Examining the effect of land cover type on the accuracy of the precipitation products suggests that both IMERG-E and IMERG-L show large and highly unrealistic overestimations over inland water bodies and permanent wetlands. The results of the current study highlight the potential of IMERG-F as a valuable source of data for precipitation monitoring in the region. Full article
(This article belongs to the Special Issue Remote Sensing of Precipitation Extremes)
Show Figures

Figure 1

Figure 1
<p>General location of the study region in Northern Hemisphere (<b>a</b>) and topography of the region (<b>b</b>).</p>
Full article ">Figure 2
<p>Spatial values of (POD), (FAR), and (CSI) in the study area for IMERG-E (<b>a</b>,<b>d</b>,<b>g</b>), IMERG-L (<b>b</b>,<b>e</b>,<b>h</b>), and IMERG-F (<b>c</b>,<b>f</b>,<b>i</b>).</p>
Full article ">Figure 3
<p>Monthly amount of precipitation in the study area averaged over 2001–2020 derived from Asfezari, IMERG-E, IMERG-L, and IMERG-F.</p>
Full article ">Figure 4
<p>Scatterplots of precipitation for IMERG-E, IMERG-L, and IMERG-F compared to Asfezari for winter (<b>a</b>–<b>c</b>), spring (<b>d</b>–<b>f</b>), summer (<b>g</b>–<b>i</b>), and fall (<b>j</b>–<b>l</b>).</p>
Full article ">Figure 5
<p>Seasonal values of precipitation for winter (<b>a</b>–<b>d</b>), spring (<b>e</b>–<b>h</b>), summer (<b>i</b>–<b>l</b>), and fall (<b>m</b>–<b>p</b>) derived from Asfezari, IMERG-E, IMERG-L, and IMERG-F.</p>
Full article ">Figure 6
<p>Density-colored scatterplots of IMERG-E, IMERG-L, and IMERG-F against Asfezari, for winter (<b>a</b>–<b>c</b>), spring (<b>d</b>–<b>f</b>), summer (<b>g</b>–<b>i</b>), and fall (<b>j</b>–<b>l</b>) over the study region. The color represents the occurrence frequency.</p>
Full article ">Figure 7
<p>Annual values of precipitation derived from (<b>a</b>) IMERG-E, (<b>b</b>) IMERG-L, and (<b>c</b>) IMERG-F, and Asfezari (<b>d</b>) averaged over 2001–2020 and annual density-colored scatterplots of IMERG-E (<b>e</b>), IMERG-L (<b>f</b>), and IMERG-F (<b>g</b>) against Asfezari over the study region. The color represents the occurrence frequency.</p>
Full article ">Figure 8
<p>Land cover types of the study area derived from MCD12Q1.061 (<b>a</b>) along with focus on the permanent wet lands and inland water bodies (<b>b</b>) and Map of Google Earth image depicting earth surface features (<b>c</b>).</p>
Full article ">Figure 9
<p>The bias (%) of IMERG-E (<b>a</b>), IMERGE-L (<b>b</b>), and IMERG-F (<b>c</b>) against Asfezari for each of the elevation levels.</p>
Full article ">Figure 10
<p>Annual values of precipitation derived from (<b>a</b>) IMERG-E V07, (<b>b</b>) IMERG-L V07, and (<b>c</b>) IMERG-F V07, and Asfezari (<b>d</b>) averaged over 2001–2020 and annual density-colored scatterplots of IMERG-E V07 (<b>e</b>), IMERG-L V07 (<b>f</b>), and IMERG-F V07 (<b>g</b>) against Asfezari over the study region. The color represents the occurrence frequency.</p>
Full article ">Figure 11
<p>The bias (%) of IMERG-E V07 (<b>a</b>), IMERGE-L V07 (<b>b</b>), and IMERG-F V07 (<b>c</b>) against Asfezari for each of the elevation levels.</p>
Full article ">Figure 12
<p>Long-term means of fixed threshold extreme precipitation indices (R10, R20, CWD, CDD) for Asfezari (<b>a</b>,<b>e</b>,<b>i</b>,<b>m</b>), IMERG-E (<b>b</b>,<b>f</b>,<b>j</b>,<b>n</b>), IMERG-L(<b>c</b>,<b>g</b>,<b>k</b>,<b>o</b>), and IMERG-F (<b>d</b>,<b>h</b>,<b>l</b>,<b>p</b>).</p>
Full article ">Figure 13
<p>Density-colored scatterplots of extreme precipitation indices (R10, R20, CWD, CDD) for IMERG-E (<b>a</b>,<b>d</b>,<b>g</b>,<b>j</b>), IMERG-L (<b>b</b>,<b>e</b>,<b>h</b>,<b>k</b>), and IMERG-F (<b>c</b>,<b>f</b>,<b>i</b>,<b>l</b>) against the Asfezari dataset. The color represents the occurrence frequency.</p>
Full article ">Figure 14
<p>Long-term means of grid-related extreme precipitation indices (R95p, R99p, R95pTOT, and R99pTOT) for Asfezari (<b>a</b>,<b>e</b>,<b>i</b>,<b>m</b>), IMERG-E (<b>b</b>,<b>f</b>,<b>j</b>,<b>n</b>), IMERG-L(<b>c</b>,<b>g</b>,<b>k</b>,<b>o</b>), and IMERG-F (<b>d</b>,<b>h</b>,<b>l</b>,<b>p</b>).</p>
Full article ">Figure 15
<p>Density-colored scatterplots of extreme precipitation indices (R95p, R99p, R95pTOT, and R99pTOT) for IMERG-E (<b>a</b>,<b>d</b>,<b>g</b>,<b>j</b>), IMERG-L (<b>b</b>,<b>e</b>,<b>h</b>,<b>k</b>), and IMERG-F (<b>c</b>,<b>f</b>,<b>i</b>,<b>l</b>) against the Asfezari dataset. The color represents the occurrence frequency.</p>
Full article ">Figure 16
<p>Long-term means of non-threshold indices extreme precipitation indices (Rx1day (mm), SDII (mm), and PRCPTOT (mm)) for Asfezari (<b>a</b>,<b>e</b>,<b>i</b>), IMERG-E (<b>b</b>,<b>f</b>,<b>j</b>), IMERG-L(<b>c</b>,<b>g</b>,<b>k</b>), and IMERG-F (<b>d</b>,<b>h</b>,<b>l</b>).</p>
Full article ">Figure 17
<p>Density-colored scatterplots of extreme precipitation indices (Rx1day, SDII, PRCPTOT) for IMERG-E (<b>a</b>,<b>d</b>,<b>g</b>), IMERG-L (<b>b</b>,<b>e</b>,<b>h</b>), and IMERG-F (<b>c</b>,<b>f</b>,<b>i</b>) against the Asfezari dataset. The color represents the occurrence frequency.</p>
Full article ">Figure 18
<p>Temporal variation in the fixed threshold indices generated from Asfezari, IMERG-E, IMERG-L, and IMERG-F over 2001–2020.</p>
Full article ">Figure 19
<p>Temporal variation in the grid-related threshold indices generated from Asfezari, IMERG-E, IMERG-L, and IMERG-F over 2001–2020.</p>
Full article ">Figure 20
<p>Temporal variation in the non-threshold indices generated from Asfezari, IMERG-E, IMERG-L, and IMERG-F over 2001–2020.</p>
Full article ">Figure 21
<p>Relevant error box plots for Asfezari, IMERG-E, IMERG-L, and IMERG-F for fixed threshold indices (<b>a</b>–<b>d</b>), grid-related threshold indices (<b>e</b>–<b>h</b>), and non-threshold indices (<b>i</b>–<b>k</b>). The whiskers denote the maximum and minimum values in the data. The boxes extending from Q1 to Q3 show the median, while the red + symbols show outliers.</p>
Full article ">
18 pages, 2212 KiB  
Article
Improvements to a Crucial Budyko-Fu Parameter and Evapotranspiration Estimates via Vegetation Optical Depth over the Yellow River Basin
by Xingyi Wang and Jiaxin Jin
Remote Sens. 2024, 16(15), 2777; https://doi.org/10.3390/rs16152777 - 29 Jul 2024
Viewed by 366
Abstract
Against the backdrop of global warming and vegetation restoration, research on the evapotranspiration mechanism of the Yellow River basin has become a hot topic. The Budyko-Fu model is widely used to estimate basin-scale evapotranspiration, and its crucial parameter ω is used to characterize [...] Read more.
Against the backdrop of global warming and vegetation restoration, research on the evapotranspiration mechanism of the Yellow River basin has become a hot topic. The Budyko-Fu model is widely used to estimate basin-scale evapotranspiration, and its crucial parameter ω is used to characterize the underlying surface and climate characteristics of different basins. However, most studies only use factors such as the normalized difference vegetation index (NDVI), which represents the greenness of vegetation, to quantify the relationship between ω and the underlying surface, thereby neglecting richer vegetation information. In this study, we used long time-series multi-source remote sensing data from 1988 to 2015 and stepwise regression to establish dynamic estimation models of parameter ω for three subwatersheds of the upper Yellow River and quantify the contribution of underlying surface factors and climate factors to this parameter. In particular, vegetation optical depth (VOD) was introduced to represent plant biomass to improve the applicability of the model. The results showed that the dynamic estimation models of parameter ω established for the three subwatersheds were reasonable (R2 = 0.60, 0.80, and 0.40), and parameter ω was significantly correlated with the VOD and standardized precipitation evapotranspiration index (SPEI) in all watersheds. The dominant factors affecting the parameter in the different subwatersheds also differed, with underlying surface factors mainly affecting the parameter in the watershed before Longyang Gorge (BLG) (contributing 64% to 76%) and the watershed from Lanzhou to Hekou Town (LHT) (contributing 63% to 83%) and climate factors mainly affecting the parameter in the watershed from Longyang Gorge to Lanzhou (LGL) (contributing 75% to 93%). The results of this study reveal the changing mechanism of evapotranspiration in the Yellow River watershed and provide an important scientific basis for regional water balance assessment, global change response, and sustainable development. Full article
(This article belongs to the Special Issue Remote Sensing for Terrestrial Hydrologic Variables)
Show Figures

Figure 1

Figure 1
<p>Three subwatersheds of the upper Yellow River basin.</p>
Full article ">Figure 2
<p>Technology roadmap for the study.</p>
Full article ">Figure 3
<p>Trend charts of evapotranspiration (<b>a</b>), potential evapotranspiration (<b>b</b>), and precipitation (<b>c</b>) in the three subwatersheds of the upper Yellow River from 1988 to 2015 as well as the trend chart of parameter <math display="inline"><semantics> <mrow> <mi mathvariant="sans-serif">ω</mi> </mrow> </semantics></math> (<b>d</b>) after the moving average treatment.</p>
Full article ">Figure 4
<p>Distribution of parameter <math display="inline"><semantics> <mrow> <mi mathvariant="sans-serif">ω</mi> </mrow> </semantics></math> of the three subwatersheds on the Budyko curve.</p>
Full article ">Figure 5
<p>The variation trends of underlying surface factors VOD (<b>a</b>) and NDVI (<b>b</b>), as well as climate factors SPEI (<b>c</b>) and TMP (<b>d</b>) in the three subwatersheds from 1988 to 2015.</p>
Full article ">Figure 6
<p>Spearman correlation analysis heatmap between parameter <math display="inline"><semantics> <mrow> <mi mathvariant="sans-serif">ω</mi> </mrow> </semantics></math> and the respective variable factors in the BLG (<b>a</b>), LGL (<b>b</b>), and LHT (<b>c</b>), * represents significant correlation between variables.</p>
Full article ">Figure 7
<p>Residual plot of true and predicted values for parameter <math display="inline"><semantics> <mrow> <mi mathvariant="sans-serif">ω</mi> </mrow> </semantics></math> (<b>a</b>–<b>c</b>) and watershed evapotranspiration (<b>d</b>–<b>f</b>).</p>
Full article ">Figure 8
<p>Quantification of the contribution of factors to parameter <math display="inline"><semantics> <mrow> <mi mathvariant="sans-serif">ω</mi> </mrow> </semantics></math> using the standardized coefficient method (<b>a</b>) and R<sup>2</sup> decomposition method (<b>b</b>).</p>
Full article ">
16 pages, 5066 KiB  
Article
Analysis of a Rainstorm Process in Nanjing Based on Multi-Source Observational Data and Lagrangian Method
by Yuqing Mao, Youshan Jiang, Cong Li, Yi Shi and Daili Qian
Atmosphere 2024, 15(8), 904; https://doi.org/10.3390/atmos15080904 - 29 Jul 2024
Viewed by 300
Abstract
Using multi-source observation data including automatic stations, radar, satellite, new detection equipment, and the Fifth Generation European Centre for Medium-Range Weather Forecasts Reanalysis (ERA-5) data, along with the Hybrid Single-Particle Lagrangian Integrated Trajectory (HYSPLIT) platform, an analysis was conducted on a rainstorm process [...] Read more.
Using multi-source observation data including automatic stations, radar, satellite, new detection equipment, and the Fifth Generation European Centre for Medium-Range Weather Forecasts Reanalysis (ERA-5) data, along with the Hybrid Single-Particle Lagrangian Integrated Trajectory (HYSPLIT) platform, an analysis was conducted on a rainstorm process that occurred in Nanjing on 15 June 2020, with the aim of providing reference for future urban flood control planning and heavy rainfall forecasting and early warning. The results showed that this rainstorm process was generated under the background of an eastward-moving northeast cold vortex and a southward retreat of the Western Pacific Subtropical High. Intense precipitation occurred near the region of large top brightness temperature (TBB) gradient values or the center of low TBB values on the northern side of the convective cloud cluster. During the heavy precipitation period, the differential propagation phase shift rate (KDP), differential reflectivity factor (ZDR), and zero-lag correlation coefficient (ρHV) detected by the S-band dual-polarization radar all increased significantly. The vertical structure of the wind field detected by the wind profile radar provided a good indication of changes in precipitation intensity, showing a strong correspondence between the timing of maximum precipitation and the intrusion of upper-level cold air. The abrupt increase in the integrated liquid water content observed by the microwave radiometer can serve as an important indicator of the onset of stronger precipitation. During the Meiyu season in Nanjing, convective precipitation was mainly composed of small to medium raindrops with diameters less than 3 mm, with falling velocities of raindrops mainly clustering between 2 and 6 m·s−1. The rainstorm process featured four water vapor transport channels: the mid-latitude westerly channel, the Indian Ocean channel, the South China Sea channel, and the Pacific Ocean channel. During heavy rainfall, the Pacific Ocean water vapor channel was the main channel at the middle and lower levels, while the South China Sea water vapor channel was the main channel at the upper level, both accounting for a trajectory proportion of 34.2%. Full article
Show Figures

Figure 1

Figure 1
<p>(<b>a</b>) Distribution of precipitation in Nanjing from 20:00 on 14 June to 20:00 on 15 June 2020. (<b>b</b>) Variation of hourly rainfall of Jiangning Lulang Station from 00:00 to 20:00 on June 15 2020.</p>
Full article ">Figure 2
<p>Height field (unit: dagpm) and wind field (unit: m·s<sup>−1</sup>) at 08:00 on 15 June 2020. (<b>a</b>) 200 hPa, (<b>b</b>) 500 hPa, (<b>c</b>) 850 hPa. (The red shadow is the wind speed area of &gt;30 m·s<sup>−1</sup>, the blue shadow is the relative humidity area of &gt;60%, and the black dot is the location of Nanjing).</p>
Full article ">Figure 3
<p>TBB feature lines and heavy precipitation areas at 05:00, 06:00, 07:00, and 10:00 on 15 June 2020 (yellow line: −32 °C, orange line: −52 °C, purple line: −62 °C, red line: −72 °C; Blue circle: hourly rainfall ≥20 mm, red circle: hourly rainfall ≥40 mm).</p>
Full article ">Figure 4
<p>The time-height profile of Z (unit: dBz), K<sub>DP</sub> (unit: °/km), Z<sub>DR</sub> (unit: dB), and <math display="inline"><semantics> <mrow> <msub> <mrow> <mi mathvariant="sans-serif">ρ</mi> </mrow> <mrow> <mi mathvariant="normal">H</mi> <mi mathvariant="normal">V</mi> </mrow> </msub> </mrow> </semantics></math> of Jiangning Lulang Station from 06:00 to 12:00 on 15 June 2020.</p>
Full article ">Figure 5
<p>The skew T–logp diagram and vertical wind profile of the Nanjing radiosonde at 20:00 14 June 2020 (<b>a</b>) and 08:00 15 June 2020 (<b>b</b>). The dark, black, solid lines indicate the state profiles. The red solid lines represent the temperature stratification profiles. The green solid lines are the dew point stratification profiles. The blue dashed lines represent the wet adiabatic profiles. The orange dashed lines represent the dry adiabatic profiles. The green dashed lines represent the iso-saturated specific humidity profiles. The gray-sloping solid lines represent the isothermal profiles. The cyan-blue dashed lines represent the 0 °C isotherm profiles. The red area represents CAPE.</p>
Full article ">Figure 6
<p>The evolution of the wind profile at Nanjing Station from 03:00 to 15:00 on 15 June 2020.</p>
Full article ">Figure 7
<p>The profiles of water vapor density (unit: g·m<sup>−3</sup>) and liquid water content (unit: g·m<sup>−3</sup>) of the microwave radiometer and hourly rainfall (unit: mm) at Nanjing Station from 00:00 on 15 June to 00:00 on 16 June 2020.</p>
Full article ">Figure 8
<p>The variation over time of raindrop concentration (lgN(D)), raindrop diameter (unit: mm), and falling velocity (unit: m·s<sup>−1</sup>) at Nanjing Station from 00:00 on 15 June to 00:00 on 16 June 2020.</p>
Full article ">Figure 9
<p>The trace distribution of the whole layer and the upper, middle and lower levels of the air masses on the 5th day of the backward tracking at 20:00 on 14 June, 08:00 and 20:00 on 15 June 2020.</p>
Full article ">Figure 10
<p>The trajectories from four water vapor sources at 20:00 on June 14, 08:00, and 20:00 on 15 June 2020. The color of the track indicates the specific humidity (g·kg<sup>−1</sup>), and the thickness of the track indicates the track intensity.</p>
Full article ">Figure 11
<p>Waterbody distribution in Nanjing.</p>
Full article ">
19 pages, 4131 KiB  
Article
Emission Control in Expressway Systems: Vehicle Emission Inventory and Policy Scenario Analysis
by Jingxu Chen, Junyi Chen, Dawei Chen and Xiuyu Shen
Systems 2024, 12(8), 273; https://doi.org/10.3390/systems12080273 - 29 Jul 2024
Viewed by 493
Abstract
Expressway systems play a vital role in facilitating intercity travels for both passengers and freights, which are also a significant source of vehicle emissions within the transportation sector. This study investigates vehicle emissions from expressway systems using the COPERT model to develop multi-year [...] Read more.
Expressway systems play a vital role in facilitating intercity travels for both passengers and freights, which are also a significant source of vehicle emissions within the transportation sector. This study investigates vehicle emissions from expressway systems using the COPERT model to develop multi-year emission inventories for different pollutants, covering the past and future trends from 2005 to 2030. Thereinto, an integrated SARIMA-SVR method is designed to portray the temporal variation of vehicle population, and the possible future trends of expressway vehicle emissions are predicted through policy scenario analysis. The Jiang–Zhe–Hu Region of China is taken as the case study to analyze emission control in expressway systems. The results indicate that (1) carbon monoxide (CO) and volatile organic compounds (VOCs) present a general upward trend primarily originating from passenger vehicles, while nitrogen oxides (NOx) and inhalable particles (PM) display a slowing upward trend with fluctuations mainly sourcing from freight vehicles; (2) vehicle population constraint is an effective emission control policy, but upgrading the medium- and long-haul transportation structure is necessary to meet the continuous growth of intercity trips. Expressway vehicle emission reduction effectiveness can be further enhanced by curtailing the update frequency of emission standards, along with the scrapping of high-emission vehicles. Full article
Show Figures

Figure 1

Figure 1
<p>The expressway network of the JZH region.</p>
Full article ">Figure 2
<p>Compositions of various vehicle emission standards from 2005 to 2020: (<b>a</b>) passenger vehicles; and (<b>b</b>) freight vehicles.</p>
Full article ">Figure 3
<p>Flowchart of the integrated SARIMA-SVR method.</p>
Full article ">Figure 4
<p>Research framework.</p>
Full article ">Figure 5
<p>Expressway VP variation in the JZH region from 2005 to 2020.</p>
Full article ">Figure 6
<p>Total emission trends by vehicle types in the JZH region from 2005 to 2020.</p>
Full article ">Figure 7
<p>Emission compositions of eight vehicle types in the JZH region between 2005 and 2020: (<b>a</b>) CO; (<b>b</b>) VOC; (<b>c</b>) NOx; and (<b>d</b>) PM2.5.</p>
Full article ">Figure 8
<p>Expressway vehicle emission trends of the JZH region between 2021 and 2030 under nine scenarios: (<b>a</b>) CO; (<b>b</b>) VOC; (<b>c</b>) NOx; and (<b>d</b>) PM2.5.</p>
Full article ">
17 pages, 11358 KiB  
Article
Fiduciary-Free Frame Alignment for Robust Time-Lapse Drift Correction Estimation in Multi-Sample Cell Microscopy
by Stefan Baar, Masahiro Kuragano, Naoki Nishishita, Kiyotaka Tokuraku and Shinya Watanabe
J. Imaging 2024, 10(8), 181; https://doi.org/10.3390/jimaging10080181 - 29 Jul 2024
Viewed by 500
Abstract
When analyzing microscopic time-lapse observations, frame alignment is an essential task to visually understand the morphological and translation dynamics of cells and tissue. While in traditional single-sample microscopy, the region of interest (RoI) is fixed, multi-sample microscopy often uses a single microscope that [...] Read more.
When analyzing microscopic time-lapse observations, frame alignment is an essential task to visually understand the morphological and translation dynamics of cells and tissue. While in traditional single-sample microscopy, the region of interest (RoI) is fixed, multi-sample microscopy often uses a single microscope that scans multiple samples over a long period of time by laterally relocating the sample stage. Hence, the relocation of the optics induces a statistical RoI offset and can introduce jitter as well as drift, which results in a misaligned RoI for each sample’s time-lapse observation (stage drift). We introduce a robust approach to automatically align all frames within a time-lapse observation and compensate for frame drift. In this study, we present a sub-pixel precise alignment approach based on recurrent all-pairs field transforms (RAFT); a deep network architecture for optical flow. We show that the RAFT model pre-trained on the Sintel dataset performed with near perfect precision for registration tasks on a set of ten contextually unrelated time-lapse observations containing 250 frames each. Our approach is robust for elastically undistorted and translation displaced (x,y) microscopic time-lapse observations and was tested on multiple samples with varying cell density, obtained using different devices. The approach only performed well for registration and not for tracking of the individual image components like cells and contaminants. We provide an open-source command-line application that corrects for stage drift and jitter. Full article
Show Figures

Figure 1

Figure 1
<p><b>Image stabilization:</b> An input containing a set of sequential images, where n is the frame number. The displacement information of the two frames is computed by first using recurrent all-pairs field transforms (RAFT) [<a href="#B16-jimaging-10-00181" class="html-bibr">16</a>] trained on the Sintel dataset [<a href="#B17-jimaging-10-00181" class="html-bibr">17</a>] to estimate the vector field describing the apparent motion (translation) of each pixel. The transversal displacement between frames (stabilized time-lapse observation) is determined by computing the median vector of the estimated vector field (pink arrows) and therefore eliminating the group motion of a set of sparse objects (blue arrows) within the RoI.</p>
Full article ">Figure 2
<p><b>Lateral RoI displacement synthesis:</b> (<b>Left</b>): time evolution sample of the instrumental jitter exhibited by Incucyte SX1. (<b>Right</b>): Randomly generated lateral jitter, introduced to unperturbed time-lapse observations obtained with an inverted Nikon Ti-E microscope.</p>
Full article ">Figure 3
<p><b>Implementation overview:</b> Schematic of the data processing procedure used to correct frame-to-frame jitter is presented from left to right. For a set of frames, each frame and its following frame are compared and its dense vector filed is computed via optical flow. Notice that individual objects within the field of view (FOV) can move in random directions with random velocities (displacement amplitudes), independently of the underlying group motion (<math display="inline"><semantics> <mrow> <mi>d</mi> <mi>x</mi> </mrow> </semantics></math>, <math display="inline"><semantics> <mrow> <mi>d</mi> <mi>y</mi> </mrow> </semantics></math>). However, from the displacement vector histograms, one can identify a single peak that for each direction (<math display="inline"><semantics> <mrow> <mo>Δ</mo> <mi>x</mi> </mrow> </semantics></math>, <math display="inline"><semantics> <mrow> <mo>Δ</mo> <mi>y</mi> </mrow> </semantics></math>) characterized the displacement direction and amplitude, which is best characterized by the median. The median of the displacement field is then used to perform an affine transformation (in x and y) to correct each frame.</p>
Full article ">Figure 4
<p><b>Framing Methods:</b> An example set of three frames is presented. The four most plausible framing methods (<b>A</b>–<b>D</b>) are displayed. Maximum framing contains the information of all frames.</p>
Full article ">Figure 5
<p><b>Displacement dispersion:</b> Estimated with recurrent all-pairs field transforms and the median of the resulting vector field for ten sample observations without an image lock-plate (no position feedback) are presented on the (<b>Left</b>). Seven samples utilizing an image lock-plate and instrumental feedback loop are presented on the (<b>Right</b>). The horizontal and vertical displacement dispersions are colored in orange and blue, respectively.</p>
Full article ">Figure 6
<p><b>Jitter, drift, and correction amplitude:</b> for two samples 230208E3-3_RA and 230208F3-1_KNK808v2 at the top and bottom, respectively. (<b>Left</b>): NN displacement amplitude per frame. (<b>Center</b>): The accumulated displacement, characterizing drift, where the red and black lines denote the median of the x and y components, respectively. (<b>Right</b>): The correction amplitude for each frame, determined by the accumulated displacement and its median. 230208E3-3_RA exhibits a relatively low drift in both directions (<b>upper center</b>), but a high frame-to-frame displacement amplitude in the x direction. 230208F3-1_KNK808v2 exhibits strong drift (<b>lower center</b>) and a lower frame-to-frame displacement amplitude for both directions.</p>
Full article ">Figure 7
<p><b>Displacement dispersion:</b> for cross-correlation and the dense optical flow methods, based on Lucas–Kanade and RAFT. (<b>Left</b>): for the Nearest Neighbor (NN) displacement for each n and n + 1 frame. (<b>Right</b>): for all frame permutations of the entire observation.</p>
Full article ">Figure 8
<p><b>Time distance-dependent offset error:</b> Distance matrix for all frame permutations of the vertical (<b>left</b>) and horizontal (<b>right</b>) components for the RAFT-based stabilization approach.</p>
Full article ">Figure 9
<p><b>Offset error in relation to the frame distance:</b> for various sampling ranges (orange to blue), for the vertical (<b>left</b>) and horizontal (<b>right</b>) components, utilizing RAFT stabilization. The positive vertical axis presents the correspondence between the nth and n+1st frame, the negative inverts the frame order and represents the n-1st and nth frame. The green circle indicate patterns that break point symmetry.</p>
Full article ">Figure 10
<p><b>Free parameter dependence:</b> (<b>Left</b>): Estimation error in relation to the upscaling factor of the phased cross-correlation-based stabilization approach. (<b>Center</b>): Estimation error for Lucas–Kanade-based stabilization, which mainly depends on the radius. (<b>Right</b>): RAFT-based stabilization error depending on the number of iterations.</p>
Full article ">Figure 11
<p><b>Displacement analysis:</b> of sample 230208E3-3 RA. (<b>A</b>) shows the first (<b>top</b>) and second frames (<b>bottom</b>). (<b>B</b>) shows the corresponding displacement maps produced by RAFT. Regions R1 to R3 indicate cells that exhibited modality but were not registered in the (<b>B</b>). (<b>C</b>) Histograms corresponding to B for both RAFT (blue) and Lucas–Kanade (orange). The corresponding median and mode are indicated as solid and dashed lines, respectively. The black arrows near the color bars of (<b>B</b>) indicate the x-component of the median presented in the figures in (<b>C</b>).</p>
Full article ">Figure 12
<p><b>CPU-GPU precision:</b> The residual of 32 bit/16 bit precision is plotted for each frame pare.</p>
Full article ">Figure A1
<p><b>Benchmark:</b> The average processing time for a set of image sizes and each image stabilization approach is presented on the left-hand side. The right-hand side presents the associated standard deviation for one thousand samples each.</p>
Full article ">Figure A2
<p><b>Jitter, drift, and correction amplitude:</b> for all samples presented in <a href="#jimaging-10-00181-f005" class="html-fig">Figure 5</a> (<b>Left</b>) for samples not utilizing an image lock-plate. (<b>Left</b>): NN displacement amplitude per frame. (<b>Center</b>): The accumulated displacement, characterizing drift, where the red and black lines denote the median of the x and y component, respectively. (<b>Right</b>): The correction amplitude for each frame, determined by the accumulated displacement and its median, in both directions.</p>
Full article ">Figure A3
<p><b>Jitter, drift, and correction amplitude:</b> Continuation of <a href="#jimaging-10-00181-f0A2" class="html-fig">Figure A2</a>.</p>
Full article ">Figure A4
<p><b>Jitter, drift, and correction amplitude:</b> For all samples presented in <a href="#jimaging-10-00181-f005" class="html-fig">Figure 5</a> (right for samples utilizing an image lock-plate. (<b>Left</b>): NN displacement amplitude per frame. (<b>Center</b>): The accumulated displacement, characterizing drift, where the red and black lines denote the median of the x and y components, respectively. (<b>Right</b>): The correction amplitude for each frame, determined by the accumulated displacement and its median, in both directions.</p>
Full article ">Figure A5
<p><b>Jitter, drift, and correction amplitude:</b> Continuation of <a href="#jimaging-10-00181-f0A4" class="html-fig">Figure A4</a>.</p>
Full article ">
15 pages, 8591 KiB  
Article
Regional Winter Wheat Yield Prediction and Variable Importance Analysis Based on Multisource Environmental Data
by Hao Xu, Hongfei Yin, Yaohui Liu, Biao Wang, Hualu Song, Zhaowen Zheng, Xiaohu Zhang, Li Jiang and Shuai Wang
Agronomy 2024, 14(8), 1623; https://doi.org/10.3390/agronomy14081623 - 24 Jul 2024
Viewed by 584
Abstract
Timely and accurate predictions of winter wheat yields are key to ensuring food security. In this research, winter wheat yield prediction models for six provinces were established using a random forest (RF) model. Two methods were employed to analyze feature variables. RF partial [...] Read more.
Timely and accurate predictions of winter wheat yields are key to ensuring food security. In this research, winter wheat yield prediction models for six provinces were established using a random forest (RF) model. Two methods were employed to analyze feature variables. RF partial dependence plots were generated to demonstrate the nonlinear relationships between the feature variables and yield, and bivariate Moran’s I was considered to identify the spatial associations between variables. Results showed that when environmental data from key growth periods were used for prediction model establishment, the root mean square error (RMSE) varied between 200 and 700 kg/ha, and the coefficient of determination (R2) exceeded 0.5. Feature variable analysis results indicated that the longitude, latitude, topography and normalized difference vegetation index (NDVI) were important variables. Below the threshold, the yield gradually increased with increasing NDVI. Bivariate Moran’s I results showed that there was zonal distribution of meteorological elements. Within a large spatial range, the change in environmental variables due to the latitude and longitude should be accounted for in modeling, but the influence of collinearity between the feature variables should be eliminated via variable importance analysis. Full article
(This article belongs to the Section Precision and Digital Agriculture)
Show Figures

Figure 1

Figure 1
<p>Research area: (<b>a</b>) geographical location, (<b>b</b>) DEM, and (<b>c</b>) spatial distribution of the mean statistical yield values (2014–2019) in the six provinces; for each color, the basic spatial unit is the county.</p>
Full article ">Figure 2
<p>Analysis workflow. The index i indicates the months adopted in this research, with 3, 4 and 5 corresponding to March, April and May, respectively.</p>
Full article ">Figure 3
<p>Box plot of the statistical yield. The upper line of the box is the maximum yield, the midline is the median yield, the lower line is the minimum yield, and the black triangle denotes an outlier.</p>
Full article ">Figure 4
<p>Relationship between the statistical yield and RF-simulated yield.</p>
Full article ">Figure 5
<p>As the feature variables changed and the yield variation exceeded 200 kg/ha, the partial dependence plots showed nonlinear relationships between the feature variables and yield.</p>
Full article ">Figure 6
<p>Radar diagram of the bivariate Moran’s I values between the yield and the feature variables.</p>
Full article ">Figure 7
<p>Radar diagram of the bivariate Moran’s I values between the longitude and environmental variables.</p>
Full article ">Figure 8
<p>Radar diagram of the bivariate Moran’s I values between the latitude and environmental variables.</p>
Full article ">Figure 9
<p>Changes in the nonlinear relationship between the environmental variables and yield using partial dependence plots before (<b>upper</b>) and after (<b>lower</b>) the deletion of the latitude and longitude.</p>
Full article ">
Back to TopTop