[go: up one dir, main page]
More Web Proxy on the site http://driver.im/
You seem to have javascript disabled. Please note that many of the page functionalities won't work as expected without javascript enabled.
 
 
Sign in to use this feature.

Years

Between: -

Subjects

remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline

Journals

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

Search Results (6,410)

Search Parameters:
Keywords = topography

Order results
Result details
Results per page
Select all
Export citation of selected articles as:
28 pages, 20186 KiB  
Article
Long-Term Statistical Analysis of Severe Weather and Climate Events in Greece
by Vassiliki Kotroni, Antonis Bezes, Stavros Dafis, Dimitra Founda, Elisavet Galanaki, Christos Giannaros, Theodore Giannaros, Athanasios Karagiannidis, Ioannis Koletsis, George Kyros, Konstantinos Lagouvardos, Katerina Papagiannaki and Georgios Papavasileiou
Atmosphere 2025, 16(1), 105; https://doi.org/10.3390/atmos16010105 (registering DOI) - 18 Jan 2025
Abstract
The Mediterranean faces frequent heavy precipitation, deadly heatwaves, and wildfires fueled by its climate. Greece, with its complex topography, experiences severe and extreme weather events that have escalated in recent years and are projected to continue rising under future climate conditions. This paper [...] Read more.
The Mediterranean faces frequent heavy precipitation, deadly heatwaves, and wildfires fueled by its climate. Greece, with its complex topography, experiences severe and extreme weather events that have escalated in recent years and are projected to continue rising under future climate conditions. This paper analyzes severe weather events and trends in Greece from 2010 to 2023, leveraging data from an expanded network of weather stations spanning across Greece, as well as long-term meteorological data from the reference weather station in the center of Athens. The focus includes analysis of heat waves, intense rainfall and droughts, thunderstorms, hail, tornadoes, and fire weather conditions. The societal impact of severe weather events is also discussed. The paper aims to provide both long-term (1901–2023) and recent year analyses (2010–2023). The main results show that between 2010 and 2023, Greece experienced: nearly one heatwave per summer; heavy rainfall events, most common in winter and autumn, showing a significant increase, particularly in the eastern Aegean and western continental Greece; dry spells, which are longest in southern Greece; thunderstorm and hail events peaking in spring and summer; fire weather conditions and risk peaking in southern Greece. Finally, societal impacts from weather hazards have increased in Greece over the past 14 years, with flash floods being the most frequent and damaging events, while public preparedness and effective risk communication remain low. Full article
Show Figures

Figure 1

Figure 1
<p>Topographic map of Greece (elevation in m) and locations referred to in the text.</p>
Full article ">Figure 2
<p>Frequency (total number) of hot days, heat waves (HWs) and HW days at NOA for 4 consecutive climatic sub-periods (1901–1930, 1931–1960, 1961–1990, and 1991–2020).</p>
Full article ">Figure 3
<p>Frequency (number of days) of different extreme precipitation categories in Athens (NOA) for 4 consecutive climatic sub-periods (1901–2020).</p>
Full article ">Figure 4
<p>The 2010–2023 average values of observed HWN for summer.</p>
Full article ">Figure 5
<p>The 2010–2023 average values of observed HWD for summer.</p>
Full article ">Figure 6
<p>The 2010–2023 average values of observed HWF for summer.</p>
Full article ">Figure 7
<p>Thunder days in Greece over the 52 prefectures based on ZEUS lightning data.</p>
Full article ">Figure 8
<p>Hail reports in Greece during the period 2010–2023.</p>
Full article ">Figure 9
<p>Maximum hail diameter reports in Greece during the period 2010–2023.</p>
Full article ">Figure 10
<p>Tornado reports in Greece during the period 2010–2023. Blue triangles denote waterspouts, while the red triangles denote tornadoes over land.</p>
Full article ">Figure 11
<p>The 2010–2023 average values of observed R20 for (<b>a)</b> Winter, (<b>b</b>) Spring, (<b>c</b>) Summer, and (<b>d</b>) Autumn.</p>
Full article ">Figure 12
<p>The 2010–2023 average values of observed SPII for (<b>a</b>) Winter, (<b>b</b>) Spring, (<b>c</b>) Summer, and (<b>d</b>) Autumn.</p>
Full article ">Figure 13
<p>Seasonal R20 trends for (<b>a</b>) Winter, (<b>b</b>) Spring, (<b>c</b>) Summer, and (<b>d</b>) Autumn.</p>
Full article ">Figure 14
<p>Seasonal SPII trends for (<b>a</b>) Winter, (<b>b</b>) Spring, (<b>c</b>) Summer, and (<b>d</b>) Autumn.</p>
Full article ">Figure 14 Cont.
<p>Seasonal SPII trends for (<b>a</b>) Winter, (<b>b</b>) Spring, (<b>c</b>) Summer, and (<b>d</b>) Autumn.</p>
Full article ">Figure 15
<p>The 2010–2023 average values of observed CDD for: (<b>a</b>) Winter, (<b>b</b>) Spring, (<b>c</b>) Summer, and (<b>d</b>) Autumn.</p>
Full article ">Figure 16
<p>Seasonal CDD trends for (<b>a</b>) Winter, (<b>b</b>) Spring, (<b>c</b>) Summer, and (<b>d</b>) Autumn.</p>
Full article ">Figure 17
<p>Annual mean values of the period 2010–2023 for: (<b>a</b>) FWImax, (<b>b</b>) FWI30, and (<b>c</b>) FLD.</p>
Full article ">Figure 18
<p>Annual distribution of weather-related phenomena in Greece in 2000–2023 by phenomenon.</p>
Full article ">Figure 19
<p>Annual distribution of fatal flood events, associated flood fatalities in Greece and average fatalities per event in 2000–2023.</p>
Full article ">
25 pages, 25562 KiB  
Article
A Mapping Method Fusing Forward-Looking Sonar and Side-Scan Sonar
by Hong Liu, Xiufen Ye, Hanwen Zhou and Hanjie Huang
J. Mar. Sci. Eng. 2025, 13(1), 166; https://doi.org/10.3390/jmse13010166 (registering DOI) - 18 Jan 2025
Viewed by 32
Abstract
In modern ocean exploration, forward-looking sonar (FLS) provides real-time 2D imaging of the seabed ahead, but its detection range is relatively limited. Conversely, side-scan sonar (SSS) enables large-scale imaging of the seabed during movement but struggles to effectively image areas directly beneath the [...] Read more.
In modern ocean exploration, forward-looking sonar (FLS) provides real-time 2D imaging of the seabed ahead, but its detection range is relatively limited. Conversely, side-scan sonar (SSS) enables large-scale imaging of the seabed during movement but struggles to effectively image areas directly beneath the sensor. Integrating FLS and SSS offers a promising solution by leveraging their complementary strengths to achieve comprehensive seabed mapping. However, no prior research has explored this fusion approach. This paper presents a novel method for FLS and SSS fusion mapping. Firstly, a novel sonar image enhancement method based on equalization is proposed, enabling simultaneous enhancement and grayscale unification of two sonar images. Additionally, an effective area extraction approach for FLS images, grounded on the approximate erosion method, is introduced to produce high-quality FLS mapping. Furthermore, by examining the data distribution in FLS and SSS mappings, the standard deviation of these datasets is utilized to refine the grayscale distribution of FLS mapping, thereby enhancing the grayscale distribution similarity between the two mapping results. Finally, FLS map data are seamlessly integrated into the gaps of the SSS map, resulting in a fused, comprehensive seabed representation. Large-scale experiments demonstrate that the proposed method effectively combines the strengths of FLS and SSS, producing complete and detailed seabed topography maps. Simultaneously, numerous ablation experiments are conducted to evaluate the impact of various parameters on fusion mapping, providing guidelines for selecting the optimal parameters. This fusion approach, thus, holds significant practical value for ocean exploration and seabed mapping applications. Full article
(This article belongs to the Section Ocean Engineering)
Show Figures

Figure 1

