[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
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
remove_circle_outline

Search Results (3,468)

Search Parameters:
Keywords = natural river

Order results
Result details
Results per page
Select all
Export citation of selected articles as:
24 pages, 5994 KiB  
Article
Mapping Natural Populus euphratica Forests in the Mainstream of the Tarim River Using Spaceborne Imagery and Google Earth Engine
by Jiawei Zou, Hao Li, Chao Ding, Suhong Liu and Qingdong Shi
Remote Sens. 2024, 16(18), 3429; https://doi.org/10.3390/rs16183429 (registering DOI) - 15 Sep 2024
Abstract
Populus euphratica is a unique constructive tree species within riparian desert areas that is essential for maintaining oasis ecosystem stability. The Tarim River Basin contains the most densely distributed population of P. euphratica forests in the world, and obtaining accurate distribution data in [...] Read more.
Populus euphratica is a unique constructive tree species within riparian desert areas that is essential for maintaining oasis ecosystem stability. The Tarim River Basin contains the most densely distributed population of P. euphratica forests in the world, and obtaining accurate distribution data in the mainstream of the Tarim River would provide important support for its protection and restoration. We propose a new method for automatically extracting P. euphratica using Sentinel-1 and 2 and Landsat-8 images based on the Google Earth Engine cloud platform and the random forest algorithm. A mask of the potential distribution area of P. euphratica was created based on prior knowledge to save computational resources. The NDVI (Normalized Difference Vegetation Index) time series was then reconstructed using the preferred filtering method to obtain phenological parameter features, and the random forest model was input by combining the phenological parameter, spectral index, textural, and backscattering features. An active learning method was employed to optimize the model and obtain the best model for extracting P. euphratica. Finally, the map of natural P. euphratica forests with a resolution of 10 m in the mainstream of the Tarim River was obtained. The overall accuracy, producer’s accuracy, user’s accuracy, kappa coefficient, and F1-score of the map were 0.96, 0.98, 0.95, 0.93, and 0.96, respectively. The comparison experiments showed that simultaneously adding backscattering and textural features improved the P. euphratica extraction accuracy, while textural features alone resulted in a poor extraction effect. The method developed in this study fully considered the prior and posteriori information and determined the feature set suitable for the P. euphratica identification task, which can be used to quickly obtain accurate large-area distribution data of P. euphratica. The method can also provide a reference for identifying other typical desert vegetation. Full article
Show Figures

Figure 1

Figure 1
<p>Geographical location of the study area and the distribution of sample points. (<b>a</b>): location of the study area in Xinjiang province in China; (<b>b</b>): training dataset distribution; (<b>c</b>): detailed sample area showing <span class="html-italic">P. euphratica</span> and non–<span class="html-italic">P. euphratica</span> in a Sentinel-2 false-color image.</p>
Full article ">Figure 2
<p>Distribution of validation dataset. The black solid line represents the range of the study area; the red and yellow points represent <span class="html-italic">P. euphratica</span> and non–<span class="html-italic">P. euphratica</span>, respectively.</p>
Full article ">Figure 3
<p>Workflow of the research.</p>
Full article ">Figure 4
<p>Threshold segmentation effect of MNDWI and NDVI. (<b>a</b>): false color image of Jieran Lik Reservoir in Xinjiang Province; (<b>b</b>): statistical result of the corresponding frequency distribution of MNDWI values of water and other ground objects in area (<b>a</b>); (<b>c</b>): false color image of Pazili Tamu in Xinjiang; (<b>d</b>): statistical result for the corresponding frequency distribution of NDVI values of desert bare land and other ground objects in region (<b>c</b>).</p>
Full article ">Figure 5
<p>Comparison of NDVI data before and after spatiotemporal fusion: (<b>a</b>) NDVI data derived from Sentinel-2 before fusion, (<b>b</b>) NDVI data after fusion.</p>
Full article ">Figure 6
<p>Comparison of the effects of different filter functions for: (<b>a</b>) <span class="html-italic">P. euphratica</span>; (<b>b</b>) <span class="html-italic">Tamarix</span>; (<b>c</b>) allee tree; (<b>d</b>) farmland; (<b>e</b>) wetland; (<b>f</b>) urban tree.</p>
Full article ">Figure 7
<p>Comparison between phenological curves of six typical vegetation species. Phenology parameters of (<b>a</b>) <span class="html-italic">P. euphratica</span>, (<b>b</b>) <span class="html-italic">Tamarix</span>, (<b>c</b>) allee tree, (<b>d</b>) farmland, (<b>e</b>) wetland, and (<b>f</b>) urban tree.</p>
Full article ">Figure 8
<p>Importance of different features in the RF classification.</p>
Full article ">Figure 9
<p>Natural <span class="html-italic">P. euphratica</span> forest maps extracted using four feature combinations: (<b>a</b>) PS, (<b>b</b>) PSB, (<b>c</b>) PST, and (<b>d</b>) PSBT.</p>
Full article ">Figure 10
<p>Comparison of <span class="html-italic">P. euphratica</span> extraction results using different feature combinations on Sentinel-2 standard false color images. Rows 1 to 4 show the identification of <span class="html-italic">P. euphratica</span> in desert areas, <span class="html-italic">P. euphratica</span>-dense areas, agricultural areas, and large river areas, respectively. The green area represents the classification result of <span class="html-italic">P. euphratica</span>. The yellow circle corresponding to each row is the area where the extraction results of different feature combinations are quite different.</p>
Full article ">Figure 11
<p>(<b>a</b>) Distribution of natural <span class="html-italic">P. euphratica</span> forest in the mainstream of the Tarim River. (<b>b</b>): UAV image of healthy <span class="html-italic">P. euphratica</span>, (<b>c</b>): classification result of healthy <span class="html-italic">P. euphratica</span>, (<b>d</b>): UAV image of unhealthy <span class="html-italic">P. euphratica</span>, (<b>e</b>): classification result of unhealthy <span class="html-italic">P. euphratica</span>, (<b>f</b>): UAV image of dense <span class="html-italic">P. euphratica</span>, (<b>g</b>): classification result of dense <span class="html-italic">P. euphratica</span>, (<b>h</b>): UAV image of sparse <span class="html-italic">P. euphratica</span>, (<b>i</b>): classification result of sparse <span class="html-italic">P. euphratica</span>. The green area represents the classification results of <span class="html-italic">P. euphratica</span>.</p>
Full article ">Figure 12
<p>Mixed pixel problems associated with <span class="html-italic">P. euphratica</span>: (<b>a</b>) <span class="html-italic">P. euphratica</span> occupying less than one pixel; (<b>b</b>) sandy soil interfering with the reflected signal of <span class="html-italic">P. euphratica</span>. The red box represents a pixel on the images for clearer observation. Basemaps of row 1-2 are UAV images while row 3 are Sentinel-2 standard false color images.</p>
Full article ">
29 pages, 5760 KiB  
Review
Concentrations of Organochlorine, Organophosphorus, and Pyrethroid Pesticides in Rivers Worldwide (2014–2024): A Review
by Acela López-Benítez, Alfredo Guevara-Lara, Miguel A. Domínguez-Crespo, José A. Andraca-Adame and Aidé M. Torres-Huerta
Sustainability 2024, 16(18), 8066; https://doi.org/10.3390/su16188066 (registering DOI) - 15 Sep 2024
Viewed by 18
Abstract
The extensive use of pesticides has led to the contamination of natural resources, sometimes causing significant and irreversible damage to the environment and human health. Even though the use of many pesticides is banned, these compounds are still being found in rivers worldwide. [...] Read more.
The extensive use of pesticides has led to the contamination of natural resources, sometimes causing significant and irreversible damage to the environment and human health. Even though the use of many pesticides is banned, these compounds are still being found in rivers worldwide. In this review, 205 documents have been selected to provide an overview of pesticide contamination in rivers over the last 10 years (2014–2024). After these documents were examined, information of 47 river systems was organized according to the types of pesticides most frequently detected, including organochloride, organophosphorus, and pyrethroid compounds. A total of 156 compounds were classified, showing that 46% of these rivers contain organochlorine compounds, while 40% exhibit organophosphorus pesticides. Aldrin, hexachlorocyclohexane, and endosulfan were the predominant organochlorine pesticides with concentration values between 0.4 and 37 × 105 ng L−1. Chlorpyrifos, malathion, and diazinon were the main organophosphorus pesticides with concentrations between 1 and 11 × 105 ng L−1. Comparing the pesticide concentrations with standard guidelines, we found that the Ganga River in India (90 ng L−1), the Owan and Okura Rivers in Nigeria (210 and 9 × 103 ng L−1), and the Dong Nai River in Vietnam (68 ng L−1) exceed the permissible levels of aldrin (30 ng L−1). Full article
(This article belongs to the Section Hazards and Sustainability)
Show Figures

Figure 1

Figure 1
<p>Research methodology followed in this review.</p>
Full article ">Figure 2
<p>Classification of the documents mentioned in this review.</p>
Full article ">Figure 3
<p>Chemical structures and LD<sub>50</sub> values of some organochlorine pesticides.</p>
Full article ">Figure 4
<p>Concentrations of organochlorine pesticides detected in rivers around the world: (<b>a</b>) highest concentration values; (<b>b</b>) intermediate concentration values.</p>
Full article ">Figure 5
<p>Chemical structures and LD<sub>50</sub> values of some organophosphorus pesticides.</p>
Full article ">Figure 6
<p>Concentrations of organophosphorus pesticides detected in rivers around the world: (<b>a</b>) highest concentration values; (<b>b</b>) intermediate concentration values.</p>
Full article ">Figure 7
<p>Chemical structures and LD<sub>50</sub> values of some pyrethroid pesticides.</p>
Full article ">Figure 8
<p>Concentrations of pyrethroid pesticides detected in rivers around the world: (<b>a</b>) highest concentration values; (<b>b</b>) intermediate concentration values.</p>
Full article ">Figure 9
<p>Occurrence of pesticides in rivers worldwide (2014–2024).</p>
Full article ">Figure 10
<p>Geographical distribution of rivers across continents (2014–2024).</p>
Full article ">Figure 11
<p>Summary of the pesticides detected in the waters of 47 rivers from 2014 to 2024.</p>
Full article ">Figure 12
<p>Pesticide concentrations in rivers, considering the reported sampling year. (A: aldrin, AC: acephate, BHC: benzene hexachloride, CP: chlorpyrifos, CY: cypermethrin, D: deltamethrin, DVC: dichlorvos, DN: dieldrin, DT: dichlorodiphenyltrichloroethane, DZ: diazinon, E: endosulfan, EN: endrin, G: glyphosate, H: heptachlor, HCH: hexachlorocyclohexane, M: malathion, MT: metamidophos, P: parathion, PF: profenofos, Q: quinalphos, and T: triazophos).</p>
Full article ">
22 pages, 3621 KiB  
Article
Estimating Non-Stationary Extreme-Value Probability Distribution Shifts and Their Parameters Under Climate Change Using L-Moments and L-Moment Ratio Diagrams: A Case Study of Hydrologic Drought in the Goat River Near Creston, British Columbia
by Isaac Dekker, Kristian L. Dubrawski, Pearce Jones and Ryan MacDonald
Hydrology 2024, 11(9), 154; https://doi.org/10.3390/hydrology11090154 (registering DOI) - 14 Sep 2024
Viewed by 168
Abstract
Here, we investigate the use of rolling-windowed L-moments (RWLMs) and L-moment ratio diagrams (LMRDs) combined with a Multiple Linear Regression (MLR) machine learning algorithm to model non-stationary low-flow hydrological extremes with the potential to simultaneously understand time-variant shape, scale, location, and probability distribution [...] Read more.
Here, we investigate the use of rolling-windowed L-moments (RWLMs) and L-moment ratio diagrams (LMRDs) combined with a Multiple Linear Regression (MLR) machine learning algorithm to model non-stationary low-flow hydrological extremes with the potential to simultaneously understand time-variant shape, scale, location, and probability distribution (PD) shifts under climate change. By employing LMRDs, we analyse changes in PDs and their parameters over time, identifying key environmental predictors such as lagged precipitation for September 5-day low-flows. Our findings indicate a significant relationship between total August precipitation L-moment ratios (LMRs) and September 5-day low-flow LMRs (τ2-Precipitation and τ2-Discharge: R2 = 0.675, p-values < 0.001; τ3-Precipitation and τ3-Discharge: R2 = 0.925, p-value for slope < 0.001, intercept not significant with p = 0.451, assuming α = 0.05 and a 31-year RWLM), which we later refine and use for prediction within our MLR algorithm. The methodology, applied to the Goat River near Creston, British Columbia, aids in understanding the implications of climate change on water resources, particularly for the yaqan nuʔkiy First Nation. We find that future low-flows under climate change will be outside the Natural Range of Variability (NROV) simulated from historical records (assuming a constant PD). This study provides insights that may help in adaptive water management strategies necessary to help preserve Indigenous cultural rights and practices and to help sustain fish and fish habitat into the future. Full article
(This article belongs to the Section Hydrological and Hydrodynamic Processes and Modelling)
Show Figures

Figure 1

Figure 1
<p>L-moment ratio diagrams (LMRDs) for: (<b>a</b>) August total precipitation (mm) and (<b>b</b>) September 5-day low-flow (<math display="inline"><semantics> <mrow> <msup> <mi mathvariant="normal">m</mi> <mn>3</mn> </msup> <mo>/</mo> <mi mathvariant="normal">s</mi> </mrow> </semantics></math>). Each panel includes the L-Coefficient of Variation (L-CV;<math display="inline"><semantics> <msub> <mi>τ</mi> <mn>2</mn> </msub> </semantics></math>) versus L-skewness (<math display="inline"><semantics> <msub> <mi>τ</mi> <mn>3</mn> </msub> </semantics></math>) and L-kurtosis versus L-skewness ratios.</p>
Full article ">Figure 2
<p>Relationship between L-moment ratios (LMRs) of August total precipitation (mm) and September 5-day low-flow (<math display="inline"><semantics> <mrow> <msup> <mi mathvariant="normal">m</mi> <mn>3</mn> </msup> <mo>/</mo> <mi mathvariant="normal">s</mi> </mrow> </semantics></math>): (<b>a</b>) <math display="inline"><semantics> <msub> <mi>τ</mi> <mn>2</mn> </msub> </semantics></math> and (<b>b</b>) <math display="inline"><semantics> <msub> <mi>τ</mi> <mn>3</mn> </msub> </semantics></math>. <math display="inline"><semantics> <msub> <mi>τ</mi> <mn>2</mn> </msub> </semantics></math>-Precipitation and <math display="inline"><semantics> <msub> <mi>τ</mi> <mn>2</mn> </msub> </semantics></math>-Discharge: <span class="html-italic">p</span>-values &lt; 0.001; <math display="inline"><semantics> <msub> <mi>τ</mi> <mn>3</mn> </msub> </semantics></math>-Precipitation and <math display="inline"><semantics> <msub> <mi>τ</mi> <mn>3</mn> </msub> </semantics></math>-Discharge: <span class="html-italic">p</span>-value for slope &lt; 0.001, intercept not significant with <span class="html-italic">p</span> = 0.451, assuming <math display="inline"><semantics> <mi>α</mi> </semantics></math> = 0.05 and a 31-year rolling-windowed L-moments (RDLMs).</p>
Full article ">Figure 3
<p>Comparison of predicted and observed L-moments (LMs; testing during training) using a 31-year rolling-window. The plots display the predicted (green) and observed (blue) values for the first (<b>a</b>), second (<b>b</b>), third (<b>c</b>), and fourth (<b>d</b>) LMs. Each subplot includes the Overall Mean Squared Error (MSE) between the predicted and observed values computed by summing and averaging the best-fit model Squared Error (SE) for each step in the forward chaining process. The equations plotted alongside the model are derived from the final (best-fit) iteration (index 38), which demonstrated the lowest SE.</p>
Full article ">Figure 4
<p>Location, scale, and shape parameters estimated using [<a href="#B45-hydrology-11-00154" class="html-bibr">45</a>]’s method of L-moments (LMs) for the Generalized Extreme Value (GEV) probability distribution (PD) for the observed (blue) and predicted (testing during training; dashed red) LMs under a 31-year rolling-window.</p>
Full article ">Figure 5
<p>LMRDs using 31-year windows under two Representative Concentration Pathways (RCP) scenarios. Panels (<b>a</b>–<b>c</b>) correspond to the RCP 4.5 scenario, while panels (<b>d</b>–<b>f</b>) correspond to the RCP 8.5 scenario. Diagrams show: (<b>a</b>,<b>d</b>) L-CV/L-skewness, (<b>b</b>,<b>e</b>) L-kurtosis/L-skewness, with theoretical PDs described in [<a href="#B45-hydrology-11-00154" class="html-bibr">45</a>] (distributions include Generalized Logistic (GLO), Generalized Extreme Value (GEV), Generalized Pareto (GPA), Generalized Normal (GNO), Pearson Type III (PE3), Wakeby Lower Bound (WAK_LB), and All Distribution Lower Bound (ALL_LB)). Plots (<b>c</b>,<b>f</b>) show the distribution count for each window. The observed LMRs for the 5-day September low-flow at the Water Survey of Canada (WSC) Goat River Near Erikson Hydrometric Gauge Station (<a href="https://wateroffice.ec.gc.ca/report/data_availability_e.html?type=historical&amp;station=08NH004&amp;parameter_type=Flow&amp;wbdisable=true" target="_blank">08NH004</a>) are plotted alongside simulated future data derived from Multiple Linear Regression (MLR) driven with total August precipitation LMs. Future data are generated using a splice of six Coupled Model Intercomparison Project Phase 5 (CMIP5) series downscaled climate models (median of “ACCESS1-0”, “CanESM2”, “CCSM4”, “CNRM-CM5”, “HadGEM2-ES”, and “MPI-ESM-LR” from 2018 to 2100) downloaded using the single cell extraction tool from the Pacific Climate Impacts Consortium (<a href="https://pacificclimate.org/data/gridded-hydrologic-model-output" target="_blank">PCIC</a>). Historical climate data are downloaded from Historical Climate Data Online (HCDO) repository for the Creston station (Climate ID <a href="https://climate.weather.gc.ca/climate_data/daily_data_e.html?hlyRange=%7C&amp;dlyRange=1912-06-01%7C2017-12-31&amp;mlyRange=1912-01-01%7C2007-02-01&amp;StationID=1111&amp;Prov=BC&amp;urlExtension=_e.html&amp;searchType=stnName&amp;optLimit=yearRange&amp;StartYear=1840&amp;EndYear=2024&amp;selRowPerPage=25&amp;Line=0&amp;searchMethod=contains&amp;Month=12&amp;Day=2&amp;txtStationName=Creston&amp;timeframe=2&amp;Year=2017" target="_blank">1142160</a>; available from 1996 to 2017).</p>
Full article ">Figure 6
<p>Results of the L-moments (LMs) derived from Multiple Linear Regression (MLR) fit to a Generalized Extreme Value (GEV) probability distribution (PD) to produce shape, scale, and location parameters: (<b>a</b>) GEV parameters (shape, scale, and location) over 144 rolling windowed time units under Representative Concentration Pathway (RCP) 4.5 and (<b>b</b>) RCP 8.5.</p>
Full article ">Figure 7
<p>LMs derived from MLR fit to a GEV PD to produce shape, scale, and location parameters to derive median flows (<math display="inline"><semantics> <mrow> <msup> <mi mathvariant="normal">m</mi> <mn>3</mn> </msup> <mo>/</mo> <mi mathvariant="normal">s</mi> </mrow> </semantics></math>) and confidence intervals (CIs) estimated from percentiles: (<b>a</b>) simulated September 5-day low-flows driven by MLR regression using August total precipitation LMs over a 31-year rolling time window under (<b>a</b>) RCP 4.5 and (<b>b</b>) RCP 8.5. The dashed red line denotes a flow of &lt;1 m<sup>3</sup>/s. Note: each simulation is based on <span class="html-italic">n</span> = 1000 iterations for both panels.</p>
Full article ">Figure 8
<p>Results of the LMs derived from MLR fit to the best-fit probability distribution (PD) (distributions described in [<a href="#B45-hydrology-11-00154" class="html-bibr">45</a>]) to produce shape, scale, and location parameters to derive median flows (<math display="inline"><semantics> <mrow> <msup> <mi mathvariant="normal">m</mi> <mn>3</mn> </msup> <mo>/</mo> <mi mathvariant="normal">s</mi> </mrow> </semantics></math>) and confidence intervals (CIs) estimated numerically from percentiles: (<b>a</b>) simulated September 5-day low-flows driven by MLR regression using August total precipitation LMs over a 31-year moving window under RCP 4.5 and (<b>b</b>) RCP 8.5. The dashed red line denotes a flow of &lt;1 m<sup>3</sup>/s. Note: each simulation is based on <span class="html-italic">n</span> = 1000 iterations for both panels.</p>
Full article ">Figure 9
<p>Overall standardized Mean Square Error (MSE) across different window sizes during model training.</p>
Full article ">Figure 10
<p>Sensitivity of window size on location, scale, and shape parameters for September 5-day low-flow estimated using the method of LMs described in [<a href="#B45-hydrology-11-00154" class="html-bibr">45</a>] derived from MLR driven by total August precipitation LMs for six common distributions described in [<a href="#B45-hydrology-11-00154" class="html-bibr">45</a>] (Generalized Extreme Value (GEV; (<b>a</b>–<b>c</b>)), Generalized Logistic (GLO; (<b>d</b>–<b>f</b>)), Generalized Normal (GNO; (<b>g</b>–<b>i</b>)), Pearson Type III (PE3; (<b>j</b>–<b>l</b>)), and Generalized Pareto (GPA; (<b>m</b>–<b>o</b>)). The solid line displays data under the Representative Concentration Pathway (RCP) 4.5 emission scenario, while the dashed line displays the RCP 8.5 emissions scenario.</p>
Full article ">Figure 11
<p>Low-flow exceedance and cumulative exceedance probability for the Goat River near Erikson Gauge Station, showing values less than 2.7 cubic meters per second (<math display="inline"><semantics> <mrow> <msup> <mi mathvariant="normal">m</mi> <mn>3</mn> </msup> <mo>/</mo> <mi mathvariant="normal">s</mi> </mrow> </semantics></math>) assuming <span class="html-italic">n</span> = 1000 simulations and a Generalized Extreme Value (GEV) probability distribution (PD). Future data assume a Representative Concentration Pathway (RCP) 4.5 emissions scenario.</p>
Full article ">
20 pages, 21202 KiB  
Article
Distribution Characteristics and Genesis Mechanism of Ground Fissures in Three Northern Counties of the North China Plain
by Chao Xue, Mingdong Zang, Zhongjian Zhang, Guoxiang Yang, Nengxiong Xu, Feiyong Wang, Cheng Hong, Guoqing Li and Fujiang Wang
Sustainability 2024, 16(18), 8027; https://doi.org/10.3390/su16188027 (registering DOI) - 13 Sep 2024
Viewed by 295
Abstract
The North China Plain is among the regions most afflicted by ground fissure disasters in China. Recent urbanization has accelerated ground fissure activity in the three counties of the northern North China Plain, posing significant threats to both the natural environment and socioeconomic [...] Read more.
The North China Plain is among the regions most afflicted by ground fissure disasters in China. Recent urbanization has accelerated ground fissure activity in the three counties of the northern North China Plain, posing significant threats to both the natural environment and socioeconomic sustainability. Despite the increased attention, a lack of comprehensive understanding persists due to delayed recognition and limited research. This study conducted field visits and geological surveys across 43 villages and 80 sites to elucidate the spatial distribution patterns of ground fissures in the aforementioned counties. By integrating these findings with regional geological data, we formulated a causative model to explain ground fissure formation. Our analysis reveals a concentration of ground fissures near the Niuxi and Rongxi faults, with the former exhibiting the most extensive distribution. The primary manifestations of ground fissures include linear cracks and patch-shaped collapse pits, predominantly oriented in east-west and north-south directions, indicating tensile failure with minimal vertical displacement. Various factors contribute to ground fissure development, including fault activity, ancient river channel distribution, bedrock undulations, rainfall, and ground settlement. Fault activity establishes a concealed fracture system in shallow geotechnical layers, laying the groundwork for ground fissure formation. Additionally, the distribution of ancient river channels and bedrock undulations modifies regional stress fields, further facilitating ground fissure emergence. Rainfall and differential ground settlement serve as triggering mechanisms, exposing ground fissures at the surface. This research offers new insights into the causes of ground fissures in the northern North China Plain, providing crucial scientific evidence for sustaining both the natural environment and the socio-economic stability of the region. Full article
(This article belongs to the Topic Natural Hazards and Disaster Risks Reduction, 2nd Volume)
Show Figures

Figure 1

Figure 1
<p>Distribution of the ground fissures and the faults in the study area. F1—Rongxi fault, F2—Rongdong fault, F3—Rongcheng fault, F4—Niuxi fault, F5—Niudong fault; A1—Xushui depression, A2—Rongcheng uplift, A3—Langgu depression, A4—Niutuo Town uplift, A5—Baxian depression; A-A′—section line; and D<sub>1</sub> and D<sub>2</sub>—drilling wells.</p>
Full article ">Figure 2
<p>Tectonic profile A-A′ in the study area [<a href="#B37-sustainability-16-08027" class="html-bibr">37</a>]. Location of the profile line is shown in <a href="#sustainability-16-08027-f001" class="html-fig">Figure 1</a>. F1—Rongxi fault, F2—Rongdong fault, F3—Rongcheng fault, F4—Niuxi fault, F5—Niudong fault.</p>
Full article ">Figure 3
<p>Groundwater level change in the study area [<a href="#B38-sustainability-16-08027" class="html-bibr">38</a>].</p>
Full article ">Figure 4
<p>Stratigraphic column profile in the study area based on drilling data [<a href="#B34-sustainability-16-08027" class="html-bibr">34</a>,<a href="#B35-sustainability-16-08027" class="html-bibr">35</a>]. The sites of drilling wells D<sub>1</sub> and D<sub>2</sub> are shown in <a href="#sustainability-16-08027-f001" class="html-fig">Figure 1</a>.</p>
Full article ">Figure 5
<p>Distribution of the ground fissures in Zhangweizhuangtou village and surrounding areas: (<b>a</b>–<b>i</b>) typical photos of the ground fissures; f1—Zhangweizhuangtou ground fissure.</p>
Full article ">Figure 6
<p>Distribution of the ground fissures in Beihoutai village, Nanhoutai village, Jiaguang village and surrounding areas: (<b>a</b>–<b>c</b>,<b>e</b>) typical photos of the wall fissures; (<b>d</b>) typical photos of the house ground subsidence; and (<b>f</b>) typical photos of the house floor fissures; F1—Rongxi fault; f2—Beihoutai ground fissure.</p>
Full article ">Figure 7
<p>Distribution of the ground fissures in Longwanxi village and surrounding areas: (<b>a</b>–<b>h</b>) typical photos of the ground fissures; f3–f5: Longwanxi ground fissures.</p>
Full article ">Figure 8
<p>Distribution of the ground fissures in Beizhang village and surrounding areas: (<b>a</b>–<b>g</b>) typical photos of the wall fissures; f6—Beizhang ground fissure.</p>
Full article ">Figure 9
<p>Distribution of the ground fissures in Dongangezhuang village and surrounding areas: (<b>a</b>–<b>e</b>): typical photos of the ground fissures; f7–f9: Dongangezhuang ground fissures.</p>
Full article ">Figure 10
<p>Formation process of rainfall-induced ground fissures under the combined influence of fault activity and rainfall erosion: (<b>a</b>) fault activity initiates the formation of concealed fissures near the surface; (<b>b</b>) infiltration of surface water leads to erosion of the soil layer, migration of soil particles, widening of cracks, and the creation of cavities; (<b>c</b>) fissures propagate upward, causing surface soil to collapse into linear fissures or collapse pits.</p>
Full article ">Figure 11
<p>Contour map of land subsidence rate in the study area (2016) [<a href="#B39-sustainability-16-08027" class="html-bibr">39</a>]. f1–f9: typical ground fissures in the study area and the details are shown in <a href="#sustainability-16-08027-t002" class="html-table">Table 2</a>.</p>
Full article ">Figure 12
<p>Pre-Cenozoic bedrock buried depth contour map and paleochannel distribution map in the study area [<a href="#B40-sustainability-16-08027" class="html-bibr">40</a>,<a href="#B41-sustainability-16-08027" class="html-bibr">41</a>]. f1–f9: typical ground fissures in the study area and the details are shown in <a href="#sustainability-16-08027-t002" class="html-table">Table 2</a>.</p>
Full article ">Figure 13
<p>Formation process of palaeochannel-type ground fissures: (<b>a</b>) original formation state; (<b>b</b>) the initial pumping resulted in uneven settlement of the strata, resulting in a tensile stress concentration area at the shoulder of the palaeochannel and forming hidden cracks; (<b>c</b>) further pumping causes uneven ground settlement to intensify, and hidden cracks develop and then appear on the surface; and (<b>d</b>) stereogram of genetic mechanism of palaeochannel type ground fissures.</p>
Full article ">Figure 14
<p>Formation process of bedrock ridge-type ground fissures: (<b>a</b>) original formation state; (<b>b</b>) the initial pumping results in uneven settlement of the strata, resulting in a tensile stress concentration area at the bedrock ridge and forming hidden cracks; (<b>c</b>) further pumping causes uneven ground settlement to intensify, and hidden cracks develop and then appear on the surface; and (<b>d</b>) stereogram of genetic mechanism of bedrock ridge-type ground fissures.</p>
Full article ">Figure 15
<p>Formation process of bedrock step-type ground fissures: (<b>a</b>) original formation state; (<b>b</b>) the initial pumping results in uneven formation settlement, and the tension stress concentration area is generated in the sudden change of terrain, forming hidden cracks; (<b>c</b>) further pumping causes uneven ground settlement to intensify, and hidden cracks develop and then appear on the surface; and (<b>d</b>) stereogram of genetic mechanism of bedrock step-type ground fissures.</p>
Full article ">
15 pages, 4826 KiB  
Article
Assessing Evapotranspiration Changes in Response to Cropland Expansion in Tropical Climates
by Leonardo Laipelt, Julia Brusso Rossi, Bruno Comini de Andrade, Morris Scherer-Warren and Anderson Ruhoff
Remote Sens. 2024, 16(18), 3404; https://doi.org/10.3390/rs16183404 - 13 Sep 2024
Viewed by 186
Abstract
The expansion of cropland in tropical regions has significantly accelerated in recent decades, triggering an escalation in water demand and changing the total water loss to the atmosphere (evapotranspiration). Additionally, the increase in areas dedicated to agriculture in tropical climates coincides with an [...] Read more.
The expansion of cropland in tropical regions has significantly accelerated in recent decades, triggering an escalation in water demand and changing the total water loss to the atmosphere (evapotranspiration). Additionally, the increase in areas dedicated to agriculture in tropical climates coincides with an increased frequency of drought events, leading to a series of conflicts among water users. However, detailed studies on the impacts of changes in water use due to agriculture expansion, including irrigation, are still lacking. Furthermore, the higher presence of clouds in tropical environments poses challenges for the availability of high-resolution data for vegetation monitoring via satellite images. This study aims to analyze 37 years of agricultural expansion using the Landsat collection and a satellite-based model (geeSEBAL) to assess changes in evapotranspiration resulting from cropland expansion in tropical climates, focusing on the São Marcos River Basin in Brazil. It also used a methodology for estimating daily evapotranspiration on days without satellite images. The results showed a 34% increase in evapotranspiration from rainfed areas, mainly driven by soybean cultivation. In addition, irrigated areas increased their water use, despite not significantly changing water use at the basin scale. Conversely, natural vegetation areas decreased their evapotranspiration rates by 22%, suggesting possible further implications with advancing changes in land use and land cover. Thus, this study underscores the importance of using satellite-based evapotranspiration estimates to enhance our understanding of water use across different land use types and scales, thereby improving water management strategies on a large scale. Full article
Show Figures

Figure 1

Figure 1
<p>São Marcos River Basin: location in Brazil (<b>a</b>), climate zones according to Köppen–Geiger classification and irrigation pivots (<b>b</b>), and Normalized Difference Vegetation Index (NDVI) values computed using average composition of Landsat 8 for 2021 (<b>c</b>).</p>
Full article ">Figure 2
<p>Daily average of <span class="html-italic">ET</span> illustrated as boxplot for each land cover and use (<b>a</b>) for Landsat scenes between 1986 and 2022. We also illustrated the seasonal monthly average of <span class="html-italic">ET</span> (<b>b</b>), and trends of annual average <span class="html-italic">ET</span> for different land types, with natural vegetation (forest and savanna) demonstrating positive trends over the years, as well as irrigated areas (<b>c</b>).</p>
Full article ">Figure 3
<p>Changes in the <span class="html-italic">ET</span> spatial patterns for the São Marcos River Basin from 1986 to 2022 (<b>a</b>). The contribution of the water usage for each land cover and use between 1986 and 2022 is shown in (<b>b</b>), whereas (<b>c</b>) illustrates changes in land cover and use.</p>
Full article ">Figure 4
<p>Annual composition ET (mm day<sup>−1</sup>) in the São Marcus River Basin between 1986 (<b>a</b>) and 2021 (<b>b</b>). Highlighted plots showed the expressive number of pivot irrigation systems over the basin for specific locations.</p>
Full article ">Figure 5
<p>Monthly <span class="html-italic">ET</span> in the São Marcus River basin was analyzed for each month of one water year (2019 and 2020). During the dry season (May to September), precipitation is limited and radiation availability is high, being a water-limited environment. Consequently, lower <span class="html-italic">ET</span> values are observed during the dry season, while the wet season increases <span class="html-italic">ET</span> rates due to higher precipitation availability.</p>
Full article ">Figure 6
<p>Seasonal differences in daily <span class="html-italic">ET</span> for irrigated and rainfed croplands in the São Marcus River Basin (<b>a</b>), and the difference between both estimations (<b>b</b>). We used a simplified method to fill the gap between Landsat scenes by interpolating <span class="html-italic">EF</span> over time and multiplying with the respective reference <span class="html-italic">ET</span>.</p>
Full article ">
15 pages, 11451 KiB  
Article
Impact of Climate Change on Distribution of Suitable Niches for Black Locust (Robinia pseudoacacia L.) Plantation in China
by Shanchao Zhao, Hesong Wang and Yang Liu
Forests 2024, 15(9), 1616; https://doi.org/10.3390/f15091616 - 13 Sep 2024
Viewed by 178
Abstract
Black locust (Robinia pseudoacacia L.), one of the major afforestation species adopted in vegetation restoration, is notable for its rapid root growth and drought resistance. It plays a vital role in improving the natural environment and soil fertility, contributing significantly to soil [...] Read more.
Black locust (Robinia pseudoacacia L.), one of the major afforestation species adopted in vegetation restoration, is notable for its rapid root growth and drought resistance. It plays a vital role in improving the natural environment and soil fertility, contributing significantly to soil and water conservation and biodiversity protection. However, compared with natural forests, due to the low diversity, simple structure and poor stability, planted forests including Robinia pseudoacacia L. are more sensitive to the changing climate, especially in the aspects of growth trend and adaptive range. Studying the ecological characteristics and geographical boundaries of Robinia pseudoacacia L. is therefore important to explore the adaptation of suitable niches to climate change. Here, based on 162 effective distribution records in China and 22 environmental variables, the potential distribution of suitable niches for Robinia pseudoacacia L. plantations in past, present and future climates was simulated by using a Maximum Entropy (MaxEnt) model. The results showed that the accuracy of the MaxEnt model was excellent and the area under the curve (AUC) value reached 0.937. Key environmental factors constraining the distribution and suitable intervals were identified, and the geographical distribution and area changes of Robinia pseudoacacia L. plantations in future climate scenarios were also predicted. The results showed that the current suitable niches for Robinia pseudoacacia L. plantations covered 9.2 × 105 km2, mainly distributed in the Loess Plateau, Huai River Basin, Sichuan Basin, eastern part of the Yunnan–Guizhou Plateau, Shandong Peninsula, and Liaodong Peninsula. The main environmental variables constraining the distribution included the mean temperature of the driest quarter, precipitation of driest the quarter, temperature seasonality and altitude. Among them, the temperature of the driest quarter was the most important factor. Over the past 90 years, the suitable niches in the Sichuan Basin and Yunnan–Guizhou Plateau have not changed significantly, while the suitable niches north of the Qinling Mountains have expanded northward by 2° and the eastern area of Liaoning Province has expanded northward by 1.2°. In future climate scenarios, the potential suitable niches for Robinia pseudoacacia L. are expected to expand significantly in both the periods 2041–2060 and 2061–2080, with a notable increase in highly suitable niches, widely distributed in southern China. A warning was issued for the native vegetation in the above-mentioned areas. This work will be beneficial for developing reasonable afforestation strategies and understanding the adaptability of planted forests to climate change. Full article
(This article belongs to the Section Forest Meteorology and Climate Change)
Show Figures

Figure 1

Figure 1
<p>The distribution records of <span class="html-italic">Robinia pseudoacacia</span> L. and the approximate range of the Loess Plateau and the Yunnan–Guizhou Plateau.</p>
Full article ">Figure 2
<p>Receiver operating characteristic (ROC) curve of the MaxEnt model used in this study.</p>
Full article ">Figure 3
<p>Response curves of <span class="html-italic">Robinia pseudoacacia</span> L. plantations to the main environmental factors.</p>
Full article ">Figure 4
<p>Potential distribution areas of <span class="html-italic">Robinia pseudoacacia</span> L. plantations in the current climate of China.</p>
Full article ">Figure 5
<p>Potential distribution areas of <span class="html-italic">Robinia pseudoacacia</span> L. plantations for the period from 1931 to 1960 in China.</p>
Full article ">Figure 6
<p>Potential distribution areas of <span class="html-italic">Robinia pseudoacacia</span> L. plantations in the period from 1961 to 1990 in China.</p>
Full article ">Figure 7
<p>Potential distribution areas of <span class="html-italic">Robinia pseudoacacia</span> L. plantations in future climate change scenarios (2041–2060).</p>
Full article ">Figure 7 Cont.
<p>Potential distribution areas of <span class="html-italic">Robinia pseudoacacia</span> L. plantations in future climate change scenarios (2041–2060).</p>
Full article ">Figure 8
<p>Potential distribution areas of <span class="html-italic">Robinia pseudoacacia</span> L. plantations in future climate change scenarios (2061–2080).</p>
Full article ">
21 pages, 13186 KiB  
Article
Variations in Microstructure and Collapsibility Mechanisms of Malan Loess across the Henan Area of the Middle and Lower Reaches of the Yellow River
by Yi Wei and Zhiquan Huang
Appl. Sci. 2024, 14(18), 8220; https://doi.org/10.3390/app14188220 - 12 Sep 2024
Viewed by 219
Abstract
The Henan area of the middle and lower reaches of the Yellow River is situated within the third sedimentary loess area, positioned as the southeasternmost segment within the transitional belt connecting the Loess Plateau with the North China Plain. Addressing concerns related to [...] Read more.
The Henan area of the middle and lower reaches of the Yellow River is situated within the third sedimentary loess area, positioned as the southeasternmost segment within the transitional belt connecting the Loess Plateau with the North China Plain. Addressing concerns related to loess collapse, landslides, and subgrade settlement across various regions attributable to the collapsible nature of Malan loess in western Henan, this study undertook collapsibility testing of undisturbed Malan loess in the province. The different mechanisms of loess collapsibility in different regions were explained from the microstructure by using the indoor immersion-compression test double-line method, scanning electron microscope (SEM), and particles and cracks analysis system (PCAS). The relationship between quantitative factors of microstructure and collapsibility of loess was analyzed by linear regression analysis. The findings indicate that under identical overburden pressure and immersion conditions, the collapsibility of Malan loess in western Henan diminishes progressively from west to east. Microstructural tests were conducted on various loess specimens using scanning electron microscopy, revealing that the distribution of loess particles is notably concentrated in the Xingyang and Gongyi areas, leading to a reduction in pore area compared to the Shanzhou and Mianchi areas. While the Mianchi and Shanzhou areas exhibit a loose arrangement of loess particles, those in Xingyang and Gongyi are comparatively denser. Analysis of microstructural images through the particles and cracks analysis system elucidated that the pore arrangement in the Gongyi and Xingyang areas is more stable than in the Mianchi and Shanzhou areas. Additionally, there is a gradual concentration in particle distribution, accompanied by an increase in agglomeration degree. According to the analysis and comparison of microstructure and quantitative parameters of four groups of loess samples before and after collapsibility, it is revealed that the change mechanism underlying loess collapsibility in various regions of western Henan primarily stems from the external factors influencing the microstructural alterations within the loess. The microstructural determinants contributing to collapsibility changes in different regions encompass three principal aspects: Firstly, modifications in the grain morphology of the Malan loess skeleton in western Henan are notable. Secondly, variations in the internal pore characteristics of loess microstructure are observed. Thirdly, disparities exist in the interconnections between soil particles. The findings of this research hold significant worth for improving construction safety and geological hazard prevention within the Loess region of western Henan. Full article
Show Figures

Figure 1

Figure 1
<p>Location of the sampling sites.</p>
Full article ">Figure 2
<p>Double-line method.</p>
Full article ">Figure 3
<p>Preservation of soil samples.</p>
Full article ">Figure 4
<p>The curve of collapsibility coefficient versus pressure in different areas.</p>
Full article ">Figure 5
<p>The curve of the collapsibility coefficient versus pressure under different initial water content ((<b>a</b>) Shanzhou, (<b>b</b>) Mianchi, (<b>c</b>) Gongyi, (<b>d</b>) Xingyang).</p>
Full article ">Figure 6
<p>The curve of collapsibility coefficient versus pressure in different areas (<span class="html-italic">ω</span> = 6.79%).</p>
Full article ">Figure 7
<p>Sem images illustrating the overall structure of loess in the western Henan area. ((<b>a</b>) Shanzhou Microstructure before soaking, (<b>b</b>) Shanzhou Microstructure after soaking, (<b>c</b>) Mianchi Microstructure before soaking, (<b>d</b>) Mianchi Microstructure after soaking, (<b>e</b>) Gongyi Microstructure before soaking, (<b>f</b>) Gongyi Microstructure after soaking, (<b>g</b>) Xingyang Microstructure before soaking, (<b>h</b>) Xingyang Microstructure after soaking).</p>
Full article ">Figure 7 Cont.
<p>Sem images illustrating the overall structure of loess in the western Henan area. ((<b>a</b>) Shanzhou Microstructure before soaking, (<b>b</b>) Shanzhou Microstructure after soaking, (<b>c</b>) Mianchi Microstructure before soaking, (<b>d</b>) Mianchi Microstructure after soaking, (<b>e</b>) Gongyi Microstructure before soaking, (<b>f</b>) Gongyi Microstructure after soaking, (<b>g</b>) Xingyang Microstructure before soaking, (<b>h</b>) Xingyang Microstructure after soaking).</p>
Full article ">Figure 8
<p>Sem images depicting the overall structure of loess in the western Henan area (<span class="html-italic">ω</span> = 6.79%). ((<b>a</b>) Shanzhou microstructure before soaking, (<b>b</b>) Shanzhou microstructure after soaking, (<b>c</b>) Mianchi microstructure before soaking, (<b>d</b>) Mianchi microstructure after soaking, (<b>e</b>) Gongyi microstructure before soaking, (<b>f</b>) Gongyi microstructure after soaking, (<b>g</b>) Xingyang microstructure before soaking, (<b>h</b>) Xingyang microstructure after soaking).</p>
Full article ">Figure 8 Cont.
<p>Sem images depicting the overall structure of loess in the western Henan area (<span class="html-italic">ω</span> = 6.79%). ((<b>a</b>) Shanzhou microstructure before soaking, (<b>b</b>) Shanzhou microstructure after soaking, (<b>c</b>) Mianchi microstructure before soaking, (<b>d</b>) Mianchi microstructure after soaking, (<b>e</b>) Gongyi microstructure before soaking, (<b>f</b>) Gongyi microstructure after soaking, (<b>g</b>) Xingyang microstructure before soaking, (<b>h</b>) Xingyang microstructure after soaking).</p>
Full article ">Figure 9
<p>Comparison between the original image and the binary image.</p>
Full article ">Figure 10
<p>Variation curves of quantitative microstructure parameters in different regions. ((<b>a</b>) average form factor, (<b>b</b>) probability entropy, (<b>c</b>) porosity distribution dimension, (<b>d</b>) pore area ratio).</p>
Full article ">Figure 11
<p>Relationship between the quantitative factors of microstructure and collapsibility coefficient. ((<b>a</b>) average form factor, (<b>b</b>) probability entropy, (<b>c</b>) porosity distribution dimension, (<b>d</b>) pore area ratio).</p>
Full article ">Figure 11 Cont.
<p>Relationship between the quantitative factors of microstructure and collapsibility coefficient. ((<b>a</b>) average form factor, (<b>b</b>) probability entropy, (<b>c</b>) porosity distribution dimension, (<b>d</b>) pore area ratio).</p>
Full article ">Figure 12
<p>Variation curves of quantitative microstructural parameters in different regions. ((<b>a</b>) Average form factor, (<b>b</b>) Probability entropy, (<b>c</b>) Porosity distribution dimension, (<b>d</b>) Pore area ratio).</p>
Full article ">Figure 13
<p>Relationship between the quantitative factors of microstructure and collapsibility coefficient. ((<b>a</b>) average form factor, (<b>b</b>) probability entropy, (<b>c</b>) porosity distribution dimension, (<b>d</b>) pore area ratio).</p>
Full article ">Figure 13 Cont.
<p>Relationship between the quantitative factors of microstructure and collapsibility coefficient. ((<b>a</b>) average form factor, (<b>b</b>) probability entropy, (<b>c</b>) porosity distribution dimension, (<b>d</b>) pore area ratio).</p>
Full article ">
16 pages, 3006 KiB  
Article
Biomonitoring of Waters and Tambacu (Colossoma macropomum × Piaractus mesopotamicus) from the Amazônia Legal, Brazil
by Karuane Saturnino da Silva Araújo, Thiago Machado da Silva Acioly, Ivaneide Oliveira Nascimento, Francisca Neide Costa, Fabiano Corrêa, Ana Maria Gagneten and Diego Carvalho Viana
Water 2024, 16(18), 2588; https://doi.org/10.3390/w16182588 - 12 Sep 2024
Viewed by 355
Abstract
Fish farming is increasingly important globally and nationally, playing a crucial role in fish production for human consumption. Monitoring microbiological and chemical contaminants from water discharge is essential to mitigate the risk of contaminating water and fish for human consumption. This study analyzes [...] Read more.
Fish farming is increasingly important globally and nationally, playing a crucial role in fish production for human consumption. Monitoring microbiological and chemical contaminants from water discharge is essential to mitigate the risk of contaminating water and fish for human consumption. This study analyzes the physicochemical and E. coli parameters of water and tambacu fish muscles (Colossoma macropomum × Piaractus mesopotamicus) in Western Maranhão, Brazil. It also includes a qualitative characterization of zooplankton in the ponds. Samples were collected from tambacu ponds in a dam system fed by natural watercourses from the Tocantins River tributaries, located at the connection of the Brazilian savanna and Amazon biomes. The physicochemical and E. coli parameters of water did not meet national standards. The zooplankton community included Rotifera, Cladocera, Copepoda, and Protozoa representatives, with no prior studies on zooplankton in the region, making these findings unprecedented. The biological quality of freshwater is crucial in fish farming, as poor quality can lead to decreased productivity and fish mortality, raising significant food safety concerns. The water quality studied is related to the potential influence of untreated wastewater as a source of contamination, leaving the studied region still far from safe water reuse practices. The findings on chemical and E. coli contamination of fish farming waters concern human health and emphasize the need for appropriate regulations. Full article
Show Figures

Figure 1

Figure 1
<p>Location and key characteristics of the study area, Maranhão, Brazil.</p>
Full article ">Figure 2
<p>Specimen of tambacu (<span class="html-italic">Colossoma macropomum × Piaractus mesopotamicus</span>) from the Amazônia legal, Brazil.</p>
Full article ">Figure 3
<p>Zooplankton species sampled in the water of tambacu fish farming tanks.</p>
Full article ">
20 pages, 10622 KiB  
Article
Machine Learning Model for River Discharge Forecast: A Case Study of the Ottawa River in Canada
by M. Almetwally Ahmed and S. Samuel Li
Hydrology 2024, 11(9), 151; https://doi.org/10.3390/hydrology11090151 - 12 Sep 2024
Viewed by 202
Abstract
River discharge is an essential input to hydrosystem projects. This paper aimed to modify the group method of data handling (GMDH) to create a new artificial intelligent forecast model (abbreviated as MGMDH) for predicting discharges at river cross-sections (CSs). The basic idea was [...] Read more.
River discharge is an essential input to hydrosystem projects. This paper aimed to modify the group method of data handling (GMDH) to create a new artificial intelligent forecast model (abbreviated as MGMDH) for predicting discharges at river cross-sections (CSs). The basic idea was to optimise the weights for selected hydrometric and meteorological predictors. One novelty of this study was that MGMDH could take the discharge observed from a neighbouring CS as a predictor when observations from the CS of interest had ceased. Another novelty was that MGMDH could include meteorological parameters as extra predictors. The model was validated using data from natural rivers. For given lead times, MGMDH automatically determined the best forecast equations, consistent with physical river hydraulics laws. This automation minimised computing time while improving accuracy. The model gave reliable forecasts, with a coefficient of determination greater than 0.978. For lead times close to the advection time from upstream to the CS of interest, the forecast had the highest reliability. MGMDH results compared well with some other machine learning models, like neural networks and the adaptive structure of the group method of data handling. It has potential applications for efficiently forecasting discharge and offers a tool to support flood management. Full article
(This article belongs to the Section Water Resources and Risk Management)
Show Figures

Figure 1

Figure 1
<p>(<b>a</b>) Close-up view of the Ottawa River between hydrometric stations 02KF009 (CS upstream or CSU) and 02KF005 (CS downstream or CSD); (<b>b</b>) broad view of the stream network, watershed boundaries, and outlet points, illustrating how these watersheds flow into and connect with the Ottawa River.</p>
Full article ">Figure 2
<p>Schematic time series of continuous observations of <span class="html-italic">q</span> (solid black curve), discontinued observations of <span class="html-italic">q</span> (solid red curve), and predicted future values of <math display="inline"><semantics> <mrow> <mover accent="true"> <mrow> <mi>q</mi> </mrow> <mo stretchy="false">^</mo> </mover> </mrow> </semantics></math> (dashed curves).</p>
Full article ">Figure 3
<p>Definition diagram of river flow: (<b>a</b>) top view of the river channel; (<b>b</b>) CS at upstream (CSU) with discharge <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>q</mi> </mrow> <mrow> <mi>U</mi> </mrow> </msub> </mrow> </semantics></math> and water level <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>η</mi> </mrow> <mrow> <mi>U</mi> </mrow> </msub> </mrow> </semantics></math> (above a certain reference datum); (<b>c</b>) CS at downstream (CSD) with discharge <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>q</mi> </mrow> <mrow> <mi>D</mi> </mrow> </msub> </mrow> </semantics></math> and water level <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>η</mi> </mrow> <mrow> <mi>D</mi> </mrow> </msub> </mrow> </semantics></math>.</p>
Full article ">Figure 4
<p>Flowchart of the methods for river discharge forecast.</p>
Full article ">Figure 5
<p>Time series of hourly averaged variable: (<b>a</b>) discharge <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>q</mi> </mrow> <mrow> <mi>D</mi> </mrow> </msub> </mrow> </semantics></math>; (<b>b</b>) water level <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>η</mi> </mrow> <mrow> <mi>D</mi> </mrow> </msub> </mrow> </semantics></math>, observed from the CS of interest (02KF005); (<b>c</b>) discharge <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>q</mi> </mrow> <mrow> <mi>U</mi> </mrow> </msub> </mrow> </semantics></math>; (<b>d</b>) water level <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>η</mi> </mrow> <mrow> <mi>U</mi> </mrow> </msub> </mrow> </semantics></math>, observed from 02KF009, covering a period of 180 days (1 January–30 June 2023). The dotted lines divide the time series into two parts: one for model training, and the other for model testing (the same in subsequent figures).</p>
Full article ">Figure 6
<p>Time series of observed hourly averaged variable: (<b>a</b>) <math display="inline"><semantics> <mrow> <mi>T</mi> </mrow> </semantics></math>; (<b>b</b>) <math display="inline"><semantics> <mrow> <mi>θ</mi> </mrow> </semantics></math>; (<b>c</b>) <math display="inline"><semantics> <mrow> <mi>ϕ</mi> </mrow> </semantics></math>; (<b>d</b>) <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>P</mi> </mrow> <mrow> <mi>a</mi> <mi>t</mi> <mi>m</mi> </mrow> </msub> </mrow> </semantics></math>; (<b>e</b>) <math display="inline"><semantics> <mrow> <mi>P</mi> </mrow> </semantics></math> for the period of 1 January–30 June 2023 at WMO station (ID: 71063) located at 45°23′00″ N, 75°43′00″ W.</p>
Full article ">Figure 7
<p>Values of <math display="inline"><semantics> <mrow> <mover accent="true"> <mrow> <mi>q</mi> </mrow> <mo stretchy="false">^</mo> </mover> </mrow> </semantics></math> predicted from Equation (14) for lead times: (<b>a</b>) <math display="inline"><semantics> <mrow> <mi>δ</mi> <mi>t</mi> </mrow> </semantics></math> = 2, (<b>b</b>) 4, (<b>c</b>) 6, (<b>d</b>) 8, (<b>e</b>) 10, (<b>f</b>) 12, (<b>g</b>) 16, and (<b>h</b>) 18 h, in comparison with observed <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>q</mi> </mrow> <mrow> <mi>D</mi> </mrow> </msub> </mrow> </semantics></math> (test data points in <a href="#hydrology-11-00151-f005" class="html-fig">Figure 5</a>a).</p>
Full article ">Figure 8
<p>Time series of hourly-averaged discharges, observed at CSD (black curve), and forecasted using Equation (14) for the training period (blue curve) and for the testing period (red curve). The forecast is for lead times: (<b>a</b>) <math display="inline"><semantics> <mrow> <mi>δ</mi> <mi>t</mi> </mrow> </semantics></math> = 2; (<b>b</b>) 4; (<b>c</b>) 6; (<b>d</b>) 8; (<b>e</b>) 10; (<b>f</b>) 12; (<b>g</b>) 16; and (<b>h</b>) 18 h.</p>
Full article ">Figure 9
<p>Performance of the forecast model: (<b>a</b>) AIC <span class="html-italic">c</span>; (<b>b</b>) normalised RMSE <math display="inline"><semantics> <mrow> <mover accent="true"> <mrow> <mi>ε</mi> </mrow> <mo stretchy="false">~</mo> </mover> </mrow> </semantics></math>; (<b>c</b>) the coefficient of determination <math display="inline"><semantics> <mrow> <msup> <mrow> <mi>R</mi> </mrow> <mrow> <mn>2</mn> </mrow> </msup> </mrow> </semantics></math>; (<b>d</b>) mean absolute relative error <math display="inline"><semantics> <mrow> <mfenced open="|" close="|" separators="|"> <mrow> <msub> <mrow> <mi>ε</mi> </mrow> <mrow> <mi>a</mi> </mrow> </msub> </mrow> </mfenced> </mrow> </semantics></math>.</p>
Full article ">Figure 10
<p>Reliability (Equation (7)) of the best model functions for: (<b>a</b>) the case of discontinued discharge observation; and (<b>b</b>) the case of continuous discharge observation.</p>
Full article ">Figure 11
<p>Time series of 15-min-averaged discharges, observed at CSD of the Boise River (black curve), and forecasted for the training period (blue curve) and for the testing period (red curve). The forecast is for lead times: (<b>a</b>) <math display="inline"><semantics> <mrow> <mi>δ</mi> <mi>t</mi> </mrow> </semantics></math> = 2; (<b>b</b>) 4; (<b>c</b>) 6; (<b>d</b>) 8; (<b>e</b>) 10; (<b>f</b>) 12; (<b>g</b>) 18; and (<b>h</b>) 24 h.</p>
Full article ">Figure 12
<p>Comparison of performance between MGMDH and other MLMs. The lead time is: (<b>a</b>) <math display="inline"><semantics> <mrow> <mi>δ</mi> <mi>t</mi> </mrow> </semantics></math> = 12 h; and (<b>b</b>) <math display="inline"><semantics> <mrow> <mi>δ</mi> <mi>t</mi> </mrow> </semantics></math> = 2 h.</p>
Full article ">Figure 13
<p>Time series of hourly averaged discharges, observed at CSD of the Missouri River (black curve), and forecasted for the training period (blue curve) and for the testing period (red curve). The forecast is for lead times: (<b>a</b>) <math display="inline"><semantics> <mrow> <mi>δ</mi> <mi>t</mi> </mrow> </semantics></math> = 2; (<b>b</b>) 4; (<b>c</b>) 6; (<b>d</b>) 8; (<b>e</b>) 10; (<b>f</b>) 12; (<b>g</b>) 18; and (<b>h</b>) 24 h.</p>
Full article ">
20 pages, 27493 KiB  
Article
Development and Application of an Environmental Vulnerability Index (EVI) for Identifying Priority Restoration Areas in the São Francisco River Basin, Brazil
by Clívia Dias Coelho, Demetrius David da Silva, Ricardo Santos Silva Amorim, Bruno Nery Fernandes Vasconcelos, Ernani Lopes Possato, Elpídio Inácio Fernandes Filho, Pedro Christo Brandão, José Ambrósio Ferreira Neto and Lucas Vieira Silva
Land 2024, 13(9), 1475; https://doi.org/10.3390/land13091475 - 12 Sep 2024
Viewed by 210
Abstract
The environmental vulnerability diagnosis of a river basin depends on a holistic analysis of its environmental aspects and degradation factors. Based on this diagnosis, the definition of priority areas where interventions for environmental recovery should be carried out is fundamental, since financial and [...] Read more.
The environmental vulnerability diagnosis of a river basin depends on a holistic analysis of its environmental aspects and degradation factors. Based on this diagnosis, the definition of priority areas where interventions for environmental recovery should be carried out is fundamental, since financial and natural resources are limited. In this study, we developed a methodology to assess these fragilities using an environmental vulnerability index (EVI) that combines physical and environmental indicators related to the natural sensitivity of ecosystems and their exposure to anthropogenic factors. The developed EVI was applied to the headwater region of the São Francisco River Basin (SFRB), Brazil. The proposed index was based on the AHP multicriteria analysis and was adapted to include four variables representative of the study area: Land Use Adequacy, Burned Area, Erosion Susceptibility, and quantitative water balance. The EVI analysis highlighted that the presence of easily erodible soils, associated with sloping areas and land use above their capacity, generate the most vulnerable areas in the headwaters of the SFRB. The highest EVI values are primarily linked to regions with shallow, easily erodible soils like Leptosols and Cambisols, found in steep areas predominantly used for pasture. In the SFBR, the greatest vulnerability was observed within a 5 km buffer around conservation units, covering approximately 32.4% of the total area. The results of this study indicate where resources should be applied for environmental preservation in the basin under study, directing the allocation of efforts to areas with lower resilience to maintain ecosystem services. Full article
Show Figures

Figure 1

Figure 1
<p>Study area: selected sub-basins in the headwaters of the São Francisco River Basin.</p>
Full article ">Figure 2
<p>Steps to create the environmental vulnerability index (EVI).</p>
Full article ">Figure 3
<p>(<b>a</b>) Land use capability and (<b>b</b>) land use intensity in the headwaters of São Francisco.</p>
Full article ">Figure 4
<p>Number of exceeding classes in the headwaters of the São Francisco River Basin.</p>
Full article ">Figure 5
<p>Recurrence of fires in the headwaters of the São Francisco River Basin.</p>
Full article ">Figure 6
<p>(<b>a</b>) Rainfall erosivity, (<b>b</b>) soil erodibility, and (<b>c</b>) slope of the headwaters of the São Francisco River Basin.</p>
Full article ">Figure 7
<p>Erosion susceptibility in the headwaters of the São Francisco River Basin.</p>
Full article ">Figure 8
<p>Quantitative water balance in the headwaters of the São Francisco River Basin.</p>
Full article ">Figure 9
<p>Environmental Vulnerability Index (EVI) in the headwaters of the São Francisco River Basin.</p>
Full article ">
20 pages, 5895 KiB  
Article
Comprehensive Zoning Strategies for Flood Disasters in China
by Huipan Li, Yuan Wang, Liying Ping, Na Li and Peng Zhao
Water 2024, 16(17), 2546; https://doi.org/10.3390/w16172546 - 9 Sep 2024
Viewed by 331
Abstract
The frequency of global floods has increased, posing significant threats to economic development and human safety. Existing flood risk zoning studies in disaster prevention lack integration of the natural–economic–social chain and urban resilience factors. This study addresses this gap by constructing flood disaster [...] Read more.
The frequency of global floods has increased, posing significant threats to economic development and human safety. Existing flood risk zoning studies in disaster prevention lack integration of the natural–economic–social chain and urban resilience factors. This study addresses this gap by constructing flood disaster risk and intensity indices using data from 31 provinces and 295 prefectural-level cities in China from 2011 to 2022. These indices incorporate natural (rainfall), economic (GDP), and social (population, built-up area) indicators to assess the flood likelihood and loss degree, providing comprehensive risk and intensity ratings. The study also examines the impact of resilience factors—environmental (green space), infrastructural (rainwater pipeline density), and natural resource (watershed areas)—on flood intensity. Findings reveal that high-risk regions are mainly in the Yangtze River Basin and southern regions, while high-intensity regions are primarily in the middle and lower Yangtze River and certain northwestern cities. Increasing rainwater pipeline density mitigates flood impacts in high-risk, high-intensity areas, while expanding green spaces and pipelines are effective in high-risk, low-intensity regions. This paper proposes a comprehensive flood hazard zoning mechanism integrating natural, economic, and social factors with urban resilience, offering insights and a scientific basis for urban flood management. Full article
(This article belongs to the Section Urban Water Management)
Show Figures

Figure 1

Figure 1
<p>Research framework diagram of the comprehensive zoning strategy for flood disasters in China.</p>
Full article ">Figure 2
<p>Risk and intensity index change plots.</p>
Full article ">Figure 3
<p>Heat map of flood disaster risk index and flood disaster intensity index of all provinces in 2011–2022.</p>
Full article ">Figure 4
<p>Flood disaster risk-intensity grade diagram. In which, (<b>a</b>) is a map of the national flood risk scale and (<b>b</b>) is a map of the national flood intensity scale. BJ—Beijing, TJ—Tianjin, HE—Hebei, SX—Shanxi, IM—Inner Mongolia; HL—Heilongjiang, JL—Jilin, LN—Liaoning; SH—Shanghai, JS—Jiangsu, ZJ—Zhejiang, AH—Anhui, FJ—Fujian, JX—Jiangxi, SD—Shandong; HA—Henan, HB—Hubei, HN—Hunan; GD—Guangdong, GX—Guangxi, HI—Hainan; CQ—Chongqing, SC—Sichuan, GZ—Guizhou, YN—Yunnan; SN—Shaanxi, GS—Gansu, QH—Qinghai, NX—Ningxia, XJ—Xinjiang; XZ—Xizang, HK—Hong Kang, MO—Macao, TW—Taiwan.</p>
Full article ">Figure 5
<p>Flood disaster risk-intensity comprehensive zoning map. In which, (<b>a</b>) is a flood exposure risk-intensity quadrant map, and (<b>b</b>) is a flood exposure risk-intensity quadrant map.</p>
Full article ">Figure 6
<p>Change plot of the risk and intensity index of each flood zone.</p>
Full article ">Figure 7
<p>Stage changes of disaster zoning in each province.</p>
Full article ">Figure 8
<p>Focus on the regional flood disaster level diagram. In which, (<b>a</b>) is the flooding intensity level map, (<b>b</b>) is the built-up area rainwater pipe percentage level map, (<b>c</b>) is the built-up area green space percentage level map, and (<b>d</b>) is the share of water bodies level map. The part of the purple line is the high-risk, low-intensity area, and the area beyond the purple line is the high-risk, high-intensity area.</p>
Full article ">Figure 9
<p>Focus on the proportion of green space in regional built-up areas and the density of rainwater pipe networks in regional built-up areas.</p>
Full article ">Figure 10
<p>Changes in urban resilience index and flood disaster evaluation index.</p>
Full article ">
43 pages, 24204 KiB  
Article
Support Vector Machine Algorithm for Mapping Land Cover Dynamics in Senegal, West Africa, Using Earth Observation Data
by Polina Lemenkova
Earth 2024, 5(3), 420-462; https://doi.org/10.3390/earth5030024 - 6 Sep 2024
Viewed by 494
Abstract
This paper addresses the problem of mapping land cover types in Senegal and recognition of vegetation systems in the Saloum River Delta on the satellite images. Multi-seasonal landscape dynamics were analyzed using Landsat 8-9 OLI/TIRS images from 2015 to 2023. Two image classification [...] Read more.
This paper addresses the problem of mapping land cover types in Senegal and recognition of vegetation systems in the Saloum River Delta on the satellite images. Multi-seasonal landscape dynamics were analyzed using Landsat 8-9 OLI/TIRS images from 2015 to 2023. Two image classification methods were compared, and their performance was evaluated in the GRASS GIS software (version 8.4.0, creator: GRASS Development Team, original location: Champaign, Illinois, USA, currently multinational project) by means of unsupervised classification using the k-means clustering algorithm and supervised classification using the Support Vector Machine (SVM) algorithm. The land cover types were identified using machine learning (ML)-based analysis of the spectral reflectance of the multispectral images. The results based on the processed multispectral images indicated a decrease in savannas, an increase in croplands and agricultural lands, a decline in forests, and changes to coastal wetlands, including mangroves with high biodiversity. The practical aim is to describe a novel method of creating land cover maps using RS data for each class and to improve accuracy. We accomplish this by calculating the areas occupied by 10 land cover classes within the target area for six consecutive years. Our results indicate that, in comparing the performance of the algorithms, the SVM classification approach increased the accuracy, with 98% of pixels being stable, which shows qualitative improvements in image classification. This paper contributes to the natural resource management and environmental monitoring of Senegal, West Africa, through advanced cartographic methods applied to remote sensing of Earth observation data. Full article
Show Figures

Figure 1

Figure 1
<p>Study area with segments of the Landsat images shown on a topographic map of Senegal. Software: GMT. Map source: Author.</p>
Full article ">Figure 2
<p>Data capture of Landsat images from the USGS EarthExplorer repository.</p>
Full article ">Figure 3
<p>Landsat images in RGB colors covering the Cape Verde Peninsula region and Saloum River Delta, West Senegal, in February: (<b>a</b>) 2015; (<b>b</b>) 2018; (<b>c</b>) 2020; (<b>d</b>) 2021; (<b>e</b>) 2022; (<b>f</b>) 2023.</p>
Full article ">Figure 4
<p>Workflow scheme illustrating the data and the main methodological steps. Software: R version 4.3.3, library DiagrammeR version 1.0.11. Diagram source: Author.</p>
Full article ">Figure 5
<p>False color composites of the Landsat 8-9 OLI/TIRS images with vegetation colored red, using a combination of spectral bands 5 (Near Infrared (NIR)), 4 (Red), and 3 (Green) of the Landsat OLI sensor covering the study area in the Cape Verde Peninsula region and Saloum River Delta, West Senegal, using February scenes: (<b>a</b>) 2015; (<b>b</b>) 2018; (<b>c</b>) 2020; (<b>d</b>) 2021; (<b>e</b>) 2022; (<b>f</b>) 2023.</p>
Full article ">Figure 6
<p>Land cover types in Senegal according to the FAO classification scheme. Software: QGIS v. 3.22. Map source: Author.</p>
Full article ">Figure 7
<p>Classification of the Landsat images from 2020 covering the Cape Verde Peninsula region and the Saloum River Delta, West Senegal: (<b>a</b>) 2015; (<b>b</b>) 2018; (<b>c</b>) 2020; (<b>d</b>) 2021; (<b>e</b>) 2022; (<b>f</b>) 2023.</p>
Full article ">Figure 8
<p>Results of the Support Vector Machine (SVM)-based classification of the Landsat images covering the Cape Verde Peninsula region and Saloum River Delta, West Senegal: (<b>a</b>) February 2015; (<b>b</b>) February 2018; (<b>c</b>) February 2020; (<b>d</b>) February 2021; (<b>e</b>) February 2022; (<b>f</b>) February 2023.</p>
Full article ">Figure 9
<p>Accuracy evaluated based on the pixel confidence levels with rejection probability values for the Landsat images covering the Cape Verde Peninsula region and Saloum River Delta, West Senegal: (<b>a</b>) 2015; (<b>b</b>) 2018; (<b>c</b>) 2020; (<b>d</b>) 2021; (<b>e</b>) 2022;(<b>f</b>) 2023.</p>
Full article ">
16 pages, 1401 KiB  
Article
Impact of Climate Change on Agricultural Production and Food Security: A Case Study in the Mekong River Delta of Vietnam
by Tran Trong Phuong, Tran Duc Vien, Cao Truong Son, Doan Thanh Thuy and Stefan Greiving
Sustainability 2024, 16(17), 7776; https://doi.org/10.3390/su16177776 - 6 Sep 2024
Viewed by 722
Abstract
Vietnam is a country highly vulnerable to climate change. Specifically, climate change has seriously impacted all aspects of Vietnam’s economic and social life, especially agricultural production. In this article, we focus on analyzing the current situation and the impacts of climate change on [...] Read more.
Vietnam is a country highly vulnerable to climate change. Specifically, climate change has seriously impacted all aspects of Vietnam’s economic and social life, especially agricultural production. In this article, we focus on analyzing the current situation and the impacts of climate change on agricultural production and food security in Vietnam, especially in the Mekong River Delta (MRD) region. Vietnam’s climate change scenarios (RCP4.5 and RCP 8.5) have warned of serious increases in temperature, rainfall, and sea level rises for the MRD in coming times. This will lead to a risk of flooding in nearly 50% of the region’s area and will seriously affect agricultural production in many aspects such as soil quality degradation, scarcity of water resources, increased droughts and floods, reduced crop productivity, and so on. These impacts will reduce Vietnam’s food supply capacity, but do not compromise national food security from a short-term perspective. Faced with this situation, the Government of Vietnam has implemented many comprehensive measures to transform agriculture towards ecology, sustainability, and low carbon emissions, with the goal of green growth and neutral carbon emissions by 2050. In particular, the focus is on combining nature-based solutions with the application of modern science and technology in agricultural production, raising awareness and the response capacity of domestic people, with international cooperation in addressing climate change issues. Full article
(This article belongs to the Special Issue Sustainable Agriculture and Food Security)
Show Figures

Figure 1

Figure 1
<p>Location of the Mekong Delta region on the map of Vietnam.</p>
Full article ">Figure 2
<p>Framework for studying the impact of climate change on agricultural production and food security in the Mekong Delta.</p>
Full article ">Figure 3
<p>Agricultural land area of the MRD affected by saltwater intrusion in 2015–2016.</p>
Full article ">Figure 4
<p>Estimated flooded land area in the MRD under the two scenarios RCP4.5 and RCP8.5.</p>
Full article ">
21 pages, 10428 KiB  
Article
Three Decades of Inundation Dynamics in an Australian Dryland Wetland: An Eco-Hydrological Perspective
by Indishe P. Senanayake, In-Young Yeo and George A. Kuczera
Remote Sens. 2024, 16(17), 3310; https://doi.org/10.3390/rs16173310 - 6 Sep 2024
Viewed by 600
Abstract
Wetland ecosystems are experiencing rapid degradation due to human activities, particularly the diversion of natural flows for various purposes, leading to significant alterations in wetland hydrology and their ecological functions. However, understanding and quantifying these eco-hydrological changes, especially concerning inundation dynamics, presents a [...] Read more.
Wetland ecosystems are experiencing rapid degradation due to human activities, particularly the diversion of natural flows for various purposes, leading to significant alterations in wetland hydrology and their ecological functions. However, understanding and quantifying these eco-hydrological changes, especially concerning inundation dynamics, presents a formidable challenge due to the lack of long-term, observation-based spatiotemporal inundation information. In this study, we classified wetland areas into ten equal-interval classes based on inundation probability derived from a dense, 30-year time series of Landsat-based inundation maps over an Australian dryland riparian wetland, Macquarie Marshes. These maps were then compared with three simplified vegetation patches in the area: river red gum forest, river red gum woodland, and shrubland. Our findings reveal a higher inundation probability over a small area covered by river red gum forest, exhibiting persistent inundation over time. In contrast, river red gum woodland and shrubland areas show fluctuating inundation patterns. When comparing percentage inundation with the Normalized Difference Vegetation Index (NDVI), we observed a notable agreement in peaks, with a lag time in NDVI response. A strong correlation between NDVI and the percentage of inundated area was found in the river red gum woodland patch. During dry, wet, and intermediate years, the shrubland patch consistently demonstrated similar inundation probabilities, while river red gum patches exhibited variable probabilities. During drying events, the shrubland patch dried faster, likely due to higher evaporation rates driven by exposure to solar radiation. However, long-term inundation probability exhibited agreement with the SAGA wetness index, highlighting the influence of topography on inundation probability. These findings provide crucial insights into the complex interactions between hydrological processes and vegetation dynamics in wetland ecosystems, underscoring the need for comprehensive monitoring and management strategies to mitigate degradation and preserve these vital ecosystems. Full article
(This article belongs to the Special Issue Remote Sensing for Land Degradation and Drought Monitoring II)
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>(<b>a</b>–<b>c</b>) location of the study area, Macquarie Marshes. (<b>d</b>) LiDAR-derived 1 m digital elevation model (DEM) of Macquarie Marshes. (<b>e</b>) Temporal average NDVI values over the Marshes as captured by the Landsat 8 collections from 2013 to 2020.</p>
Full article ">Figure 2
<p>Inundated area in Northern Marshes as captured by the inundation maps developed with Landsat 5, 7, and 8 datasets using the RaFMIC approach [<a href="#B31-remotesensing-16-03310" class="html-bibr">31</a>].</p>
Full article ">Figure 3
<p>Annual probability of inundation over the Marshes from 1988 to 2020, as captured by the Landsat 5-, 7-, and 8-based inundation maps [<a href="#B31-remotesensing-16-03310" class="html-bibr">31</a>]. Landsat 5-, 7-, and 8-derived maps are demarcated using red, blue, and green map borders, respectively.</p>
Full article ">Figure 4
<p>(<b>a</b>) Annual stream discharge to the Marshes from Marebone Weir (#421090) and Marebone Break (#421088), and (<b>b</b>) annual rainfall captured by the PERSIANN data and rain gauges, #051042 and #051057.</p>
Full article ">Figure 5
<p>(<b>a</b>–<b>c</b>) Inundation probability maps derived from Landsat 5-, 7-, and 8-based inundation maps classified into ten inundation probability classes. (<b>d</b>) Classified inundation probability map based on all the Landsat-derived inundation maps (i.e., Landsat 5, 7, and 8), collectively. <span class="html-italic">n</span> is the number of inundation maps used to develop each probability of inundation map. (<b>e</b>) Classification of vegetation over the Northern Marshes in 2013 based on Bowen et al. [<a href="#B75-remotesensing-16-03310" class="html-bibr">75</a>].</p>
Full article ">Figure 6
<p>Time record of inundated area in each probability of inundation class over the Northern Marshes as captured collectively by the inundation maps derived from Landsat 5, 7, and 8 image collections. Orange dots indicate each data point.</p>
Full article ">Figure 7
<p>(<b>a</b>) The three simplified vegetation patches dominated by river red gum (RRG) forest, RRG woodland, and shrubland (which encompasses common reed, mixed marsh/water couch, and terrestrial vegetation) over the Northern Marshes. (<b>b</b>–<b>d</b>) Time series of percentage areal inundation over the three vegetation patches as collectively captured by the Landsat 5-, 7-, and 8-based inundation maps.</p>
Full article ">Figure 8
<p>Time series between the percentage areal inundation and areal average NDVI values over the three vegetation patches: (<b>a</b>) river red gum forest, (<b>b</b>) river red gum woodland, and (<b>c</b>) shrubland in the Northern Marshes as captured by the Landsat 7-based inundation and NDVI products. (<b>d</b>–<b>f</b>) Linear regressions between Percentage areal inundation and NDVI over the three vegetation patches in the Northern Marshes as captured collectively by Landsat 5,- 7-, and 8-based products.</p>
Full article ">Figure 9
<p>Inundation probability of generally (<b>a</b>) dry, (<b>b</b>) normal, and (<b>c</b>) wet years as captured by Landsat 8-based inundation maps over the Northern Marshes with the three vegetation patches. (<b>d</b>) Average rainfall over the area from gauges #051042 and #051057 with the total discharge of Marebone Weir and Marebone Break from 2013 to 2019.</p>
Full article ">Figure 10
<p>(<b>a</b>) Two drying events as captured by the Landsat 5-based inundation maps from 8 September 1990 to 30 January 1991 (six inundation maps) and from 29 August 1998 to 5 February 1999 (five inundation maps). (<b>b</b>,<b>c</b>) Probability of inundation captured by the Landsat 5-based inundation maps during these two drying events. (<b>d</b>) SAGA Wetness Index (SWI) over the area derived from the 1 m LiDAR-derived DEM.</p>
Full article ">
20 pages, 14326 KiB  
Article
The Impact of Sandbars on Bank Protection Structures in Low-Land Reaches: Case of Ganges and Brahmaputra-Jamuna
by Shampa, Hussain Muhammad Muktadir, Israt Jahan Nejhum, A. K. M. Saiful Islam, Md. Munsur Rahman and G. M. Tarekul Islam
Water 2024, 16(17), 2523; https://doi.org/10.3390/w16172523 - 5 Sep 2024
Viewed by 584
Abstract
Sandbars are an integral part of the alluvial river’s geophysical system due to these rivers’ wide sediment availability and varied transport capacity. The sandbars’ evolution and translation considerably influence the stability of the riverbank. However, while designing the riverbank protection structures (RBPS), the [...] Read more.
Sandbars are an integral part of the alluvial river’s geophysical system due to these rivers’ wide sediment availability and varied transport capacity. The sandbars’ evolution and translation considerably influence the stability of the riverbank. However, while designing the riverbank protection structures (RBPS), the impact of such sandbars is often overlooked, as the evolution of such bars is quite uncertain in terms of location, amplitude, and translation. This study evaluates the localized impact of sandbars on bank protection structures in two types of alluvial rivers: meandering (Ganges) and braided (Brahmaputra-Jamuna), utilizing time series satellite images, hydraulic characteristics, and numerical modeling. We found that sandbar development initiates width adjustment in both meandering and braided rivers when the ratio of width to depth surpasses 90. In the case of meandering rivers, riverbank erosion mostly occurs as a result of the presence of alternate bars or point bars. Sandbars in a meandering river (Ganges) can lead to an approximate 18% increase in flow depth. The depth-averaged velocity is anticipated to rise by approximately 29%, and the tractive force may increase by a factor of 1.6. On the other hand, the braided river (the Brahmaputra-Jamuna) underwent significant bank erosion due to the presence of both free unit and hybrid types of bars. In such rivers, the depth of the flow may experience a notable increase of 18%, while the depth-averaged velocity undergoes an approximate increase of 50%, and the tractive force has the potential to grow by a factor of 5.3. Consequently, we recommend allowing the natural evolution of sandbars while preserving the riverbank (where needed only) through RBPS, considering these additional loads. Full article
(This article belongs to the Section Hydrology)
Show Figures

Figure 1

Figure 1
<p>Study areas for (<b>A</b>) Ganges and (<b>B</b>) Brahmaputra-Jamuna. The bottom figure shows the longitudinal profile of the river’s valley from source to sink.</p>
Full article ">Figure 2
<p>Definition sketch of channel width and bar width.</p>
Full article ">Figure 3
<p>Grid and bathymetry were used in the models. BWDB commonly refers levels to the Public Works Datum (PWD), which is 0.46 m below the MSL.</p>
Full article ">Figure 4
<p>The boundary conditions used in the models. (<b>a</b>) Ganges model and (<b>b</b>) Brahmaputra-Jamuna model.</p>
Full article ">Figure 5
<p>The comparison of simulated and measured bed levels. (<b>a</b>) Ganges and (<b>b</b>) Brahmaputra-Jamuna.</p>
Full article ">Figure 5 Cont.
<p>The comparison of simulated and measured bed levels. (<b>a</b>) Ganges and (<b>b</b>) Brahmaputra-Jamuna.</p>
Full article ">Figure 6
<p>Planform changes of Ganges around the study area for the last five decades, derived from Landsat images.</p>
Full article ">Figure 7
<p>The changes in bar alignment for the last two decades. The background satellite image is the Landsat image, dated February 2023.</p>
Full article ">Figure 8
<p>The relationship between channel width, bar width, and bank erosion of Ganges reach.</p>
Full article ">Figure 9
<p>The relationship between the bank erosion and bar–channel properties of Ganges reach. (<b>a</b>) relationship between bar-channel width ratio to bank erosion (<b>b</b>) relationship between bar-channel width ratio to bar amplitude.</p>
Full article ">Figure 10
<p>Planform changes of Brahmaputra-Jamuna around the study area for the last five decades, derived from Landsat images.</p>
Full article ">Figure 11
<p>The changes in bar alignment in the study have reached off the Brahmaputra-Jamuna River over the last two decades. The background satellite image is the Landsat image, dated February 2023.</p>
Full article ">Figure 12
<p>The relationship between channel width, bar width, and bank erosion of Brahmaputra-Jamuna.</p>
Full article ">Figure 13
<p>The relationship between bank erosion and bar–channel properties of Brahmaputra-Jamuna reach. (<b>a</b>) relationship between bar-channel width ratio to bank erosion (<b>b</b>) relationship between bar-channel width ratio to bar amplitude.</p>
Full article ">Figure 14
<p>The distribution of (<b>top</b>) flow depth, the (<b>middle</b>) velocity at the peak discharge, and the (<b>bottom</b>) bed level change under the M<sub>BS</sub> and M<sub>EX</sub> scenarios at the Ganges.</p>
Full article ">Figure 15
<p>The distribution of (<b>top</b>) flow depth, the (<b>middle</b>) velocity at the peak discharge, and the (<b>bottom</b>) bed level change under the B<sub>bs</sub> and B<sub>ex</sub> scenarios at Brahmaputra-Jamuna.</p>
Full article ">Figure 16
<p>The generation of vortex flow generated near the sandbar: (<b>a</b>) Ganges and (<b>b</b>) Brahmaputra-Jamuna.</p>
Full article ">
Back to TopTop