[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

Article Types

Countries / Regions

Search Results (252)

Search Parameters:
Keywords = real topography

Order results
Result details
Results per page
Select all
Export citation of selected articles as:
18 pages, 5416 KiB  
Article
Study of the Impact of Surface Topography on Wear Resistance
by Ben Wang, Wei Zhang and Zhongxun Liu
Coatings 2024, 14(9), 1128; https://doi.org/10.3390/coatings14091128 - 2 Sep 2024
Viewed by 435
Abstract
The surface texture parameter is a real reflection of the surface topography, which is closely related to the tribological properties of the surface, and the study of the correlation between the surface texture parameter and wear resistance is of great significance in revealing [...] Read more.
The surface texture parameter is a real reflection of the surface topography, which is closely related to the tribological properties of the surface, and the study of the correlation between the surface texture parameter and wear resistance is of great significance in revealing the tribological influence mechanism of the surface and realising the functional manufacturing of the surface. This paper takes the ball-end milling surface as the research object, establishes the three-dimensional simulation model of the surface topography, and analyses the surface topography and the surface texture parameter change rule. Based on the improved correlation analysis model, the correlation characteristics between the surface texture parameters, and between the surface texture parameters and the relative wear rate of per unit sliding distance KV, were investigated, and the prediction model of KV was established based on the surface texture parameters Sku, Sa, Sxp, Sp, and Ssk, and the correctness of the model was verified by experiments. The study in this paper provides a new idea to further reveal the relationship between surface topographical features and wear resistance and to guide the functional manufacturing of surfaces. Full article
Show Figures

Figure 1

Figure 1
<p>Surface topography after ball-end milling machining.</p>
Full article ">Figure 2
<p>Surface topography of ball-end milling at different z-values: (<b>a</b>) Group A; (<b>b</b>) Group B; (<b>c</b>) Group C; (<b>d</b>) Group D; (<b>e</b>) Group E; (<b>f</b>) Group F; (<b>g</b>) Group G; (h) Group H; (<b>i</b>) Group I.</p>
Full article ">Figure 2 Cont.
<p>Surface topography of ball-end milling at different z-values: (<b>a</b>) Group A; (<b>b</b>) Group B; (<b>c</b>) Group C; (<b>d</b>) Group D; (<b>e</b>) Group E; (<b>f</b>) Group F; (<b>g</b>) Group G; (h) Group H; (<b>i</b>) Group I.</p>
Full article ">Figure 3
<p>Finite element simulation model of sliding wear process: (<b>a</b>) geometrical modelling of the sliding wear process; (<b>b</b>) meshing of the contact area.</p>
Full article ">Figure 4
<p>Characteristics of variation of surface texture parameters: (<b>a</b>) <span class="html-italic">S</span>p and <span class="html-italic">S</span>v; (<b>b</b>) <span class="html-italic">S</span>z and <span class="html-italic">S</span>xp; (<b>c</b>) <span class="html-italic">S</span>q and <span class="html-italic">S</span>a; (<b>d</b>) <span class="html-italic">S</span>sk ; (<b>e</b>) <span class="html-italic">S</span>ku ; (<b>f</b>) <span class="html-italic">S</span>dq ; (<b>g</b>) <span class="html-italic">S</span>dr.</p>
Full article ">Figure 5
<p>Change characteristics of contact pressure during sliding wear: (<b>a</b>) Group A; (<b>b</b>) Group B; (<b>c</b>) Group C; (<b>d</b>) Group D; (<b>e</b>) Group E; (<b>f</b>) Group F; (<b>g</b>) Group G; (h) Group H; (<b>i</b>) Group I.</p>
Full article ">Figure 6
<p>Variation characteristics of wear deformation: (<b>a</b>) Group A; (<b>b</b>) Group B; (<b>c</b>) Group C; (<b>d</b>) Group D; (<b>e</b>) Group E; (<b>f</b>) Group F; (<b>g</b>) Group G; (h) Group H; (<b>i</b>) Group I.</p>
Full article ">Figure 7
<p>Ball-end milling experiment equipment: (<b>a</b>) 5-axis vertical machining centre; (<b>b</b>) ball-end milling machining process.</p>
Full article ">Figure 8
<p>Cutting tools for ball-end milling: (<b>a</b>) cutting tool shank; (<b>b</b>) blade.</p>
Full article ">Figure 9
<p>Surface after ball-end milling: (<b>a</b>) Groups M; (<b>b</b>) Groups N.</p>
Full article ">Figure 10
<p>Surface topography after ball-end milling process: (<b>a</b>) Groups M; (<b>b</b>) Groups N.</p>
Full article ">Figure 11
<p>Reciprocating sliding friction and wear test instrument: (<b>a</b>) MFT-5000 multifunctional friction and wear tester; (<b>b</b>) reciprocating sliding friction wear test.</p>
Full article ">Figure 12
<p>Surface after sliding wear: (<b>a</b>) Groups M; (<b>b</b>) Groups N.</p>
Full article ">
26 pages, 4486 KiB  
Article
Bleached Hair as Standard Template to Insight the Performance of Commercial Hair Repair Products
by Eva Martins, Pedro Castro, Alessandra B. Ribeiro, Carla F. Pereira, Francisca Casanova, Rui Vilarinho, Joaquim Moreira and Óscar L. Ramos
Cosmetics 2024, 11(5), 150; https://doi.org/10.3390/cosmetics11050150 - 28 Aug 2024
Viewed by 1218
Abstract
The increasing demand for effective hair care products has highlighted the necessity for rigorous claims substantiation methods, particularly for products that target specific hair types. This is essential because the effectiveness of a product can vary significantly based on the hair’s condition and [...] Read more.
The increasing demand for effective hair care products has highlighted the necessity for rigorous claims substantiation methods, particularly for products that target specific hair types. This is essential because the effectiveness of a product can vary significantly based on the hair’s condition and characteristics. A well-defined bleaching protocol is crucial for creating a standardized method to assess product efficacy, especially for products designed to repair damaged hair. The objective of this study was to create a practical bleaching protocol that mimics real-world consumer experiences, ensuring that hair samples exhibit sufficient damage for testing. This approach allows for a reliable assessment of how well various products can repair hair. The protocol serves as a framework for evaluating hair properties and the specific effects of each product on hair structure. Color, brightness, lightness, morphology, and topography were primarily used to understand the big differences in the hair fiber when treated with two repair benchmark products, K18® and Olaplex®, in relation to the Bleached hair. The devised bleaching protocol proved to be a fitting framework for assessing the properties of hair and the unique characteristics of each tested product within the hair fiber. This protocol offers valuable insights and tools for substantiating consumer claims, with morphological and mechanical methods serving as indispensable tools for recognizing and validating claims related to hair. The addition of K18® and Olaplex® demonstrated an increase in hair brightness (Y) and lightness (L* and a*) in relation to the Bleached samples, which were considered relevant characteristics for consumers. Olaplex®’s water-based nature creates a visible inner sheet, effectively filling empty spaces and improving the disulfide linkage network. This enhancement was corroborated by the increased number of disulfide bonds and evident changes in the FTIR profile. In contrast, K18®, owing to the lipophilic nature of its constituents, resulted in the formation of an external layer above the fiber. The composition of each of the products had a discrete impact on the fiber distribution, which was an outcome relevant to the determination of spreadability by consumers. Full article
(This article belongs to the Special Issue 10th Anniversary of Cosmetics—Recent Advances and Perspectives)
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Analytical tools employed in this study for hair fiber characterization.</p>
Full article ">Figure 2
<p>Photograph of hair tress appearance (<b>A</b>); hair tress brightness (Y) (<b>B</b>), lightness (<b>C</b>), a* (<b>D</b>), and b* (<b>E</b>) for K18<sup>®</sup> and Olaplex<sup>®</sup>-treated samples in comparison to bleached samples. The I, II, and III Roman numerals showed the variation of change for Y, L*, and a*, respectively, for hair tresses after the application of K18<sup>®</sup> or Olaplex<sup>®</sup> in relation to the bleached hair. Asterisks (*) denote a statistically significant difference (<span class="html-italic">p</span> &lt; 0.05). Values represent the mean ± SD.</p>
Full article ">Figure 3
<p>SEM micrographs of three distinct hair monofilament regions: (<b>A</b>–<b>C</b>) root, (<b>D</b>–<b>F</b>) middle, and end (<b>G</b>–<b>I</b>) of the bleached hair samples before bleaching optimization. SEM imaging to observe hair samples: (<b>J</b>–<b>L</b>) Bleached (<b>M</b>–<b>O</b>) K18<sup>®</sup>, and (<b>P</b>–<b>R</b>) Olaplex<sup>®</sup> of the surface morphology at 1000×, 2500×, and 5000× magnifications.</p>
Full article ">Figure 4
<p>Atomic Force Microscopy (AFM) images and calculated surface roughness of Bleached, K18<sup>®</sup> and Olaplex<sup>®</sup>-treated hair samples. 2D and 3D topography images (20 μm × 20 μm) for the (<b>A</b>,<b>B</b>) Bleached hair sample; (<b>C</b>,<b>D</b>) K18<sup>®</sup> and (<b>E</b>,<b>F</b>) Olaplex<sup>®</sup>. Quantitative parameters of the topography for hair samples as (<b>G</b>) Average Roughness—Ra (nm); (<b>H</b>) RMS Roughness—Rq (nm); and (<b>I</b>) Peak to valley Roughness Rt (nm). The asterisks (*) denote a statistically significant difference (<span class="html-italic">p</span> &lt; 0.05). Values represent the mean ± SD of at least two individual measurements.</p>
Full article ">Figure 5
<p>Density (<b>A</b>) and protein (<b>B</b>) content (%) results obtained for Bleached, K18<sup>®,</sup> and Olaplex<sup>®</sup> hair-treated samples. The asterisks (*) denote a statistically significant difference (<span class="html-italic">p</span> &lt; 0.05). Values represent the mean ± SD.</p>
Full article ">Figure 6
<p>FTIR profiles of hair samples. (<b>A</b>) Overlapped FTIR spectra of Bleached, K18<sup>®</sup>, and Olaplex<sup>®</sup> hair samples. (<b>B</b>–<b>D</b>) Amide I deconvolution peaks of Bleached, K18<sup>®</sup>, and Olaplex<sup>®</sup>, respectively.</p>
Full article ">Figure 7
<p>Disulfide bonds assessed by Raman spectroscopy (<b>A</b>) and S-R bonds were evaluated using Ellman’s chemical procedure for the tested samples (<b>B</b>). Statistical significance between the samples was verified by an asterisk (*) symbol at <span class="html-italic">p</span> &lt; 0.05.</p>
Full article ">Figure 8
<p>Thermal analysis of Bleached, K18<sup>®</sup>, and Olaplex<sup>®</sup> hair samples. (<b>A</b>) DSC thermograms of hair samples heated from 20 °C to 270 °C at a heating rate of 10 °C/min. The DSC curves for water evaporation and keratin denaturation (°C) are presented as I and II, respectively. (<b>B</b>–<b>F</b>) DSC results of water evaporation, water removal enthalpy (mW/mg), keratin peak and denaturation enthalpy of keratin (mW/mg), and α-helix of keratin peak for statistical analysis determination. III and IV symbols showed the variation of keratin peak in relation to Bleached. The asterisks (*) denote a statistically significant difference (<span class="html-italic">p</span> &lt; 0.05). Values represent the mean ± SD (<span class="html-italic">n</span> = 6).</p>
Full article ">Figure 9
<p>TGA analysis of K18<sup>®</sup> and Olapex<sup>®</sup> hair-treated samples in relation to the bleached sample (used as a control). (<b>A</b>) Weight (%) and (<b>B</b>) DTG (%, °C) curves of the TGA experiments.</p>
Full article ">Figure 10
<p>Physical and mechanical properties of hair fibers. (<b>A</b>–<b>C</b>) Thickness (μm) of Bleached, K18<sup>®</sup>, and Olaplex<sup>®</sup>, respectively. (<b>D</b>) Mean and Standard deviation of hair samples (n = 25). (<b>E</b>) Tensile strength (N), I symbol shows the variation of tensile strength in relation to Bleached, and (<b>F</b>) Extensibility (mm) of the hair strands obtained by Tensile tests. The asterisks (*) denote a statistically significant difference (<span class="html-italic">p</span> &lt; 0.05).</p>
Full article ">Figure 11
<p>Schematic representation of Bleached, K18<sup>®</sup> and Olaplex<sup>®</sup> after full characterization of hair fiber.</p>
Full article ">
17 pages, 9708 KiB  
Article
Analysing Temporal Evolution of OpenStreetMap Waterways Completeness in a Mountain Region of Portugal
by Elisabete S. Veiga Monteiro and Glória Rodrigues Patrício
Remote Sens. 2024, 16(17), 3159; https://doi.org/10.3390/rs16173159 - 27 Aug 2024
Viewed by 289
Abstract
In recent decades, the creation and availability of Voluntary Geographic Information (VGI) have changed the paradigm associated with the production of Geospatial Information (GI), since, due to its free access, citizens can view, analyse, process, and validate this type of data. One of [...] Read more.
In recent decades, the creation and availability of Voluntary Geographic Information (VGI) have changed the paradigm associated with the production of Geospatial Information (GI), since, due to its free access, citizens can view, analyse, process, and validate this type of data. One of the most popular examples of VGI is the collaborative OpenStreetMap (OSM) project which covers a wide range of themes or characteristics associated with the real world. One of these themes is the feature “waterway” that represents watercourses. The quality of OSM data characteristics is a topic that has been published by many authors in recent years, particularly on the analysis of the completeness indicator. However, few references are found in the literature about studies that analyse the completeness of OSM watercourses or even watercourses obtained by other sources. All this motivated the authors to develop a study that aims to analyse the completeness of these specific lines that have so much relevance to hydrologists. The study presents an analysis of the variation over time in completeness/coverage of the OSM “waterway” feature in the period between 2014 and 2023 in a mountainous region included in the Mondego River basin, located in the Inland of Portugal. The methodology applied is supported by classical methods of measuring the completeness of lines that may be found in the literature. The total length of the watercourses was calculated and compared in percentage terms with the total length of the reference watercourses for dates under analysis. The watercourses of the military official hydrography of the 1/25,000 scale were used as a reference. The relation of the OSM completeness with some indicators related to terrain surface (altitude, slope, and location/proximity settlements) was also analysed. The choice of these indicators was motivated by the fact that the study area has strong mountain characteristics and is crossed by the main Portuguese river. The analysis was performed using the Shuttle Radar Topography Mission Digital Elevation Model (SRTM DEM) data and satellite image of Geographic Information System software. The results show that the completeness of this OSM feature (waterway) has a slight increase, considering the amplitude of the studied period (nine years) and the fact that, nowadays, digital mobile devices enable easy access to satellite images, allowing the digitalization of geographic entities or objects of the real world remotely. Regarding the indicator altitude, slope, and location/proximity of the settlements, we believe that there is no influence of these indicators on the evolution of the completeness of the OSM waterways in the study area. Full article
Show Figures

Figure 1

Figure 1
<p>Location of the study area: rectangular region within the Mondego River basin, which corresponds to the area of sheet number 201 of the M888 cartographic series of the Geospatial Information Center of the Portuguese Army.</p>
Full article ">Figure 2
<p>Reference data—Army Geospatial Information Center Hydrography at 1/25,000 scale.</p>
Full article ">Figure 3
<p>OSM waterways extracted in the study area in 2014 (orange colour) and in 2023 (green colour), over the satellite image.</p>
Full article ">Figure 4
<p>Z1 (orange colour) and Z2 (blue colour) zones were randomly defined and included in the study area.</p>
Full article ">Figure 5
<p>OSM watercourses extracted in 2014 (orange colour) and 2023 (green colour) and the reference watercourses.</p>
Full article ">Figure 6
<p>Total length of watercourses of the reference drainage network and of OSM watercourses in 2014 and 2023 (<b>left</b>) and completeness values of OSM watercourses in 2014 and 2023 (<b>right</b>).</p>
Full article ">Figure 7
<p>SRTM Digital Elevation Model of the study area and OSM watercourses extracted in 2014 (orange colour) and 2023 (green colour).</p>
Full article ">Figure 8
<p>Slope map of the study area with the representation of OSM watercourses in 2014 (orange colour) and in 2023 (green colour).</p>
Full article ">Figure 9
<p>News’s OSM waterways (extracted in 2023—green colour) crossing some settlements in the north bank of Mondego River (<b>a</b>) and a detailed image exemplifying the western settlement crossed (<b>b</b>).</p>
Full article ">Figure 10
<p>Reference watercourses and OSM watercourses in the Z1 zone in 2014 and 2023.</p>
Full article ">Figure 11
<p>Reference watercourses and OSM watercourses in the Z2 zone in 2014 and 2023.</p>
Full article ">Figure 12
<p>Completeness of OSM waterways (in 2014 and 2023) versus terrain altitude (DEM) in Z1 (left) and Z2 (right) zones.</p>
Full article ">Figure 13
<p>Completeness of OSM waterways (in 2014 and 2023) versus terrain slope in Z1 (left) and Z2 (right) zones.</p>
Full article ">Figure 14
<p>New’s OSM waterways (2023) in northeast and east part of Z1 zone.</p>
Full article ">Figure 15
<p>New’s OSM waterways (2023) in western and south part of Z2 zone.</p>
Full article ">
82 pages, 17098 KiB  
Review
Statistical Dynamics and Subgrid Modelling of Turbulence: From Isotropic to Inhomogeneous
by Jorgen S. Frederiksen, Vassili Kitsios and Terence J. O’Kane
Atmosphere 2024, 15(8), 921; https://doi.org/10.3390/atmos15080921 - 31 Jul 2024
Cited by 1 | Viewed by 543
Abstract
Turbulence is the most important, ubiquitous, and difficult problem of classical physics. Feynman viewed it as essentially unsolved, without a rigorous mathematical basis to describe the statistical dynamics of this most complex of fluid motion. However, the paradigm shift came in 1959, with [...] Read more.
Turbulence is the most important, ubiquitous, and difficult problem of classical physics. Feynman viewed it as essentially unsolved, without a rigorous mathematical basis to describe the statistical dynamics of this most complex of fluid motion. However, the paradigm shift came in 1959, with the formulation of the Eulerian direct interaction approximation (DIA) closure by Kraichnan. It was based on renormalized perturbation theory, like quantum electrodynamics, and is a bare vertex theory that is manifestly realizable. Here, we review some of the subsequent exciting achievements in closure theory and subgrid modelling. We also document in some detail the progress that has been made in extending statistical dynamical turbulence theory to the real world of interactions with mean flows, waves and inhomogeneities such as topography. This includes numerically efficient inhomogeneous closures, like the realizable quasi-diagonal direct interaction approximation (QDIA), and even more efficient Markovian Inhomogeneous Closures (MICs). Recent developments include the formulation and testing of an eddy-damped Markovian anisotropic closure (EDMAC) that is realizable in interactions with transient waves but is as efficient as the eddy-damped quasi-normal Markovian (EDQNM). As well, a similarly efficient closure, the realizable eddy-damped Markovian inhomogeneous closure (EDMIC) has been developed. Moreover, we present subgrid models that cater for the complex interactions that occur in geophysical flows. Recent progress includes the determination of complete sets of subgrid terms for skilful large-eddy simulations of baroclinic inhomogeneous turbulent atmospheric and oceanic flows interacting with Rossby waves and topography. The success of these inhomogeneous closures has also led to further applications in data assimilation and ensemble prediction and generalization to quantum fields. Full article
(This article belongs to the Special Issue Isotropic Turbulence: Recent Advances and Current Challenges)
Show Figures

Figure 1

Figure 1
<p>Comparison of transient kinetic energy spectra <inline-formula><mml:math id="mm1012"><mml:semantics><mml:mrow><mml:mi>E</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mi>k</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:semantics></mml:math></inline-formula> for the DIA, SCFT and LET closures compared with DNS models for initial spectrum B at initial (<inline-formula><mml:math id="mm230"><mml:semantics><mml:mrow><mml:mi>t</mml:mi><mml:mo>=</mml:mo><mml:mn>0</mml:mn></mml:mrow></mml:semantics></mml:math></inline-formula>) and final times (<inline-formula><mml:math id="mm231"><mml:semantics><mml:mrow><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mi>max</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mn>0.8</mml:mn></mml:mrow></mml:semantics></mml:math></inline-formula>). Adapted from Frederiksen and Davies (2000) [<xref ref-type="bibr" rid="B31-atmosphere-15-00921">31</xref>].</p>
Full article ">Figure 2
<p>Comparison of transient kinetic energy spectra <inline-formula><mml:math id="mm1013"><mml:semantics><mml:mrow><mml:mi>E</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mi>k</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:semantics></mml:math></inline-formula> for the DIA and RDIA closures compared with DNS models for initial spectrum B at initial (<inline-formula><mml:math id="mm252"><mml:semantics><mml:mrow><mml:mi>t</mml:mi><mml:mo>=</mml:mo><mml:mn>0</mml:mn></mml:mrow></mml:semantics></mml:math></inline-formula>) and final times (<inline-formula><mml:math id="mm253"><mml:semantics><mml:mrow><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mi>max</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mn>0.8</mml:mn></mml:mrow></mml:semantics></mml:math></inline-formula>). Adapted from Frederiksen and Davies (2004) [<xref ref-type="bibr" rid="B33-atmosphere-15-00921">33</xref>].</p>
Full article ">Figure 3
<p>Comparison of transient kinetic energy spectra <inline-formula><mml:math id="mm1014"><mml:semantics><mml:mrow><mml:mi>E</mml:mi><mml:mo stretchy="false">(</mml:mo><mml:mi>k</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:semantics></mml:math></inline-formula> for the EDQNM and EDMAC models for the runs in <xref ref-type="table" rid="atmosphere-15-00921-t002">Table 2</xref> at the initial (<inline-formula><mml:math id="mm350"><mml:semantics><mml:mrow><mml:mi>t</mml:mi><mml:mo>=</mml:mo><mml:mn>0</mml:mn></mml:mrow></mml:semantics></mml:math></inline-formula>) and final times (<inline-formula><mml:math id="mm351"><mml:semantics><mml:mrow><mml:msub><mml:mi>t</mml:mi><mml:mrow><mml:mi>max</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mn>0.4</mml:mn></mml:mrow></mml:semantics></mml:math></inline-formula>). Adapted from Frederiksen and O’Kane (2024) [<xref ref-type="bibr" rid="B65-atmosphere-15-00921">65</xref>].</p>
Full article ">Figure 4
<p>The diagrammatic representation of the integral terms in the QDIA mean-field <inline-formula><mml:math id="mm1015"><mml:semantics><mml:mrow><mml:mo>&lt;</mml:mo><mml:msub><mml:mi>ζ</mml:mi><mml:mstyle mathvariant="bold" mathsize="normal"><mml:mi>k</mml:mi></mml:mstyle></mml:msub><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy="false">)</mml:mo><mml:mo>&gt;</mml:mo></mml:mrow></mml:semantics></mml:math></inline-formula> tendency equation. Note that the Jacobian requires no diagrammatic representation. For the DIA, the mean-field is zero.</p>
Full article ">Figure 5
<p>The diagrammatic representation of the integral terms in the QDIA diagonal two-point two-time cumulant <inline-formula><mml:math id="mm1016"><mml:semantics><mml:mrow><mml:msub><mml:mi>C</mml:mi><mml:mstyle mathvariant="bold" mathsize="normal"><mml:mi>k</mml:mi></mml:mstyle></mml:msub><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo>,</mml:mo><mml:msup><mml:mi>t</mml:mi><mml:mo>′</mml:mo></mml:msup><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:semantics></mml:math></inline-formula> tendency equation. The terms that the QDIA has in common with the DIA are within the red rectangle.</p>
Full article ">Figure 6
<p>The diagrammatic representation of the integral terms in the QDIA diagonal response function <inline-formula><mml:math id="mm1017"><mml:semantics><mml:mrow><mml:msub><mml:mi>R</mml:mi><mml:mstyle mathvariant="bold" mathsize="normal"><mml:mi>k</mml:mi></mml:mstyle></mml:msub><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo>,</mml:mo><mml:msup><mml:mi>t</mml:mi><mml:mo>′</mml:mo></mml:msup><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:semantics></mml:math></inline-formula> tendency equation. The terms that the QDIA has in common with the DIA are within the red rectangle.</p>
Full article ">Figure 7
<p>Initial (<bold>a</bold>) and evolved day 10 (<bold>b</bold>) <inline-formula><mml:math id="mm1018"><mml:semantics><mml:mrow><mml:mi>Q</mml:mi><mml:mi>D</mml:mi><mml:mi>I</mml:mi><mml:mi>A</mml:mi><mml:mo stretchy="false">[</mml:mo><mml:mover accent="true"><mml:mi>ζ</mml:mi><mml:mo>¯</mml:mo></mml:mover><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy="false">)</mml:mo><mml:mo stretchy="false">]</mml:mo></mml:mrow></mml:semantics></mml:math></inline-formula>, (<bold>c</bold>) <inline-formula><mml:math id="mm495"><mml:semantics><mml:mrow><mml:mi>M</mml:mi><mml:mi>I</mml:mi><mml:msup><mml:mi>C</mml:mi><mml:mrow><mml:mi>X</mml:mi><mml:mo>=</mml:mo><mml:mn>0</mml:mn></mml:mrow></mml:msup><mml:mo stretchy="false">[</mml:mo><mml:mover accent="true"><mml:mi>ζ</mml:mi><mml:mo>¯</mml:mo></mml:mover><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy="false">)</mml:mo><mml:mo stretchy="false">]</mml:mo></mml:mrow></mml:semantics></mml:math></inline-formula>, (<bold>d</bold>) <inline-formula><mml:math id="mm496"><mml:semantics><mml:mrow><mml:mi>M</mml:mi><mml:mi>I</mml:mi><mml:msup><mml:mi>C</mml:mi><mml:mrow><mml:mi>X</mml:mi><mml:mo>=</mml:mo><mml:mstyle scriptlevel="+1"><mml:mfrac bevelled="true"><mml:mn>1</mml:mn><mml:mn>2</mml:mn></mml:mfrac></mml:mstyle></mml:mrow></mml:msup><mml:mo stretchy="false">[</mml:mo><mml:mover accent="true"><mml:mi>ζ</mml:mi><mml:mo>¯</mml:mo></mml:mover><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy="false">)</mml:mo><mml:mo stretchy="false">]</mml:mo></mml:mrow></mml:semantics></mml:math></inline-formula>, (<bold>e</bold>) <inline-formula><mml:math id="mm497"><mml:semantics><mml:mrow><mml:mi>M</mml:mi><mml:mi>I</mml:mi><mml:msup><mml:mi>C</mml:mi><mml:mrow><mml:mi>X</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msup><mml:mo stretchy="false">[</mml:mo><mml:mover accent="true"><mml:mi>ζ</mml:mi><mml:mo>¯</mml:mo></mml:mover><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy="false">)</mml:mo><mml:mo stretchy="false">]</mml:mo></mml:mrow></mml:semantics></mml:math></inline-formula>, and (<bold>f</bold>) ensemble of DNS Rossby wave streamfunction in non-dimensional units. Multiply values by <inline-formula><mml:math id="mm498"><mml:semantics><mml:mrow><mml:mstyle scriptlevel="+1"><mml:mfrac><mml:mn>1</mml:mn><mml:mn>4</mml:mn></mml:mfrac></mml:mstyle><mml:msubsup><mml:mi>a</mml:mi><mml:mi>e</mml:mi><mml:mn>2</mml:mn></mml:msubsup><mml:msup><mml:mo mathvariant="bold">Ω</mml:mo><mml:mrow><mml:mo>−</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msup><mml:mo>=</mml:mo><mml:mn>740</mml:mn><mml:msup><mml:mrow><mml:mrow><mml:mo> </mml:mo><mml:mi>km</mml:mi></mml:mrow></mml:mrow><mml:mn>2</mml:mn></mml:msup><mml:msup><mml:mi mathvariant="normal">s</mml:mi><mml:mrow><mml:mrow><mml:mo>-</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:mrow></mml:msup></mml:mrow></mml:semantics></mml:math></inline-formula> to convert to units of <inline-formula><mml:math id="mm499"><mml:semantics><mml:mrow><mml:msup><mml:mrow><mml:mi>km</mml:mi></mml:mrow><mml:mn>2</mml:mn></mml:msup><mml:msup><mml:mi mathvariant="normal">s</mml:mi><mml:mrow><mml:mrow><mml:mo>-</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:mrow></mml:msup></mml:mrow></mml:semantics></mml:math></inline-formula>. From Frederiksen and O’Kane (2022) [<xref ref-type="bibr" rid="B75-atmosphere-15-00921">75</xref>].</p>
Full article ">Figure 8
<p>The initial and evolved non-dimensional kinetic energy spectra for the four inhomogeneous abridged closures in panels (<bold>a</bold>) <inline-formula><mml:math id="mm1019"><mml:semantics><mml:mrow><mml:mi>Q</mml:mi><mml:mi>D</mml:mi><mml:mi>I</mml:mi><mml:mi>A</mml:mi><mml:mo stretchy="false">[</mml:mo><mml:mover accent="true"><mml:mi>ζ</mml:mi><mml:mo>¯</mml:mo></mml:mover><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy="false">)</mml:mo><mml:mo stretchy="false">]</mml:mo></mml:mrow></mml:semantics></mml:math></inline-formula>, (<bold>b</bold>) <inline-formula><mml:math id="mm528"><mml:semantics><mml:mrow><mml:mi>M</mml:mi><mml:mi>I</mml:mi><mml:msup><mml:mi>C</mml:mi><mml:mrow><mml:mi>X</mml:mi><mml:mo>=</mml:mo><mml:mn>0</mml:mn></mml:mrow></mml:msup><mml:mo stretchy="false">[</mml:mo><mml:mover accent="true"><mml:mi>ζ</mml:mi><mml:mo>¯</mml:mo></mml:mover><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy="false">)</mml:mo><mml:mo stretchy="false">]</mml:mo></mml:mrow></mml:semantics></mml:math></inline-formula>, (<bold>c</bold>) <inline-formula><mml:math id="mm529"><mml:semantics><mml:mrow><mml:mi>M</mml:mi><mml:mi>I</mml:mi><mml:msup><mml:mi>C</mml:mi><mml:mrow><mml:mstyle scriptlevel="+1"><mml:mfrac bevelled="true"><mml:mn>1</mml:mn><mml:mn>2</mml:mn></mml:mfrac></mml:mstyle></mml:mrow></mml:msup><mml:mo stretchy="false">[</mml:mo><mml:mover accent="true"><mml:mi>ζ</mml:mi><mml:mo>¯</mml:mo></mml:mover><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy="false">)</mml:mo><mml:mo stretchy="false">]</mml:mo></mml:mrow></mml:semantics></mml:math></inline-formula>, (<bold>d</bold>) <inline-formula><mml:math id="mm530"><mml:semantics><mml:mrow><mml:mi>M</mml:mi><mml:mi>I</mml:mi><mml:msup><mml:mi>C</mml:mi><mml:mrow><mml:mi>X</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msup><mml:mo stretchy="false">[</mml:mo><mml:mover accent="true"><mml:mi>ζ</mml:mi><mml:mo>¯</mml:mo></mml:mover><mml:mo stretchy="false">(</mml:mo><mml:mi>t</mml:mi><mml:mo stretchy="false">)</mml:mo><mml:mo stretchy="false">]</mml:mo></mml:mrow></mml:semantics></mml:math></inline-formula>, as compared to ensemble mean of 1800 DNSs. In each panel we show the initial mean energy (solid blue), initial transient energy (dashed blue), evolved DNS mean energy (solid black), evolved DNS transient energy (dashed black), evolved closure mean energy (solid red) and evolved closure transient energy (dashed red). Multiplication by <inline-formula><mml:math id="mm531"><mml:semantics><mml:mrow><mml:mstyle scriptlevel="+1"><mml:mfrac><mml:mn>1</mml:mn><mml:mn>4</mml:mn></mml:mfrac></mml:mstyle><mml:msubsup><mml:mi>a</mml:mi><mml:mi>e</mml:mi><mml:mn>2</mml:mn></mml:msubsup><mml:msup><mml:mo mathvariant="bold">Ω</mml:mo><mml:mrow><mml:mo>−</mml:mo><mml:mn>2</mml:mn></mml:mrow></mml:msup><mml:mo>=</mml:mo><mml:mn>5.4</mml:mn><mml:mo>×</mml:mo><mml:msup><mml:mrow><mml:mn>10</mml:mn></mml:mrow><mml:mn>4</mml:mn></mml:msup><mml:msup><mml:mrow><mml:mrow><mml:mo> </mml:mo><mml:mi mathvariant="normal">m</mml:mi></mml:mrow></mml:mrow><mml:mn>2</mml:mn></mml:msup><mml:msup><mml:mi mathvariant="normal">s</mml:mi><mml:mrow><mml:mrow><mml:mo>-</mml:mo><mml:mn>2</mml:mn></mml:mrow></mml:mrow></mml:msup></mml:mrow></mml:semantics></mml:math></inline-formula> is required to convert to units of <inline-formula><mml:math id="mm532"><mml:semantics><mml:mrow><mml:msup><mml:mi mathvariant="normal">m</mml:mi><mml:mn>2</mml:mn></mml:msup><mml:msup><mml:mi mathvariant="normal">s</mml:mi><mml:mrow><mml:mrow><mml:mo>-</mml:mo><mml:mn>2</mml:mn></mml:mrow></mml:mrow></mml:msup></mml:mrow></mml:semantics></mml:math></inline-formula>. From Frederiksen and O’Kane (2022) [<xref ref-type="bibr" rid="B75-atmosphere-15-00921">75</xref>].</p>
Full article ">Figure 9
<p>The physical space plots of the eddy streamfunction (<inline-formula><mml:math id="mm1020"><mml:semantics><mml:mrow><mml:msup><mml:mi mathvariant="normal">m</mml:mi><mml:mn>2</mml:mn></mml:msup><mml:msup><mml:mi mathvariant="normal">s</mml:mi><mml:mrow><mml:mo>−</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:semantics></mml:math></inline-formula>) for both DNS and QDIA models at the initial day (29 October) and day 5 (3 November) of the breeding cycle, and after five forecast days (8 November). Adapted from O’Kane and Frederiksen (2008) [<xref ref-type="bibr" rid="B68-atmosphere-15-00921">68</xref>].</p>
Full article ">Figure 10
<p>(<bold>a</bold>) EDQNM-based nondimensional subgrid-scale parameterizations for T31 from Frederiksen and Davies (1997) [<xref ref-type="bibr" rid="B126-atmosphere-15-00921">126</xref>], for drain (line A), backscatter (line B) and net (line C) eddy viscosities. (<bold>b</bold>) as in (<bold>a</bold>) for DNS-based subgrid terms calculated using the stochastic modelling approach of Frederiksen and Kepert (2006) [<xref ref-type="bibr" rid="B127-atmosphere-15-00921">127</xref>]. Adapted from Frederiksen and Davies (1997) [<xref ref-type="bibr" rid="B126-atmosphere-15-00921">126</xref>] and Frederiksen and Kepert (2006) [<xref ref-type="bibr" rid="B127-atmosphere-15-00921">127</xref>].</p>
Full article ">Figure 11
<p>(<bold>a</bold>) Kinetic energy spectra from Frederiksen and Davies (1997) [<xref ref-type="bibr" rid="B126-atmosphere-15-00921">126</xref>] for isotropised observations from January 1979 (line A), DNS at T63 (line B), LES at T31 (line C) and kinetic energy plus and minus one standard deviation (line D). LES with EDQNM-based stochastic backscatter and drain eddy viscosity is scaled by 10<sup>−1</sup>, and LES with EDQNM-based net eddy viscosity is scaled by 10<sup>−2</sup>. (<bold>b</bold>) as in (<bold>a</bold>) for DNS-based subgrid terms calculated using the stochastic modelling approach of Frederiksen and Kepert (2006) [<xref ref-type="bibr" rid="B127-atmosphere-15-00921">127</xref>]. Adapted from Frederiksen and Davies (1997) [<xref ref-type="bibr" rid="B126-atmosphere-15-00921">126</xref>] and Frederiksen and Kepert (2006) [<xref ref-type="bibr" rid="B127-atmosphere-15-00921">127</xref>].</p>
Full article ">Figure 12
<p>Kinetic energy spectra at 250 hPa for LES at T63 (solid black line) for DNS at T126 truncated back to T63 (black dashed line) and DNS plus and minus one standard deviation (red dotted lines). Top set of spectra are for LES with only bare dissipation. Middle set, shifted down one decade from true value, are for LES with deterministic net dissipation. Bottom set, shifted down two decades, are for LES with drain dissipation and stochastic backscatter. Adapted from Zidikheri and Frederiksen (2009) [<xref ref-type="bibr" rid="B156-atmosphere-15-00921">156</xref>].</p>
Full article ">Figure 13
<p>QG baroclinic atmospheric flow: (<bold>a</bold>) anisotropic drain eddy viscosity at <inline-formula><mml:math id="mm1021"><mml:semantics><mml:mrow><mml:msub><mml:mi>T</mml:mi><mml:mi>R</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mn>63</mml:mn></mml:mrow></mml:semantics></mml:math></inline-formula>; (<bold>b</bold>) anisotropic backscatter eddy viscosity at <inline-formula><mml:math id="mm631"><mml:semantics><mml:mrow><mml:msub><mml:mi>T</mml:mi><mml:mi>R</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mn>63</mml:mn></mml:mrow></mml:semantics></mml:math></inline-formula>; (<bold>c</bold>) isotropised drain, backscatter and net eddy viscosities for multiple truncation levels; (<bold>d</bold>) value of the drain eddy viscosity at <inline-formula><mml:math id="mm632"><mml:semantics><mml:mrow><mml:mi>n</mml:mi><mml:mo>=</mml:mo><mml:msub><mml:mi>T</mml:mi><mml:mi>R</mml:mi></mml:msub></mml:mrow></mml:semantics></mml:math></inline-formula> for various truncations, levels and atmospheres, with solid line the scaling law for the drain eddy viscosity, and dashed line for the net; (<bold>e</bold>) as in (<bold>d</bold>) but for the spectral slope; and (<bold>f</bold>) comparison between DNS (dashed black) and LES (solid magenta) using the parameterization variants of anisotropic stochastic (AS), anisotropic deterministic (AD), isotropic stochastic (IS), isotropic deterministic (ID), scaling law stochastic (LS) and scaling law deterministic (LD). Adapted from Kitsios et al. (2012) [<xref ref-type="bibr" rid="B157-atmosphere-15-00921">157</xref>].</p>
Full article ">Figure 14
<p>QG baroclinic oceanic flow simulations: (<bold>a</bold>) anisotropic drain eddy viscosity at <inline-formula><mml:math id="mm1022"><mml:semantics><mml:mrow><mml:msub><mml:mi>T</mml:mi><mml:mi>R</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mn>63</mml:mn></mml:mrow></mml:semantics></mml:math></inline-formula>; (<bold>b</bold>) anisotropic backscatter eddy viscosity at <inline-formula><mml:math id="mm650"><mml:semantics><mml:mrow><mml:msub><mml:mi>T</mml:mi><mml:mi>R</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mn>63</mml:mn></mml:mrow></mml:semantics></mml:math></inline-formula>; (<bold>c</bold>) isotropised drain, backscatter and net eddy viscosities for multiple truncation levels; (<bold>d</bold>) value of the drain eddy viscosity at <inline-formula><mml:math id="mm651"><mml:semantics><mml:mrow><mml:mi>n</mml:mi><mml:mo>=</mml:mo><mml:msub><mml:mi>T</mml:mi><mml:mi>R</mml:mi></mml:msub></mml:mrow></mml:semantics></mml:math></inline-formula> for various truncations, levels and oceans; (<bold>e</bold>) as in (<bold>d</bold>) but for the spectral slope; and (<bold>f</bold>) comparison between DNS (dashed black) and LES (solid magenta using the parameterization variants of anisotropic stochastic (AS), anisotropic deterministic (AD), isotropic stochastic (IS), isotropic deterministic (ID), scaling law stochastic (LS) and scaling law deterministic (LD). Adapted from Kitsios et al. (2013) [<xref ref-type="bibr" rid="B160-atmosphere-15-00921">160</xref>].</p>
Full article ">Figure 15
<p>Primitive equation model-based subgrid-scale parameterizations as a function of wavenumber <inline-formula><mml:math id="mm1023"><mml:semantics><mml:mi>n</mml:mi></mml:semantics></mml:math></inline-formula> and vertical <inline-formula><mml:math id="mm657"><mml:semantics><mml:mrow><mml:mi>σ</mml:mi><mml:mo>−</mml:mo></mml:mrow></mml:semantics></mml:math></inline-formula>level: (<bold>a</bold>) drain dissipation for vorticity; (<bold>b</bold>) backscatter dissipation for vorticity; (<bold>c</bold>) net dissipation for vorticity; (<bold>d</bold>) net dissipation for divergence; (<bold>e</bold>) net dissipation for temperature; and (<bold>f</bold>) net dissipation for surface pressure. Adapted from Frederiksen et al. (2015) [<xref ref-type="bibr" rid="B138-atmosphere-15-00921">138</xref>].</p>
Full article ">Figure 16
<p>Upper diagonal component of net subgrid dissipation matrices from Kitsios et al. (2017) [<xref ref-type="bibr" rid="B163-atmosphere-15-00921">163</xref>] illustrated in the horizontal wavenumber plane (<inline-formula><mml:math id="mm1024"><mml:semantics><mml:mrow><mml:msub><mml:mi>k</mml:mi><mml:mi>x</mml:mi></mml:msub><mml:mo>≥</mml:mo><mml:mn>0</mml:mn><mml:mo>,</mml:mo><mml:msub><mml:mi>k</mml:mi><mml:mi>z</mml:mi></mml:msub><mml:mo>≥</mml:mo><mml:mn>0</mml:mn></mml:mrow></mml:semantics></mml:math></inline-formula>) using rectangular truncation <inline-formula><mml:math id="mm691"><mml:semantics><mml:mrow><mml:msub><mml:mi>T</mml:mi><mml:mi>R</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mn>63</mml:mn></mml:mrow></mml:semantics></mml:math></inline-formula> for wall normal scales: (<bold>a</bold>) <inline-formula><mml:math id="mm692"><mml:semantics><mml:mrow><mml:mrow><mml:mrow><mml:msub><mml:mi>λ</mml:mi><mml:mi>y</mml:mi></mml:msub><mml:msub><mml:mi>u</mml:mi><mml:mi>τ</mml:mi></mml:msub></mml:mrow><mml:mo>/</mml:mo><mml:mrow><mml:mi>ν</mml:mi><mml:mo>=</mml:mo><mml:mn>100</mml:mn></mml:mrow></mml:mrow></mml:mrow></mml:semantics></mml:math></inline-formula>; (<bold>b</bold>) <inline-formula><mml:math id="mm693"><mml:semantics><mml:mrow><mml:mrow><mml:mrow><mml:msub><mml:mi>λ</mml:mi><mml:mi>y</mml:mi></mml:msub><mml:msub><mml:mi>u</mml:mi><mml:mi>τ</mml:mi></mml:msub></mml:mrow><mml:mo>/</mml:mo><mml:mi>ν</mml:mi></mml:mrow><mml:mo>=</mml:mo><mml:mn>10</mml:mn></mml:mrow></mml:semantics></mml:math></inline-formula>; (<bold>c</bold>) <inline-formula><mml:math id="mm694"><mml:semantics><mml:mrow><mml:mrow><mml:mrow><mml:msub><mml:mi>λ</mml:mi><mml:mi>y</mml:mi></mml:msub><mml:msub><mml:mi>u</mml:mi><mml:mi>τ</mml:mi></mml:msub></mml:mrow><mml:mo>/</mml:mo><mml:mi>ν</mml:mi></mml:mrow><mml:mo>=</mml:mo><mml:mn>2</mml:mn></mml:mrow></mml:semantics></mml:math></inline-formula>; and (<bold>d</bold>) <inline-formula><mml:math id="mm695"><mml:semantics><mml:mrow><mml:mrow><mml:mrow><mml:msub><mml:mi>λ</mml:mi><mml:mi>y</mml:mi></mml:msub><mml:msub><mml:mi>u</mml:mi><mml:mi>τ</mml:mi></mml:msub></mml:mrow><mml:mo>/</mml:mo><mml:mi>ν</mml:mi></mml:mrow><mml:mo>=</mml:mo><mml:mn>0.4</mml:mn></mml:mrow></mml:semantics></mml:math></inline-formula> where <inline-formula><mml:math id="mm696"><mml:semantics><mml:mrow><mml:msub><mml:mi>λ</mml:mi><mml:mi>y</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mn>1</mml:mn><mml:mo>−</mml:mo><mml:mi>cos</mml:mi><mml:mrow><mml:mo>(</mml:mo><mml:mrow><mml:mrow><mml:mrow><mml:mn>2</mml:mn><mml:mi>π</mml:mi></mml:mrow><mml:mo>/</mml:mo><mml:mrow><mml:msub><mml:mi>k</mml:mi><mml:mi>y</mml:mi></mml:msub></mml:mrow></mml:mrow></mml:mrow><mml:mo>)</mml:mo></mml:mrow></mml:mrow></mml:semantics></mml:math></inline-formula> with <inline-formula><mml:math id="mm697"><mml:semantics><mml:mrow><mml:msub><mml:mi>k</mml:mi><mml:mi>y</mml:mi></mml:msub></mml:mrow></mml:semantics></mml:math></inline-formula>, the Chebyshev polynomial index. All coefficients are symmetric about the line <inline-formula><mml:math id="mm698"><mml:semantics><mml:mrow><mml:msub><mml:mi>k</mml:mi><mml:mi>x</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mn>0</mml:mn></mml:mrow></mml:semantics></mml:math></inline-formula>. Adapted from Kitsios et al. (2017) [<xref ref-type="bibr" rid="B163-atmosphere-15-00921">163</xref>].</p>
Full article ">Figure 17
<p>Comparison of the DNS (black dotted line) and LES (magenta solid line) on the basis of time-averaged kinetic energy spectra. The two-dimensional spectrum is compared at <inline-formula><mml:math id="mm1025"><mml:semantics><mml:mrow><mml:msup><mml:mi>y</mml:mi><mml:mo>+</mml:mo></mml:msup><mml:mo>=</mml:mo><mml:mrow><mml:mrow><mml:mi>y</mml:mi><mml:msub><mml:mi>u</mml:mi><mml:mi>τ</mml:mi></mml:msub></mml:mrow><mml:mo>/</mml:mo><mml:mrow><mml:mi>ν</mml:mi><mml:mo>=</mml:mo></mml:mrow></mml:mrow><mml:mn>133</mml:mn></mml:mrow></mml:semantics></mml:math></inline-formula> for a LES with: (<bold>a</bold>) no subgrid parameterization; (<bold>b</bold>) net dissipation; and (<bold>c</bold>) drain dissipation with stochastic backscatter. Contour levels radiating out from the origin are <inline-formula><mml:math id="mm701"><mml:semantics><mml:mrow><mml:msup><mml:mrow><mml:mn>10</mml:mn></mml:mrow><mml:mrow><mml:mo>−</mml:mo><mml:mn>2.5</mml:mn></mml:mrow></mml:msup><mml:mo>,</mml:mo><mml:mo> </mml:mo><mml:msup><mml:mrow><mml:mn>10</mml:mn></mml:mrow><mml:mrow><mml:mo>−</mml:mo><mml:mn>3</mml:mn></mml:mrow></mml:msup><mml:mo>,</mml:mo><mml:mo> </mml:mo><mml:msup><mml:mrow><mml:mn>10</mml:mn></mml:mrow><mml:mrow><mml:mo>−</mml:mo><mml:mn>3.5</mml:mn></mml:mrow></mml:msup><mml:mo>,</mml:mo><mml:mo> </mml:mo><mml:msup><mml:mrow><mml:mn>10</mml:mn></mml:mrow><mml:mrow><mml:mo>−</mml:mo><mml:mn>4</mml:mn></mml:mrow></mml:msup><mml:mo>,</mml:mo><mml:mo> </mml:mo><mml:msup><mml:mrow><mml:mn>10</mml:mn></mml:mrow><mml:mrow><mml:mo>−</mml:mo><mml:mn>4.5</mml:mn></mml:mrow></mml:msup><mml:mo>,</mml:mo><mml:mo> </mml:mo><mml:msup><mml:mrow><mml:mn>10</mml:mn></mml:mrow><mml:mrow><mml:mo>−</mml:mo><mml:mn>5</mml:mn></mml:mrow></mml:msup><mml:mo>,</mml:mo><mml:mo> </mml:mo><mml:msup><mml:mrow><mml:mn>10</mml:mn></mml:mrow><mml:mrow><mml:mo>−</mml:mo><mml:mn>5.5</mml:mn></mml:mrow></mml:msup></mml:mrow></mml:semantics></mml:math></inline-formula>. The vertical axis in (<bold>a</bold>) is applicable to (<bold>b</bold>,<bold>c</bold>). Adapted from Kitsios et al. (2017) [<xref ref-type="bibr" rid="B163-atmosphere-15-00921">163</xref>].</p>
Full article ">Figure 18
<p>Subgrid-scale QDIA closure integral terms plotted in (<inline-formula><mml:math id="mm1026"><mml:semantics><mml:mrow><mml:msub><mml:mstyle mathvariant="bold" mathsize="normal"><mml:mi>k</mml:mi></mml:mstyle><mml:mi>x</mml:mi></mml:msub><mml:mo>,</mml:mo><mml:msub><mml:mstyle mathvariant="bold" mathsize="normal"><mml:mi>k</mml:mi></mml:mstyle><mml:mi>y</mml:mi></mml:msub></mml:mrow></mml:semantics></mml:math></inline-formula>)-space for the general case. Adapted from O’Kane and Frederiksen (2008) [<xref ref-type="bibr" rid="B71-atmosphere-15-00921">71</xref>].</p>
Full article ">Figure 19
<p>Subgrid-scale residual Jacobian terms <inline-formula><mml:math id="mm1027"><mml:semantics><mml:mrow><mml:msup><mml:mo>∇</mml:mo><mml:mrow><mml:mo>−</mml:mo><mml:mn>2</mml:mn></mml:mrow></mml:msup><mml:msub><mml:mi>J</mml:mi><mml:mi>S</mml:mi></mml:msub><mml:mo stretchy="false">(</mml:mo><mml:mi>ψ</mml:mi><mml:mo>,</mml:mo><mml:mi>ζ</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:semantics></mml:math></inline-formula> and <inline-formula><mml:math id="mm738"><mml:semantics><mml:mrow><mml:msup><mml:mo>∇</mml:mo><mml:mrow><mml:mo>−</mml:mo><mml:mn>2</mml:mn></mml:mrow></mml:msup><mml:msub><mml:mi>J</mml:mi><mml:mi>S</mml:mi></mml:msub><mml:mo stretchy="false">(</mml:mo><mml:mi>ψ</mml:mi><mml:mo>,</mml:mo><mml:mi>h</mml:mi><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:semantics></mml:math></inline-formula>, in terms of spectral components, with the zonal wind <inline-formula><mml:math id="mm739"><mml:semantics><mml:mi>U</mml:mi></mml:semantics></mml:math></inline-formula> and streamfunction <inline-formula><mml:math id="mm740"><mml:semantics><mml:mi>ψ</mml:mi></mml:semantics></mml:math></inline-formula> overlaid. Adapted from O’Kane and Frederiksen (2008) [<xref ref-type="bibr" rid="B71-atmosphere-15-00921">71</xref>].</p>
Full article ">Figure 20
<p>The comparison of DNS of an atmospheric flow of <inline-formula><mml:math id="mm1028"><mml:semantics><mml:mrow><mml:mi>T</mml:mi><mml:mo>=</mml:mo><mml:mn>63</mml:mn></mml:mrow></mml:semantics></mml:math></inline-formula> resolution to an LES of <inline-formula><mml:math id="mm759"><mml:semantics><mml:mrow><mml:msub><mml:mi>T</mml:mi><mml:mi>R</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mn>15</mml:mn></mml:mrow></mml:semantics></mml:math></inline-formula>. Base flow properties: (<bold>a</bold>) top level mean zonal velocity; and (<bold>b</bold>) topography. At the bottom level are the meridional acceleration of the (<bold>c</bold>) mean-field Jacobian; (<bold>d</bold>) eddy-mean-field; and (<bold>e</bold>) eddy–topographic terms. (<bold>f</bold>) The top-level kinetic energy spectra of DNS (dashed black) and stochastic LES (magenta). Each set of spectra has LES with the labelled terms included in the mean subgrid tendency term. Adapted from Kitsios and Frederiksen (2019) [<xref ref-type="bibr" rid="B137-atmosphere-15-00921">137</xref>].</p>
Full article ">Figure 21
<p>The comparison of the DNS of <inline-formula><mml:math id="mm1029"><mml:semantics><mml:mrow><mml:mi>T</mml:mi><mml:mo>=</mml:mo><mml:mn>504</mml:mn></mml:mrow></mml:semantics></mml:math></inline-formula> resolution and LES of an oceanic flow. Base flow properties: (<bold>a</bold>) top level mean zonal velocity; and (<bold>b</bold>) absolute value of the topographic gradient. For a <inline-formula><mml:math id="mm762"><mml:semantics><mml:mrow><mml:msub><mml:mi>T</mml:mi><mml:mi>R</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mn>252</mml:mn></mml:mrow></mml:semantics></mml:math></inline-formula> truncation, the bottom level meridional acceleration of the: (<bold>c</bold>) mean-field Jacobian; (<bold>d</bold>) eddy-mean-field; and (<bold>e</bold>) eddy–topographic terms. (<bold>f</bold>) The top-level kinetic energy spectra of DNS (dashed black) and stochastic LES (magenta). The top three spectra are for LES with <inline-formula><mml:math id="mm763"><mml:semantics><mml:mrow><mml:msub><mml:mi>T</mml:mi><mml:mi>R</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mn>252</mml:mn></mml:mrow></mml:semantics></mml:math></inline-formula> and bottom three with <inline-formula><mml:math id="mm764"><mml:semantics><mml:mrow><mml:msub><mml:mi>T</mml:mi><mml:mi>R</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mn>378</mml:mn></mml:mrow></mml:semantics></mml:math></inline-formula>. The spectra are labelled with LES eddy–eddy parameterizations of no subgrid model (none), anisotropic stochastic (AS), or anisotropic deterministic (AD). The mean subgrid tendency component contains all the terms. Adapted from Kitsios et al. (2023) [<xref ref-type="bibr" rid="B155-atmosphere-15-00921">155</xref>].</p>
Full article ">
16 pages, 9394 KiB  
Article
Analysis of Different Height Correction Models for Tropospheric Delay Grid Products over the Yunnan Mountains
by Fangrong Zhou, Luohong Li, Yifan Wang, Zelin Dai, Chenchen Ding, Hui Li and Yunbin Yuan
Atmosphere 2024, 15(8), 872; https://doi.org/10.3390/atmos15080872 - 23 Jul 2024
Viewed by 415
Abstract
Accurate tropospheric delays are of great importance for both Global Navigation Satellite System (GNSS)-based positioning and precipitable water vapor monitoring. The gridded tropospheric delay products, including zenith hydrostatic delays (ZHD) and zenith wet delays (ZWD), are the most ideal method for accessing accurate [...] Read more.
Accurate tropospheric delays are of great importance for both Global Navigation Satellite System (GNSS)-based positioning and precipitable water vapor monitoring. The gridded tropospheric delay products, including zenith hydrostatic delays (ZHD) and zenith wet delays (ZWD), are the most ideal method for accessing accurate tropospheric delays. The vertical adjustment method is critical for implementing the gridded tropospheric products. In this work, we consider the different models used for grid products and assess their performance over Yunnan mountains with complex topography. We summarize the main results as follows: (1) The products can provide accurate ZHD with mean biases of −2.6 mm and mean Standard Deviation (STD) of 1.5 mm while the ZWD results from grid products show a performance with biases of −0.4 mm and STD of 1.3 cm over the Yunnan area. (2) The Tv-based model shows a better performance than the T0-based model and IGPZWD in rugged areas with large height differences. The grid products can provide hourly ZHD with biases of 3 mm and wet delay with mean biases of within 2 cm and mean STD of below 3 cm in the Yunnan mountains, which exhibit a large height difference of around 1.5 km. (3) The radiosondes results confirm that the Tv-based model has an obvious advantage in calculating ZHD height corrections for differences within 2 km while the T0-model suffers from a loss in accuracy in the case of large height differences. If the site is located more than 1 km below the reference height, the IGPZWD model can provide a better ZWD with a mean bias of 1.5 cm and a mean STD of 1.7 cm. With vertical reduction models, the grid products can provide accurate ZHD and ZWD in real time, even if in complex area. Full article
(This article belongs to the Section Upper Atmosphere)
Show Figures

Figure 1

Figure 1
<p>The distribution of grid points (purple dots), GNSS stations (red triangles), meteorological stations (grey dots), and radiosonde (blue pentagram) over Yunnan and the surrounding area.</p>
Full article ">Figure 2
<p>The difference in ZHD (<b>a</b>) and ZWD (<b>b</b>) at grid point (26.5° N, 99.5° E, 3118.39 m) from the grid products compared with ERA5. The text from left to right labeled in the figure is Bias, RMS, and STD, respectively.</p>
Full article ">Figure 3
<p>Mean Bias and STD of ZHD derived from grid products with respect to ERA5 over the Yunnan mountains.</p>
Full article ">Figure 4
<p>Mean Bias and STD of ZWD derived from products with respect to ERA5 over the Yunnan mountains.</p>
Full article ">Figure 5
<p>The performances in ZHD (mm) from grid products with IGPZWD, <span class="html-italic">T</span><sub>0</sub>-Based model, and <span class="html-italic">T<sub>v</sub></span>-based model at station 56586 (<b>a</b>) and station 48803 (<b>b</b>). The text from left to right labeled in the figure is Bias and STD for each model, respectively. The height differences between the site and grid reference height in meters are listed.</p>
Full article ">Figure 6
<p>Bias and STD for ZHD with IGPZWD and THCN model implemented at KMIN (<b>a</b>) and THCN (<b>b</b>). The listed height differences between site and grid reference height are in meters.</p>
Full article ">Figure 7
<p>The performances in two vertical reduction models for ZWD (cm) in summer (<b>upper</b>) and autumn (<b>lower</b>) at station YNYL.</p>
Full article ">Figure 8
<p>The performance of three vertical reduction models for ZHD (mm) at different height differences at radiosondes station 56778 (25.0° N, 102.3° E), 56964 (22.8° N, 101.0° E), and 56985 (23.4° N, 103.3° E). The zero height refers to the average height of four grid points near the RS sites. The bold text labeled in the figure is the Station Number.</p>
Full article ">Figure 9
<p>The performance of IGPZWD and CWSH models for ZWD (cm) vertical reduction at different height differences at six radiosondes stations. The zero height refers to the average height of four grid points near the radiosondes sites. The bold text labeled in the figure is the Station Number.</p>
Full article ">Figure 10
<p>The ZWD (mm) of IGPZWD, CWSH models, and RS-ZWD with different height differences at radiosondes station 56571. The text from left to right labeled in the figure is Bias and STD for each model, respectively. The height differences between the site and grid reference height in meters are listed.</p>
Full article ">Figure 11
<p>The ZWD errors (mm) after vertical adjustment with the IGPZWD and CWSH models at radiosondes station 56571. The RS-ZWD is used as the reference.</p>
Full article ">
29 pages, 7421 KiB  
Article
Continuous Online Semantic Implicit Representation for Autonomous Ground Robot Navigation in Unstructured Environments
by Quentin Serdel, Julien Marzat and Julien Moras
Robotics 2024, 13(7), 108; https://doi.org/10.3390/robotics13070108 - 18 Jul 2024
Viewed by 780
Abstract
While mobile ground robots have now the physical capacity of travelling in unstructured challenging environments such as extraterrestrial surfaces or devastated terrains, their safe and efficient autonomous navigation has yet to be improved before entrusting them with complex unsupervised missions in such conditions. [...] Read more.
While mobile ground robots have now the physical capacity of travelling in unstructured challenging environments such as extraterrestrial surfaces or devastated terrains, their safe and efficient autonomous navigation has yet to be improved before entrusting them with complex unsupervised missions in such conditions. Recent advances in machine learning applied to semantic scene understanding and environment representations, coupled with modern embedded computational means and sensors hold promising potential in this matter. This paper therefore introduces the combination of semantic understanding, continuous implicit environment representation and smooth informed path-planning in a new method named COSMAu-Nav. It is specifically dedicated to autonomous ground robot navigation in unstructured environments and adaptable for embedded, real-time usage without requiring any form of telecommunication. Data clustering and Gaussian processes are employed to perform online regression of the environment topography, occupancy and terrain traversability from 3D semantic point clouds while providing an uncertainty modeling. The continuous and differentiable properties of Gaussian processes allow gradient based optimisation to be used for smooth local path-planning with respect to the terrain properties. The proposed pipeline has been evaluated and compared with two reference 3D semantic mapping methods in terms of quality of representation under localisation and semantic segmentation uncertainty using a Gazebo simulation, derived from the 3DRMS dataset. Its computational requirements have been evaluated using the Rellis-3D real world dataset. It has been implemented on a real ground robot and successfully employed for its autonomous navigation in a previously unknown outdoor environment. Full article
Show Figures

Figure 1

Figure 1
<p>Example visualisation of the proposed data treatment pipeline producing the simultaneous regression of variables of interest for robot navigation, with an intermediate data compression step.</p>
Full article ">Figure 2
<p>Scheme of the GPR environment representation pipeline. <span class="html-italic">p</span> denotes 3D points coordinates, <span class="html-italic">l</span> semantic labels, <math display="inline"><semantics> <mi>κ</mi> </semantics></math> network confidences, <span class="html-italic">T</span> terrain traversability values, <span class="html-italic">d</span> distances to closest obstacles, <span class="html-italic">s</span> the terrain slope and <math display="inline"><semantics> <mi>σ</mi> </semantics></math> the GPR standard deviation. <math display="inline"><semantics> <mrow> <mi>x</mi> <mo>,</mo> <mi>y</mi> </mrow> </semantics></math> are 2D coordinates. Variables marked with * are products of the compression and fusion process, <math display="inline"><semantics> <msub> <mrow> <mi>H</mi> <mo>|</mo> </mrow> <mrow> <mi>R</mi> <mo>→</mo> <mi>W</mi> </mrow> </msub> </semantics></math> is the transform between the robot sensor frame and the world frame.</p>
Full article ">Figure 3
<p>Scheme of the COSMAu-Nav architecture performing online smooth local path-planning from semantic point clouds and robot localisation. The local planner goal can be given by an external global navigation module such as the one described in <a href="#sec6dot2dot3-robotics-13-00108" class="html-sec">Section 6.2.3</a>.</p>
Full article ">Figure 4
<p>Visualisation of the Gazebo simulation (reproduced from [<a href="#B2-robotics-13-00108" class="html-bibr">2</a>]) showing the 3D environment mesh derived from the 3DRMS dataset, the mobile robot model and its camera module outputs.</p>
Full article ">Figure 5
<p>Comparative evaluation of the clustering compression methods and their parameters on the ratio between GPR environment representation accuracy and computational load.</p>
Full article ">Figure 6
<p>Comparative evaluation of the impact of the process variance threshold <math display="inline"><semantics> <msub> <mi>t</mi> <msup> <mi>σ</mi> <mn>2</mn> </msup> </msub> </semantics></math> on the GPR environment representation quality for both compression methods.</p>
Full article ">Figure 7
<p>Diagram of the baseline SMaNA architecture [<a href="#B2-robotics-13-00108" class="html-bibr">2</a>] with different candidate environment representation methods.</p>
Full article ">Figure 8
<p>Visualisation of the environment representations produced by the compared methods from sensor data gathered along a trajectory inside the described simulation environment with no input data degradation and nominal parameters. Robot trajectory is displayed in white in all 4 figures. Top left displays the reference navigation grid built from the simulator ground-truth. Top right and bottom left display the navigation graphs built by SMaNA respectively using Semantic Octomap and Kimera-Semantics. Bottom left displays the inference of the GPR environment representation on a 2D uniform grid of same resolution.</p>
Full article ">Figure 9
<p>Results of the comparative evaluation of COSMAu-Nav environment representation process with the application of Semantic Octomap and Kimera Semantics in SMaNA in terms of environment representation quality in function of the mapping and inference resolution.</p>
Full article ">Figure 10
<p>Results of the comparative evaluation of COSMAu-Nav environment representation process with the application of Semantic Octomap and Kimera Semantics in SMaNA in terms of environment representation quality in function of the robot localisation noise.</p>
Full article ">Figure 11
<p>Results of the comparative evaluation of COSMAu-Nav environment representation process with the application of Semantic Octomap and Kimera Semantics in SMaNA in terms of environment representation quality in function of the semantic segmentation error level.</p>
Full article ">Figure 12
<p>Results of the comparative evaluation of COSMAu-Nav environment representation process against the integration of Semantic Octomap and Kimera Semantics in SMaNA in terms of computational performances when processing real-time data from the sequence <span class="html-italic">00</span> of the Rellis-3D dataset.</p>
Full article ">Figure 13
<p>Scheme of the sensor data treatment pipeline generating the robot localisation and 3D point clouds enriched with semantic labels from data produced by a 3D LiDAR and an RGB camera mounted on a ground robot.</p>
Full article ">Figure 14
<p>Example data produced by the online pre-treatment pipeline embedded on the robot.</p>
Full article ">Figure 15
<p>Visualisation of COSMAu-Nav local GPR environment representation and its different planning outputs during the robot autonomous navigation mission.</p>
Full article ">Figure 16
<p>Visual result of the presented navigation task outcome. <b>Left</b>: aerial view of the reference map cloud with the trajectory followed by the robot and the given successive way-points. <b>Right</b>: global navigation environment reconstructed by COSMAu-Nav with the recorded robot key-poses.</p>
Full article ">Figure A1
<p>Detailed results of the ablation study presented in <a href="#sec5dot2-robotics-13-00108" class="html-sec">Section 5.2</a> displaying the influence of the compared compression methods parameters over the map quality metrics and computational characteristics.</p>
Full article ">Figure A1 Cont.
<p>Detailed results of the ablation study presented in <a href="#sec5dot2-robotics-13-00108" class="html-sec">Section 5.2</a> displaying the influence of the compared compression methods parameters over the map quality metrics and computational characteristics.</p>
Full article ">
27 pages, 4322 KiB  
Article
Adaptive Filtering with Fitted Noise Estimate (AFFiNE): Blink Artifact Correction in Simulated and Real P300 Data
by Kevin E. Alexander, Justin R. Estepp and Sherif M. Elbasiouny
Bioengineering 2024, 11(7), 707; https://doi.org/10.3390/bioengineering11070707 - 12 Jul 2024
Viewed by 580
Abstract
(1) Background: The electroencephalogram (EEG) is frequently corrupted by ocular artifacts such as saccades and blinks. Methods for correcting these artifacts include independent component analysis (ICA) and recursive-least-squares (RLS) adaptive filtering (-AF). Here, we introduce a new method, AFFiNE, that applies Bayesian adaptive [...] Read more.
(1) Background: The electroencephalogram (EEG) is frequently corrupted by ocular artifacts such as saccades and blinks. Methods for correcting these artifacts include independent component analysis (ICA) and recursive-least-squares (RLS) adaptive filtering (-AF). Here, we introduce a new method, AFFiNE, that applies Bayesian adaptive regression spline (BARS) fitting to the adaptive filter’s reference noise input to address the known limitations of both ICA and RLS-AF, and then compare the performance of all three methods. (2) Methods: Artifact-corrected P300 morphologies, topographies, and measurements were compared between the three methods, and to known truth conditions, where possible, using real and simulated blink-corrupted event-related potential (ERP) datasets. (3) Results: In both simulated and real datasets, AFFiNE was successful at removing the blink artifact while preserving the underlying P300 signal in all situations where RLS-AF failed. Compared to ICA, AFFiNE resulted in either a practically or an observably comparable error. (4) Conclusions: AFFiNE is an ocular artifact correction technique that is implementable in online analyses; it can adapt to being non-stationarity and is independent of channel density and recording duration. AFFiNE can be utilized for the removal of blink artifacts in situations where ICA may not be practically or theoretically useful. Full article
(This article belongs to the Section Biosignal Processing)
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Block diagram of the AFFiNE method. The reference noise input, r<sub>v</sub>[<span class="html-italic">n</span>], is conditioned by BARS to be free of bidirectional contamination from EEG sources close to the recording VEOG electrodes, thus resulting in a more ideal reference noise input signal for RLS-AF. The resulting adaptive filter is a FIR of length M, which operates on the conditioned reference noise input signal and is then linearly subtracted from the recorded EEG signal, s[<span class="html-italic">n</span>], to produce an EEG signal that is free of blink artifacts, e[<span class="html-italic">n</span>].</p>
Full article ">Figure 2
<p>This flowchart graphically outlines the main components of the current study with the relevant section numbers provided in parenthesis. Those sections primarily pertaining to the simulated data are shown in grey, and those pertaining primarily to the real data are in black; sections common to both datasets are shown in white.</p>
Full article ">Figure 3
<p>Stimulus timing sequence for the visual oddball paradigm. Each stimulus appeared on the screen for 100 ms, with inter-stimulus intervals varying on a uniform distribution between 2.0 and 2.5 s. During the inter-stimulus interval, a fixation cross appeared on the screen.</p>
Full article ">Figure 4
<p>Example signals used in the creation of the simulated dataset. (<b>A</b>) Example standard and oddball ERPs (grey and black, respectively) at electrode Pz, and (<b>B</b>) the spatial weights used to project the signal to all simulated electrodes. (<b>C</b>) An example simulated blink signal at FPz, and (<b>D</b>) its spatial weights used to project the signal across the scalp.</p>
Full article ">Figure 5
<p>Example time segment of the simulated signals and the summed (‘Total’) simulated signal at electrode Pz. The total simulated signal consisted of summed ‘Noise’, ‘ERP’, and ‘Blink’ signals. The blink samples in the blink signal were placed such that they coincided with half of the ERP samples.</p>
Full article ">Figure 6
<p>Histogram showing trial counts for each participant contained in the sample of real data used to compare the performance of the correction methods. A total of 17 participants were included with trial counts ranging from 18 (2 participants) to 33 (1 participant).</p>
Full article ">Figure 7
<p>Standard, oddball, and difference wave ERPs calculated from the simulated dataset for all Artifact by Correction Method combinations at FPz and Pz. Because these are simulated data, ‘Truth’ is the actual ERP calculated in each condition using only the simulated cortical source (i.e., task-relevant and -irrelevant) data. Data are shown with y-axis breaks in FPz/Blink Present to visualize the Uncorrected blink artifact maximum amplitude without detracting from the lower-amplitude features of the Truth and Corrected (ICA, RLS-AF, and AFFiNE) waveforms.</p>
Full article ">Figure 8
<p>Topographic MAE from each of the standard, oddball, and difference wave ERPs, both Blink Present and Blink Absent Artifact conditions, and the ICA, RLS, and AFFiNE levels of the Correction Method. The ‘Uncorrected’ level of the Correction Method is omitted here because no correction is reasonably attempted in the Unfiltered condition.</p>
Full article ">Figure 9
<p>Simulated data MAE values at FPz and Pz electrodes over 30 oddball and 30 standard trials in the Blink Present (<b>A</b>) and Blink Absent (<b>B</b>) Artifact conditions. Individual MAE values are plotted as a ‘swarmplot’ to aid in the visualization of the discrete data by offsetting data points in the x-dimension; this offset in the x-dimension is for visualization only and does not encode any information about the distribution of the presented data. Grey circles indicate distribution means.</p>
Full article ">Figure 10
<p>Real data, original and corrected, grand-averaged ERP waveforms at electrodes FPz and Pz. Pz Oddball and difference wave ERP measurement windows are indicated as shaded portions of the x-axis. Data are shown with y-axis breaks in FPz/Blink Present to visualize the Uncorrected blink artifact maximum amplitude without detracting from the lower-amplitude features of the Truth and Corrected (ICA, RLS-AF, and AFFiNE) waveforms.</p>
Full article ">Figure 11
<p>Real data mean participant MAEs. Errors were calculated using the Uncorrected data in the Blink Absent Artifact condition as the “truth data”.</p>
Full article ">Figure 12
<p>Oddball and difference wave ERP amplitude (<b>A</b>) and latency (<b>B</b>) measure absolute errors after performing each correction method on the Blink Absent trials. Individual absolute errors are plotted as a ‘swarmplot’ to aid in the visualization of the discrete data by offsetting data points in the x-dimension; this offset in the x-dimension is for visualization only and does not encode any information about the distribution of the presented data. Grey circles indicate distribution means. Practically speaking, latency error was either zero or one sample point (i.e., 3.9 ms), with rare exceptions, across all Correction Methods.</p>
Full article ">Figure 13
<p>Oddball and Difference wave ERP amplitude (<b>A</b>) and latency (<b>B</b>) measure absolute errors after performing each correction method on the Blink Present trials. Individual absolute errors are plotted as a ‘swarmplot’ to aid in the visualization of the discrete data by offsetting data points in the x-dimension; this offset in the x-dimension is for visualization only and does not encode any information about the distribution of the presented data. Grey circles indicate distribution means.</p>
Full article ">Figure 14
<p>Real data from a single participant; uncorrected and corrected ERP waveforms at electrodes FPz and Pz. Pz Oddball and difference wave ERP measurement windows are indicated as shaded portions of the x-axis. Data are shown with y-axis breaks in FPz/Blink Present to visualize Uncorrected blink artifact maximum amplitudes without detracting from the lower-amplitude features of the Truth and Corrected (ICA, RLS-AF, and AFFiNE) waveforms.</p>
Full article ">Figure 15
<p>Real data from an additional, single participant; uncorrected and corrected ERP waveforms at electrodes FPz and Pz. Pz Oddball and difference wave ERP measurement windows are indicated as shaded portions of the x-axis. Data are shown with y-axis breaks in FPz/Blink Present to visualize Uncorrected blink artifact maximum amplitudes without detracting from the lower-amplitude features of the Truth and Corrected (ICA, RLS-AF, and AFFiNE) waveforms.</p>
Full article ">
17 pages, 1824 KiB  
Systematic Review
The Complex Role Played by the Default Mode Network during Sexual Stimulation: A Cluster-Based fMRI Meta-Analysis
by Joana Pinto, Camila Comprido, Vanessa Moreira, Marica Tina Maccarone, Carlotta Cogoni, Ricardo Faustino, Duarte Pignatelli and Nicoletta Cera
Behav. Sci. 2024, 14(7), 570; https://doi.org/10.3390/bs14070570 - 5 Jul 2024
Viewed by 772
Abstract
The default mode network (DMN) is a complex network that plays a significant and active role during naturalistic stimulation. Previous studies that have used naturalistic stimuli, such as real-life stories or silent or sonorous films, have found that the information processing involved a [...] Read more.
The default mode network (DMN) is a complex network that plays a significant and active role during naturalistic stimulation. Previous studies that have used naturalistic stimuli, such as real-life stories or silent or sonorous films, have found that the information processing involved a complex hierarchical set of brain regions, including the DMN nodes. The DMN is not involved in low-level features and is only associated with high-level content-related incoming information. The human sexual experience involves a complex set of processes related to both external context and inner processes. Since the DMN plays an active role in the integration of naturalistic stimuli and aesthetic perception with beliefs, thoughts, and episodic autobiographical memories, we aimed at quantifying the involvement of the nodes of the DMN during visual sexual stimulation. After a systematic search in the principal electronic databases, we selected 83 fMRI studies, and an ALE meta-analysis was calculated. We performed conjunction analyses to assess differences in the DMN related to stimulus modalities, sex differences, and sexual orientation. The results show that sexual stimulation alters the topography of the DMN and highlights the DMN’s active role in the integration of sexual stimuli with sexual schemas and beliefs. Full article
(This article belongs to the Special Issue Neural Correlates of Cognitive and Affective Processing)
Show Figures

Figure 1

Figure 1
<p>Flowchart of the selection process. ** From Page MJ, et al. [<a href="#B25-behavsci-14-00570" class="html-bibr">25</a>].</p>
Full article ">Figure 2
<p>Characteristics of the included studies. (<b>A</b>). Four doughnut charts of the number of participants for each gender group (men, women, transsexuals, and nonbinary) and each sexual orientation subgroup for men and women (heterosexual, homosexual, bisexual, or healthy) and for transsexuals (transsexual male to female and transsexual female to male) (<b>B</b>). Boxplot of the age distribution of men, women, and both transsexual subgroups (<b>C</b>). Histograms of the age distribution of the male and female participants (<b>D</b>). Doughnut charts related to the number of types of MRI machine magnetic field strengths and tasks used in the studies (<b>E</b>). Bar plot of the number of studies produced by each country. All the graphs were created using MATLAB (vers. 2022b).</p>
Full article ">Figure 3
<p>ALE map of the resulting DMN nodes during sexual stimulation. Maps are overimposed on a 2 × 2 × 2 mm MNI template according to neurological convention. The colored bar denotes the corresponding ALE value ranges indicated on the maps (FDR pN &lt; 0.01).</p>
Full article ">Figure 4
<p>Conjunction ALE map results. Activation likelihood maps resulting from the contrast between videos and pictures (TOP), between men and women (bottom left), and between heterosexual and homosexual subjects (bottom right). The contrast map “video vs. pictures” (FDR pN &lt; 0.01) showed convergent activity in the ventral parahippocampal regions of the DMN. The two contrast maps at the bottom (<span class="html-italic">p</span> &lt; 0.001) show the involvement of both subcortical and cortical DMN nodes in sex differences and sexual preference. The maps are superimposed on an MNI 2 × 2 × 2 mm template according to neurological convention. The colored bars under each map indicate the ALE value ranges.</p>
Full article ">
10 pages, 2856 KiB  
Article
Real-Time Plasmonic Strain Sensors Based on Surface Relief Diffraction Gratings
by Yazan Bdour and Ribal Georges Sabat
Micromachines 2024, 15(7), 863; https://doi.org/10.3390/mi15070863 - 30 Jun 2024
Viewed by 644
Abstract
Large-scale diffraction gratings were fabricated in surface relief on azobenzene thin films and transferred to flexible PDMS substrates using soft lift-off lithography. The PDMS gratings were strained along the grating vector axis and the resulting surface topography was analyzed using diffraction angle measurements, [...] Read more.
Large-scale diffraction gratings were fabricated in surface relief on azobenzene thin films and transferred to flexible PDMS substrates using soft lift-off lithography. The PDMS gratings were strained along the grating vector axis and the resulting surface topography was analyzed using diffraction angle measurements, AFM imagery and surface plasmon resonance (SPR) spectra. All measurement methods exhibited a linear response in strain indicating the useability of these sensors in real-world applications. For SPR-based strain sensing, an increasing pitch and a decreasing modulation depth were observed with increasing strain. The SPR peak shifted by ~1.0 nm wavelength and the SPR intensity decreased by ~0.3 a.u. per percentage of applied strain. The tested PDMS samples retained their integrity even after multiple cycles of stretching and relaxation, making them a suitable strain sensor. Full article
(This article belongs to the Special Issue Plasmonic Sensors and Their Applications)
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Image the in-house-built fixture used to apply linear strain to the nanostructured PDMS sample. The fixture features a side-mounted knob for controlled movement of the clamps, enabling precise incremental stretching. A linear metric scale with a Vernier scale is used to accurately measure the separation distance as the sample is strained. Additionally, the fixture includes an optical mount handle on the opposite side of the knob, facilitating easier setup during experiments.</p>
Full article ">Figure 2
<p>Plot of the change in pitch of the gratings as the PDMS sample undergoes strain and relaxation obtained optically through the change in the non-zero diffraction orders angles. During stretching, the pitch increases by 3.16 nm per percentage of strain applied, while upon relaxation, the pitch alteration is observed to be 3.38 nm per percentage of strain relieved.</p>
Full article ">Figure 3
<p>(<b>a</b>) Morphological analysis of the gratings using AFM, displaying the average pitch and modulation depth at varying levels of applied strain. (<b>b</b>) AFM image obtained at approximately 15% applied strain, illustrating the topographical image obtained of the gratings under mechanical stress.</p>
Full article ">Figure 4
<p>(<b>a</b>) Normalized transmission spectra illustrating the plasmonic response of a 40 nm Au-coated PDMS grating as strain increases. The spectra highlight the shift in the SPR peak with applied strain. (<b>b</b>) Corresponding shift in SPR wavelength peak as the sample undergoes both strain and relaxation, indicating the sensitivity of the gratings to mechanical stress.</p>
Full article ">Figure 5
<p>(<b>a</b>) Graph illustrating the correlation between the modulation depth of the grating and the intensity of the SPR peak. As the grating depth decreases under strain, the SPR peak intensity exhibits a corresponding reduction. (<b>b</b>) Plot demonstrating the decrease in intensity of the SPR peak as a function of strain.</p>
Full article ">
24 pages, 9897 KiB  
Review
Surface Deformation of Biocompatible Materials: Recent Advances in Biological Applications
by Sunhee Yoon, Ahmed Fuwad, Seorin Jeong, Hyeran Cho, Tae-Joon Jeon and Sun Min Kim
Biomimetics 2024, 9(7), 395; https://doi.org/10.3390/biomimetics9070395 - 28 Jun 2024
Cited by 1 | Viewed by 659
Abstract
The surface topography of substrates is a crucial factor that determines the interaction with biological materials in bioengineering research. Therefore, it is important to appropriately modify the surface topography according to the research purpose. Surface topography can be fabricated in various forms, such [...] Read more.
The surface topography of substrates is a crucial factor that determines the interaction with biological materials in bioengineering research. Therefore, it is important to appropriately modify the surface topography according to the research purpose. Surface topography can be fabricated in various forms, such as wrinkles, creases, and ridges using surface deformation techniques, which can contribute to the performance enhancement of cell chips, organ chips, and biosensors. This review provides a comprehensive overview of the characteristics of soft, hard, and hybrid substrates used in the bioengineering field and the surface deformation techniques applied to the substrates. Furthermore, this review summarizes the cases of cell-based research and other applications, such as biosensor research, that utilize surface deformation techniques. In cell-based research, various studies have reported optimized cell behavior and differentiation through surface deformation, while, in the biosensor and biofilm fields, performance improvement cases due to surface deformation have been reported. Through these studies, we confirm the contribution of surface deformation techniques to the advancement of the bioengineering field. In the future, it is expected that the application of surface deformation techniques to the real-time interaction analysis between biological materials and dynamically deformable substrates will increase the utilization and importance of these techniques in various fields, including cell research and biosensors. Full article
(This article belongs to the Special Issue Dynamical Response of Biological System and Biomaterial 2024)
Show Figures

Figure 1

Figure 1
<p>Types of various surface deformations. (<b>A</b>) A phase diagram for predicting the formation of various deformation patterns. The circular markers indicate fold condition, the square marker indicates wrinkle condition, and the triangular markers indicate crease condition. Reprinted from [<a href="#B113-biomimetics-09-00395" class="html-bibr">113</a>] under the Creative Commons CC BY-NC-ND license, (<b>B</b>) examples of the different surface deformations discussed in this review. Since double (or period double) is rarely used in the biological applications covered in this review, it is considered a type of fold and is not described separately.</p>
Full article ">Figure 2
<p>Wrinkling formation by stretching and releasing. (<b>A</b>) Process, (<b>B</b>) real wrinkling surface. Reprinted with permission from [<a href="#B118-biomimetics-09-00395" class="html-bibr">118</a>].</p>
Full article ">Figure 3
<p>Creasing formation by swelling. (<b>A</b>) Process, (<b>B</b>) real creasing surface. Reprinted with permission from [<a href="#B122-biomimetics-09-00395" class="html-bibr">122</a>].</p>
Full article ">Figure 4
<p>Folding and ridge formation by compression based on the modulus ratio (μ<sub>f</sub>/μ<sub>s</sub>). (<b>A</b>) Process, (<b>B</b>) real fold surface (The red arrow indicates the valley, a characteristic point that marks the bifurcation of the wrinkle and fold), (<b>C</b>) real ridge surface. Reprinted from [<a href="#B22-biomimetics-09-00395" class="html-bibr">22</a>] under the Creative Commons CC-BY license.</p>
Full article ">Figure 5
<p>Image (<b>A</b>) shows that the hBM-MSCs cultured on the winkled surface could be stimulated into osteogenic differentiation. Image (<b>B</b>) shows the analysis of different wrinkle sizes on the morphology of the hBM-MSCs. The microscopic images illustrate the cells cultured on the different surfaces with different surface topography. Cell cytoskeleton and nucleus were stained using TRITC-labeled phalloidin (red) and DAPI (blue). The white arrow shows the direction of the wrinkles (Scale bar for all images is 100 μm). Image (<b>C</b>) shows the analysis of cell orientation, cell area, and cell aspect ratio (CAR) for different surface topographies (Data are shown as mean ± standard deviation (SD), and * <span class="html-italic">p</span> &lt; 0.05, ** <span class="html-italic">p</span> &lt; 0.01, *** <span class="html-italic">p</span> &lt; 0.001). Reprinted from [<a href="#B141-biomimetics-09-00395" class="html-bibr">141</a>] under the Creative Commons CC-BY license.</p>
Full article ">Figure 6
<p>Image (<b>A</b>) shows the fabrication of a nanoscale ridge pattern array using UV-assisted capillary force lithography. Image (<b>B</b>) shows the SEM image of the PUA ridge-patterned surface. Image (<b>C</b>) shows a cross-sectional SEM image of the substrate. Image (<b>D</b>) shows the neural differentiation of hESCs cultured on the flat surface and on the ridge surface. The cells were immunostained with DAPI, tuj1, brachyury, PDX1, and nestin. Images (<b>a</b>–<b>c</b>) show the cell differentiation on flat surface for a five-day culture, images (<b>d</b>–<b>f</b>) show the cell culture on the ridge surface over 5 days, and images (<b>g</b>–<b>i</b>) show the culture over 10 days on the ridge-patterned surface (Scale bar for all images is 100 μm). Reprinted with permission from [<a href="#B142-biomimetics-09-00395" class="html-bibr">142</a>].</p>
Full article ">Figure 7
<p>Images (<b>A</b>–<b>C</b>) show representative parallel and vertical topographies together with the flat control (no pattern on the PUA surface). Images (<b>D</b>–<b>F</b>) show time-lapse microscopic images of in vitro wound healing with respect to the orientation and densities of ridge/groove topography. The yellow double arrow indicate orientation of the patterns. Image (<b>G</b>) shows fibroblast coverage, which was quantified for each topographical orientation and density. The results indicated that the perpendicular topography exhibited significantly higher coverage rates in comparison to the smooth control surface and the topography with parallel alignment. Reprinted with permission from [<a href="#B18-biomimetics-09-00395" class="html-bibr">18</a>].</p>
Full article ">Figure 8
<p>Image (<b>A</b>) shows the schematic process of fabricating the winkled gradient device using the plasma oxidation process and then air exposure under ambient conditions. (<b>B</b>) Focal length images of the cells after two days of culture (Scale bars = 22 μm). (<b>C</b>) AFM images of surface topography with different metal oxides (Scale bars = 4 μm). (<b>D</b>) Fluorescence images of the cells along the full length of the device (Scale bars = 1 mm). Reprinted from [<a href="#B147-biomimetics-09-00395" class="html-bibr">147</a>] under the Creative Commons CC-BY-NC-ND license.</p>
Full article ">Figure 9
<p>Image (<b>A</b>) shows the schematic illustration of hierarchically wrinkled mesoporous ceramics platform fabrication using sol-gel reaction for cell alignment. The sol-gel precursor is coated on the cleaned PET substrate and, subsequently, PDMS casting is performed; finally, the dynamic interference release and flipping of PDMS/SiO<sub>2</sub> generates the hierarchically wrinkled surface. Image (<b>B</b>) shows the optical images of the surface synthesized using sol-gel precursor mixtures with different weights of MTES:TEOS. Image (<b>C</b>) shows cell images and analysis results. (<b>a</b>) Optical images of blood vessel cells on the wrinkled surface. The white arrows indicate wrinkle valley. (<b>b</b>,<b>c</b>) The fluorescence images show the comparison of cell alignment on the flat surface and the wrinkled surface. (<b>d</b>) Illustration of cell growth angles (θ°) and wrinkle orientation. (<b>e</b>) The angular histogram of cells on the surface was measured. (<b>f</b>) The FWHM of cell alignment distribution. Reprinted with permission from [<a href="#B148-biomimetics-09-00395" class="html-bibr">148</a>].</p>
Full article ">Figure 10
<p>(<b>A</b>) WST-1 assay analysis showing the cell metabolic activity on the biointerfaces compared to the control surface; after proliferation over 6 days, it increased by 25%. Reprinted with permission from [<a href="#B150-biomimetics-09-00395" class="html-bibr">150</a>]. (<b>B</b>) Effect of the wrinkled surface on cell alignment. The cells were analyzed with phalloidin (cytoskeleton) and DAPI (nuclear) assays, with red showing the cytoskeleton and blue showing the nuclear material. The figures show that cells on the wrinkled surfaces (<b>b</b>) exhibited alignment compared to the cells on the flat surface (<b>a</b>). The white arrow indicates wrinkle orientation (Scale bar = 200 μm). Reprinted with permission from [<a href="#B151-biomimetics-09-00395" class="html-bibr">151</a>]. (<b>C</b>) Alignment of hiPSC-CM cells on the static SMP-PEM substrate. (<b>a</b>,<b>b</b>) The figures show the AFM image of flat and wrinkled surface(Scale bars = 5 μm). (<b>c</b>,<b>d</b>) The figures show SEM images of a wrinkled surface and cell (Scale bars = 10 μm). (<b>e</b>,<b>f</b>) The fluorescent images show randomly orientated cells and highly aligned cell on a flat surface and wrinkled surface, respectively (Scale bars = 50 μm). The white arrow indicates the wrinkle direction. (<b>g</b>,<b>h</b>) Further, the images also show the well-aligned myofibril structure on the wrinkled surface. The white arrow indicates wrinkle direction (Scale bars = 10 μm). Reprinted with permission from [<a href="#B152-biomimetics-09-00395" class="html-bibr">152</a>].</p>
Full article ">Figure 11
<p>(<b>A</b>) Schematic figures for the fabrication of Ag/CNTs-PDMS film sensor. A thin CNT film is transferred to PDMS, and the pre-strain and release strain generate the wrinkle structure. Subsequently, the Ag film sputters and the electrodes are assembled for flexible and wearable sensor. (<b>B</b>) SEM image showing the morphology of a wrinkled structure according to the strain strength. (<b>a</b>–<b>e</b>) Original CNTs, (<b>f</b>–<b>j</b>) Ag/CNTs-PDMS composite films at tensile strain 0%, 5%, 10%, 5%, and 0%, respectively. (<b>C</b>) Film sensor monitoring of weak human body signals. (<b>a</b>) Sensor’s location overview. Signals from (<b>b</b>) finger bending, (<b>c</b>) the wrist, (<b>d</b>) the upper lip, and (<b>e</b>) the chest. Reprinted from [<a href="#B153-biomimetics-09-00395" class="html-bibr">153</a>] under the Creative Commons CC-BY license.</p>
Full article ">Figure 12
<p>(<b>A</b>) Schematics for the fabrication of a nanoparticle array at the top of the ridged crack. Image (<b>a</b>) shows ordered nanoparticles at the tip of wrinkled structure. Image (<b>b</b>) shows electromagnetic field simulation of two nanoparticles. Image (<b>c</b>) shows manufacturing process. A silver nanofilm is sputtered on the pre-stretched PDMS. By releasing the films, triangle cross-section wrinkles are formed. Colloid drops are located on the substrate and only gold nanoparticles remain after evaporation. (<b>d</b>) AFM and (<b>e</b>) SEM images show the aligned nanoparticles at the top of the wrinkled structure. (<b>B</b>) The gap between two nanoparticles can be tuned by (<b>a</b>) mechanical stretch and release force. (<b>b</b>) XFDTD simulations show that the electromagnetic effect is controlled by the nanogap. Image (<b>c</b>) shows SEM image of a nano-gap under different stretching forces. (<b>d</b>) The nanogap is determined using applied stress and SERS intensity is maximized at 2% stretching condition. (<b>C</b>) Electromagnetic field is maximized (<b>a</b>) between gold nanoparticle and silver substrate. (<b>b</b>) The SERS effect enhanced the wrinkled structure under the nanoparticle condition. Reprinted with permission from [<a href="#B156-biomimetics-09-00395" class="html-bibr">156</a>].</p>
Full article ">Figure 13
<p>(<b>A</b>) Schematics for the fabrication of a micro pattern on the PDMS substrate. (<b>a</b>) Stretching, UV/ozone treatment, and the releasing process can generate surface deformation patterns. (<b>b</b>,<b>c</b>) The SEM image shows a constant wrinkle pattern. (<b>B</b>) Water contact angle measured under various surface condition. (<b>a</b>) On the flat surface, the contact angle is 110° and, (<b>b</b>,<b>c</b>) as the stress increases, the contact angle also increases. (<b>d</b>) Finally, in the cassie state, the contact angle reaches close to 180°. The blue line indicates increased height of wrinkled structure. (<b>C</b>) Black points show the experimental results of surface deformation corresponding to the theoretical graphs. Reprinted from [<a href="#B160-biomimetics-09-00395" class="html-bibr">160</a>] under the Creative Commons CC-BY license.</p>
Full article ">Figure 14
<p>(<b>A</b>) Schematics for the fabrication of wrinkled PDMS. Image (<b>a</b>) shows plasma treatment without applied strain. In contrast, to form a uniform wrinkle structure, (<b>b</b>) a20% pre-stretch, oxygen plasma, and releasing sequence were applied. (<b>c</b>) Schematic image of highly controlled and regular buckling of the surface. Image (<b>d</b>) shows strain-induced dynamic wrinkles. (<b>e</b>) The SEM image shows that the wrinkle patterns are well formed on the PDMS surface. (<b>B</b>) Fluorescence microscopy image showing (<b>a</b>) the attached <span class="html-italic">P. aeruginosa</span> PA-14 bacteria forming a pattern at the valley of the wrinkle, (<b>b</b>) while bacteria appear to be randomly distributed on the flat PDMS. (<b>c</b>) SEM image of bacteria on the surface features. (<b>C</b>) Compared to a static culture, a dynamic vibration condition drastically inhibits bacterial attachment. However, introduction of dwell time weakens the effect of the vibration. (<b>a</b>) Continuous strain condition, with no dwell time. (<b>b</b>–<b>d</b>) Introduction of dwell time of 5, 18, 60 min between the strain cycles. Reprinted from [<a href="#B163-biomimetics-09-00395" class="html-bibr">163</a>] under the Creative Commons CC-BY license.</p>
Full article ">
14 pages, 29541 KiB  
Article
A Type of Scale-Oriented Terrain Pattern Derived from Normalized Topographic Relief Layers and Its Interpretation
by Xi Nan, Ainong Li, Zhengwei He and Jinhu Bian
ISPRS Int. J. Geo-Inf. 2024, 13(6), 209; https://doi.org/10.3390/ijgi13060209 - 17 Jun 2024
Viewed by 629
Abstract
Topographic scale characteristics contain valuable information for interpreting landform structures, which is crucial for understanding the spatial differentiation of landforms across large areas. However, the absence of parameters that specifically describe the topographic scale characteristics hinders the quantitative representation of regional topography from [...] Read more.
Topographic scale characteristics contain valuable information for interpreting landform structures, which is crucial for understanding the spatial differentiation of landforms across large areas. However, the absence of parameters that specifically describe the topographic scale characteristics hinders the quantitative representation of regional topography from the perspective of spatial scales. In this study, false-color composite images were generated using normalized topographic relief data, showing a type of scale-oriented terrain pattern. Subsequent analysis indicated a direct correlation between the luminance of the patterns and the normalized topographic relief. Additionally, a linear correlation exists between the color of the patterns and the change rate in normalized topographic relief. Based on the analysis results, the issue of characterizing topographic scale effects was transformed into a problem of interpreting terrain patterns. The introduction of two parameters, flux and curl of topographic field, allowed for the interpretation of the terrain patterns. The assessment indicated that the calculated values of topographic field flux are equivalent to the luminance of the terrain patterns and the variations in the topographic field curl correspond with the spatial differentiation of colors in the terrain patterns. This study introduced a new approach to analyzing topographic scale characteristics, providing a pathway for quantitatively describing scale effects and automatically classifying landforms at a regional scale. Through exploratory analysis on artificially constructed simple DEMs and verification in four typical geomorphological regions of real terrain, it was shown that the terrain pattern method has better intuitiveness than the scale signature approach. It can reflect the scale characteristics of terrain in continuous space. Compared to the MTPCC image, the terrain parameters derived from the terrain pattern method further quantitatively describe the scale effects of the terrain. Full article
Show Figures

Figure 1

Figure 1
<p>CIE Lab color mode. The ‘L’ channel represents luminance. The ‘a’ channel represents the range from red to green and the ‘b’ channel represents the range from yellow to blue.</p>
Full article ">Figure 2
<p>Sample data. (<b>a</b>) Three-dimensional view of the data. (<b>b</b>) Contours with topographic feature lines on the DEM.</p>
Full article ">Figure 3
<p>The DEM of validation area covers the southeastern part of the Qinghai–Tibet Plateau, Hengduan Mts., Sichuan Basin, and the hilly regions in Fujian and Zhejiang Provinces, China, exhibiting diverse landforms.</p>
Full article ">Figure 4
<p>False-color images with different normalized topographic relief layers. (<b>a</b>) Composite image created by blending R5, R7, and R9, with ridge and saddle lines overlayed. (<b>b</b>) Composite image of R5, R9, and R13, with the contours overlayed as a reference of topography and the contours acting as a reference for the terrain’s morphology, corresponding to the pattern colors. (<b>c</b>) Composite image of R5, R11, and R17, with six points marked off for further analysis.</p>
Full article ">Figure 5
<p>Image showing terrain patterns of the actual land surface, created using R9, R19, and R29. (<b>a</b>) Stretching branch of Gangdise, fault-block mountains; (<b>b</b>) Qionglai Mts., fold mountains; (<b>c</b>) Karst landform area (southern part); and (<b>d</b>) southeast hills in China, granite hilly area.</p>
Full article ">Figure 6
<p>Luminance of the patterns. The boxplots span the 25th to 75th percentiles. (<b>a</b>) Violin plots of the pattern luminance for three sections in the sample DEM, with the red dots in the plot indicating mean values. (<b>b</b>) Correlation between luminance and normalized topographic relief.</p>
Full article ">Figure 6 Cont.
<p>Luminance of the patterns. The boxplots span the 25th to 75th percentiles. (<b>a</b>) Violin plots of the pattern luminance for three sections in the sample DEM, with the red dots in the plot indicating mean values. (<b>b</b>) Correlation between luminance and normalized topographic relief.</p>
Full article ">Figure 7
<p>Scale signature of the normalized topographic relief at observation sites.</p>
Full article ">Figure 8
<p>(<b>a</b>) Color vectors of the grid cell in <a href="#ijgi-13-00209-f004" class="html-fig">Figure 4</a>b. (<b>b</b>) Correlation between the change rate of the normalized relief and the rotation angle of the color vectors (only blue-hued grid cells selected). (<b>c</b>) Correlation between the change rate of the normalized relief and the rotation angle of the color vectors (only red-hued grid cells selected).</p>
Full article ">Figure 9
<p>Comparison between the luminance of the pattern (L) and the flux of the topographic field <math display="inline"><semantics> <mrow> <mover accent="true"> <mi>φ</mi> <mo stretchy="true">¯</mo> </mover> </mrow> </semantics></math>. (<b>a</b>) Luminance of the sample image. (<b>b</b>) Topographic flux in the sample area. (<b>c</b>) Scatter plot showing L and <math display="inline"><semantics> <mrow> <mover accent="true"> <mi>φ</mi> <mo stretchy="true">¯</mo> </mover> </mrow> </semantics></math>.</p>
Full article ">Figure 10
<p>Comparison between <span class="html-italic">β</span><sub>curl</sub>, the characteristic angle of the topographic field curl, and <span class="html-italic">β</span><sub>color</sub>, the characteristic angle of the pattern color vector. (<b>a</b>,<b>a-1</b>–<b>a-4</b>) Overall distribution of <span class="html-italic">β</span><sub>color</sub> in the study area and local magnifications in four geomorphic sample areas. (<b>b-1</b>–<b>b-4</b>) Distribution of <span class="html-italic">β</span><sub>color</sub> in four geomorphic sample areas. (<b>c-1</b>–<b>c-4</b>) Histogram of <span class="html-italic">β</span><sub>curl</sub>, with the red dashed lines indicating the positions of the extreme values. (<b>d-1</b>–<b>d-4</b>) Histogram of <span class="html-italic">β</span><sub>color</sub>, with the red dashed lines indicating the positions of the extreme values.</p>
Full article ">Figure 11
<p>Topographic field flux in different geomorphic areas. (<b>a</b>) Sub-area i is dominated by the Qionglai Mts., including parts of the Longmen Mts., with large undulating mountains and high mountains as the main features. Sub-area ii is dominated by small undulating low mountains. Sub-area iii is the Chengdu Plain. (<b>b</b>) Topographic field flux for the three sub-areas in the window size range of 9 to 29 pixels, with an observation baseline scale of 5 km.</p>
Full article ">Figure 12
<p>Changes in the topographic field curl of the Gangdise Mts. branch. (<b>a</b>) Spatial distribution of <span class="html-italic">β</span><sub>curl</sub>; both sub-areas i and ii are typical distribution of fault-block mountains. (<b>b</b>) The numerical distribution of <span class="html-italic">β</span><sub>curl</sub> in sub-areas i and ii, with the red dashed lines indicating the positions of the extreme values.</p>
Full article ">
18 pages, 6597 KiB  
Article
A Performance Comparison of 3D Survey Instruments for Their Application in the Cultural Heritage Field
by Irene Lunghi, Emma Vannini, Alice Dal Fovo, Valentina Di Sarno, Alessandra Rocco and Raffaella Fontana
Sensors 2024, 24(12), 3876; https://doi.org/10.3390/s24123876 - 15 Jun 2024
Viewed by 651
Abstract
Thanks to the recent development of innovative instruments and software with high accuracy and resolution, 3D modelling provides useful insights in several sectors (from industrial metrology to cultural heritage). Moreover, the 3D reconstruction of objects of artistic interest is becoming mandatory, not only [...] Read more.
Thanks to the recent development of innovative instruments and software with high accuracy and resolution, 3D modelling provides useful insights in several sectors (from industrial metrology to cultural heritage). Moreover, the 3D reconstruction of objects of artistic interest is becoming mandatory, not only because of the risks to which works of art are increasingly exposed (e.g., wars and climatic disasters) but also because of the leading role that the virtual fruition of art is taking. In this work, we compared the performance of four 3D instruments based on different working principles and techniques (laser micro-profilometry, structured-light topography and the phase-shifting method) by measuring four samples of different sizes, dimensions and surface characteristics. We aimed to assess the capabilities and limitations of these instruments to verify their accuracy and the technical specifications given in the suppliers’ data sheets. To this end, we calculated the point densities and extracted several profiles from the models to evaluate both their lateral (XY) and axial (Z) resolution. A comparison between the nominal resolution values and those calculated on samples representative of cultural artefacts was used to predict the performance of the instruments in real case studies. Overall, the purpose of this comparison is to provide a quantitative assessment of the performance of the instruments that allows for their correct application to works of art according to their specific characteristics. Full article
(This article belongs to the Special Issue Stereo Vision Sensing and Image Processing)
Show Figures

Figure 1

Figure 1
<p>Samples used for 3D reconstruction: (<b>a</b>) acrylic paint mock-up: for the ten acrylic colours on the white base, the five areas of increasing thickness are numbered from 1 to 5 and delimited by vertical dotted lines, (<b>b</b>) canvases 1 and 3, (<b>c</b>) canvas 2 and (<b>d</b>) plaster statue.</p>
Full article ">Figure 2
<p>Sample canvas 2: 3D outcomes generated by (<b>a</b>) MP with 50 mm focal length and 20 µm sampling step; (<b>b</b>) MICRON3D scanner; (<b>c</b>) EinScan scanner; (<b>d</b>) FARO Focus scanner.</p>
Full article ">Figure 3
<p>Point cloud details of sample canvas 2 generated by (<b>a</b>) 100 mm focal length and 200 µm sampling step (24 ± 2 pp/mm<sup>2</sup>); (<b>b</b>) EinScan scanner (27 ± 1 pp/mm<sup>2</sup>); (<b>c</b>) FARO Focus scanner (30 ± 3 pp/mm<sup>2</sup>); (<b>d</b>) 50 mm focal length and 100 µm sampling step (95 ± 7 pp/mm<sup>2</sup>); (<b>e</b>) MICRON3D scanner (145 ± 10 pp/mm<sup>2</sup>); (<b>f</b>) 50 mm focal length and 50 µm sampling step (375 ± 25 pp/mm<sup>2</sup>); (<b>g</b>) 50 mm focal length and 20 µm sampling step (2448 ± 178 pp/mm<sup>2</sup>).</p>
Full article ">Figure 4
<p>Details of the meshes of sample canvas 2 generated by micro-profilometry with 50 mm focal length and (<b>a</b>) 20 µm and (<b>b</b>) 50 µm sampling steps and 100 µm sampling steps and (<b>c</b>) 50 mm, (<b>d</b>) 100 mm and (<b>e</b>) 200 mm focal lengths.</p>
Full article ">Figure 4 Cont.
<p>Details of the meshes of sample canvas 2 generated by micro-profilometry with 50 mm focal length and (<b>a</b>) 20 µm and (<b>b</b>) 50 µm sampling steps and 100 µm sampling steps and (<b>c</b>) 50 mm, (<b>d</b>) 100 mm and (<b>e</b>) 200 mm focal lengths.</p>
Full article ">Figure 5
<p>(<b>a</b>) Plaster statue, with the dotted arrow indicating the acquisition direction of the profiles shown in (<b>b</b>) and in detail in (<b>c</b>).</p>
Full article ">Figure 6
<p>(<b>a</b>–<b>c</b>) Detail of canvases 1, 2 and 3, respectively. Top: colour map, with the arrow indicating the profile direction, and simulated raking light image for MP_f50_s20; bottom: profiles as measured by the different instruments.</p>
Full article ">Figure 7
<p>(<b>a</b>) Acrylic paint mock-up with the red rectangle highlighting the analysed area (light blue pigment). Colour map of (<b>b</b>) the MP_f50_s50 point cloud, (<b>c</b>) the selected area on the raw point cloud and (<b>d</b>) its conditioned surface. The white dotted arrows in (<b>c</b>,<b>d</b>) show the direction along which the profiles were extracted, which are plotted in blue and green, respectively (<b>e</b>). The five areas of increasing thickness are numbered from 1 to 5 in (<b>a</b>,<b>d</b>,<b>e</b>).</p>
Full article ">Figure 7 Cont.
<p>(<b>a</b>) Acrylic paint mock-up with the red rectangle highlighting the analysed area (light blue pigment). Colour map of (<b>b</b>) the MP_f50_s50 point cloud, (<b>c</b>) the selected area on the raw point cloud and (<b>d</b>) its conditioned surface. The white dotted arrows in (<b>c</b>,<b>d</b>) show the direction along which the profiles were extracted, which are plotted in blue and green, respectively (<b>e</b>). The five areas of increasing thickness are numbered from 1 to 5 in (<b>a</b>,<b>d</b>,<b>e</b>).</p>
Full article ">Figure 8
<p>Conditioned surfaces of (<b>a</b>) MP_f50_s20, (<b>b</b>) MP_f200_s100, (<b>e</b>) EinScan, (<b>f</b>) MICRON3D and (<b>g</b>) MP_f200_s100 point clouds; (<b>c</b>) raw colour map of FARO Focus point cloud. The dotted arrows in (<b>a</b>), (<b>b</b>,<b>c</b>) indicate the direction of the extracted profiles plotted in (<b>d</b>), while the white dotted arrows in (<b>a</b>,<b>e</b>–<b>g</b>) indicate the direction of the extracted profiles plotted in (<b>h</b>).</p>
Full article ">Figure 8 Cont.
<p>Conditioned surfaces of (<b>a</b>) MP_f50_s20, (<b>b</b>) MP_f200_s100, (<b>e</b>) EinScan, (<b>f</b>) MICRON3D and (<b>g</b>) MP_f200_s100 point clouds; (<b>c</b>) raw colour map of FARO Focus point cloud. The dotted arrows in (<b>a</b>), (<b>b</b>,<b>c</b>) indicate the direction of the extracted profiles plotted in (<b>d</b>), while the white dotted arrows in (<b>a</b>,<b>e</b>–<b>g</b>) indicate the direction of the extracted profiles plotted in (<b>h</b>).</p>
Full article ">Figure 9
<p>(<b>a</b>) Conditioned surface of the MP_f50_20 cloud with the areas selected for the thickness computation highlighted in light green on the pigment and in purple at the sides on the preparation layer. (<b>b</b>) Mean thickness and standard deviation for each painted area, shown on the <span class="html-italic">x</span>-axis in a colour scale.</p>
Full article ">
22 pages, 6626 KiB  
Article
An Efficient GPS Algorithm for Maximizing Electric Vehicle Range
by Karim Aboelsoud, Hatem Y. Diab, Mahmoud Abdelsalam and Moutaz M. Hegaze
Appl. Sci. 2024, 14(11), 4858; https://doi.org/10.3390/app14114858 - 4 Jun 2024
Viewed by 671
Abstract
Although the main purpose of conventional geographical positioning systems (GPSs) is to determine either the fastest path or the shortest distance to a destination, this function may not be enough for electric vehicles (EVs). This is simply because the fastest/shortest path may consume [...] Read more.
Although the main purpose of conventional geographical positioning systems (GPSs) is to determine either the fastest path or the shortest distance to a destination, this function may not be enough for electric vehicles (EVs). This is simply because the fastest/shortest path may consume relatively higher energy when compared to other paths depending on the nature, speed limit, and topography of the road. This means that the driving range of the EV per charge decreases dramatically. This paper aims to develop a new algorithm and model dedicated for EV GPS which not only selects shortest/fastest routes, but also focuses on the most energy efficient route. This is achieved by considering many factors including aerodynamics, wind speed, topology of roads, with a clear objective of reducing the energy consumed from the battery. A MATLAB Simulink model is developed and validated with real-life case studies to ensure the results are realistic and accurate. Full article
Show Figures

Figure 1

Figure 1
<p>Global EV and ICE market share forecasting as (%).</p>
Full article ">Figure 2
<p>EV and ICE journey time plus refueling/recharging time comparison.</p>
Full article ">Figure 3
<p>(<b>a</b>) Longitudinal motion and resistances of the vehicle on a flat terrain. (<b>b</b>) Longitudinal motion and resistances of the vehicle on an inclined road (θ is the gradient under consideration).</p>
Full article ">Figure 4
<p>Flowchart for comparing two suggested routes and their energy calculation method.</p>
Full article ">Figure 5
<p>MATLAB Simulink model interface [<a href="#B37-applsci-14-04858" class="html-bibr">37</a>].</p>
Full article ">Figure 6
<p>The FTP75 Driving Cycle (speed vs. time) used for Scenario 1.</p>
Full article ">Figure 7
<p>(<b>a</b>) SOC% for Scenario 1 using different vehicle models. (<b>b</b>) Driving Cycle used in Scenario 2 for a constant 120 km/h velocity.</p>
Full article ">Figure 8
<p>Driving Cycle used in Case 2.</p>
Full article ">Figure 9
<p>SOC (%) difference for three scenarios in Case 2 at different wind speeds and directions.</p>
Full article ">Figure 10
<p>The routes suggested by Google Maps with distance and time.</p>
Full article ">Figure 11
<p>Different topologies of routes in Case 3 [<a href="#B46-applsci-14-04858" class="html-bibr">46</a>].</p>
Full article ">Figure 12
<p>(<b>a</b>) Route 1 Driving Cycle for ID3 Volkswagen. (<b>b</b>) SOC% for the energy consumed in Route 1.</p>
Full article ">Figure 13
<p>(<b>a</b>) Route 2 Driving Cycle for ID3 Volkswagen. (<b>b</b>) SOC% for the energy consumed in Route 2.</p>
Full article ">Figure 14
<p>Google Maps suggestion for Case 4.</p>
Full article ">Figure 15
<p>(<b>a</b>) Driving Cycle for Case 4. (<b>b</b>) Inclinations in the route selected for Case 4 [<a href="#B46-applsci-14-04858" class="html-bibr">46</a>].</p>
Full article ">Figure 16
<p>(<b>a</b>) SOC% for the energy consumed in Route 1 by Volkswagen ID3. (<b>b</b>) SOC% for the energy consumed in Route 1 by Volkswagen E-Lavida.</p>
Full article ">Figure 16 Cont.
<p>(<b>a</b>) SOC% for the energy consumed in Route 1 by Volkswagen ID3. (<b>b</b>) SOC% for the energy consumed in Route 1 by Volkswagen E-Lavida.</p>
Full article ">Figure 17
<p>(<b>a</b>) Battery SOC equals 1 at the start of the journey by E-Lavida. (<b>b</b>) Battery SOC decreased to 0.9 at the end of the journey.</p>
Full article ">Figure 18
<p>(<b>a</b>) Battery SOC equals 100% at the start of the journey by ID3. (<b>b</b>) Battery SOC decreased to 92% at the end of the journey. (<b>c</b>) Some real data from the vehicle.</p>
Full article ">
17 pages, 4006 KiB  
Article
Energy Assessment of Alternative City Bus Lines: A Case Study in Gijón, Spain
by Jose Diaz, Borja Pérez and Francisco J. Fernández
Sustainability 2024, 16(10), 4101; https://doi.org/10.3390/su16104101 - 14 May 2024
Viewed by 701
Abstract
The progressive substitution of traditional buses with electric ones in urban public transport constitutes a fundamental challenge towards sustainable mobility. This paper presents a methodological approach to assess energy consumption in urban bus networks, focusing on a city with varied topography and examining [...] Read more.
The progressive substitution of traditional buses with electric ones in urban public transport constitutes a fundamental challenge towards sustainable mobility. This paper presents a methodological approach to assess energy consumption in urban bus networks, focusing on a city with varied topography and examining alternative bus lines for similar trips. Utilizing a quasi-static longitudinal model, real GPS data, and open access terrain models, the analysis aims to provide a nuanced understanding of energy performance, considering factors such as stop characteristics, gradients, and driving styles. The influence of driving style on commercial speed is observed to be modest, yet significant, in terms of energy efficiency. This research identifies the most advantageous line for transitioning from Internal Combustion Engine (ICE) to Electric Motor (EM) powertrains, resulting in a significant 68% reduction in energy consumption. Beyond specific line details, this methodology offers insights applicable to medium-sized cities, emphasizing efficient route prioritization and providing enhanced user information for informed decisions in the context of sustainable transportation solutions. Full article
(This article belongs to the Collection Sustainability in Urban Transportation Planning)
Show Figures

Figure 1

Figure 1
<p>Aerial photograph of Gijón indicating the relief contours and bus lines: L02 (white), L10 (grey), and L15 (black). Terminus stops (white dots) and a sample user location (red dot) are shown.</p>
Full article ">Figure 2
<p>Block diagram showing sub-systems for internal combustion engine powertrain (<b>top</b>) and electric motor powertrain (<b>bottom</b>) models.</p>
Full article ">Figure 3
<p>Velocity profile comparison for line L10. Passenger stops marked with red dots.</p>
Full article ">Figure 4
<p>Velocity and elevation profiles for the main driving cycles considered in this study Passenger stops marked with red dots.</p>
Full article ">Figure 5
<p>Stop frequency vs. commercial speed (<b>a</b>) and vs. standstill time (<b>b</b>). Small dots represent statistical data of worldwide operators from [<a href="#B18-sustainability-16-04101" class="html-bibr">18</a>].</p>
Full article ">Figure 6
<p>Energy consumption at the motor (solid dots for ICE; open dots for EM) and energy demand at wheel (bars) indicating rolling (r), aerodynamic (a), and braking (b) contributions. The elements represented in blue correspond to the semi-synthetic cycles with <span class="html-italic">v<sub>max</sub></span> = 30 km/h.</p>
Full article ">Figure 7
<p>Relation between braking energy and driving cycle characteristics (gradient <span class="html-italic">α</span> vs. acceleration <span class="html-italic">a</span>) for braking periods (<span class="html-italic">e<sub>t</sub></span> &lt; 0) in normal cycles (<b>left</b>) and semi-synthetic cycles with <span class="html-italic">v<sub>max</sub></span> = 30 km/h (<b>right</b>). Color contours represent the frequency of given conditions during braking periods. The constant <span class="html-italic">e<sub>t</sub></span> curves (white) are computed for a 13 t vehicle with 37 kWh/100 km resistance due to aerodynamic and rolling resistances. Bottom plots represent the average condition for each cycle.</p>
Full article ">
17 pages, 5075 KiB  
Article
CNN-BiLSTM: A Novel Deep Learning Model for Near-Real-Time Daily Wildfire Spread Prediction
by Mohammad Marjani, Masoud Mahdianpari and Fariba Mohammadimanesh
Remote Sens. 2024, 16(8), 1467; https://doi.org/10.3390/rs16081467 - 20 Apr 2024
Cited by 3 | Viewed by 2103
Abstract
Wildfires significantly threaten ecosystems and human lives, necessitating effective prediction models for the management of this destructive phenomenon. This study integrates Convolutional Neural Network (CNN) and Bidirectional Long Short-Term Memory (BiLSTM) modules to develop a novel deep learning model called CNN-BiLSTM for near-real-time [...] Read more.
Wildfires significantly threaten ecosystems and human lives, necessitating effective prediction models for the management of this destructive phenomenon. This study integrates Convolutional Neural Network (CNN) and Bidirectional Long Short-Term Memory (BiLSTM) modules to develop a novel deep learning model called CNN-BiLSTM for near-real-time wildfire spread prediction to capture spatial and temporal patterns. This study uses the Visible Infrared Imaging Radiometer Suite (VIIRS) active fire product and a wide range of environmental variables, including topography, land cover, temperature, NDVI, wind informaiton, precipitation, soil moisture, and runoff to train the CNN-BiLSTM model. A comprehensive exploration of parameter configurations and settings was conducted to optimize the model’s performance. The evaluation results and their comparison with benchmark models, such as a Long Short-Term Memory (LSTM) and CNN-LSTM models, demonstrate the effectiveness of the CNN-BiLSTM model with IoU of F1 Score of 0.58 and 0.73 for validation and training sets, respectively. This innovative approach offers a promising avenue for enhancing wildfire management efforts through its capacity for near-real-time prediction, marking a significant step forward in mitigating the impact of wildfires. Full article
Show Figures

Figure 1

Figure 1
<p>The geographical location of the study area, highlighting the distribution of active fires obtained from the VIIRS dataset. The daily temperature information was extracted from ERA-5, and precipitation data was sourced from the Australian Landscape Water Balance.</p>
Full article ">Figure 2
<p>Patch extraction process and data structure preparation: section (<b>A</b>) shows initial data structure before patch extraction with daily temporal resolution; section (<b>B</b>) provides information about the data shapes; section (<b>C</b>) indicates the 50% overlap for two consecutive patches (N and N + 1) in x direction.</p>
Full article ">Figure 3
<p>The overall framework of TD-CNNs. In this visualization, the orange cubes represent the input tensors, and the purple cubes signify the extracted features using CNNs. Additionally, the variable T denotes the time dimension.</p>
Full article ">Figure 4
<p>The architecture of the proposed CNN-BiLSTM model.</p>
Full article ">Figure 5
<p>Schematic definition of TP, FN, FP, and TN.</p>
Full article ">Figure 6
<p>Schematic definition of precision and recall.</p>
Full article ">Figure 7
<p>Histograms of the IoU, recall, precision, and F1 of the CNN-BiLSTM model for the validation set (<b>top</b>) and training set (<b>bottom</b>).</p>
Full article ">Figure 8
<p>CNN-BiLSTM predictions for 20 validation samples. The limitation of the model in predictions is highlighted with an orange arrow.</p>
Full article ">Figure 9
<p>The IoU score of validation set for different threshold values. The vertical red line shows where the threshold value is optimal and the green cross marker indicates the IoU based on the optimal threshold value.</p>
Full article ">Figure 10
<p>Histograms of the IoU for different time steps.</p>
Full article ">Figure 11
<p>Daily correlation between predictions and environment variables.</p>
Full article ">
Back to TopTop