Figure 1
<p>Typical SSS image. The (<b>left</b>) panel depicts a rock pile, while the (<b>right</b>) panel shows a sand pattern on the seafloor.</p>
Full article ">Figure 2
<p>Typical FLS images. (<b>a</b>) The FLS imaging result of a rock pile on the seabed. (<b>b</b>) Corresponding FLS fan-shaped image. The areas circled by red lines represent the effective imaging regions of the FLS, while the remaining areas constitute the background.</p>
Full article ">Figure 3
<p>Schematic diagram of the unmanned ship carrying FLS and SSS for seabed scanning. The red lines indicate the square-cone imaging area of the FLS, while the yellow lines represent the fan-shaped imaging area of the SSS. The center line connects the centers of the FLS and SSS imaging areas and is parallel to the center lines of the FLS and SSS devices.</p>
Full article ">Figure 4
<p>Data processing flowchart of the proposed method.</p>
Full article ">Figure 5
<p>Flowchart of grayscale equalization and geocoding method for SSS. (<b>a</b>) Original SSS image. (<b>b</b>) Filtered SSS image. (<b>c</b>) Grayscale equalized SSS image. (<b>d</b>) Resulting SSS map.</p>
Full article ">Figure 6
<p>Schematic diagram of the coordinate system definition for SSS geocoding. <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>P</mi> </mrow> <mrow> <mi>S</mi> <mi>S</mi> <mi>S</mi> </mrow> </msub> </mrow> </semantics></math> represents the physical center coordinates of the SSS device in WCS. <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>P</mi> </mrow> <mrow> <mi>r</mi> </mrow> </msub> </mrow> </semantics></math> is the coordinate of a point on the seafloor in the WCS.</p>
Full article ">Figure 7
<p>Flowchart of grayscale equalization and geocoding method for FLS. (<b>a</b>) Original FLS image. (<b>b</b>) FLS image processed with the mean filtering method. (<b>c</b>) Binary image distinguishing the effective imaging area from the noise background, calculated using the approximate erosion method. The white area represents the effective imaging area (pixel value is 1), while the black area represents the noise background (pixel value is 0). (<b>d</b>) Grayscale-balanced image. (<b>e</b>) Effective imaging area map after grayscale equalization and removal of the noise background. (<b>f</b>) FLS mapping result after geocoding the grayscale-balanced image.</p>
Full article ">Figure 8
<p>Schematic diagram for effective region extraction based on the approximate erosion method. (<b>a</b>) illustration of the calculation process on an FLS Image. (<b>b</b>) Binary image representing the calculated effective imaging region. <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>x</mi> </mrow> <mrow> <mi>u</mi> <mi>p</mi> </mrow> </msub> </mrow> </semantics></math> keeps increasing and <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>x</mi> </mrow> <mrow> <mi>d</mi> <mi>o</mi> <mi>w</mi> <mi>n</mi> </mrow> </msub> </mrow> </semantics></math> keeps decreasing, approaching the effective imaging area in the middle step by step. <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>W</mi> </mrow> <mrow> <mi>w</mi> <mi>i</mi> <mi>d</mi> <mi>t</mi> <mi>h</mi> </mrow> </msub> </mrow> </semantics></math> is the width of FLS image. <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>H</mi> </mrow> <mrow> <mi>h</mi> <mi>e</mi> <mi>i</mi> <mi>g</mi> <mi>h</mi> <mi>t</mi> </mrow> </msub> </mrow> </semantics></math> is the height of FLS image.</p>
Full article ">Figure 9
<p>Schematic diagram of FLS imaging coordinate system. <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>P</mi> </mrow> <mrow> <mi>F</mi> <mi>L</mi> <mi>S</mi> </mrow> </msub> </mrow> </semantics></math> is the coordinate of the FLS physical center in the WCS. <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>P</mi> </mrow> <mrow> <mi>n</mi> </mrow> </msub> </mrow> </semantics></math> is the coordinate of a point on the seabed in the WCS. (<b>a</b>) Three-dimensional diagram. (<b>b</b>) Bottom view. (<b>c</b>) Side view.</p>
Full article ">Figure 10
<p>A schematic diagram of data sampling in the overlapping region. The blue square represents the overlapping area used to calculate the average standard deviation, corresponding to a flat seabed. The green lines represent the positions where sampling is performed within the overlapping region.</p>
Full article ">Figure 11
<p>The pixel data of the two sonar maps corresponding to the green line in <a href="#jmse-13-00166-f010" class="html-fig">Figure 10</a>. (<b>a</b>) The values of all the pixels along the line. (<b>b</b>) The statistical histogram of these data. The horizontal axis represents the original pixel values divided by 4 due to the small sample size.</p>
Full article ">Figure 12
<p>The pixel data of the two sonar maps after adjustment. (<b>a</b>) The values of all the pixels along the line. (<b>b</b>) The statistical histogram of these data. The horizontal axis represents the original pixel values divided by 4 due to the small sample size.</p>
Full article ">Figure 13
<p>The FLS and SSS maps after grayscale adjustment. (<b>a</b>) The result of directly overlaying the FLS map onto the SSS map. (<b>b</b>) The result of filling the gaps in the SSS map with the FLS map data.</p>
Full article ">Figure 14
<p>Unmanned surface vessel and equipment used to collect data.</p>
Full article ">Figure 15
<p>(<b>a</b>) ES900 side scan sonar. (<b>b</b>) M750d forward-looking sonar.</p>
Full article ">Figure 16
<p>Schematic diagram of data sampling in the overlapping region. The green lines represent the positions where single ping sampling is conducted within the overlapping region. The blue square represents the overlapping area used to calculate the average standard deviation.</p>
Full article ">Figure 17
<p>The sampling data generated by different filter parameters, <math display="inline"><semantics> <mrow> <msup> <mrow> <mi>M</mi> </mrow> <mrow> <mi>F</mi> <mi>L</mi> <mi>S</mi> </mrow> </msup> </mrow> </semantics></math> in Ping 1. (<b>a</b>) <math display="inline"><semantics> <mrow> <msup> <mrow> <mi>M</mi> </mrow> <mrow> <mi>F</mi> <mi>L</mi> <mi>S</mi> </mrow> </msup> <mo>=</mo> <mn>10</mn> </mrow> </semantics></math>. (<b>b</b>) <math display="inline"><semantics> <mrow> <msup> <mrow> <mi>M</mi> </mrow> <mrow> <mi>F</mi> <mi>L</mi> <mi>S</mi> </mrow> </msup> <mo>=</mo> <mn>50</mn> </mrow> </semantics></math>. (<b>c</b>) <math display="inline"><semantics> <mrow> <msup> <mrow> <mi>M</mi> </mrow> <mrow> <mi>F</mi> <mi>L</mi> <mi>S</mi> </mrow> </msup> <mo>=</mo> <mn>100</mn> </mrow> </semantics></math>. (<b>d</b>) <math display="inline"><semantics> <mrow> <msup> <mrow> <mi>M</mi> </mrow> <mrow> <mi>F</mi> <mi>L</mi> <mi>S</mi> </mrow> </msup> <mo>=</mo> <mn>200</mn> </mrow> </semantics></math>.</p>
Full article ">Figure 18
<p>Sonar maps generated by different filter parameters, <math display="inline"><semantics> <mrow> <msup> <mrow> <mi>M</mi> </mrow> <mrow> <mi>F</mi> <mi>L</mi> <mi>S</mi> </mrow> </msup> </mrow> </semantics></math>. (<b>a</b>) <math display="inline"><semantics> <mrow> <msup> <mrow> <mi>M</mi> </mrow> <mrow> <mi>F</mi> <mi>L</mi> <mi>S</mi> </mrow> </msup> <mo>=</mo> <mn>10</mn> </mrow> </semantics></math>. (<b>b</b>) <math display="inline"><semantics> <mrow> <msup> <mrow> <mi>M</mi> </mrow> <mrow> <mi>F</mi> <mi>L</mi> <mi>S</mi> </mrow> </msup> <mo>=</mo> <mn>50</mn> </mrow> </semantics></math>. (<b>c</b>) <math display="inline"><semantics> <mrow> <msup> <mrow> <mi>M</mi> </mrow> <mrow> <mi>F</mi> <mi>L</mi> <mi>S</mi> </mrow> </msup> <mo>=</mo> <mn>100</mn> </mrow> </semantics></math>. (<b>d</b>) <math display="inline"><semantics> <mrow> <msup> <mrow> <mi>M</mi> </mrow> <mrow> <mi>F</mi> <mi>L</mi> <mi>S</mi> </mrow> </msup> <mo>=</mo> <mn>200</mn> </mrow> </semantics></math>.</p>
Full article ">Figure 19
<p>The sampling data generated by different filter parameters <math display="inline"><semantics> <mrow> <msup> <mrow> <mi>M</mi> </mrow> <mrow> <mi>S</mi> <mi>S</mi> <mi>S</mi> </mrow> </msup> </mrow> </semantics></math> in Ping 1. (<b>a</b>) <math display="inline"><semantics> <mrow> <msup> <mrow> <mi>M</mi> </mrow> <mrow> <mi>S</mi> <mi>S</mi> <mi>S</mi> </mrow> </msup> <mo>=</mo> <mn>5</mn> </mrow> </semantics></math>. (<b>b</b>) <math display="inline"><semantics> <mrow> <msup> <mrow> <mi>M</mi> </mrow> <mrow> <mi>S</mi> <mi>S</mi> <mi>S</mi> </mrow> </msup> <mo>=</mo> <mn>20</mn> </mrow> </semantics></math>. (<b>c</b>) <math display="inline"><semantics> <mrow> <msup> <mrow> <mi>M</mi> </mrow> <mrow> <mi>S</mi> <mi>S</mi> <mi>S</mi> </mrow> </msup> <mo>=</mo> <mn>40</mn> </mrow> </semantics></math>. (<b>d</b>) <math display="inline"><semantics> <mrow> <msup> <mrow> <mi>M</mi> </mrow> <mrow> <mi>S</mi> <mi>S</mi> <mi>S</mi> </mrow> </msup> <mo>=</mo> <mn>60</mn> </mrow> </semantics></math>.</p>
Full article ">Figure 20
<p>Sonar maps generated by different filter parameters, <math display="inline"><semantics> <mrow> <msup> <mrow> <mi>M</mi> </mrow> <mrow> <mi>S</mi> <mi>S</mi> <mi>S</mi> </mrow> </msup> </mrow> </semantics></math>. (<b>a</b>) <math display="inline"><semantics> <mrow> <msup> <mrow> <mi>M</mi> </mrow> <mrow> <mi>S</mi> <mi>S</mi> <mi>S</mi> </mrow> </msup> <mo>=</mo> <mn>5</mn> </mrow> </semantics></math>. (<b>b</b>) <math display="inline"><semantics> <mrow> <msup> <mrow> <mi>M</mi> </mrow> <mrow> <mi>S</mi> <mi>S</mi> <mi>S</mi> </mrow> </msup> <mo>=</mo> <mn>20</mn> </mrow> </semantics></math>. (<b>c</b>) <math display="inline"><semantics> <mrow> <msup> <mrow> <mi>M</mi> </mrow> <mrow> <mi>S</mi> <mi>S</mi> <mi>S</mi> </mrow> </msup> <mo>=</mo> <mn>40</mn> </mrow> </semantics></math>. (<b>d</b>) <math display="inline"><semantics> <mrow> <msup> <mrow> <mi>M</mi> </mrow> <mrow> <mi>S</mi> <mi>S</mi> <mi>S</mi> </mrow> </msup> <mo>=</mo> <mn>60</mn> </mrow> </semantics></math>.</p>
Full article ">Figure 21
<p>Sonar map results generated by different parameters, <math display="inline"><semantics> <mrow> <mi>κ</mi> </mrow> </semantics></math>. (<b>a</b>) <math display="inline"><semantics> <mrow> <mi>κ</mi> <mo>=</mo> <mn>30</mn> </mrow> </semantics></math>. (<b>b</b>) <math display="inline"><semantics> <mrow> <mi>κ</mi> <mo>=</mo> <mn>40</mn> </mrow> </semantics></math>. (<b>c</b>) <math display="inline"><semantics> <mrow> <mi>κ</mi> <mo>=</mo> <mn>5</mn> </mrow> </semantics></math>0. (<b>d</b>) <math display="inline"><semantics> <mrow> <mi>κ</mi> <mo>=</mo> <mn>60</mn> </mrow> </semantics></math>.</p>
Full article ">Figure 22
<p>Comparison of large-scale fusion map and SSS map results (distance of more than 800 m). The above image represents the result of seabed mapping using SSS alone. The center images provide a zoomed-in view of a specific local area. The image below shows the seabed mapping results generated by the proposed method.</p>
Full article ">Figure 23
<p>Comparison of large-scale fusion map and SSS map results (data in Sanya, China). The above image represents the result of seabed mapping using SSS alone. The center images provide a zoomed-in view of a specific local area. The image below shows the seabed mapping results generated by the proposed method.</p>
Full article ">Figure 24
<p>Fusion maps generated using different image filtering methods. (<b>a</b>) Bilateral filtering. (<b>b</b>) Mean filtering. (<b>c</b>) Median filtering. (<b>d</b>) Gaussian filtering.</p>
Full article ">Figure 25
<p>Fusion maps generated by different fusion methods. (<b>a</b>) Fixed-weight fusion method. (<b>b</b>) This study. (<b>c</b>) Distance-weighted fusion method. (<b>d</b>) Maximum-value fusion method.</p>
Full article ">Figure 25 Cont.
<p>Fusion maps generated by different fusion methods. (<b>a</b>) Fixed-weight fusion method. (<b>b</b>) This study. (<b>c</b>) Distance-weighted fusion method. (<b>d</b>) Maximum-value fusion method.</p>
Full article ">
24 pages, 5440 KiB  
Article
Environmental Covariates for Sampling Optimization and Pest Prediction in Soybean Crops
by Cenneya Lopes Martins, Maiara Pusch, Wesley Augusto Conde Godoy and Lucas Rios do Amaral
AgriEngineering 2025, 7(1), 21; https://doi.org/10.3390/agriengineering7010021 (registering DOI) - 18 Jan 2025
Viewed by 71
Abstract
Insect pest infestations can vary due to spatial differences in microclimates and food availability within agroecosystems. Covariates can reflect these environmental conditions. This study tested whether using environmental covariates in two-phase sample optimization improved the spatial predictions for soybean insect pests. During the [...] Read more.
Insect pest infestations can vary due to spatial differences in microclimates and food availability within agroecosystems. Covariates can reflect these environmental conditions. This study tested whether using environmental covariates in two-phase sample optimization improved the spatial predictions for soybean insect pests. During the 2021–2022 crop season, insect pest samples were collected at 50 georeferenced points in a commercial soybean field in Brazil, alongside data on environmental covariates such as vegetation indices, soil properties, terrain topography, and distances from riparian areas. Three covariates were selected using correlation and principal component analysis (PCA). In the 2022–2023 crop season, sample designs were optimized using the iterative algorithm optimization of sample configurations using spatial simulated annealing (SPSANN) using the selected covariates, resulting in two optimized designs that were compared to a regular grid. Data from the three sampling designs comprising 50 points were evaluated using geostatistical methods, regression analysis (pest abundance), and classification (pest presence or absence) via the random forest algorithm. The data showed no spatial dependence, making using geostatistical interpolators inappropriate. However, a multi-objective optimized sampling design, tailored to refine configurations for identifying and estimating variograms and spatial trends essential for spatial interpolation, produced the most accurate predictions. Therefore, a two-phase sample optimization with prior in situ selection of environmental covariates improves pest predictions in agricultural systems, contributing to more efficient and sustainable agricultural management. Full article
Show Figures

Figure 1

Figure 1
<p>Diagram of the two-phase sample optimization research using environmental covariates.</p>
Full article ">Figure 2
<p>Location of the experimental area showing the two fields. Cartographic base: IBGE, 2023. Basemap: Google Satellite.</p>
Full article ">Figure 3
<p>Environmental covariates. Vegetation indices (cycle image): EVI, NDVI, NDRE, SFDVI, and DVI ((<b>A</b>–<b>E</b>), respectively); soil clay content (<b>F</b>); slope (<b>G</b>); river distance (<b>H</b>), and riparian forest distance (<b>I</b>).</p>
Full article ">Figure 4
<p>Sampling designs in Phase 1 (<b>A</b>) and Phase 2 (<b>B</b>–<b>E</b>). (<b>A</b>) Phase 1: 28 optimized MSSD sampling points combined with 22 random points (50 points). Phase 2: Regular grid (<b>B</b>), optimized CORR design (<b>C</b>), optimized SPAN design (<b>D</b>), each with 50 points, and external dataset (20 points).</p>
Full article ">Figure 5
<p>PCA results using the median of pests and the mean of the VIs in the soybean cycle.</p>
Full article ">Figure 6
<p>Scatter plots and metrics of RF regression modeling using environmental covariates in predicting the total pests using the regular (squares), CORR (stars), and SPAN (triangles) sampling designs ((<b>A</b>–<b>C</b>), respectively), and the prediction of <span class="html-italic">E. heros</span> in the same sampling designs (<b>D</b>–<b>F</b>). The red line represents the 1:1 ideal relationship (observed = predicted), while the dashed line indicates the regression line of the model predictions.</p>
Full article ">Figure 7
<p>Scatter plots and external validation metrics for the total pest predictions using the regular (squares), CORR (stars), and SPAN (triangles) sampling designs ((<b>A</b>–<b>C</b>), respectively). The red line represents the 1:1 ideal relationship (observed = predicted), while the dashed line indicates the regression line of the model predictions.</p>
Full article ">Figure 8
<p>Environmental covariates, pest sampling points, and prediction maps. (<b>A</b>–<b>C</b>) Environmental covariates selected in phase 1: soil clay content, NDVI, and distance from river. (<b>D</b>–<b>F</b>) Predicted maps using the RF regression algorithm, with environmental covariates as predictors of total pests abundance in the regular, CORR, and SPAN sampling designs, respectively. (<b>G</b>–<b>I</b>) Predicted maps using the RF classifier algorithm, with environmental covariates as predictors of the presence and absence of <span class="html-italic">E. heros</span> in the regular, CORR, and SPAN sampling designs, respectively.</p>
Full article ">
7 pages, 1757 KiB  
Case Report
Combined Multilayered Amniotic Membrane Graft and Fibrin Glue as a Surgical Management of Limbal Dermoid Cyst
by Maria Poddi, Vito Romano, Alfredo Borgia, Floriana Porcaro, Carlo Cagini and Marco Messina
J. Clin. Med. 2025, 14(2), 607; https://doi.org/10.3390/jcm14020607 (registering DOI) - 18 Jan 2025
Viewed by 116
Abstract
Background/Objectives: To report the cosmetic, clinical, and visual outcomes of a combined surgical approach for treating a corneal/limbal dermoid using excision and a three-layered amniotic membrane graft with fibrin glue. Methods: An 18-year-old female presented with impaired vision and ocular discomfort caused by [...] Read more.
Background/Objectives: To report the cosmetic, clinical, and visual outcomes of a combined surgical approach for treating a corneal/limbal dermoid using excision and a three-layered amniotic membrane graft with fibrin glue. Methods: An 18-year-old female presented with impaired vision and ocular discomfort caused by a prominent dome-shaped limbal congenital dermoid on the inferotemporal cornea, resulting in a significant aesthetic concern. A full assessment, including refraction, best-corrected visual acuity (BCVA), corneal topography, aberrometry and anterior segment OCT (AS-OCT) was conducted to plan the surgical approach. The dermoid was excised under peribulbar anaesthesia using manual lamellar dissection, followed by the application of 0.02% Mitomycin C and a multilayered amniotic membrane graft with fibrin glue. A bandage contact lens was applied and removed after three weeks, with postoperative treatment including topical antibiotics and steroids. Follow-ups were conducted on day 1, at 1 week, 3 weeks, 2 months, 6 months, 1 year, and 2 years. Results: Histopathological examination confirmed the mesoblastic nature of the lesion. Significant improvements in BCVA and ocular symptoms were observed. Corneal topography showed ocular surface regularization with reduction of high order aberrations and point spread function. AS-OCT showed complete integration of the amniotic membrane, with full epithelial coverage of the defect. The healing process was uneventful and the ocular surface remained stable throughout the entire follow-up, without complications or recurrence. Conclusions: This approach of dermoid excision, multilayered amniotic membrane and fibrin glue restored vision effectively, with notable improvements in ocular surface and cosmetic outcomes, without recurrence over two years. Full article
(This article belongs to the Special Issue Keratitis and Keratopathy: New Insights into Diagnosis and Treatment)
Show Figures

Figure 1

Figure 1
<p>Preoperative slit lamp photograph of inferotemporal limbal dermoid (<b>A</b>). Sirius topography image of anterior tangential curvature revealing a significant irregularity of the ocular surface (<b>B</b>). Aberrometric assessments showing Point Spread Function (PSF) with Strehl Ratio (<b>C</b>) and whole corneal aberrations (<b>D</b>,<b>F</b>). Heidelberg Spectralis anterior segment optical coherence tomography (AS − OCT cornea module) scan showing the hyperreflective lesion without a clear evidence of its depth (<b>E</b>).</p>
Full article ">Figure 2
<p>1 week postoperative slit lamp image showing the multilayered amniotic membranes graft perfectly covering the defect after the dermoid removal (<b>A</b>). Sirius topography images revealing a significant improvement in the regularity of the ocular surface (<b>B</b>). Aberrometric assessment with a substantial reduction in total corneal aberrations including both LOAs and HOAs, along with a notable improvement in PSF and Strehl ratio (<b>C</b>,<b>D</b>,<b>F</b>). Heidelberg Spectralis OCT images demonstrating a complete integration (graft) of the three-layered amniotic membrane into the corneoscleral surface (<b>E</b>).</p>
Full article ">Figure 3
<p>Slit lamp images taken two years post-surgery showing no signs of recurrence (<b>A</b>,<b>B</b>), proper corneal thickness (<b>C</b>), and no evidence of neovascularization or significant scarring (<b>A</b>,<b>B</b>). The epithelial map (<b>D</b>) and OCT scans (<b>E</b>) show a compensatory epithelial thickening in the inferotemporal cornea (between 63 µm and 68 µm).</p>
Full article ">
14 pages, 1021 KiB  
Article
Temperature: An Influencing Factor on the Rheological and Energetic Parameters of Acid Pressure Technology Operations
by Gabriel Hernández-Ramírez, Antonio Bernardo-Sánchez, Aristides Alejandro Legrá-Lobaina, Laura Álvarez de Prado, Rodney Martínez-Rojas, Liudmila Pérez-García, Leonel Garcell-Puyáns, Jose Fernández-Ordás and Javier Menéndez
Minerals 2025, 15(1), 86; https://doi.org/10.3390/min15010086 (registering DOI) - 17 Jan 2025
Viewed by 174
Abstract
In this work, a study was carried out on the effect of temperature as the main influencing factor on the rheological behavior of lateritic suspensions, raw material for the operations of pressurized acid technology (HPAL) used to obtain nickel and cobalt from of [...] Read more.
In this work, a study was carried out on the effect of temperature as the main influencing factor on the rheological behavior of lateritic suspensions, raw material for the operations of pressurized acid technology (HPAL) used to obtain nickel and cobalt from of oxidized ores. From studies of X-ray diffraction, X-ray fluorescence, particle size analyzer, and mathematical modeling, the behavior of the interactions and rheological characteristics of the analyzed samples were obtained. In this study, it was concluded that the use of mathematical models that relate the temperature up to 90 °C and the energy parameters of the pumping system of flows, loads, hydraulic losses, power, and efficiency would allow finding ways to increase and stabilize the flow of fed hydromixture with a flow rate of 1600 m3/h and a solids concentration of 48% (w/w) and guarantee the efficiency of the technological process. Full article
(This article belongs to the Section Mineral Processing and Extractive Metallurgy)
23 pages, 10808 KiB  
Article
Disaster-Pregnant Environment Stability Evaluation of Geohazards in the Yellow River–Huangshui River Valley, China
by Tengyue Zhang, Qiang Zhou, Weidong Ma, Yuan Gao, Hanmei Li and Qiuyang Zhang
Sustainability 2025, 17(2), 732; https://doi.org/10.3390/su17020732 (registering DOI) - 17 Jan 2025
Viewed by 207
Abstract
This study aims to identify the key factors contributing to the destabilization of the geohazard disaster-pregnant environment in the Yellow River–Huangshui River Valley and provide a robust scientific basis for proactive disaster prevention, management of disaster chains, and mitigation of multi-hazard clusters in [...] Read more.
This study aims to identify the key factors contributing to the destabilization of the geohazard disaster-pregnant environment in the Yellow River–Huangshui River Valley and provide a robust scientific basis for proactive disaster prevention, management of disaster chains, and mitigation of multi-hazard clusters in unstable regions. The research focuses on the Yellow River–Huangshui River Valley, evaluating the stability of its geohazard disaster-pregnant environment. The disaster-pregnant environment is classified into static and dynamic categories. The static disaster-pregnant environment encompasses factors such as lithology, fracture density, topography, slope, river network density, and vegetation cover. The dynamic disaster-pregnant environment incorporates variables such as extreme rainfall, consecutive rainy days, annual rainfall averages, monthly high temperatures, monthly maximum temperature variations, average annual air temperatures, and human activities. A random forest model was employed to quantitatively assess the stability of the geohazard disaster-pregnant environment in the Yellow River–Huangshui River Valley. The findings indicated that (1) extreme indicators were the primary contributors to the destabilization of the disaster-pregnant environment, with very heavy rainfall contributing 28% and consecutive rainy days contributing 27%. Human activities ranked next, accounting for 15%. (2) Unstable regions for static, dynamic, and integrated disaster-pregnant environments accounted for 44%, 45%, and 44% of the study area, respectively, with all unstable areas concentrated in river valley regions. (3) The overall trend of stability in the disaster-pregnant environment was characterized by widespread instability. Extremely unstable areas were predominantly located in river valley regions, largely influenced by human activities. Conversely, only 0.1% of the region exhibited signs of stability, and 2.1% showed a tendency toward extreme stability. Full article
15 pages, 8746 KiB  
Article
Self-Assembly Strategy for Synthesis of WO3@TCN Heterojunction: Efficient for Photocatalytic Degradation and Hydrogen Production via Water Splitting
by Li Zhou, Wenjie Zhang, Zezhao Huang, Feng Hu, Peng Li and Xiaoquan Yao
Molecules 2025, 30(2), 379; https://doi.org/10.3390/molecules30020379 - 17 Jan 2025
Viewed by 201
Abstract
Herein, a WO3@TCN photocatalyst was successfully synthesized using a self-assembly method, which demonstrated effectiveness in degrading organic dyestuffs and photocatalytic evolution of H2. The synergistic effect between WO3 and TCN, along with the porous structure of TCN, facilitated [...] Read more.
Herein, a WO3@TCN photocatalyst was successfully synthesized using a self-assembly method, which demonstrated effectiveness in degrading organic dyestuffs and photocatalytic evolution of H2. The synergistic effect between WO3 and TCN, along with the porous structure of TCN, facilitated the formation of a heterojunction that promoted the absorption of visible light, accelerated the interfacial charge transfer, and inhibited the recombination of photogenerated electron–hole pairs. This led to excellent photocatalytic performance of 3%WO3@TCN in degrading TC and catalyzing H2 evolution from water splitting under visible-light irradiation. After modulation, the optimal 3%WO3@TCN exhibited a maximal degradation rate constant that was twofold higher than that of TCN alone and showed continuous H2 generation in the photocatalytic hydrogen evolution. Mechanistic studies revealed that •O2 constituted the major active species for the photocatalytic degradation of tetracycline. Experimental and DFT results verified the electronic transmission direction of WO3@TCN heterojunction. Overall, this study facilitates the structural design of green TCN-based heterojunction photocatalysts and expands the application of TCN in the diverse photocatalytic processes. Additionally, this study offers valuable insights into strategically employing acid regulation modulation to enhance the performance of carbon nitride-based photocatalysts by altering the topography of WO3@TCN composite material dramatically. Full article
Show Figures

Figure 1

Figure 1
<p>(<b>a</b>,<b>b</b>) SEM images of TCN and 3%WO<sub>3</sub>@TCN; (<b>c</b>,<b>d</b>) TEM and HRTEM images of TCN and 3%WO<sub>3</sub>@TCN; (<b>e</b>,<b>f</b>) EDS elemental mappings of 3%WO<sub>3</sub>@TCN material and the corresponding elemental mappings of C, N, O, and W.</p>
Full article ">Figure 2
<p>XPS spectra of (<b>a</b>) TCN and 3%WO<sub>3</sub>@TCN; (<b>b</b>) W 4f spectra of 3%WO<sub>3</sub>@TCN photocatalysts; (<b>c</b>) C1s and (<b>d</b>) N1s spectra of TCN and 3%WO<sub>3</sub>@TCN.</p>
Full article ">Figure 3
<p>(<b>a</b>) The FT-IR spectra of TCN and 3%WO<sub>3</sub>@TCN catalysts; (<b>b</b>) XRD patterns; (<b>c</b>) N<sub>2</sub> adsorption–desorption isotherms; and (<b>d</b>) pore size distribution of the WO<sub>3</sub>, TCN, and 3%WO<sub>3</sub>@TCN catalysts.</p>
Full article ">Figure 4
<p>(<b>a</b>) UV-Vis DRS spectra; (<b>b</b>) band gap energy; (<b>c</b>) VB-XPS spectra; and (<b>d</b>) schematic diagram of energy band levels of WO<sub>3</sub>, TCN, and 3%WO<sub>3</sub>@TCN.</p>
Full article ">Figure 5
<p>The work functions of (<b>a</b>) TCN; (<b>b</b>) WO<sub>3</sub>; and (<b>c</b>) WO<sub>3</sub>@TCN. The electron distributions of (<b>d</b>) TCN; (<b>e</b>) WO<sub>3</sub>; and (<b>f</b>) WO<sub>3</sub>@TCN.</p>
Full article ">Figure 6
<p>(<b>a</b>) Photocatalytic degradation of TC with different catalysts; (<b>b</b>) pseudo-first-order kinetic fitting curves; (<b>c</b>) effects of different free radical inhibitor on TC degradation. Reaction conditions: [catalysts]<sub>0</sub> = 0.2 mg·mL<sup>−1</sup>, [TC]<sub>0</sub> = 20 mg·L<sup>−1</sup>, T = 25 °C; (<b>d</b>) photocatalytic efficiency of 3%WO<sub>3</sub>@TCN toward TC removal in the presence of different scavengers.</p>
Full article ">Figure 7
<p>Schematic illustration of the proposed photocatalytic mechanism by 3%WO<sub>3</sub>@TCN under visible light irradiation.</p>
Full article ">Scheme 1
<p>Synthetic procedure of TCN and WO<sub>3</sub>@TCN.</p>
Full article ">
26 pages, 7227 KiB  
Article
Uncertainty-Based Scale Identification and Process–Topography Interaction Analysis via Bootstrap: Application to Grit Blasting
by François Berkmans, Julie Lemesle, Robin Guibert, Michal Wieczorowski, Christopher Brown and Maxence Bigerelle
Fractal Fract. 2025, 9(1), 48; https://doi.org/10.3390/fractalfract9010048 - 17 Jan 2025
Viewed by 276
Abstract
Finding the relevant scale to observe the influence of a process is one of the most important purposes of multiscale surface characterization. This study investigates various methods to determine a pertinent scale for evaluating the relationship between the relative area and grit blasting [...] Read more.
Finding the relevant scale to observe the influence of a process is one of the most important purposes of multiscale surface characterization. This study investigates various methods to determine a pertinent scale for evaluating the relationship between the relative area and grit blasting pressure. Several media types were tested alongside two different methods for calculating the relative area and three bootstrapping approaches for scale determination through regression. Comparison with the existing literature highlights innovations in roughness parameter characterization, particularly the advantages of relative area over traditional parameters like Sa. This study also discusses the relevance of different media types in influencing surface topography. Additionally, insights from a similar study on the multiscale Sdq parameter and blasting pressure correlation are integrated, emphasizing a scale relevance akin to our Sdr method’s 120 µm cut-off length. Overall, our findings suggest a pertinent scale of 10,000 µm2 for the Patchwork method and a 120 µm cut-off length for the Sdr method, derived from bootstrapping on residual regression across all media. At the relevant scale, every value of R2 inferior to 0.83 is not significant with the threshold of 5% for the two methods of calculation of the relative area. This study enhances the understanding of how media types and blasting pressures impact surface topography, offering insights for refining material processing and surface treatment strategies. Full article
(This article belongs to the Special Issue Fractal Analysis and Its Applications in Materials Science)
Show Figures

Figure 1

Figure 1
<p>Comparison of two methods, Sdr (ISO 25178-2) and Patchwork, for calculating relative areas of surface topographies created by blasting with glass beads. The points represent the median of the relative area values, categorized by calculation method and pressure. Blue symbols indicate the median points for the Patchwork method, while red symbols correspond to the Sdr method. The scale refers to the cut-off length of the low-pass Gaussian filter applied in the Sdr calculation. For the Patchwork method, the tile size in µm<sup>2</sup> is equal to half the square of the cut-off length.</p>
Full article ">Figure 2
<p>Surface topographies of TA6V surfaces grit-blasted at 2 bar (<b>a</b>), 4 bar (<b>b</b>), and 8 bar (<b>c</b>) with the C300 medium. The aggressiveness of the medium can make it difficult to assess visually the gradation in blasting intensity. More surface topographies are shown in <a href="#app1-fractalfract-09-00048" class="html-app">Appendix A</a>.</p>
Full article ">Figure 3
<p>Diagram of the two calculation methods used in this study, shown in terms of relative length. The blue continuous line represents a real surface. The green line, a linear interpolation between measured height points, represents our measured profile (the Sdr method calculates the relative length at the sampling scale). The red line illustrates the profile obtained by the Patchwork method.</p>
Full article ">Figure 4
<p>Results of the linear regressions of the relative area as a function of pressure for the two calculation methods. Simulations from 0 to 9 are obtained from bootstrapping replication of the real data and then averaged. The results come from measurements performed on surfaces blasted with the C 300 medium (corundum). Each simulation corresponds to an R<sup>2</sup> value, which is then averaged.</p>
Full article ">Figure 5
<p>Analysis of the R<sup>2</sup> distributions according to the scale of calculation for relative area under hypotheses H1 (<b>a</b>) and H0 (<b>b</b>) for the three bootstrapping methods: simple bootstrap (<b>i</b>), bootstrap based on pairs (<b>ii</b>), and bootstrap based on residuals (<b>iii</b>). The tile size of the Patchwork method (in µm<sup>2</sup>) is equal to half the square of the cut-off length of the Sdr method. Two plots are proposed for each bootstrapping method: the first one based on the media (<b>c</b>,<b>e</b>,<b>g</b>) and the second one based on the method of the relative area calculation, Sdr or Patchwork (<b>d</b>,<b>f</b>,<b>h</b>).</p>
Full article ">Figure 5 Cont.
<p>Analysis of the R<sup>2</sup> distributions according to the scale of calculation for relative area under hypotheses H1 (<b>a</b>) and H0 (<b>b</b>) for the three bootstrapping methods: simple bootstrap (<b>i</b>), bootstrap based on pairs (<b>ii</b>), and bootstrap based on residuals (<b>iii</b>). The tile size of the Patchwork method (in µm<sup>2</sup>) is equal to half the square of the cut-off length of the Sdr method. Two plots are proposed for each bootstrapping method: the first one based on the media (<b>c</b>,<b>e</b>,<b>g</b>) and the second one based on the method of the relative area calculation, Sdr or Patchwork (<b>d</b>,<b>f</b>,<b>h</b>).</p>
Full article ">Figure 6
<p>Surface topographies of TA6V samples grit-blasted at 2 bar (<b>a</b>), 4 bar (<b>b</b>), and 8 bar (<b>c</b>) with the C300 medium. The range of height varies significantly. The surfaces are the same as those presented in <a href="#fractalfract-09-00048-f002" class="html-fig">Figure 2</a> but this time filtered with a low-pass Gaussian filter at a 120 µm cut off (the relevance scale).</p>
Full article ">Figure 7
<p>Distributions of the R<sup>2</sup> values at all scales under H1 (<b>a</b>) and H0 (<b>b</b>) for every method of bootstrapping computation: simple bootstrap (<b>i</b>), paired bootstrap (<b>ii</b>), and bootstrap based on residuals (<b>iii</b>). The black lines on the H0 plots are the threshold value at 95% of the R<sup>2</sup> distribution: 0.59 (<b>bi</b>), 0.91 (<b>bii</b>), and 0.83 (<b>biii</b>).</p>
Full article ">Figure 8
<p>Evolution of the slope (<b>i</b>) and intercept (<b>ii</b>) as a function of scale for H1 (<b>a</b>) and H0 (<b>b</b>) using bootstrap based on residuals.</p>
Full article ">Figure 9
<p>Distribution of the R<sup>2</sup> values by medium at the relevant scale for the Patchwork (<b>i</b>) and Sdr (<b>ii</b>) methods and for H1 (<b>a</b>) and H0 (<b>b</b>). The digits after 250 indicate the blasting series (e.g., G 250-1 = first series of the G250 medium).</p>
Full article ">Figure 10
<p>Box plots of the relative area values by pressure at the relevance scale (tile size between 10,000 µm<sup>2</sup> and 14,000 µm<sup>2</sup> for the Patchwork method and cut-off length of 120 µm for the Sdr method). The results are presented by medium (<b>a</b>–<b>e</b>) and calculation method (<b>i</b>,<b>ii</b>).</p>
Full article ">Figure 10 Cont.
<p>Box plots of the relative area values by pressure at the relevance scale (tile size between 10,000 µm<sup>2</sup> and 14,000 µm<sup>2</sup> for the Patchwork method and cut-off length of 120 µm for the Sdr method). The results are presented by medium (<b>a</b>–<b>e</b>) and calculation method (<b>i</b>,<b>ii</b>).</p>
Full article ">Figure 11
<p>Bivariate density (intercept, slope) of the linear regression at the relevant scale between relative area for the three media of grit blasting and the two methods of relative area calculation (Patchwork, Sdr) obtained by bootstrap on residuals. The red frame is a zoom with ellipses of confidence at 95%.</p>
Full article ">Figure A1
<p>Surface topographies of blasted surface using the medium G 100 at (<b>a</b>) 2 bar of pressure, (<b>b</b>) 4 bar of pressure, and (<b>c</b>) 8 bar of pressure.</p>
Full article ">Figure A2
<p>Surface topographies of blasted surface using the medium G 250 at (<b>a</b>) 2 bar of pressure, (<b>b</b>) 4 bar of pressure, and (<b>c</b>) 8 bar of pressure.</p>
Full article ">Figure A3
<p>Surface topographies of blasted surface using the medium C 300 at (<b>a</b>) 2 bar of pressure, (<b>b</b>) 4 bar of pressure, and (<b>c</b>) 8 bar of pressure.</p>
Full article ">
23 pages, 12001 KiB  
Article
Enhancing Off-Road Topography Estimation by Fusing LIDAR and Stereo Camera Data with Interpolated Ground Plane
by Gustav Sten, Lei Feng and Björn Möller
Sensors 2025, 25(2), 509; https://doi.org/10.3390/s25020509 - 16 Jan 2025
Viewed by 254
Abstract
Topography estimation is essential for autonomous off-road navigation. Common methods rely on point cloud data from, e.g., Light Detection and Ranging sensors (LIDARs) and stereo cameras. Stereo cameras produce dense point clouds with larger coverage but lower accuracy. LIDARs, on the other hand, [...] Read more.
Topography estimation is essential for autonomous off-road navigation. Common methods rely on point cloud data from, e.g., Light Detection and Ranging sensors (LIDARs) and stereo cameras. Stereo cameras produce dense point clouds with larger coverage but lower accuracy. LIDARs, on the other hand, have higher accuracy and longer range but much less coverage. LIDARs are also more expensive. The research question examines whether incorporating LIDARs can significantly improve stereo camera accuracy. Current sensor fusion methods use LIDARs’ raw measurements directly; thus, the improvement in estimation accuracy is limited to only LIDAR-scanned locations The main contribution of our new method is to construct a reference ground plane through the interpolation of LIDAR data so that the interpolated maps have similar coverage as the stereo camera’s point cloud. The interpolated maps are fused with the stereo camera point cloud via Kalman filters to improve a larger section of the topography map. The method is tested in three environments: controlled indoor, semi-controlled outdoor, and unstructured terrain. Compared to the existing method without LIDAR interpolation, the proposed approach reduces average error by 40% in the controlled environment and 67% in the semi-controlled environment, while maintaining large coverage. The unstructured environment evaluation confirms its corrective impact. Full article
Show Figures

Figure 1

Figure 1
<p>Beam distribution dependant on distance from sensor.</p>
Full article ">Figure 2
<p>Sensor and software setup. (<b>a</b>) Visualizing of how the LIDAR and stereo camera were mounted; (<b>b</b>) software setup for recording data.</p>
Full article ">Figure 3
<p>Process of mapping point clouds to elevation map. (<b>a</b>) Single point cloud mapping to elevation map; (<b>b</b>) multiple point clouds mapping to elevation map.</p>
Full article ">Figure 4
<p>Example of point clouds from stereo camera (<b>a</b>), LIDAR (<b>b</b>), and the actual ground truth at the center of the point cloud (<b>c</b>). Note that (<b>c</b>) is zoomed in.</p>
Full article ">Figure 5
<p>Interpolation methodology. (<b>a</b>) Description of the interpolation direction in the gridmap. (<b>b</b>) Trapezoid function for variance between two measured points, <math display="inline"><semantics> <msub> <mi>p</mi> <mn>1</mn> </msub> </semantics></math> and <math display="inline"><semantics> <msub> <mi>p</mi> <mn>2</mn> </msub> </semantics></math>.</p>
Full article ">Figure 6
<p>Interpolated map and its corresponding variance. x and y are grid cell indexes.</p>
Full article ">Figure 7
<p>Single sensor elevation map.</p>
Full article ">Figure 8
<p>Estimation Errors of the Two Sensors.</p>
Full article ">Figure 9
<p>Fused elevation maps.</p>
Full article ">Figure 10
<p>Estimation errors of the two fusion methods.</p>
Full article ">Figure 11
<p>Photograph of the test area, with critical measurement points marked.</p>
Full article ">Figure 12
<p>Example of raw point clouds with the objects highlighted.</p>
Full article ">Figure 13
<p>Stereo camera and LIDAR maps with their resulting variance.</p>
Full article ">Figure 14
<p>Fused maps with their resulting variance.</p>
Full article ">Figure 15
<p>Photograph of the test area.</p>
Full article ">Figure 16
<p>Example of raw point clouds.</p>
Full article ">Figure 17
<p>Stereo camera and LIDAR maps with their resulting variance.</p>
Full article ">Figure 18
<p>Fused maps with their resulting variance.</p>
Full article ">Figure 19
<p>Estimation of both fusion methods along Y = 34.</p>
Full article ">
15 pages, 2514 KiB  
Article
Impact of Spatial Configuration of Bioretention Cells on Catchment Hydrological Performance Under Extreme Rainfall Conditions with Different Stormwater Flow Paths
by Xu Liu, Jun Huang, Sicheng Zheng, Li Wang, Yimin Huang and Zebin Yu
Water 2025, 17(2), 233; https://doi.org/10.3390/w17020233 - 16 Jan 2025
Viewed by 276
Abstract
Bioretention cells (BCs) are widely used to manage urban runoff due to their positive impact on runoff control. Current research primarily focuses on optimizing the internal structural design of bioretention cells, while studies on the interactions between their spatial configuration, topography, and land [...] Read more.
Bioretention cells (BCs) are widely used to manage urban runoff due to their positive impact on runoff control. Current research primarily focuses on optimizing the internal structural design of bioretention cells, while studies on the interactions between their spatial configuration, topography, and land use types are limited. This study employs the Storm Water Management Model (SWMM) and uses extreme rainfall to analyze the influence of typical stormwater flow paths, determined by various land use types and topography, as well as the spatial configurations of bioretention cells on catchment hydrological performance. The results show the following: (1) Different stormwater flow paths significantly affect catchment hydrological performance, with series-type pathways performing best. (2) The spatial configuration of bioretention cells significantly influences catchment hydrological performance. Decentralized BCs under series-type pathways showed better performance for reducing total outflow and peak runoff, with reduction rates increasing by 7.1% and 8.8%, while centralized BCs better delayed peak times. (3) Stormwater flow paths affect BC efficiency in catchment hydrological performance. Decentralized BCs under a series-type stormwater flow path are recommended for priority use. This study provides a novel perspective for optimizing the spatial arrangement of BCs and urban stormwater management, thereby contributing to flood risk mitigation. Full article
(This article belongs to the Special Issue Watershed Hydrology and Management under Changing Climate)
Show Figures

Figure 1

Figure 1
<p>(<b>a</b>) The location of the study area. (<b>b</b>) The elevation of the study area. (<b>c</b>) The land use types in the study area.</p>
Full article ">Figure 2
<p>Schematic diagram of three stormwater flow paths: ① series flow path; ② parallel flow path; ③ hybrid flow path.</p>
Full article ">Figure 3
<p>Schematic diagram of two spatial distributions. (<b>a</b>) Centralized BC treatment. (<b>b</b>) Decentralized BC treatment.</p>
Full article ">Figure 4
<p>Hydrological differences in catchment areas under three types of stormwater flow paths. The dots on the right side of the boxplot represent individual rainfall events. The curve illustrates the normal distribution of these points, reflecting the central tendency and dispersion of the data.</p>
Full article ">Figure 5
<p>The impact of centralized BCs on catchment hydrology under three types of flow paths. The dots on the right side of the boxplot represent individual rainfall events. The curve illustrates the normal distribution of these points, reflecting the central tendency and dispersion of the data.</p>
Full article ">Figure 6
<p>The impact of decentralized BCs on catchment hydrology under three types of flow paths. The dots on the right side of the boxplot represent individual rainfall events. The curve illustrates the normal distribution of these points, reflecting the central tendency and dispersion of the data.</p>
Full article ">Figure 7
<p>The comparison of the impact of centralized and decentralized BCs on catchment hydrology under three types of flow paths.</p>
Full article ">
11 pages, 4078 KiB  
Article
Biomimetic Silicone Surfaces for Antibacterial Applications
by Marie Barshutina, Dmitry Yakubovsky, Aleksey Arsenin, Valentyn Volkov, Sergey Barshutin, Anastasiya Vladimirova and Andrei Baymiev
Polymers 2025, 17(2), 213; https://doi.org/10.3390/polym17020213 - 16 Jan 2025
Viewed by 219
Abstract
Biomimetic patterning emerges as a promising antibiotic-free approach to protect medical devices from bacterial adhesion and biofilm formation. The main advantage of this approach lies in its simplicity and scalability for industrial applications. In this study, we employ it to produce antibacterial coatings [...] Read more.
Biomimetic patterning emerges as a promising antibiotic-free approach to protect medical devices from bacterial adhesion and biofilm formation. The main advantage of this approach lies in its simplicity and scalability for industrial applications. In this study, we employ it to produce antibacterial coatings based on silicone materials, widely used in the healthcare industry. In doing so, we patterned silicone substrates with a topography of various flower petals (rose, chamomile, pansy, and magnolia) and studied the relationship between the antibacterial properties of the obtained biomimetic substrates and their surface topography. To study the surface topography of biomimetic surfaces, we used the fractal analysis of their SEM images. In particular, as a measure of surface complexity and heterogeneity, we used the values of the developed interfacial area ratio (Sdr) and lacunarity coefficient (β). In the result, we found that the bacterial area coverage of biomimetic substrates decreased exponentially with the increase in their surface complexity and heterogeneity, and prominent antibacterial properties were observed at β > 1.6 and Sdr > 50. The results of this study can be used to identify biomimetic materials with superior antibacterial properties and produce efficient antibacterial silicone coatings for biomedical and healthcare applications. Full article
(This article belongs to the Special Issue The Application of Polymers in Biomimetics)
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Schematic of the fabrication process to obtain silicone replicas of flower petals.</p>
Full article ">Figure 2
<p>SEM images illustrating the silicone replicas of rose, chamomile, pansy, and magnolia petals: (<b>a</b>,<b>d</b>,<b>g</b>,<b>j</b>) Photograph of a flower; (<b>b</b>,<b>e</b>,<b>h</b>,<b>k</b>) silicone replica of a flower petal with a 100 µm scale bar; (<b>c</b>,<b>f</b>,<b>i</b>,<b>l</b>) silicone replica of a flower petal with a 10 µm scale bar.</p>
Full article ">Figure 3
<p>Fractal analysis of SEM image topography of flower petal replicas: (<b>a</b>) lacunarity analysis; (<b>b</b>) roughness ratio analysis. The relative standard deviation (RSD) of all measurements was less than 5%.</p>
Full article ">Figure 4
<p>The fluorescent images of biomimetic and flat substrates populated with bacteria (at 20× magnification): (<b>a</b>) flat; (<b>b</b>) magnolia replica; (<b>c</b>) pansy replica; (<b>d</b>) chamomile replica; (<b>e</b>) rose replica; (<b>f</b>) quantitative analysis.</p>
Full article ">Figure 5
<p>Graphical representation of the E. coli area coverage of biomimetic surfaces as a function of: (<b>a</b>) lacunarity coefficient; (<b>b</b>) developed interfacial area ratio. The stars on both plots indicate experimental data, while the line indicates the fitted data using second order polynomials.</p>
Full article ">
24 pages, 2596 KiB  
Article
Spatiotemporal Evolution and Driving Mechanisms of kNDVI in Different Sections of the Yangtze River Basin Using Multiple Statistical Methods and the PLSPM Model
by Zhenjiang Wu, Fengmei Yao, Adeel Ahmad, Fan Deng and Jun Fang
Remote Sens. 2025, 17(2), 299; https://doi.org/10.3390/rs17020299 - 16 Jan 2025
Viewed by 375
Abstract
Spatiotemporal vegetation changes serve as a key indicator of regional ecological environmental quality and provide crucial guidance for developing strategies for regional ecological protection and sustainable development. Currently, vegetation change studies in the Yangtze River Basin primarily rely on the Normalized Difference Vegetation [...] Read more.
Spatiotemporal vegetation changes serve as a key indicator of regional ecological environmental quality and provide crucial guidance for developing strategies for regional ecological protection and sustainable development. Currently, vegetation change studies in the Yangtze River Basin primarily rely on the Normalized Difference Vegetation Index (NDVI). However, the NDVI is susceptible to atmospheric and soil conditions and exhibits saturation phenomena in areas with high vegetation coverage. In contrast, the kernel NDVI (kNDVI) demonstrates significant advantages in suppressing background noise and improving saturation thresholds through nonlinear kernel transformation, thereby enhancing sensitivity to vegetation changes. To elucidate the spatiotemporal characteristics and driving mechanisms of vegetation changes in the Yangtze River Basin, this study constructed a temporal kNDVI using MOD09GA data from 2000 to 2022. Considering sectional heterogeneity, rather than analyzing the entire region as a whole as in previous studies, this research examined spatiotemporal evolution characteristics by sections using four statistical metrics. Subsequently, Partial Least Squares Path Modeling (PLSPM) was innovatively introduced to quantitatively analyze the influence mechanisms of topographic, climatic, pedological, and socioeconomic factors. Compared to traditional correlation analysis and the geographical detector method, PLSPM, as a theoretically driven statistical method, can simultaneously process path relationships among multiple latent variables, effectively revealing the intensity and pathways of driving factors’ influences, while providing more credible and interpretable explanations for kNDVI variation mechanisms. Results indicate that the overall kNDVI in the Yangtze River Basin exhibited an upward trend, with the midstream demonstrating the most significant improvement with minimal interannual fluctuations, the upstream displaying an east-increasing and west-stable spatial pattern, and the downstream demonstrating coexisting improvement and degradation characteristics, with these trends expected to persist. Driving mechanism analysis reveals that the upstream was predominantly influenced by the climatic factor, the midstream was dominated by terrain, and the downstream displayed terrain–soil coupling effects. Based on these findings, it is recommended that the upstream focus on enhancing vegetation adaptation management to climate change, the midstream need to coordinate the relationship between topography and human activities, and the downstream should concentrate on controlling the negative impacts of urban expansion on vegetation. Full article
Show Figures

Figure 1

Figure 1
<p>Sectional map (<b>a</b>), land cover map (<b>b</b>), and elevation map (<b>c</b>) of the Yangtze River Basin.</p>
Full article ">Figure 2
<p>Line charts of kNDVI changes from 2000 to 2022 in the entire Yangtze River Basin (<b>a</b>), upstream (<b>b</b>), midstream (<b>c</b>), and downstream (<b>d</b>).</p>
Full article ">Figure 3
<p>Spatial distribution of the kNDVI in the Yangtze River Basin in (<b>a</b>) 2000 and (<b>b</b>) 2022, and (<b>c</b>) the percentage distribution of different kNDVI categories from 2000 to 2022.</p>
Full article ">Figure 4
<p>Temporal trends of kNDVI in the Yangtze River Basin from 2000 to 2022.</p>
Full article ">Figure 5
<p>Stability of kNDVI variations in the Yangtze River Basin from 2000 to 2022.</p>
Full article ">Figure 6
<p>Future trends in the kNDVI for the Yangtze River Basin.</p>
Full article ">Figure 7
<p>Framework of the structural model.</p>
Full article ">Figure 8
<p>Path coefficients diagram of the structural model in the upstream, midstream, and downstream of the Yangtze River.</p>
Full article ">
18 pages, 6529 KiB  
Article
A Novel Algorithm for Estimating the Sand Dune Density of the Taklimakan Desert Based on Remote Sensing Data
by Mingyu Wang, Yongqiang Liu, Huoqing Li, Minzhong Wang, Wen Huo and Zonghui Liu
Remote Sens. 2025, 17(2), 297; https://doi.org/10.3390/rs17020297 - 16 Jan 2025
Viewed by 340
Abstract
The dune density is an important parameter for representing the characteristics of desert geomorphology, providing a precise depiction of the undulating topography of the desert. Owing to the limitations of estimation methods and data availability, accurately quantifying dune density has posed a significant [...] Read more.
The dune density is an important parameter for representing the characteristics of desert geomorphology, providing a precise depiction of the undulating topography of the desert. Owing to the limitations of estimation methods and data availability, accurately quantifying dune density has posed a significant challenge; in response to this issue, we propose an innovative model to estimate dune density using a dune vertex search combined with four-directional orographic spectral decomposition. This study reveals several key insights: (1) Taklimakan Desert distributes approximately 5.31 × 107 dunes, with a linear regression fit R2 of 0.79 between the estimated and observed values. The average absolute error and root mean square error are calculated as 25.61 n/km2 and 30.48 n/km2, respectively. (2) The distribution of dune density across the eastern, northeastern, southern, and western parts of the Taklimakan Desert is relatively lower, while there is higher dune density in the central and northern areas. (3) The observation data constructed using the improved YOLOv8s algorithm and remote sensing imagery effectively validate the estimation results of dune density. The new algorithm demonstrates a high level of accuracy in estimating sand dune density, thereby providing crucial parameters for sub-grid orographic parameterization in desert regions. Additionally, its application potential in dust modeling appears promising. Full article
Show Figures

Figure 1

Figure 1
<p>Research area overview, (<b>a</b>) research area location, (<b>b</b>) types of dunes and locations of sampling areas, and (<b>c1</b>–<b>c13</b>) sample area Google Earth image. The desert map is provided by the National Cryosphere Desert Data Center. (<a href="http://www.ncdc.ac.cn" target="_blank">http://www.ncdc.ac.cn</a> (accessed on 20 March 2024)).</p>
Full article ">Figure 2
<p>Dunes in the TD. (<b>a</b>,<b>b</b>) Dunes in the Hade (83°42′E, 40°45′N), date: 04/2024. (<b>c</b>–<b>e</b>) Dunes in the Xiaotang (84°18′E, 40°49′N), date: 05/2024. (<b>f</b>,<b>g</b>) Dunes on both sides of the desert highway, (<b>f</b>): (84°19′E, 40°34′N), date: 10/2023 (<b>g</b>): (83°44′E, 39°15′N), date: 05/2024. (<b>h</b>–<b>j</b>) Dunes in the Tazhong (83°38′E, 38°59′N), date: 05/2024. (<b>k</b>,<b>l</b>) Yutian oasis–desert ecotone (81°28′E, 36°56′N), date: 05/2023.</p>
Full article ">Figure 3
<p>The workflow of this study.</p>
Full article ">Figure 4
<p>Extraction algorithm based on dune vertices. (<b>a</b>) DEM. (<b>b</b>) Distribution of dune vertices and areas.</p>
Full article ">Figure 5
<p>The orographic spectral decomposition of the four-directional method. (<b>a</b>) 0°, (<b>b</b>) 45°, (<b>c</b>) 90°, and (<b>d</b>) 135°. Different colors represent different terrain spectral lines.</p>
Full article ">Figure 6
<p>The framework for the acquisition of dune density observation data.</p>
Full article ">Figure 7
<p>Identification results of some sample areas. (<b>a</b>–<b>f</b>) Six randomly selected sample areas; the part circled by the orange yellow box represents the sand dunes recognized by the model.</p>
Full article ">Figure 8
<p>Dune density distribution in TD. (<b>a</b>) Dune density distribution; (<b>b</b>–<b>i</b>) Google Earth images of typical sample regions. (The display of (<b>b</b>–<b>i</b>) is not the actual size of the sample area but only a part of it).</p>
Full article ">Figure 9
<p>Accuracy of the typical sample area.</p>
Full article ">Figure 10
<p>Linear fitting results between estimated and observed values of dune density.</p>
Full article ">
29 pages, 12669 KiB  
Article
Integrated Machine Learning Approaches for Landslide Susceptibility Mapping Along the Pakistan–China Karakoram Highway
by Mohib Ullah, Haijun Qiu, Wenchao Huangfu, Dongdong Yang, Yingdong Wei and Bingzhe Tang
Land 2025, 14(1), 172; https://doi.org/10.3390/land14010172 - 15 Jan 2025
Viewed by 277
Abstract
The effectiveness of data-driven landslide susceptibility mapping relies on data integrity and advanced geospatial analysis; however, selecting the most suitable method and identifying key regional factors remains a challenging task. To address this, this study assessed the performance of six machine learning models, [...] Read more.
The effectiveness of data-driven landslide susceptibility mapping relies on data integrity and advanced geospatial analysis; however, selecting the most suitable method and identifying key regional factors remains a challenging task. To address this, this study assessed the performance of six machine learning models, including Convolutional Neural Networks (CNNs), Random Forest (RF), Categorical Boosting (CatBoost), their CNN-based hybrid models (CNN+RF and CNN+CatBoost), and a Stacking Ensemble (SE) combining CNN, RF, and CatBoost in mapping landslide susceptibility along the Karakoram Highway in northern Pakistan. Twelve geospatial factors were examined, categorized into Topography/Geomorphology, Land Cover/Vegetation, Geology, Hydrology, and Anthropogenic Influence. A detailed landslide inventory of 272 occurrences was compiled to train the models. The proposed stacking ensemble and hybrid models improve landslide susceptibility modeling, with the stacking ensemble achieving an AUC of 0.91. Hybrid modeling enhances accuracy, with CNN–RF boosting RF’s AUC from 0.85 to 0.89 and CNN–CatBoost increasing CatBoost’s AUC from 0.87 to 0.90. Chi-square (χ2) values (9.8–21.2) and p-values (<0.005) confirm statistical significance across models. This study identifies approximately 20.70% of the area as from high to very high risk, with the SE model excelling in detecting high-risk zones. Key factors influencing landslide susceptibility showed slight variations across the models, while multicollinearity among variables remained minimal. The proposed modeling approach reduces uncertainties, enhances prediction accuracy, and supports decision-makers in implementing effective landslide mitigation strategies. Full article
Show Figures

Figure 1

Figure 1
<p>Geospatial analysis of landslide risks in northern Pakistan featuring (<b>a</b>) a regional map locating the study area in Asia and (<b>b</b>) a detailed topographic map of the Karakorum Highway, landslide locations, significant earthquakes, and key settlements. Earthquake data for the study area (1970–2015) were obtained from the China Earthquake Networks Center, with data processed and clipped to the Pakistan–China Economic Corridor region by Northwest University’s College of Urban and Environmental Sciences under the National International Science and Technology Cooperation Project (Grant No. 2018YFE0100100).</p>
Full article ">Figure 2
<p>Lithological map of the study region.</p>
Full article ">Figure 3
<p>Diagram outlining the steps involved in creating a landslide susceptibility map.</p>
Full article ">Figure 4
<p>A comprehensive overview of the landslide field investigations conducted within the study area. Subfigure (<b>a</b>) shows the entire study area, with landslide polygons highlighted in red against a grayscale elevation map. A red square on this map identifies the specific area examined in greater detail. Subfigure (<b>b</b>) provides an elevation map of the region within the red square from (<b>a</b>), featuring landslide polygons also marked in red. Subfigures (<b>c</b>–<b>f</b>) delineate the boundaries of detected landslides with yellow lines in various regions: (<b>c</b>) Gilgit area, (<b>d</b>) Chilas area, (<b>e</b>) Babusar area, and (<b>f</b>) Passu area. Photos were taken during field surveys conducted on 20 July 2024, with photo credits attributed to our research team.</p>
Full article ">Figure 5
<p>Maps depicting explanatory variables in the study region: (<b>a</b>) elevation, (<b>b</b>) aspect, (<b>c</b>) curvature, (<b>d</b>) NDVI, (<b>e</b>) TWI, (<b>f</b>) slope, (<b>g</b>) rainfall, (<b>h</b>) landcover, (<b>i</b>) proximity to roads, (<b>j</b>) proximity to streams, (<b>k</b>) proximity to faults, (<b>l</b>) lithology. Additional notes for (<b>l</b>): Jms—Jurassic metamorphic and sedimentary rocks, Ks—Cretaceous sedimentary rocks, Mi—Middle Jurassic rocks, pC—Undivided Precambrian rocks, Pz—Undifferentiated Paleozoic rocks, Pzl—Lower Paleozoic rocks, PzpC—Permian—Precambrian rocks, Ti—Tertiary igneous rocks, TrCs—Carboniferous-Triassic sedimentary rocks, Trms—Triassic metamorphic and sedimentary rocks.</p>
Full article ">Figure 6
<p>(<b>a</b>) Architecture of the CNN model. (<b>b</b>) Overview of the hybrid modeling workflow.</p>
Full article ">Figure 7
<p>The loss and error curves for three models: (<b>a</b>) the CatBoost training and validation loss decreasing as the number of iterations increases, (<b>b</b>) the CNN training and validation loss decreasing as the number of epochs increases, and (<b>c</b>) the Random Forest Out-of-Bag (OOB) error and validation error reducing as the number of trees increases.</p>
Full article ">Figure 8
<p>Compares IGR and q values across selected conditioning factors. Factors such as roads, streams, and faults are noted for their proximity effects on these values.</p>
Full article ">Figure 9
<p>Landslide susceptibility maps derived using various modeling approaches: (<b>a</b>) Random Forest (RF), (<b>b</b>) Categorical Boosting (CatBoost), (<b>c</b>) Convolutional Neural Network (CNN), (<b>d</b>) Convolutional Neural Network–Random Forest (CNN–RF), (<b>e</b>) Convolutional Neural Network–Categorical Boosting (CNN–CatBoost), (<b>f</b>) Stacking Ensemble (SE).</p>
Full article ">Figure 10
<p>ROC (Receiver Operating Characteristic) curves for different machine learning models assessing landslide susceptibility.</p>
Full article ">Figure 11
<p>Bar charts displaying the importance of various features across different models for landslide susceptibility in the study area.</p>
Full article ">Figure 12
<p>The relationships between various factors and landslide occurrences using the Weight of Evidence (WoE) model. Each subplot is labeled according to the specific factor being analyzed, including (<b>a</b>) aspect, (<b>b</b>) curvature, (<b>c</b>) elevation, (<b>d</b>) proximity to faults, (<b>e</b>) land cover, (<b>f</b>) lithology, (<b>g</b>) NDVI, (<b>h</b>) rainfall, (<b>i</b>) proximity to roads, (<b>j</b>) slope, (<b>k</b>) proximity to streams, and (<b>l</b>) TWI (Topographic Wetness Index). The land cover classes are denoted by their abbreviations: Barren Land (BL), Forest Grassland (FG), Dry Farmland (DF), Grassland/Moss (G/M), Cultivated Land (CL), Permanent Snow (PS), Waterbody (WB), and Woodland (W). Additional Notes for (<b>f</b>): Jms—Jurassic metamorphic and sedimentary rocks, Ks—Cretaceous sedimentary rocks, Mi—Middle Jurassic rocks, pC—Undivided Precambrian rocks, Pz—Undifferentiated Paleozoic rocks, Pzl—Lower Paleozoic rocks, PzpC—Permian—Precambrian rocks, Ti—Tertiary igneous rocks, TrCs—Upper Carboniferous—Lower Triassic sedimentary rocks, Trms—Triassic metamorphic and sedimentary rocks.</p>
Full article ">
30 pages, 17518 KiB  
Article
Preserving History: Assessments and Climate Adaptations at the House of the Seven Gables in Salem, Massachusetts, USA
by Paul Wright, Susan Baker and Stephen S. Young
Atmosphere 2025, 16(1), 84; https://doi.org/10.3390/atmos16010084 - 15 Jan 2025
Viewed by 347
Abstract
Salem, Massachusetts, is one of the oldest cities in the United States (1629) and its coastal location on the Atlantic helped create one of the wealthiest cities in America during the late 18th century, but today its coastal location threatens many of its [...] Read more.
Salem, Massachusetts, is one of the oldest cities in the United States (1629) and its coastal location on the Atlantic helped create one of the wealthiest cities in America during the late 18th century, but today its coastal location threatens many of its buildings due to sea level rise and increased storm activity. The House of the Seven Gables, a National Historic Landmark District, consists of five important historic buildings, the most famous being The Turner Ingersoll Mansion (1668), more commonly known as The House of the Seven Gables. Considered one of the most important houses in America, it is also one of the most threatened historic buildings due to its location on Salem’s harbor. The House of the Seven Gables conducted a two-year study funded by Massachusetts Coastal Zone Management to evaluate the risks posed by climate change. This process included the use of data from groundwater monitoring wells and a tidal gauge installed on-site, along with soil samples and a detailed survey base plan including topography and subsurface infrastructure. The project team then used the Massachusetts Coastal Flood Risk Model (MC-FRM) to assess climate change impacts on the site in 2030, 2050, and 2070, and then created a plan for adaptations that should be implemented before those risks materialize. Strategies for adapting to storm surges, increasing groundwater, and intense surface water runoff were evaluated for their effectiveness and appropriateness for the historic site. The conclusion of the study resulted in a five-phase plan ending in the managed retreat of the historic buildings to higher ground on the existing site. This article goes beyond other research that suggests coastal retreats by demonstrating how to quantitatively evaluate current and future coastal issues with predictive models and how to set viable dates for adaptive solutions and a managed retreat. Full article
(This article belongs to the Special Issue Climate Change Challenges for Heritage Architecture)
Show Figures

Figure 1

Figure 1
<p>Location map of Salem and the House of the Seven Gables.</p>
Full article ">Figure 2
<p>Aerial view of Gables campus, photo by Hugh Hou, 2023 [<a href="#B20-atmosphere-16-00084" class="html-bibr">20</a>].</p>
Full article ">Figure 3
<p>Filled tidelands, tidal shoreline, and buffer areas [<a href="#B20-atmosphere-16-00084" class="html-bibr">20</a>]. Drawing by the Horsley Witten Group. Note: this figure has a scale bar to show the size of the Gables campus, but the following campus figures do not have scale bars.</p>
Full article ">Figure 4
<p>Coastal flood plain [<a href="#B20-atmosphere-16-00084" class="html-bibr">20</a>]. Drawing by the Horsley Witten Group.</p>
Full article ">Figure 5
<p>Land cover [<a href="#B20-atmosphere-16-00084" class="html-bibr">20</a>]. Drawing by the Horsley Witten Group.</p>
Full article ">Figure 6
<p>Topography and drainage patterns [<a href="#B20-atmosphere-16-00084" class="html-bibr">20</a>]. Drawing by the Horsley Witten Group.</p>
Full article ">Figure 7
<p>Building elevations [<a href="#B20-atmosphere-16-00084" class="html-bibr">20</a>]. Drawing by the Horsley Witten Group.</p>
Full article ">Figure 8
<p>Coastal flood risk models [<a href="#B20-atmosphere-16-00084" class="html-bibr">20</a>]. Drawing by the Horsley Witten Group.</p>
Full article ">Figure 9
<p>Estimated groundwater elevations [<a href="#B20-atmosphere-16-00084" class="html-bibr">20</a>]. Drawing by the Horsley Witten Group.</p>
Full article ">Figure 10
<p>Seawall and revetment [<a href="#B20-atmosphere-16-00084" class="html-bibr">20</a>]. Drawing by the Horsley Witten Group.</p>
Full article ">Figure 11
<p>Adopted from the National Park Service flood adaptation guidelines [<a href="#B30-atmosphere-16-00084" class="html-bibr">30</a>].</p>
Full article ">Figure 12
<p>Proposed stormwater improvements in phase I [<a href="#B20-atmosphere-16-00084" class="html-bibr">20</a>]. Phasing details are described in the master planning and phasing section. Drawing by the Horsley Witten Group.</p>
Full article ">Figure 13
<p>Proposed stormwater improvements in phase V [<a href="#B20-atmosphere-16-00084" class="html-bibr">20</a>]. Previous improvements from prior phases are shown grayed out but are still meant to be used. Drawing by the Horsley Witten Group.</p>
Full article ">Figure 14
<p>Building subsurface drainage systems [<a href="#B20-atmosphere-16-00084" class="html-bibr">20</a>].</p>
Full article ">Figure 15
<p>Dry floodproofing [<a href="#B20-atmosphere-16-00084" class="html-bibr">20</a>].</p>
Full article ">Figure 16
<p>Wet floodproofing [<a href="#B20-atmosphere-16-00084" class="html-bibr">20</a>].</p>
Full article ">Figure 17
<p>Elevating the structure [<a href="#B20-atmosphere-16-00084" class="html-bibr">20</a>].</p>
Full article ">Figure 18
<p>Building relocation [<a href="#B20-atmosphere-16-00084" class="html-bibr">20</a>].</p>
Full article ">Figure 19
<p>Strategies evaluation chart developed by the climate adaptation team [<a href="#B20-atmosphere-16-00084" class="html-bibr">20</a>].</p>
Full article ">Figure 20
<p>Timeline triggers map [<a href="#B20-atmosphere-16-00084" class="html-bibr">20</a>].</p>
Full article ">Figure 21
<p>Current House of the Seven Gables campus map [<a href="#B20-atmosphere-16-00084" class="html-bibr">20</a>]. Drawing by the Horsley Witten Group.</p>
Full article ">Figure 22
<p>Triggers and actions for phase I [<a href="#B20-atmosphere-16-00084" class="html-bibr">20</a>].</p>
Full article ">Figure 23
<p>Phase I map for the House of the Seven Gables campus [<a href="#B20-atmosphere-16-00084" class="html-bibr">20</a>]. Numbers refer to locations where various actions are to take place. Drawing by the Horsley Witten Group.</p>
Full article ">Figure 24
<p>Triggers and actions for phase II [<a href="#B20-atmosphere-16-00084" class="html-bibr">20</a>].</p>
Full article ">Figure 25
<p>Phase II map for the House of the Seven Gables campus [<a href="#B20-atmosphere-16-00084" class="html-bibr">20</a>]. Numbers refer to locations where various actions are to take place. Drawing by the Horsley Witten Group.</p>
Full article ">Figure 26
<p>Triggers and actions for phase III [<a href="#B20-atmosphere-16-00084" class="html-bibr">20</a>].</p>
Full article ">Figure 27
<p>Phase III map for the House of the Seven Gables campus [<a href="#B20-atmosphere-16-00084" class="html-bibr">20</a>]. Numbers refer to locations where various actions are to take place. Drawing by the Horsley Witten Group.</p>
Full article ">Figure 28
<p>Triggers and actions for phase IV [<a href="#B20-atmosphere-16-00084" class="html-bibr">20</a>].</p>
Full article ">Figure 29
<p>Phase IV map for the House of the Seven Gables campus [<a href="#B20-atmosphere-16-00084" class="html-bibr">20</a>]. Numbers refer to locations where various actions are to take place. Drawing by the Horsley Witten Group.</p>
Full article ">Figure 30
<p>Triggers and actions for phase V [<a href="#B20-atmosphere-16-00084" class="html-bibr">20</a>].</p>
Full article ">Figure 31
<p>Phase V map for the House of the Seven Gables campus [<a href="#B20-atmosphere-16-00084" class="html-bibr">20</a>]. Numbers refer to locations where various actions are to take place. Drawing by the Horsley Witten Group.</p>
Full article ">
Back to TopTop