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

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

Years

Between: -

Subjects

remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline

Journals

Article Types

Countries / Regions

Search Results (151)

Search Parameters:
Keywords = spring onset

Order results
Result details
Results per page
Select all
Export citation of selected articles as:
20 pages, 7366 KiB  
Article
How Severe Was the 2022 Flash Drought in the Yangtze River Basin?
by Liyan Yang and Jia Wei
Remote Sens. 2024, 16(22), 4122; https://doi.org/10.3390/rs16224122 - 5 Nov 2024
Viewed by 405
Abstract
Flash droughts, characterized by their rapid onset and severe impacts, have critical implications for the ecological environment and water resource security. However, inconsistent definitions of flash droughts have hindered scientific assessments of drought severity, limiting efforts in disaster prevention and mitigation. In this [...] Read more.
Flash droughts, characterized by their rapid onset and severe impacts, have critical implications for the ecological environment and water resource security. However, inconsistent definitions of flash droughts have hindered scientific assessments of drought severity, limiting efforts in disaster prevention and mitigation. In this study, we propose a new method for explicitly characterizing flash drought events, with particular emphasis on the process of soil moisture recovery. The temporal and spatial evolution of flash droughts over the Yangtze River Basin was analyzed, and the severity of the extreme flash drought in 2022 was assessed by comparing its characteristics and impacts with those of three typical dry years. Additionally, the driving factors of the 2022 flash drought were evaluated from multiple perspectives. Results indicate that the new identification method for flash droughts is reasonable and reliable. In recent years, the frequency and duration of flash droughts have significantly increased, with the Dongting Lake and Poyang Lake basins being particularly affected. Spring and summer were identified as peak seasons for flash droughts, with the middle reaches most affected in spring, while summer droughts tend to impact the entire basin. Compared to 2006, 2011, and 2013, the flash drought in 2022 affected the largest area, with the highest number of grids experiencing two flash drought events and a development rate exceeding 15%. Moreover, the summer heat in 2022 was more extreme than in the other three years, extending from spring to fall, especially during July–August. Its evolution was driven by the Western Pacific Subtropical High, which suppressed precipitation and elevated temperatures. The divergence of water vapor flux intensified water shortages, while anomalies in latent and sensible heat fluxes increased surface evaporation and heat transfer, further disturbing the regional water cycle. This study provides valuable insights for flash drought monitoring and early warning in the context of a changing climate. Full article
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>(<b>a</b>) The location and the elevation of the Yangtze River Basin. (<b>b</b>) The subbasins and average annual soil moisture of the Yangtze River Basin.</p>
Full article ">Figure 2
<p>Schematic diagram of the definition of the onset, development and recovery stages in flash drought events used in this study. v is the development rate, defined as the average rate of decline in soil moisture percentiles.</p>
Full article ">Figure 3
<p>Spatial distribution of flash drought characteristics in the Yangtze River Basin from 1950 to 2022. (<b>a</b>) Frequency: the multi-year average number of flash drought events. (<b>b</b>) Mean duration: the average number of days each drought event lasted. (<b>c</b>) Mean development rate: the average speed at which each flash drought event developed. (<b>d</b>) Mean intensity: the average intensity of each drought event. The white grid points indicate that there was no flash drought. The numbers in the (<b>a</b>–<b>d</b>) represent the subbasins.</p>
Full article ">Figure 4
<p>Trend changes in flash drought characteristics in the Yangtze River Basin from 1950 to 2022. (<b>a</b>) Regional average frequency of flash drought events. (<b>b</b>) Regional mean duration of flash drought events. (<b>c</b>) Regional average development rate of flash drought events. (<b>d</b>) Regional mean intensity of flash drought events. The black dotted line is the trend line. The solid blue line is the pre-abrupt mean. The solid red line is the post-abrupt mean. The overall trends passed the MK significance test with <span class="html-italic">p</span> &lt; 0.05 in (<b>a</b>,<b>c</b>,<b>d</b>), and the overall trends passed the MK significance test with <span class="html-italic">p</span> &lt; 0.1 in (<b>b</b>).</p>
Full article ">Figure 5
<p>Seasonal characteristics of flash droughts in the Yangtze River Basin. (<b>a</b>) The multi-year average frequency of flash drought in different seasons. (<b>b</b>) The multi-year average total days of flash drought in different seasons. (<b>c</b>) Most frequent season of first flash drought event per grid, 1950–2022. The numbers in the figure represent the subbasins.</p>
Full article ">Figure 6
<p>(<b>a</b>) The percentage of grid points affected by flash droughts relative to the total number of grid points in the entire basin for each pentad. (<b>b</b>) Spatial distribution of flash drought frequency in the Yangtze River Basin in 2022. (<b>c</b>) Spatial distribution of flash drought duration in the Yangtze River Basin in 2022. (<b>d</b>) Spatial distribution of flash drought development rate in the Yangtze River Basin in 2022. (<b>e</b>) Spatial distribution of flash drought intensity in the Yangtze River Basin in 2022. The numbers in the figure (<b>b</b>–<b>e</b>) represent the subbasins.</p>
Full article ">Figure 7
<p>(<b>a</b>) Standardized precipitation anomalies from 2003 to 2022. (<b>b</b>) Comparison of regional average flash drought characteristics in typical drought years. (<b>c</b>) Percentage of different flash drought frequency in typical drought years. (<b>d</b>) Percentage of different flash drought durations in typical drought years. (<b>e</b>) Percentage of different flash drought development rates in typical drought years. (<b>f</b>) Percentage of different flash drought intensities in typical drought years. The red error bar in (<b>c</b>–<b>f</b>) indicates the difference from the multi-year average (2003–2022), and an upward trend indicates that the value is less than the multi-year average.</p>
Full article ">Figure 8
<p>(<b>a</b>) Time series of accumulated precipitation in drought years. (<b>b</b>) Time series of maximum temperature in drought years. The horizontal axis represents the pentad of the year, and the black solid line indicates the multi-year average for the past twenty years (2003–2022) in (<b>a</b>,<b>b</b>).</p>
Full article ">Figure 9
<p>Spatial patterns of hydrometeorological variables over the Yangtze River Basin in 2022. (<b>a</b>) Precipitation in the Yangtze River Basin in June 2022. (<b>b</b>) Precipitation in the Yangtze River Basin in July 2022. (<b>c</b>) Precipitation in the Yangtze River Basin in August 2022. (<b>d</b>) Temperature in the Yangtze River Basin in June 2022. (<b>e</b>) Temperature in the Yangtze River Basin in July 2022. (<b>f</b>) Temperature in the Yangtze River Basin in August 2022.</p>
Full article ">Figure 10
<p>Spatial patterns of climatic variables over the Yangtze River Basin in 2022. (<b>a</b>) Moisture divergence in the Yangtze River Basin in June 2022. (<b>b</b>) Moisture divergence in the Yangtze River Basin in July 2022. (<b>c</b>) Moisture divergence in the Yangtze River Basin in August 2022. (<b>d</b>) Geopotential height of 500 hpa in the Yangtze River Basin in June 2022. (<b>e</b>) Geopotential height of 500 hpa in the Yangtze River Basin in July 2022. (<b>f</b>) Geopotential height of 500 hpa in the Yangtze River Basin in August 2022.</p>
Full article ">Figure 11
<p>Spatial patterns of surface heat flux over the Yangtze River Basin in 2022. (<b>a</b>) Latent heat flux in the Yangtze River Basin in June 2022. (<b>b</b>) Latent heat flux in the Yangtze River Basin in July 2022. (<b>c</b>) Latent heat flux in the Yangtze River Basin in August 2022. (<b>d</b>) Sensible heat flux in the Yangtze River Basin in June 2022. (<b>e</b>) Sensible heat flux in the Yangtze River Basin in July 2022. (<b>f</b>) Sensible heat flux in the Yangtze River Basin in August 2022. The ERA5 dataset convention for vertical fluxes is positive downwards. Negative surface latent heat flux indicates evaporation and negative surface sensible heat flux indicates heat transfer from the surface to the atmosphere.</p>
Full article ">
19 pages, 5558 KiB  
Article
Convolution Neural Network Development for Identifying Damage in Vibrating Pylons with Mass Attachments
by George D. Manolis and Georgios I. Dadoulis
Sensors 2024, 24(19), 6255; https://doi.org/10.3390/s24196255 - 27 Sep 2024
Viewed by 533
Abstract
A convolution neural network (CNN) is developed in this work to detect damage in pylons by measuring their vibratory response. More specifically, damage detection through testing relies on the development of damage-sensitive indicators, which are then used to reach a decision regarding the [...] Read more.
A convolution neural network (CNN) is developed in this work to detect damage in pylons by measuring their vibratory response. More specifically, damage detection through testing relies on the development of damage-sensitive indicators, which are then used to reach a decision regarding the existence/absence of damage, provided they have been retrieved from at least two distinct structural states. Damage indicators, however, exhibit a relatively low sensitivity regarding the onset of structural damage, further exacerbated by the low amplitude response to a variety of environmentally induced loads. To this end, a mathematical model is developed to interpret the experimental data recovered from a fixed-base pylon with a top mass attachment to transverse motion. Damage is introduced in the mathematical model in the form of springs corresponding to the cracking of the beam’s lower end. Families of numerically generated acceleration records are produced at select stations along the beam’s height, which are then used for training a CNN. Once trained, it is used to identify damage from acceleration records produced from a series of experiments. Difficulties faced by CNN in correctly identifying the presence/absence of damage in the pylon are discussed, and steps taken to improve the quality of the results are proposed. Full article
Show Figures

Figure 1

Figure 1
<p>Typical <span class="html-italic">autoencoder</span> architecture.</p>
Full article ">Figure 2
<p>Flexible pylon configurations: (<b>a</b>) Intact condition; (<b>b</b>) placement of a 5 kg mass at the top; (<b>c</b>) crack development at the base during a fatigue test involving a large number of cycles induced by an actuator at the top.</p>
Full article ">Figure 3
<p>Cantilevered pylon: (<b>a</b>) cross section and length geometry; (<b>b</b>) free-body diagram; (<b>c</b>) top segment detail; (<b>d</b>) pylon base details.</p>
Full article ">Figure 4
<p>(<b>a</b>) Cantilevered pylon under cyclic loading; (<b>b</b>) details of the top mass; (<b>c</b>) details of cracking formation at the pylon’s base [<a href="#B28-sensors-24-06255" class="html-bibr">28</a>].</p>
Full article ">Figure 5
<p>Cantilevered pylon mechanical model.</p>
Full article ">Figure 6
<p>Comparison between the analytically computed and experimentally measured first flexural eigenfunction <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>φ</mi> </mrow> <mrow> <mn>1</mn> </mrow> </msub> </mrow> </semantics></math> of the tapered pylon versus normalized pylon height <math display="inline"><semantics> <mrow> <mi>ξ</mi> </mrow> </semantics></math>.</p>
Full article ">Figure 7
<p>Pylon instrumentation (left) and typical acceleration time histories recorded by the sensors for cases: (<b>a</b>) Intact pylon, <math display="inline"><semantics> <mrow> <mn>225</mn> <mtext> </mtext> <mi mathvariant="normal">s</mi> </mrow> </semantics></math> duration; (<b>b</b>) Pylon with attached mass at the top, <math display="inline"><semantics> <mrow> <mn>225</mn> <mtext> </mtext> <mi mathvariant="normal">s</mi> </mrow> </semantics></math> duration; (<b>c</b>) Pylon with a base crack and measurements in the N–S direction, <math display="inline"><semantics> <mrow> <mn>457</mn> <mtext> </mtext> <mi mathvariant="normal">s</mi> </mrow> </semantics></math> duration; (<b>d</b>) Pylon with a base crack and measurements in the E–W direction, <math display="inline"><semantics> <mrow> <mn>457</mn> <mtext> </mtext> <mi mathvariant="normal">s</mi> </mrow> </semantics></math> duration.</p>
Full article ">Figure 8
<p>The pylon’s generalized coordinates as derived from the weight functions of the Autoencoder: From top downwards are cases (<b>a</b>) intact pylon, (<b>b</b>) intact pylon with a top mass, (<b>c</b>) pylon with a base crack and measurements in the N–S direction, (<b>d</b>) pylon with a base crack and measurements in the W–E direction.</p>
Full article ">Figure 9
<p>The first two time histories out of a bin of 400 time histories, recorded by sensor <b>a805,</b> that are generated as input for the development of spectrograms.</p>
Full article ">Figure 10
<p>Representative spectrograms for pylon scenarios: (<b>a</b>) Intact, (<b>b</b>) intact pylon with a top mass, (<b>c</b>) pylon with a base crack and N–S recording, and (<b>d</b>) pylon with a base crack and W–E recording.</p>
Full article ">Figure 11
<p>CNN architecture: Input, neuron connections, and output [<a href="#B32-sensors-24-06255" class="html-bibr">32</a>]. Note: The question mark (?) at the top was set equal to one in our development.</p>
Full article ">Figure 12
<p>CNN development: (<b>a</b>) Increase in accuracy and (<b>b</b>) decrease in the value of the <span class="html-italic">loss</span> function as a function of the <span class="html-italic">epochs</span>.</p>
Full article ">Figure A1
<p>Dimensionless modal participation factor <span class="html-italic">μ<sub>j</sub></span> as a function of the mass ratio <span class="html-italic">R</span> for flexural vibrations due to a point force applied at the top of the pylon.</p>
Full article ">
21 pages, 4579 KiB  
Article
Differentiated In-Row Soil Management in a High-Density Olive Orchard: Effects on Weed Control, Tree Growth and Yield, and Economic and Environmental Sustainability
by Enrico Maria Lodolini, Nadia Palmieri, Alberto de Iudicibus, Pompea Gabriella Lucchese, Matteo Zucchini, Veronica Giorgi, Samuele Crescenzi, Kaies Mezrioui, Davide Neri, Corrado Ciaccia and Alberto Assirelli
Agronomy 2024, 14(9), 2051; https://doi.org/10.3390/agronomy14092051 - 7 Sep 2024
Viewed by 666
Abstract
Two different in-row soil management techniques were compared in the Olive Orchard Innovation Long-term experiment of the Council for Agricultural Research and Economics, Research Centre for Olive, Fruit, and Citrus Crops in Rome, Italy. Rows were managed with an in-row rotary tiller and [...] Read more.
Two different in-row soil management techniques were compared in the Olive Orchard Innovation Long-term experiment of the Council for Agricultural Research and Economics, Research Centre for Olive, Fruit, and Citrus Crops in Rome, Italy. Rows were managed with an in-row rotary tiller and with synthetic mulching using permeable polypropylene placed after cultivar Maurino olive trees planting. The effects of the two treatments were assessed through weed soil coverage and the growth of the olive trees. Results showed better agronomic performance associated with synthetic mulching. The weed control effect along the row of a young high-density olive orchard was higher with the synthetic mulching compared to hoeing. The effect of the synthetic mulching seemed to disappear when removed from the ground (spring 2023) since no significant differences were found for tree size and yield in the two tested in-row soil management systems at the end of 2023. Finally, the growth of the young olive trees (Trunk Cross Sectional Area, Height, and Canopy expansion) measured across the three years, was higher for the synthetic mulched row than the hoed one. The use of synthetic mulching along the row positively forced the vegetative growth of the young olive trees and anticipated the onset of fruit production compared to periodical hoeing: a significantly higher fruit production was registered three years after planting. Root diameter was higher under synthetic mulching one year after planting, and no differences were observed in the following sampling dates showing similar fluctuations linked to the seasonal growth pattern. The life cycle assessment and costing highlighted that the application of mulching had a higher eco- and economic-efficiency than the periodical in-row soil hoeing. Full article
(This article belongs to the Special Issue The Impact of Mulching on Crop Production and Farmland Environment)
Show Figures

Figure 1

Figure 1
<p>Trend of the monthly mean air temperature and total precipitation during the trial.</p>
Full article ">Figure 2
<p>Compared in-row soil management treatments: synthetic mulching (<b>a</b>) and hoeing (<b>b</b>).</p>
Full article ">Figure 3
<p>Weed cover (%) in the intra-row area of the hoeing treatment in September 2020 (<b>a</b>) and November 2020 (<b>b</b>).</p>
Full article ">Figure 4
<p>Weed soil coverage status along the synthetic mulching row and along hoeing one at (<b>A</b>,<b>D</b>) 35 DAT, (<b>B</b>,<b>E</b>) 69 DAT, and (<b>C</b>,<b>F</b>) 147 DAT.</p>
Full article ">Figure 5
<p>Principal Component Analysis for the February (<b>a</b>,<b>b</b>) and the March (<b>c</b>,<b>d</b>) samplings. TOT = Total coverage (%); R = Richness; e = Evenness; PRN (%) = life form (% of perennials); CH = Canopy heights (m); SW = Seeds weight (g); SLA= Surface Leaf Area; Grass-like % = graminoid species incidence; N-fixing % = leguminous species incidence.</p>
Full article ">Figure 6
<p>Trunk Cross Sectional Area (TCSA) measured at the planting and at the end of each growing season. Data are shown as mean ± standard error of 30 replicates. Different letters indicate significant differences between treatments within each measurement date according to the Tukey–Kramer HSD test (<span class="html-italic">p</span> &lt; 0.05).</p>
Full article ">Figure 7
<p>Canopy volume (m<sup>3</sup>) calculated at the end of the 2nd and 3rd growing season. Data are shown as mean ± standard error of 30 replicates. Different letters indicate significant differences between treatments within each measurement date according to the Tukey–Kramer HSD test (<span class="html-italic">p</span> &lt;0.05).</p>
Full article ">Figure 8
<p>Root diameter (mm) measured in November 2020, December 2021, and October 2022. Data are shown as mean ± standard error of 8 replicates. Different letters indicate significant differences between the two tested treatments within the same sampling date according to the Tukey–Kramer HSD test (<span class="html-italic">p</span> &lt; 0.05). When the letters are not shown, it indicates no significant differences between treatments within the same sampling date.</p>
Full article ">Figure 9
<p>TCSA/RLD ratio (cm g) evolution from November 2020 to October 2022 for mulching and hoeing. Each point represents the ratio between the average values of TCSA and RLD within each measurement date.</p>
Full article ">Figure 10
<p>The impact of environmental impact categories per each row management strategy (FU = 1 hectare). The values are expressed as percentages in relation to the field management adopted technique with the highest environmental impact, which is expressed as 100%. Legend: AD: Abiotic depletion; Adf: Abiotic depletion-fossil fuels; GWP: Global Warming (GWP100a); ODP: Ozone Layer Depletion; HT: Human Toxicity; FW: Fresh water aquatic ecotox.; ME: Marine Aquatic Eco-Toxicity; TE: Terrestrial Eco-Toxicity; PO: Photochemical Oxidation; AC: Acidification; EU: Eutrophication.</p>
Full article ">Figure 11
<p>The sensitivity results: the impact of environmental impact categories per each row management strategy (FU = 1 Mg of product). The values are expressed as percentages in relation to the field management with the highest environmental impact, which is expressed as 100%. Legend: AD: Abiotic depletion; Adf: Abiotic depletion-fossil fuels; GWP: Global Warming (GWP100a); ODP: Ozone Layer Depletion; HT: Human Toxicity; FW: Fresh water aquatic ecotox; ME: Marine Aquatic Eco-Toxicity; TE: Terrestrial Eco-Toxicity; PO: Photochemical Oxidation; AC: Acidification; EU: Eutrophication.</p>
Full article ">
15 pages, 1434 KiB  
Article
Differential Responses of Bilberry (Vaccinium myrtillus) Phenology and Density to a Changing Environment: A Study from Western Carpathians
by Martin Kubov, Peter Fleischer, Jakub Tomes, Mohammad Mukarram, Rastislav Janík, Benson Turyasingura, Peter Fleischer and Branislav Schieber
Plants 2024, 13(17), 2406; https://doi.org/10.3390/plants13172406 - 28 Aug 2024
Viewed by 867
Abstract
Environmental factors regulate the regeneration of mountain spruce forests, with drought, wind, and bark beetles causing the maximum damage. How these factors minimise spruce regeneration is still poorly understood. We conducted this study to investigate how the phenology and population dynamics of bilberry [...] Read more.
Environmental factors regulate the regeneration of mountain spruce forests, with drought, wind, and bark beetles causing the maximum damage. How these factors minimise spruce regeneration is still poorly understood. We conducted this study to investigate how the phenology and population dynamics of bilberry (Vaccinium myrtillus L.), a dominant understory species of mountain spruce forests, are related to selected environmental factors that are modified by natural disturbances (bark beetle and wind). For this, we analysed bilberry at different sites affected by bark beetles and adjacent undisturbed forests in the Tatra National Park (TANAP) during the growing season (April–September) in 2016–2021, six years after the initial bark beetle attack. The observations were taken along an altitudinal gradient (1100–1250–1400 m a.s.l.) in two habitats (disturbed spruce forest—D, undisturbed spruce forest—U). We found that habitat and altitude influenced the onset of selected phenological phases, such as the earliest onset at low altitudes (1100 m a.s.l.) in disturbed forest stands and the latest at high altitudes (1400 m a.s.l.) in undisturbed stands. Although there were non-significant differences between habitats and altitudes, likely due to local climate conditions and the absence of a tree layer, these findings suggest that bilberry can partially thrive in disturbed forest stands. Despite temperature fluctuations during early spring, the longer growing season benefits its growth. Full article
(This article belongs to the Section Plant Ecology)
Show Figures

Figure 1

Figure 1
<p>Location of the study area in Europe and in the Slovak Republic (TANAP-Tatra National Park; 1100, 1250, 1400-altitude of research plots (m a.s.l.); D-disturbed forest stand; U-undisturbed forest stand).</p>
Full article ">Figure 2
<p>Example of one cluster plot design used in the study comprising 12 square subplots.</p>
Full article ">Figure 3
<p>Climate graphs for the three study sites for 2016-2021. (Boxplots of median values of temperature and precipitation and corresponding quartiles. Dots in some boxplots represent potential outliers).</p>
Full article ">Figure 4
<p>Changes in the density of bilberry shoots in different habitats on a vertical gradient between 2016 and 2021. (U—Undisturbed Forest; D—Disturbed Forest). (Boxplots of median values of temperature and precipitation and corresponding quartiles. Dots in some boxplots represent potential outliers).</p>
Full article ">
12 pages, 2825 KiB  
Article
An Observational Longitudinal Study on Seasonal Variations in Tourette Syndrome: Evidence for a Role of Ambient Temperature in Tic Exacerbation
by Jacopo Lamanna, Riccardo Mazzoleni, Ramona Farina, Mattia Ferro, Roberta Galentino, Mauro Porta and Antonio Malgaroli
Biomedicines 2024, 12(8), 1668; https://doi.org/10.3390/biomedicines12081668 - 26 Jul 2024
Viewed by 763
Abstract
Tourette syndrome (TS) is a high-incidence neurobehavioral disorder that generally begins in childhood. Several factors play a role in its etiology, including genetic influence and auto-immune activation by streptococcal infections. In general, symptoms subside after the end of adolescence, but, in a significant [...] Read more.
Tourette syndrome (TS) is a high-incidence neurobehavioral disorder that generally begins in childhood. Several factors play a role in its etiology, including genetic influence and auto-immune activation by streptococcal infections. In general, symptoms subside after the end of adolescence, but, in a significant number of patients, they remain in adulthood. In this study, we evaluated temporal variations in the two core clinical features of TS including tics and obsessive–compulsive disorder (OCD) symptoms. An observational longitudinal study lasting 15 months (2017–2019) was conducted on a cohort of 24 people recruited in Milan (Italy) who were diagnosed with a subtype of TS known as obsessive–compulsive tic disorder. Inclusion criteria included a global score of the Yale global tic severity scale (Y-GTSS) > 50, a Yale–Brown obsessive–compulsive scale (Y-BOCS) global score > 15, and TS onset at least one year prior. Y-GTSS and Y-BOCS data were acquired at six time points, together with local environmental data. Tics, but not OCD symptoms, were found to be more severe in spring and summer compared with winter and autumn (p < 0.001). Changes in tics displayed an appreciable oscillation pattern in the same subject and also a clear synchrony among different subjects, indicating an external orchestrating factor. Ambient temperature showed a significant correlation with Y-GTSS measurements (p < 0.001). We argue that the increase in tics observed during hot seasons can be related to increasing ambient temperature. We believe that our results can shed light on the seasonal dynamics of TS symptomatology and provide clues for preventing their worsening over the year. Full article
Show Figures

Figure 1

Figure 1
<p>Study flowchart. Criteria and numbers related to patient enrollment, non-inclusions, and drop-outs during this study.</p>
Full article ">Figure 2
<p>Y-GTSS scores along seasons. Mean values and standard deviation for Y-GTSS measurements over a period of 15 months with quarterly sampling are shown. Background colors represent colder (blue) vs. warmer (red) semesters. Horizontal bars above and below the graph represent significant differences among seasons according to post hoc multiple comparisons (Bonferroni correction). * <span class="html-italic">p</span> &lt; 0.05, ** <span class="html-italic">p</span> &lt; 0.01, *** <span class="html-italic">p</span> &lt; 0.001.</p>
Full article ">Figure 3
<p>Y-BOCS scores along seasons. Mean values and standard deviation for Y-BOCS measurements over a period of 15 months with quarterly sampling are shown (main effect of season: F<sub>(4,122)</sub> = 33.490, <span class="html-italic">p</span> &lt; 0.001, ANOVA on GLME; all pairwise contrasts with Bonferroni correction n.s.). Background colors represent colder (blue) vs. warmer (red) semesters.</p>
Full article ">Figure 4
<p>Autocorrelation analysis of seasonal Y-GTSS scores. (<b>A</b>) Five individual corrected series of Y-GTSS measurements are shown across the six seasons. The dashed line represents a spline interpolation (50 time points) of the 6 data points. (<b>B</b>) Mean (black line) and standard deviation (red area) of corrected and interpolated Y-GTSS data from the 17 subjects selected for autocorrelation analysis. (<b>C</b>) Autocorrelation of the individual series of Y-GTSS scores for each patient (random colors). Dashed and straight grey lines indicate 95% and 99% confidence intervals, respectively. (<b>D</b>) Five individual trends in Y-BOCS scores across our measuring time span. In comparison with Y-GTSS individual trends in (<b>A</b>), Y-BOCS scores do not show a clear oscillatory pattern across seasons.</p>
Full article ">Figure 5
<p>Cross-correlation analysis for detrended seasonal Y-GTSS series. (<b>A</b>) Cross-correlation between the series of Y-GTSS scores (50 points interpolation) and the associated temperature series for each subject (random colors), where the x-axis represents lag and the y-axis represents the level of correlation. Dashed and straight grey lines represent 95% and 99% confidence intervals, respectively. (<b>B</b>) Cross-correlation between each subject’s detrended series of Y-GTSS scores (50-point interpolation). The color grade represents the level of correlation (blue and red are negative and positive values, respectively). The vertical sorting (from top to bottom) was based on the lag of maximum correlation (from the largest to the smallest lag). (<b>C</b>) Maximum values of the cross-correlation between the Y-GTSS score series and temperature series are shown as a function of their lag (in months). (<b>D</b>) Negative peaks of the cross-correlation between the temperature series and Y-GTSS series are shown as a function of their lag (in months). (<b>E</b>) Negative peaks of autocorrelation in the Y-GTSS series are shown as a function of their lag (in months). Different marker shapes for plots C-E indicate three different clusters obtained using cluster analysis.</p>
Full article ">
30 pages, 10784 KiB  
Article
Phenology and Plant Functional Type Link Optical Properties of Vegetation Canopies to Patterns of Vertical Vegetation Complexity
by Duncan Jurayj, Rebecca Bowers and Jessica V. Fayne
Remote Sens. 2024, 16(14), 2577; https://doi.org/10.3390/rs16142577 - 13 Jul 2024
Viewed by 1033
Abstract
Vegetation vertical complexity influences biodiversity and ecosystem productivity. Rapid warming in the boreal region is altering patterns of vertical complexity. LiDAR sensors offer novel structural metrics for quantifying these changes, but their spatiotemporal limitations and their need for ecological context complicate their application [...] Read more.
Vegetation vertical complexity influences biodiversity and ecosystem productivity. Rapid warming in the boreal region is altering patterns of vertical complexity. LiDAR sensors offer novel structural metrics for quantifying these changes, but their spatiotemporal limitations and their need for ecological context complicate their application and interpretation. Satellite variables can estimate LiDAR metrics, but retrievals of vegetation structure using optical reflectance can lack interpretability and accuracy. We compare vertical complexity from the airborne LiDAR Land Vegetation and Ice Sensor (LVIS) in boreal Canada and Alaska to plant functional type, optical, and phenological variables. We show that spring onset and green season length from satellite phenology algorithms are more strongly correlated with vegetation vertical complexity (R = 0.43–0.63) than optical reflectance (R = 0.03–0.43). Median annual temperature explained patterns of vegetation vertical complexity (R = 0.45), but only when paired with plant functional type data. Random forest models effectively learned patterns of vegetation vertical complexity using plant functional type and phenological variables, but the validation performance depended on the validation methodology (R2 = 0.50–0.80). In correlating satellite phenology, plant functional type, and vegetation vertical complexity, we propose new methods of retrieving vertical complexity with satellite data. Full article
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>A map of Canada and Alaska showing the locations of the LiDAR flightlines used in this study. The green shading over the base map delineates the Arctic Boreal Vulnerability Experiment (ABoVE) Core Region. Blue and orange flightlines, taken from the 2017 and 2019 airborne LIDAR campaigns, respectively, were included in all analyses and model training. Red flightlines were excluded from any analysis and were used exclusively for validating random forest models.</p>
Full article ">Figure 2
<p>The ABoVE-LULC PFT composition expressed as a percentage of the total number of PFT observations for the following datasets: the entire ABoVE Domain, every LVIS flightline from 2017 and 2019 in the ABoVE Domain, the LVIS training flightlines used in this study, and the LVIS validation flightlines used in this study.</p>
Full article ">Figure 3
<p>(<b>a</b>) Box plots of plant height for each original NEON growth form, with legend colors corresponding to the reduced NEON growth groups; (<b>b</b>) the relationship between MCD and plant height, with colors denoting NEON growth groups; (<b>c</b>) the relationship between SD and plant height, with colors denoting NEON growth groups; (<b>d</b>) the relationship between stem diameter and crown height, with colors denoting specific taxa. For plots (<b>a</b>–<b>c</b>), the universal relationship is depicted by a black LOESS line and the group-specific relationships are depicted by colored LOESS lines. Gray shading indicates the 95% confidence interval of each LOESS line.</p>
Full article ">Figure 4
<p>The universal relationship between canopy height and VC is depicted by the black scatter plot and the relationships between canopy height and VC for individual PFTs are depicted by the colored LOESS lines. The gray shading around each LOESS line represents the 95 percent confidence interval. Points with canopy heights near zero but high VCs likely represent LiDAR measurements of wetlands and water that are erroneously classified.</p>
Full article ">Figure 5
<p>Violin plots and box plots of VC for each PFT. Violin width depicts VC value density while the bottom, middle, and top lines of each white box plot represent the 25th, 50th, and 75th percentiles, respectively.</p>
Full article ">Figure 6
<p>LiDAR waveforms associated with the ABoVE-LULC (<b>a</b>) coniferous forest class; (<b>b</b>) deciduous forest class; (<b>c</b>) mixed forest class; (<b>d</b>) and woodland forest class. Each PFT class is represented by 100 randomly selected waveforms. Red and blue lines represent the mean and median values from those 100 waveforms, respectively. Violins in the upper right corner of each plot show the distribution of VC values associated with each group of 100 waveforms. Refer to <a href="#remotesensing-16-02577-f005" class="html-fig">Figure 5</a> for a comparison of the distributions of VC values in this figure and the distributions of VC values associated with each PFT across the entire ABoVE study domain.</p>
Full article ">Figure 7
<p>LiDAR waveforms associated with different ranges of VC between (<b>a</b>) 0.0 and 0.2; (<b>b</b>) 0.2 and 0.4; (<b>c</b>) 0.4 and 0.6; (<b>d</b>) and 0.6 and 0.8. Each VC bin is represented by 100 randomly selected waveforms. Red and blue lines represent the mean and median values from those 100 waveforms, respectively. Refer to <a href="#remotesensing-16-02577-f004" class="html-fig">Figure 4</a> to compare these waveform structures to the relationship between canopy height and VC.</p>
Full article ">Figure 8
<p>Relationships of VC to (<b>a</b>) MAT and (<b>b</b>) PRCP. The relationship across all PFTs is depicted by the red LOESS line with gray borders and the black scatter plot. The relationships for individual PFTs are depicted by the colored LOESS lines. Light gray shading around each LOESS line represents the 95 percent confidence interval.</p>
Full article ">Figure 9
<p>The percentage of each ABoVE-LULC PFT class that is present in the different Köppen–Geiger climate classes and subclasses for every LVIS 2017 and 2019 flightline in the ABoVE Domain.</p>
Full article ">Figure 10
<p>Relationships between (<b>a</b>) LTSO, (<b>b</b>) LTGSL, (<b>c</b>) ASO, (<b>d</b>) AGSL, and VC. The relationship for all PFTs is depicted by the red LOESS line with gray borders and the black scatter plot. The relationships for individual PFTs are depicted by the colored LOESS lines. Light gray shading around each LOESS line represents the 95 percent confidence interval.</p>
Full article ">Figure 11
<p>Relationships between (<b>a</b>) SWIR2, (<b>b</b>) NBR, (<b>c</b>) NIR, (<b>d</b>) NDVI, and VC. The relationship for all PFTs is depicted by the red LOESS line with gray borders and the black scatter plot. The relationships for individual PFTs are depicted by the colored LOESS lines. Light gray shading around each LOESS line represents the 95 percent confidence interval.</p>
Full article ">Figure 12
<p>(<b>a</b>) Variable importance plot for random forest Model C; (<b>b</b>) the relationship between predicted and observed VC for Model C on the sampled validation dataset; (<b>c</b>) the relationship between predicted and observed VC for Model C on the six LVIS flightlines that were reserved from analysis, training, and testing. The red lines in plots (<b>a</b>,<b>b</b>) represent a perfect one-to-one relationship.</p>
Full article ">Figure 13
<p>(<b>a</b>) A map of Northern Canada and Alaska indicating the location of the HLS tile where random forest Model C was used to predict VC. The HLS tile in (<b>a</b>) overlaps with one of the LVIS 2017 held-out validation flightlines. A closer view of the predicted VC values for the entire HLS tile which borders Great Slave Lake in the Northwest Territories is offered in (<b>b</b>). The observed VC values for the 2017 reserved flightline are depicted in (<b>c</b>) The VC values predicted by Model C are shown in (<b>d</b>). The bright green border in (<b>d</b>) represents the boundary of the LVIS 2017 validation flightline shown in (<b>c</b>).</p>
Full article ">
24 pages, 20759 KiB  
Article
Snowmelt Onset and Caribou (Rangifer tarandus) Spring Migration
by Mariah T. Matias, Joan M. Ramage, Eliezer Gurarie and Mary J. Brodzik
Remote Sens. 2024, 16(13), 2391; https://doi.org/10.3390/rs16132391 - 29 Jun 2024
Viewed by 1349
Abstract
Caribou (Rangifer tarandus) undergo exceptionally large, annual synchronized migrations of thousands of kilometers, triggered by their shared environmental stimuli. The proximate triggers of those migrations remain mysterious, though snow characteristics play an important role due to their influence on the mechanics [...] Read more.
Caribou (Rangifer tarandus) undergo exceptionally large, annual synchronized migrations of thousands of kilometers, triggered by their shared environmental stimuli. The proximate triggers of those migrations remain mysterious, though snow characteristics play an important role due to their influence on the mechanics of locomotion. We investigate whether the snow melt–refreeze status relates to caribou movement, using previously collected Global Positioning System (GPS) caribou collar data. We analyzed 117 individual female caribou with >30,000 observations between 2007 and 2016 from the Bathurst herd in Northern Canada. We used a hierarchical model to estimate the beginning, duration, and end of spring migration and compared these statistics against snow pack melt characteristics derived from 37 GHz vertically polarized (37V GHz) Calibrated Enhanced-Resolution Brightness Temperatures (CETB) at 3.125 km resolution. The timing of migration for Bathurst caribou generally tracked the snowmelt onset. The start of migration was closely linked to the main melt onset in the wintering areas, occurring on average 2.6 days later (range −1.9 to 8.4, se 0.28, n = 10). The weighted linear regression was also highly significant (p-value = 0.002, R2=0.717). The relationship between migration arrival times and the main melt onset on the calving grounds (R2 = 0.688, p-value = 0.003), however, had a considerably more variable lag (mean 13.3 d, se 0.67, range 3.1–20.4). No migrations ended before the main melt onset at the calving grounds. Thawing conditions may provide a trigger for migration or favorable conditions that increase animal mobility, and suggest that the snow properties are more important than snow presence. Further work is needed to understand how widespread this is and why there is such a relationship. Full article
(This article belongs to the Special Issue Understanding the Movement Ecology of Wildlife on the Changing Planet)
Show Figures

Figure 1

Figure 1
<p>Study location indicating female Bathurst caribou range (2007–2016). Individual caribou locations are shown in light blue dots. Data courtesy of GNWT ECC.</p>
Full article ">Figure 2
<p>Modeled migration for 1 female caribou showing the fit of her migratory shift from <math display="inline"><semantics> <msub> <mi>m</mi> <mn>1</mn> </msub> </semantics></math> to <math display="inline"><semantics> <msub> <mi>m</mi> <mn>2</mn> </msub> </semantics></math> in 2011. (<b>A</b>) The x (longitude) and y (latitude) coordinates of the caribou between her wintering and calving grounds ranges. The light and dark blue circles show the 95th and 50th percentile for this individual’s estimated range area, and the dark blue line is showing the linear distance between the centroids of the two ranges. (<b>B</b>) Example of modeled migrations for 4 female caribou Migratory White Noise (MWN) model fits of their migratory shifts from wintering (<math display="inline"><semantics> <msub> <mi>m</mi> <mn>1</mn> </msub> </semantics></math>) to calving (<math display="inline"><semantics> <msub> <mi>m</mi> <mn>2</mn> </msub> </semantics></math>) grounds in 2011. The dotted ellipses represent the population-level ranges. (<b>C</b>,<b>D</b>) One year of longitudinal (<b>C</b>) and latitudinal (<b>D</b>) displacements of the individual. The light and dark blue lines on the time series plots are confidence intervals for the estimated means. A spring migration (significant northeastward displacement) can be observed around day 125. Based on Gurarie et al. [<a href="#B21-remotesensing-16-02391" class="html-bibr">21</a>].</p>
Full article ">Figure 3
<p>Time series of 2007 <math display="inline"><semantics> <msub> <mi>T</mi> <mi>B</mi> </msub> </semantics></math> and <math display="inline"><semantics> <mrow> <mi>D</mi> <mi>A</mi> <mi>V</mi> </mrow> </semantics></math> from a 3.125 km 37V GHz EASE-Grid (Equal-Area Scalable Earth-Grid) pixel encompassing the departure location (64.267°N, 116.213°W) of the Bathurst herd in 2007. (<b>A</b>) The brightness temperature minimum and maximum (black lines and dark gray shaded region between the high and low <math display="inline"><semantics> <msub> <mi>T</mi> <mi>B</mi> </msub> </semantics></math> values) and (<b>B</b>) extracted as the <math display="inline"><semantics> <mrow> <mi>D</mi> <mi>A</mi> <mi>V</mi> </mrow> </semantics></math> time series (blue line) shown in the bottom graph. In 2007, there were two early melt events that were detected by the algorithm, which are highlighted by pink vertical bars. The melt onset and period of high <math display="inline"><semantics> <mrow> <mi>D</mi> <mi>A</mi> <mi>V</mi> </mrow> </semantics></math> are annotated by black labeled vertical lines, and the melt–refreeze period is highlighted by a light gray vertical shaded bar. Dashed horizontal lines are thresholds used for <math display="inline"><semantics> <msub> <mi>T</mi> <mi>B</mi> </msub> </semantics></math> and <math display="inline"><semantics> <mrow> <mi>D</mi> <mi>A</mi> <mi>V</mi> </mrow> </semantics></math> in the algorithm.</p>
Full article ">Figure 4
<p>CETB observations and female caribou latitudinal movement (north–south migration) between 2007 and 2016. The top two graphs are of brightness temperature (<b>A</b>) and <math display="inline"><semantics> <mrow> <mi>D</mi> <mi>A</mi> <mi>V</mi> </mrow> </semantics></math> (<b>B</b>). The graph includes individual components of female caribou latitudinal movement through time (<b>C</b>). Individual female caribou are colored by unique collar ID. Collar ID legend has been omitted, as there are &gt;100 unique IDs within this time span. Number of individuals per year is noted below each year’s panel, where n = number of individuals.</p>
Full article ">Figure 5
<p>(<b>A</b>) Location map for Bathurst herd range. (<b>B</b>) Extent of each of the subpanels. All subsets are in geographic coordinates and have the same extent. (<b>C</b>) Time series of MODIS true color imagery with caribou locations, 2007–2016. From left to right in each panel, the date of the image is 30 March–6 April, 1–8 May, and 25 May–1 June. Caribou locations at the corresponding time are depicted by the orange circles in each image. White signifies ice or snow cover on these mostly clear composites, while green and brown show the land appearing as the snow melts away.</p>
Full article ">Figure 6
<p>Melt onset maps (early and main melt) for 2011. Panel (<b>A</b>) is an overview map of the area seen in panel (<b>B</b>). Panel (<b>B</b>) is showing early (<b>left</b>) and main (<b>right</b>) melt onset maps and female caribou departure and arrival locations in 2011. Individual departure and arrival locations are shown in orange dots, and the herd departure and arrival locations are shown as red location pins. The color bar in panel (<b>B</b>) represents the melt onset date in day of year.</p>
Full article ">Figure 7
<p>Time series of multiple datasets from the departure location (64.267°N, 116.213°W) of females in the Bathurst herd in 2007 and 2011. Upper images show snow cover from MODIS at key dates (<b>A</b>,<b>F</b>), acquired during the first early melt event, melt onset, and end of high <math display="inline"><semantics> <mrow> <mi>D</mi> <mi>A</mi> <mi>V</mi> </mrow> </semantics></math> period.The next two graphs are of brightness temperature (<b>B</b>,<b>G</b>) and <math display="inline"><semantics> <mrow> <mi>D</mi> <mi>A</mi> <mi>V</mi> </mrow> </semantics></math> (<b>C</b>,<b>H</b>). The brightness temperature minimum and maximum (black lines and dark gray shaded region) are shown in the top graph, and <math display="inline"><semantics> <mrow> <mi>D</mi> <mi>A</mi> <mi>V</mi> </mrow> </semantics></math> (blue lines) are shown in the lower graph. Early melt events are highlighted by pink vertical bars; melt onset and <math display="inline"><semantics> <mrow> <mi>D</mi> <mi>A</mi> <mi>V</mi> </mrow> </semantics></math> are annotated by black labeled vertical lines; and the melt–refreeze period is highlighted by a light gray vertical bar. The dashed horizontal lines are the thresholds used for brightness temperature and <math display="inline"><semantics> <mrow> <mi>D</mi> <mi>A</mi> <mi>V</mi> </mrow> </semantics></math> in the algorithms. The bottom two scatter plots are multiple daily readings of the displacement rate (speed) (<b>D</b>,<b>I</b>) and distance (<b>E</b>,<b>J</b>) traveled by individual female caribou in the Bathurst herd in each year. Female caribou movement events during spring migration (departure from wintering grounds and arrival at calving grounds) are annotated and highlighted by orange vertical bars. Similar plots for all years are shown in <a href="#remotesensing-16-02391-f0A1" class="html-fig">Figure A1</a>, <a href="#remotesensing-16-02391-f0A2" class="html-fig">Figure A2</a>, <a href="#remotesensing-16-02391-f0A3" class="html-fig">Figure A3</a>, <a href="#remotesensing-16-02391-f0A4" class="html-fig">Figure A4</a>, <a href="#remotesensing-16-02391-f0A5" class="html-fig">Figure A5</a>, <a href="#remotesensing-16-02391-f0A6" class="html-fig">Figure A6</a>, <a href="#remotesensing-16-02391-f0A7" class="html-fig">Figure A7</a>, <a href="#remotesensing-16-02391-f0A8" class="html-fig">Figure A8</a>, <a href="#remotesensing-16-02391-f0A9" class="html-fig">Figure A9</a> and <a href="#remotesensing-16-02391-f0A10" class="html-fig">Figure A10</a>.</p>
Full article ">Figure 8
<p>Plot of migration timing against snowmelt phenology for Bathurst caribou from 2007 to 2016. Black squares indicate the start of migration (departure from winter range), open circles indicate end of migration, and the black line indicates the duration. Triangles indicate time difference between the first early melt event (point) and melt onset (base). For the year when there is no triangle (2016), there were no early melt events. Orange and purple horizontal lines indicate melt onset at the winter and calving ground locations, respectively. Gray triangles indicate the snowmelt timing at, respectively, the wintering range: the bottom of the triangles represents the early melt onset, and the wider upper edges indicate the main melt onset date. Gray shading indicates the period of high <math display="inline"><semantics> <mrow> <mi>D</mi> <mi>A</mi> <mi>V</mi> </mrow> </semantics></math>, i.e., the melt–refreeze period on the wintering range.</p>
Full article ">Figure 9
<p>Herd departure date (<b>A</b>) and arrival date (<b>B</b>) and migration duration (<b>C</b>) vs. melt onset date, 2007–2016. Relationships between early and main melt onset dates and the female Bathurst herd’s spring migration departure (start) date. Points are colored by year. The trend line for early melt events is the dotted line and the trend line for main melt events is the solid line.</p>
Full article ">Figure A1
<p>Time series example of CETB and migration data in 2007. Please see above for detailed description of figure.</p>
Full article ">Figure A2
<p>Time series example of CETB and migration data in 2008. Please see above for detailed description of figure.</p>
Full article ">Figure A3
<p>Time series example of CETB and migration data in 2009. Please see above for detailed description of figure.</p>
Full article ">Figure A4
<p>Time series example of CETB and migration data in 2010. Please see above for detailed description of figure.</p>
Full article ">Figure A5
<p>Time series example of CETB and migration data in 2011. Please see above for detailed description of figure.</p>
Full article ">Figure A6
<p>Time series example of CETB and migration data in 2012. Please see above for detailed description of figure.</p>
Full article ">Figure A7
<p>Time series example of CETB and migration data in 2013. Please see above for detailed description of figure.</p>
Full article ">Figure A8
<p>Time series example of CETB and migration data in 2014. Please see above for detailed description of figure.</p>
Full article ">Figure A9
<p>Time series example of CETB and migration data in 2015. Please see above for detailed description of figure.</p>
Full article ">Figure A10
<p>Time series example of CETB and migration data in 2016. Please see above for detailed description of figure.</p>
Full article ">Figure A11
<p>Individual (small circles) herd-level (large circles) female Bathurst caribou migration timing by year compared to early and main melt events. Each movement event data point for departure from wintering grounds (closed circles) and arrival at calving grounds (open circles) corresponds to a melt onset date. Data points within the gray shaded region of the graph represent movement events that occurred on or after melt onset had been triggered in that region. The solid lines connecting each set of points represent one year of readings and connect the start and end dates of migration for a given year. Smaller, lighter dots are the results for individual animals. The open and closed circle data labels are the movement event dates for day of departure and arrival. (<b>A</b>) Early melt onset was computed using a requirement of one out of two consecutive days of experienced melt. (<b>B</b>) Main melt onset was computed using a melt detection requirement of five out of seven consecutive days of experienced melt. Melt onset is consistent over large areas, so many individual caribou appear distributed linearly, reflecting variation in departure or arrival, yet less variation in snow conditions experienced.</p>
Full article ">
22 pages, 33778 KiB  
Article
Synthetic Aperture Radar Monitoring of Snow in a Reindeer-Grazing Landscape
by Ida Carlsson, Gunhild Rosqvist, Jenny Marika Wennbom and Ian A. Brown
Remote Sens. 2024, 16(13), 2329; https://doi.org/10.3390/rs16132329 - 26 Jun 2024
Viewed by 1065
Abstract
Snow cover and runoff play an important role in the Arctic environment, which is increasingly affected by climate change. Over the past 30 years, winter temperatures in northern Sweden have risen by 2 °C, accompanied by an increase in precipitation. This has led [...] Read more.
Snow cover and runoff play an important role in the Arctic environment, which is increasingly affected by climate change. Over the past 30 years, winter temperatures in northern Sweden have risen by 2 °C, accompanied by an increase in precipitation. This has led to a higher incidence of thaw–freeze and rain-on-snow events. Snow properties, such as the snow depth and longevity, and the timing of snowmelt in spring significantly impact the alpine tundra vegetation. The emergent vegetation at the edge of the snow patches during spring and summer constitutes an essential nutrient supply for reindeer. We have used Sentinel-1 synthetic aperture radar (SAR) to determine the onset of the surface melt and the end of the snow cover in the core reindeer grazing area of the Laevás Sámi reindeer-herding community in northern Sweden. Using SAR data from March to August during the period 2017 to 2021, the start of the surface melt is identified by detecting the season’s backscatter minimum. The end of the snow cover is determined using a threshold approach. A comparison between the results of the analysis of the end of the snow cover from Sentinel-1 and in situ measurements, for the years 2017 to 2020, derived from an automatic weather station located in Laevásvággi reveals a 2- to 10-day difference in the snow-free ground conditions, which indicates that the method can be used to investigate when the ground is free of snow. VH data are preferred to VV data due to the former’s lower sensitivity to temporary wetting events. The outcomes from the season backscatter minimum demonstrate a distinct 25-day difference in the start of the runoff between the 5 investigated years. The backscatter minimum and threshold-based method used here serves as a valuable complement to global snowmelt monitoring. Full article
(This article belongs to the Section Ecological Remote Sensing)
Show Figures

Figure 1

Figure 1
<p>Modified illustration based on Buchelt et al. [<a href="#B36-remotesensing-16-02329" class="html-bibr">36</a>], describing the backscatter intensity during the snow melting season derived from S-1 SAR data. The start of the surface melt (SOSM) is marked by the S-1 backscatter reaching its minimum, while the end of the snowmelt (EOS) is indicated by the backscatter starting to reach a higher value after reaching the season’s minimum value [<a href="#B36-remotesensing-16-02329" class="html-bibr">36</a>].</p>
Full article ">Figure 2
<p>The area of interest for this study is the spring and summer grazing area used by reindeer of the Laevás Sámi reindeer-herding community, northern Sweden. The yellow circle marks where the automatic weather station (AWS) in Laevásvággi 18.96°E 68.04°N is located, and the area is also the calving ground for Laevás reindeers [<a href="#B42-remotesensing-16-02329" class="html-bibr">42</a>,<a href="#B44-remotesensing-16-02329" class="html-bibr">44</a>].</p>
Full article ">Figure 3
<p>The acquisitions of S-1 (blue) single-look complex data, in interferometric wide-swath mode, in ascending orbit from the Alaska Satellite Facility [<a href="#B45-remotesensing-16-02329" class="html-bibr">45</a>] downloaded on the 12th of September, 14th of December of 2022 and 6th of May 2024.</p>
Full article ">Figure 4
<p>Workflow for preprocessing the S-1 images from 2017–2021 using the Sentinel Application Platform (SNAP). This process involved obtaining the latest orbit file, splitting swaths, and debursting images, followed by merging them into a single coherent image. Calibration to the backscatter coefficient (<span class="html-italic">β</span><sup>0</sup>) was conducted according to Small (2011). Subsequent steps included multilooking, terrain flattening for pixel location rectification, and Lee Sigma speckle filtering [<a href="#B36-remotesensing-16-02329" class="html-bibr">36</a>] for noise removal. Radiometric terrain correction utilized the range-Doppler technique with a 2 m DEM [<a href="#B44-remotesensing-16-02329" class="html-bibr">44</a>].</p>
Full article ">Figure 5
<p>Simplification (the season backscatter is smoothed by the average backscattering during the season) of the seasonal backscatter in decibels from Sentinel-1 for all five years in both polarisations (<b>a</b>,<b>b</b>). VV polarisation (<b>a</b>) consistently exhibits higher backscatter values throughout the season compared to VH polarisation (<b>b</b>). VV polarisation is preferred for surface features and roughness, while VH polarisation is more suitable for detecting internal structure and volume scattering within targets like vegetation. Notably, in March 2017, there is a period characterized by lower backscatter values in the VV polarisation.</p>
Full article ">Figure 6
<p>Overview of monthly SOSM<sub><span class="html-italic">S</span>-1</sub> in the VV polarisation from 2017 to 2021 in the spring and summer grazing area of Laevás reindeer. In 2017 (<b>a</b>), a substantial SOSM<sub><span class="html-italic">S</span>-1</sub> was detected in the VV polarisation mode, covering 56% of the area in March. In 2018 (<b>b</b>), the SOSM<sub><span class="html-italic">S</span>-1</sub> started in April (red) and May (light yellow) and in May (light yellow) during 2019 (<b>c</b>). In 2020 (<b>d</b>), the SOSM<sub><span class="html-italic">S</span>-1</sub> started in May (light yellow) and June (light blue). During the year 2021 (<b>e</b>), there were large SOSM<sub><span class="html-italic">S</span>-1</sub> in May (light yellow) and in July (blue). These findings suggest that the SOSM<sub><span class="html-italic">S</span>-1</sub> varies significantly across the years, with the VV polarisation mode consistently exhibiting an earlier SOSM than the VH.</p>
Full article ">Figure 7
<p>Overview of the monthly SOSM<sub><span class="html-italic">S</span>-1</sub> in the VH polarisation over the observed years 2017 to 2021. In 2017 (<b>a</b>), the VH exhibited the largest snowmelt in May (light yellow) and June (light blue). In 2018 (<b>b</b>), a significant SOSM<sub><span class="html-italic">S</span>-1</sub> was displayed in April (red) and May (light yellow), as well in May (light yellow) in 2019 (<b>c</b>). Noteworthily, 2020 (<b>d</b>) exhibited pronounced SOSM<sub><span class="html-italic">S</span>-1</sub> peaks in May (light yellow) and June (light blue). In 2021 (<b>e</b>), the highest SOSM was recorded in May (light yellow). These findings highlight the variability in the seasonal snowmelt dynamics captured in the data.</p>
Full article ">Figure 8
<p>Monthly end of snowmelt (EOS<sub><span class="html-italic">S</span>-1</sub>) in the VV polarisation between the years 2017 and 2021 are shown in the figure. In the VV polarisation for 2017 (<b>a</b>), the deviation is evident, with the largest amount of EOS<sub><span class="html-italic">S</span>-1</sub> occurring in March, a pattern not observed in the VH polarisation (<a href="#remotesensing-16-02329-f009" class="html-fig">Figure 9</a>). The year 2018 (<b>b</b>) exhibits an early EOS<sub><span class="html-italic">S</span>-1</sub> in the VV polarisation, indicating bare ground in a significant portion of the area as early as May (light yellow). Moreover, 2019 (<b>c</b>) and 2020 (<b>d</b>) show similar EOS<sub><span class="html-italic">S</span>-1</sub> in June (light blue) and July (blue), while 2019 exhibits some earlier melting, particularly in April (red). In 2021 (<b>e</b>), the EOS<sub><span class="html-italic">S</span>-1</sub> occurrence was notable in May, with a more substantial presence observed in July. Moreover, there are areas within the region where data are not available, as evidenced across all the years.</p>
Full article ">Figure 9
<p>End of season (EOS) observations from 2017 to 2021, as depicted by the VH polarisation. The data reveal significant variations in the EOS percentages across different years and months, with certain trends standing out prominently. For instance, noticeable spikes in the snowmelt are observed in May and June across multiple years, indicating periods of accelerated melting. In 2017 (<b>a</b>), the EOS percentages remained consistently low throughout the observed months, with minimal snowmelt recorded with the largest EOS in June (light blue). In 2018 (<b>b</b>), snowmelt began to appear in April (red) and increased notably in May (light blue). In 2019 (<b>c</b>), the trend of increasing snowmelt continued into 2019, with May (light yellow) showcasing substantial melting percentages. Moreover, 2020 (<b>d</b>) witnessed a pronounced increase in snowmelt compared to previous years, particularly notable in May (light yellow) and June (light blue). In 2021 (<b>e</b>), the EOS percentages displayed a remarkable spike in May (light yellow), indicating a notably accelerated snowmelt compared to previous years.</p>
Full article ">
17 pages, 21390 KiB  
Article
Partial Oscillation Flow Control on Airfoil at Low Reynolds Numbers
by Guanxiong Li and Jingyu Wang
Appl. Sci. 2024, 14(11), 4762; https://doi.org/10.3390/app14114762 - 31 May 2024
Viewed by 542
Abstract
Among the critical factors contributing to the decline in the aerodynamic performance of near-space aircraft under low Reynolds number conditions, a significant one lies in the occurrence of laminar separation bubbles forming on the wings. Within the scope of this investigation, the primary [...] Read more.
Among the critical factors contributing to the decline in the aerodynamic performance of near-space aircraft under low Reynolds number conditions, a significant one lies in the occurrence of laminar separation bubbles forming on the wings. Within the scope of this investigation, the primary research methodology adopted involves utilizing an unsteady numerical simulation technique rooted in a spring-smoothed dynamic grid system. This study meticulously examines the aerodynamic attributes and flow patterns exhibited by an airfoil undergoing partial oscillation, thereby elucidating the underlying mechanisms through which such oscillations lead to enhanced lift and diminished drag forces. The outcomes of this research reveal that the imposition of partial oscillation engenders a noteworthy augmentation of 4.9% in the lift coefficient of the airfoil, concurrent with a substantial diminution of 15.3% in its drag coefficient when juxtaposed against the non-deforming counterpart. The oscillation frequency exerts a profound influence on both the onset location of transition and the extent of the laminar separation bubble’s development. As the oscillation frequency escalates, it follows an initial ascending trend in the lift coefficient of the airfoil, followed by a subsequent decline, whereas the drag coefficient exhibits an initial decrement prior to a rising tendency, thus indicating the existence of an optimal frequency point where the airfoil achieves its most favorable aerodynamic characteristics. It is observed that the flow control effects are optimally pronounced when the region subjected to partial oscillation is proximate to the airfoil’s leading edge or situated precisely at the centroid of the laminar separation bubble. Full article
(This article belongs to the Section Aerospace Science and Engineering)
Show Figures

Figure 1

Figure 1
<p>Grids for numerical simulation. (<b>a</b>) Grids of the calculation domain; (<b>b</b>) grids around the airfoil.</p>
Full article ">Figure 2
<p>Lift curves of NACA 0009.</p>
Full article ">Figure 3
<p>Schematic diagram of the water tunnel test system.</p>
Full article ">Figure 4
<p>Water tunnel test system and the model of airfoil. (<b>a</b>) Water tunnel test system; (<b>b</b>) test model of the airfoil.</p>
Full article ">Figure 5
<p>Calculation model of the partial oscillation on airfoil.</p>
Full article ">Figure 6
<p>Deformation of the airfoil in one cycle. (<b>a</b>) <span class="html-italic">t</span> = 0.2 <span class="html-italic">T</span>; (<b>b</b>) <span class="html-italic">t</span> = 0.4 <span class="html-italic">T</span>; (<b>c</b>) <span class="html-italic">t</span> = 0.6 <span class="html-italic">T</span>; (<b>d</b>) <span class="html-italic">t</span> = 0.8 <span class="html-italic">T</span>.</p>
Full article ">Figure 7
<p>Grid distribution near oscillation surface. (<b>a</b>) Initial grid; (<b>b</b>) refined grid.</p>
Full article ">Figure 8
<p>Structures of the laminar separation bubble.</p>
Full article ">Figure 9
<p>Laminar separation bubble obtained from the calculation and the test.</p>
Full article ">Figure 10
<p>Aerodynamic characteristics of the airfoil under different oscillation frequencies. (<b>a</b>) Lift and drag coefficient (<span class="html-italic">P</span> = 0.1<span class="html-italic">c</span>); (<b>b</b>) lift–drag ratio (<span class="html-italic">P</span> = 0.1<span class="html-italic">c</span>); (<b>c</b>) lift and drag coefficient (<span class="html-italic">P</span> = 0.6<span class="html-italic">c</span>); (<b>d</b>) lift–drag ratio (<span class="html-italic">P</span> = 0.6<span class="html-italic">c</span>).</p>
Full article ">Figure 11
<p>Pressure and viscous drag coefficients of the airfoil under different oscillation frequencies (<span class="html-italic">P</span> = 0.6<span class="html-italic">c</span>).</p>
Full article ">Figure 12
<p>Position of the transition and length of the laminar separation bubble on the airfoil under different oscillation frequencies (<span class="html-italic">P</span> = 0.6<span class="html-italic">c</span>). (<b>a</b>) Position of the transition; (<b>b</b>) length of the laminar separation bubble.</p>
Full article ">Figure 13
<p>Friction coefficient of the airfoil (<span class="html-italic">P</span> = 0.6<span class="html-italic">c</span>, <span class="html-italic">f</span> = 4 Hz).</p>
Full article ">Figure 14
<p>Shape of the LSB on the airfoil (<span class="html-italic">P</span> = 0.6<span class="html-italic">c</span>, <span class="html-italic">f</span> = 4 Hz).</p>
Full article ">Figure 15
<p>Influence of the oscillation position on aerodynamic characteristics of airfoil. (<b>a</b>) Lift and drag coefficient; (<b>b</b>) lift–drag ratio.</p>
Full article ">Figure 16
<p>Streamlines and velocity distributions of airfoil. (<b>a</b>) Rigid airfoil; (<b>b</b>) <span class="html-italic">P</span> = 0.1<span class="html-italic">c</span>; (<b>c</b>) <span class="html-italic">P</span> = 0.5<span class="html-italic">c</span>; (<b>d</b>) <span class="html-italic">P</span> = 0.6<span class="html-italic">c</span>.</p>
Full article ">Figure 17
<p>Pressure distributions on the upper surface of airfoils.</p>
Full article ">Figure 18
<p>Variation of aerodynamic force of airfoil with time. (<b>a</b>) Lift coefficient; (<b>b</b>) drag coefficient.</p>
Full article ">Figure 19
<p>Shapes of LSB at different times in one oscillation cycle (<span class="html-italic">P</span> = 0.6<span class="html-italic">c</span>). (<b>a</b>) <span class="html-italic">t</span> = 0<span class="html-italic">T</span>; (<b>b</b>) <span class="html-italic">t</span> = 0.2<span class="html-italic">T</span>; (<b>c</b>) <span class="html-italic">t</span> = 0.4<span class="html-italic">T</span>; (<b>d</b>) <span class="html-italic">t</span> = 0.6<span class="html-italic">T</span>; (<b>e</b>) <span class="html-italic">t</span> = 0.8<span class="html-italic">T</span>; (<b>f</b>) <span class="html-italic">t</span> = 1.0<span class="html-italic">T</span>.</p>
Full article ">
13 pages, 4013 KiB  
Article
Water Composition, Biomass, and Species Distribution of Vascular Plants in Lake Agmon-Hula (LAH) (1993–2023) and Nearby Surroundings: A Review
by Moshe Gophen
Water 2024, 16(10), 1450; https://doi.org/10.3390/w16101450 - 19 May 2024
Viewed by 1077
Abstract
A significant change to the land cover in the Hula Valley was carried out during the 1950s: A swampy area densely covered by aquatic vegetation and the old shallow lake Hula were drained. The natural shallow lake and swamps land cover were converted [...] Read more.
A significant change to the land cover in the Hula Valley was carried out during the 1950s: A swampy area densely covered by aquatic vegetation and the old shallow lake Hula were drained. The natural shallow lake and swamps land cover were converted into agricultural development land use in two stages: (1) Drainage that was accomplished in 1957; (2) Implementation of the renovated hydrological system structure, including the newly created shallow Lake Agmon-Hula (LAH), was completed in 2007. The long-term data record of the restored diversity of the submerged and emerged aquatic plant community, and its relation to water quality in the newly created shallow Lake Agmon-Hula LAH, was statistically evaluated. Internal interactions within the LAH ecosystem between aquatic plants and water quality, including nitrification, de-nitrification, sedimentation, photosynthetic intensity, and plant biomass and nutrient composition, were statistically evaluated. The plant community in LAH maintains a seasonal growth cycle of onset during late spring–summer and dieback accompanied by decomposed degradation during fall–early winter. The summer peak of aquatic plant biomass and consequent enhancement of photosynthetic intensity induces a pH increase during daytime and carbonate precipitation. Nevertheless, the ecosystem is aerobic and sulfate reduction and H2S concentration are negligible. The Hula reclamation project (HP) is aimed at the improvement of eco-tourism’s integration into management design. The vegetation research confirms habitat enrichment. Full article
(This article belongs to the Special Issue The Role of Vegetation in Freshwater Ecology)
Show Figures

Figure 1

Figure 1
<p>Geographical map of the Kinneret Drainage Basin (<b>left</b>) and the Hula Valley (<b>right</b>).</p>
Full article ">Figure 2
<p>Two-way quadratic prediction relationship analysis between nutrient (SO<sub>4</sub>, NO<sub>3</sub>, TN, TDS, TDN, NO<sub>2</sub>, NH<sub>4</sub>) concentrations, EC, DO, and season (month) in LAH. Results indicate high values in winter and low values in summer followed by a moderately (distinctly shaper for ammonium) increased trend in fall–early winter. All plots represent moderate, and for ammonium, “deep”, U-shaped relationships.</p>
Full article ">Figure 2 Cont.
<p>Two-way quadratic prediction relationship analysis between nutrient (SO<sub>4</sub>, NO<sub>3</sub>, TN, TDS, TDN, NO<sub>2</sub>, NH<sub>4</sub>) concentrations, EC, DO, and season (month) in LAH. Results indicate high values in winter and low values in summer followed by a moderately (distinctly shaper for ammonium) increased trend in fall–early winter. All plots represent moderate, and for ammonium, “deep”, U-shaped relationships.</p>
Full article ">Figure 3
<p>Two-way quadratic prediction relationship analysis between nutrient (TSS, TP, TDP,) concentrations, pH, NTU, ALK, and season (month) in LAH. Results indicate a high range of value variabilities and a summer increase of phosphorus (TP, TDP), pH, and TDS parameters whilst ALK and NTU represent U-shaped relationships with high variability.</p>
Full article ">Figure 4
<p>Plot of optimization of temporal (1994–2023) linear prediction (Linear Predictive Coding, LPC) w/95% confidence interval, where future values of discrete-time signals are estimated as a linear function of previously observed values of nutrients TP, TSS, TN, and TDP as well as DO and NTU during 1994–2023 in LAH. All nutrient concentrations and DO and NTU represent a trend of temporal elevation. A high range of annual variabilities of TSS is shown.</p>
Full article ">Figure 5
<p>Plot of optimization of temporal (1994–2023) linear prediction (Linear Predictive Coding, LPC) w/95% confidence interval, where future values of discrete-time signals are estimated as a linear function of previously observed values of nutrients, TDS, SO<sub>4</sub>, NO<sub>3</sub>, NH<sub>4</sub>, NO<sub>2</sub>, during 1994–2023 in LAH. All nutrient concentrations except those for NH4 and TDN represent a distinct trend of temporal decline. A high range of annual variabilities, with no change of annual temporal trend for NH4 and TDN, are indicated.</p>
Full article ">Figure 5 Cont.
<p>Plot of optimization of temporal (1994–2023) linear prediction (Linear Predictive Coding, LPC) w/95% confidence interval, where future values of discrete-time signals are estimated as a linear function of previously observed values of nutrients, TDS, SO<sub>4</sub>, NO<sub>3</sub>, NH<sub>4</sub>, NO<sub>2</sub>, during 1994–2023 in LAH. All nutrient concentrations except those for NH4 and TDN represent a distinct trend of temporal decline. A high range of annual variabilities, with no change of annual temporal trend for NH4 and TDN, are indicated.</p>
Full article ">
16 pages, 12022 KiB  
Article
Spatiotemporal Variations in Snow Cover on the Tibetan Plateau from 2003 to 2020
by Chaoxu Pu, Shuaibo Zhou, Peijun Sun, Yunchuan Luo, Siyi Li and Zhangli Sun
Water 2024, 16(10), 1364; https://doi.org/10.3390/w16101364 - 11 May 2024
Viewed by 869
Abstract
The variations in snow cover on the Tibetan Plateau play a pivotal role in comprehending climate change patterns and governing hydrological processes within the region. This study leverages daily snow cover data and the NASA Digital Elevation Model (DEM) from 2003 to 2020 [...] Read more.
The variations in snow cover on the Tibetan Plateau play a pivotal role in comprehending climate change patterns and governing hydrological processes within the region. This study leverages daily snow cover data and the NASA Digital Elevation Model (DEM) from 2003 to 2020 to analyze spatiotemporal snow cover days and assess their responsiveness to climatic shifts by integrating meteorological data. The results reveal significant spatial heterogeneity in snow cover across the Plateau, with a slight decreasing trend in annual average snow cover duration. Snow cover is predominantly observed during the spring and winter seasons, constituting approximately 32% of the total snow cover days annually. The onset and cessation of snow cover occur within a range of 120–220 days. Additionally, an increasing trend in snow cover duration below 5000 m altitude was observed, in addition to a decreasing trend above 5000 m altitude. Sub-basin analysis delineates the Tarim River Basin as exhibiting the lengthiest average annual snow cover duration of 83 days, while the Yellow River Basin records the shortest duration of 31 days. The decreasing trend in snow cover duration closely aligns with climate warming trends, characterized by a warming rate of 0.17 ± 0.54 °C per decade, coupled with a concurrent increase in precipitation at a rate of 3.09 ± 3.81 mm per year. Temperature exerts a more pronounced influence on annual snow cover duration variation compared to precipitation, as evidenced by a strong negative correlation (CC = −0.67). This study significantly augments the comprehension of hydrological cycle dynamics on the Tibetan Plateau, furnishing essential insights for informed decision-making in water resource management and ecological conservation efforts. Full article
Show Figures

Figure 1

Figure 1
<p>Study area of the Tibetan Plateau and ten major basins, the red star symbol in the upper right corner represents the location of the Tibetan Plateau.</p>
Full article ">Figure 2
<p>Spatial distribution (<b>a</b>) and temporal (<b>b</b>) characteristics of annual snow cover duration on the Tibetan Plateau (2003–2020 hydrological years). Both the left and right axes of the line chart (right) represent the number of snowfall days, and the left <span class="html-italic">y</span>-axis represents snowfall days in a particular season, while the right <span class="html-italic">y</span>-axis represents snowfall days in a hydrological year. The red dotted line represents the line trend of annual average snowfall days (non-significant).</p>
Full article ">Figure 3
<p>The spatial distribution of the SOD (<b>a</b>) and SED (<b>b</b>) across the Tibetan Plateau.</p>
Full article ">Figure 4
<p>The variations in SCD in different altitude bins of the Tibetan Plateau during the hydrological years of 2003–2020. The SOD (subfigure(<b>f</b>)) shows a significant increasing trend (<span class="html-italic">p</span> &lt; 0.05, F test), and the other trends are non-significant.</p>
Full article ">Figure 5
<p>The characteristics of SCD changes in the Tibetan Plateau from 2003 to 2020 (hydrological year). Both the left and right axes of the line charts represent the number of snowfall days, and the left <span class="html-italic">y</span>-axis represents snowfall days in a particular season, while the right <span class="html-italic">y</span>-axis represents snowfall days in a hydrological year. The trend of SCD in the Tarim Basin (Autumn) is significantly decreasing (<span class="html-italic">p</span> &lt; 0.10, F test), and trends of SCD in the Tsaidam Basin (Spring), Yellow River Basin (Spring), and Hexi Corridor Basin (Summer) are significantly increasing (<span class="html-italic">p</span> &lt; 0.10, F test). The other trends of SCD are non-significant.</p>
Full article ">Figure 6
<p>The spatial distribution of SCD across seasons from 2003 to 2020 (hydrological year) on the Tibetan Plateau.</p>
Full article ">Figure 7
<p>Spatial distribution of temperature (<b>a</b>) and precipitation (<b>b</b>) on the Tibetan Plateau from 2003 to 2020 (hydrological year), and the trend changes in snow cover days, temperature, and precipitation variation (<b>c</b>). Both the precipitation and temperature show increasing trends, with the former presenting statistically significant (<span class="html-italic">p</span> &lt; 0.05, F test), while the snowfall showed a decreasing trend (non-significant).</p>
Full article ">Figure 8
<p>The trend changes in seasonal snow cover days and meteorological factors on the Tibetan Plateau during the hydrological years from 2003 to 2020.</p>
Full article ">
13 pages, 1488 KiB  
Article
Seasonality Affects Low-Molecular-Weight Organic Acids and Phenolic Compounds’ Composition in Scots Pine Litterfall
by Anna Ilek, Monika Gąsecka, Zuzanna Magdziak, Costas Saitanis and Courtney M. Siegert
Plants 2024, 13(10), 1293; https://doi.org/10.3390/plants13101293 - 8 May 2024
Viewed by 1016
Abstract
Background and Aims: Secondary plant metabolites, including organic acids and phenolic compounds, have a significant impact on the properties of organic matter in soil, influencing its structure and function. How the production of these compounds in foliage that falls to the forest floor [...] Read more.
Background and Aims: Secondary plant metabolites, including organic acids and phenolic compounds, have a significant impact on the properties of organic matter in soil, influencing its structure and function. How the production of these compounds in foliage that falls to the forest floor as litterfall varies across tree age and seasonality are of considerable interest for advancing our understanding of organic matter dynamics. Methods: Monthly, we collected fallen needles of Scots pine (Pinus sylvestris L.) across stands of five different age classes (20, 40, 60, 80, and 100 years) for one year and measured the organic acids and phenolic compounds. Results: Seven low-molecular-weight organic acids and thirteen phenolic compounds were detected in the litterfall. No differences were observed across stand age. Significant seasonal differences were detected. Most compounds peaked during litterfall in the growing season. Succinic acid was the most prevalent organic acid in the litterfall, comprising 78% of total organic acids (351.27 ± 34.27 µg g− 1), and was 1.5 to 11.0 times greater in the summer than all other seasons. Sinapic acid was the most prevalent phenolic compound in the litterfall (42.15 µg g− 1), representing 11% of the total phenolic compounds, and was 39.8 times greater in spring and summer compared to autumn and winter. Growing season peaks in needle concentrations were observed for all thirteen phenolic compounds and two organic acids (lactic, succinic). Citric acid exhibited a definitive peak in late winter into early spring. Conclusions: Our results highlight the seasonal dynamics of the composition of secondary plant metabolites in litterfall, which is most different at the onset of the growing season. Fresh inputs of litterfall at this time of emerging biological activity likely have seasonal impacts on soil’s organic matter composition as well. Full article
(This article belongs to the Special Issue Physiological and Biochemical Responses to Abiotic Stresses in Plants)
Show Figures

Figure 1

Figure 1
<p>Organic acid content in Scots pine needles collected across four seasons (mean ± SE), where: (<b>A</b>) sum of organic acids, (<b>B</b>) succinic acid, (<b>C</b>) lactic acid, (<b>D</b>) citric acid, (<b>E</b>) malic acid, (<b>F</b>) malonic acid, (<b>G</b>) fumaric acid, (<b>H</b>) acetic acid. The same letters indicate no statistical differences in the sum of organic acid contents (red color) and individual acids between seasons (Kruskal–Wallis test, <span class="html-italic">p</span> &lt; 0.05).</p>
Full article ">Figure 2
<p>Phenolic compound contents in Scots pine needles collected across four seasons (mean ± SE), where: (<b>A</b>) 2,5-DHBA acid content, (<b>B</b>) caffeic acid, (<b>C</b>) coumaric acid, (<b>D</b>) protocatechuic acid, (<b>E</b>) ferulic acid, (<b>F</b>) t-cinnamic acid, (<b>G</b>) chlorogenic acid, (<b>H</b>) 4-HBA acid, (<b>I</b>) catechin acid, (<b>J</b>) gallic acid, (<b>K</b>) syringic acid, (<b>L</b>) vanillic acid, and (<b>M</b>) sinapic acid. The same letters indicate no statistical differences in the individual phenolic compounds between seasons (Kruskal–Wallis test, <span class="html-italic">p</span> &lt; 0.05).</p>
Full article ">Figure 3
<p>Spearman’s (non-parametric) correlation coefficient (r) matrix of the measured parameters (phenolic compounds and organic acids). Correlation coefficient values are shown in cells where the relationship is significant (<span class="html-italic">p</span> &lt; 0.05). Blank cells represent no significant correlation. The corresponding color scale denotes positive (blue-colored tones) and negative (red-colored tones) correlation levels.</p>
Full article ">
22 pages, 15666 KiB  
Article
Adaptability of Prunus cerasifera Ehrh. to Climate Changes in Multifunctional Landscape
by Djurdja Petrov, Mirjana Ocokoljić, Nevenka Galečić, Dejan Skočajić and Isidora Simović
Atmosphere 2024, 15(3), 335; https://doi.org/10.3390/atmos15030335 - 8 Mar 2024
Cited by 3 | Viewed by 1187
Abstract
Urban trees play a vital role in mitigating climate changes, maintaining the sustainability of ecosystems. This study focuses on the assessment of the resilience of cherry plums to climate changes, a fruit-bearing species that offers diverse ecosystem services within multifunctional urban and suburban [...] Read more.
Urban trees play a vital role in mitigating climate changes, maintaining the sustainability of ecosystems. This study focuses on the assessment of the resilience of cherry plums to climate changes, a fruit-bearing species that offers diverse ecosystem services within multifunctional urban and suburban landscapes. This study examines flowering and fruiting in the context of climate characteristics, expressed through the Day of the Year (DOY), Growing Degree Days (GDDs), and a yield over 17 consecutive years. The results indicate significant shifts in the DOY but not in the GDD, apart from the end of flowering. The onset of flowering was earlier and the end postponed, extending the phenophase by an average of 4 days. The cherry plum’s yield was unaffected by climate changes, including extreme events like a late-spring frost. The stability of the cherry plum was confirmed by the phenological patterns of the bullace (cherry plum and blackthorn hybrid) exhibiting repeated flowering in the warmest year of 2023. The cherry plum is an adaptive species, with a high adaptability to a changing climate and a high resistance to late-spring frosts; thus, it is a favorable choice in urban design and planning, demonstrating resilience to climate shifts and thriving in polluted urban environments. It is especially appreciated for multiple ecosystem services: biodiversity conservation in natural and semi-natural areas, yielding good provisions in challenging environments, and the preservation of ornamental values through an extended flowering phenophase. Full article
(This article belongs to the Special Issue Climate Change Impacts and Adaptation Strategies in Agriculture)
Show Figures

Figure 1

Figure 1
<p>Ecosystem services of <span class="html-italic">Prunus cerasifera</span> Ehrh., with the underlined topics covered in this paper.</p>
Full article ">Figure 2
<p>The research area— the municipality of Čukarica (Belgrade, Serbia).</p>
Full article ">Figure 3
<p>A 3D surface plot of the relative relief values derived over the study area.</p>
Full article ">Figure 4
<p>Mean seasonal air temperature and total precipitation for autumn (<b>a</b>), winter (<b>b</b>), spring (<b>c</b>), and summer (<b>d</b>) and their respective terciles for Belgrade (black lines correspond to the percentiles of the reference period), in the research period, in relation to the reference period 1991–2020.</p>
Full article ">Figure 4 Cont.
<p>Mean seasonal air temperature and total precipitation for autumn (<b>a</b>), winter (<b>b</b>), spring (<b>c</b>), and summer (<b>d</b>) and their respective terciles for Belgrade (black lines correspond to the percentiles of the reference period), in the research period, in relation to the reference period 1991–2020.</p>
Full article ">Figure 4 Cont.
<p>Mean seasonal air temperature and total precipitation for autumn (<b>a</b>), winter (<b>b</b>), spring (<b>c</b>), and summer (<b>d</b>) and their respective terciles for Belgrade (black lines correspond to the percentiles of the reference period), in the research period, in relation to the reference period 1991–2020.</p>
Full article ">Figure 5
<p>Mean monthly air temperature and sum of precipitation for February (<b>a</b>) and March (<b>b</b>) and their corresponding terciles for Belgrade (black lines correspond to the percentiles of the reference period), in the research period, in relation to the reference period 1991–2020.</p>
Full article ">Figure 6
<p>The DOY, for the period of 2007–2023, for the beginning of flowering (BF), full flowering (FF), and the end of flowering (EF) and associated linear trends.</p>
Full article ">Figure 7
<p>GDDs (°C), for the period of 2007–2023, for the beginning of flowering (BF), full flowering (FF), and the end of flowering (EF) and associated linear trends.</p>
Full article ">Figure 8
<p>Graphic representation of Spearman’s correlation coefficients for the DOY and GDDs for cherry plums, for the period of 2007–2023, in the study area.</p>
Full article ">Figure 9
<p>A phenogram of the cherry plum and bullace from January to November 2023 in the study area.</p>
Full article ">
19 pages, 1544 KiB  
Article
The Impact of a Six-Hour Light–Dark Cycle on Wheat Ear Emergence, Grain Yield, and Flour Quality in Future Plant-Growing Systems
by Helena Clauw, Hans Van de Put, Abderahman Sghaier, Trui Kerkaert, Els Debonne, Mia Eeckhout and Kathy Steppe
Foods 2024, 13(5), 750; https://doi.org/10.3390/foods13050750 - 28 Feb 2024
Cited by 3 | Viewed by 1442
Abstract
Cultivating wheat (Triticum aestivum) in a closed environment offers applications in both indoor farming and in outer-space farming. Tailoring the photoperiod holds potential to shorten the growth cycle, thereby increasing the annual number of cycles. As wheat is a long-day plant, [...] Read more.
Cultivating wheat (Triticum aestivum) in a closed environment offers applications in both indoor farming and in outer-space farming. Tailoring the photoperiod holds potential to shorten the growth cycle, thereby increasing the annual number of cycles. As wheat is a long-day plant, a night shorter than a critical length is required to induce flowering. In growth chambers, experiments were conducted to examine the impact of a 6 h light–dark cycle on the timing of wheat ear emergence, grain yield, and flour quality. Under equal daily light-integral conditions, the 6 h light–dark cycle promoted growth and development, resulting in accelerated ear emergence when compared to a 12 h cycle, additionally indicating that 12 h of darkness was excessive. To further stimulate heading and increase yield, the 6 h cycle was changed at the onset of stem elongation to a 14 h–10 h, mimicking spring conditions, and maintained until maturity. This successful transition was then combined with two levels of light intensity and nutrient solution, which did not significantly impact yield, while tillering and grain ripening did increase under higher light intensities. Moreover, it enabled manipulation of the baking quality, although lower-end falling numbers were observed. In conclusion, combining a 6 h light–dark cycle until stem elongation with a 14 h–10 h cycle presents a promising strategy for increasing future wheat production in closed environments. The observation of low falling numbers underscores the importance of factoring in flour quality when designing the wheat-growing systems of the future. Full article
(This article belongs to the Section Grain)
Show Figures

Figure 1

Figure 1
<p>Total dry weight of grains per plant in EXP1 (<b>A</b>), EXP2 (<b>B</b>), and EXP3 (<b>C</b>). HL, LL, HN, and LN refer to ‘high light intensity’, ‘low light intensity’, ‘high nutrient concentration’, and ‘low nutrient concentration’, respectively. Different letters indicate statistical differences (<span class="html-italic">p</span> &lt; 0.05) within the experiments.</p>
Full article ">Figure 2
<p>Number of shoots with grains in EXP1 (<b>A</b>), EXP2 (<b>B</b>), and EXP3 (<b>C</b>); number of grains per ear for EXP1 (<b>D</b>), EXP2 (<b>E</b>), and EXP3 (<b>F</b>); and dry weight per grain (g) for EXP1 (<b>G</b>), EXP2 (<b>H</b>), and EXP3 (<b>I</b>). HL, LL, HN, and LN refer to ‘high light intensity’, ‘low light intensity’, ‘high nutrient concentration’, and ‘low nutrient concentration’, respectively. Different letters indicate statistical differences (<span class="html-italic">p</span> &lt; 0.05) within the experiments.</p>
Full article ">Figure 3
<p>Total number of shoots for EXP2 (<b>A</b>) and EXP3 (<b>B</b>). HL, LL, HN, and LN refer to ‘high light intensity’, ‘low light intensity’, ‘high nutrient concentration’, and ‘low nutrient concentration’, respectively. Different letters indicate statistical differences (<span class="html-italic">p</span> &lt; 0.05) within the experiments.</p>
Full article ">Figure 4
<p>Average dry matter content distribution of the grains on an ear for all treatments in EXP1, EXP2, and EXP3. HL, LL, HN, and LN refer to ‘high light intensity’, ‘low light intensity’, ‘high nutrient concentration’, and ‘low nutrient concentration’, respectively. The dry matter contents, from dark to light brown, are &gt;85%, 70–85%, 50–70%, and &lt;50%.</p>
Full article ">Figure 5
<p>SPAD of flag leaves for all treatments in EXP3 on different days after sowing (DAS) during grain filling and ripening: 68 DAS (<b>A</b>) (grain filling), 97 DAS (<b>B</b>) (grain filling and ripening), and 113 DAS (<b>C</b>) (ripening). HL, LL, HN, and LN refer to ‘high light intensity’, ‘low light intensity’, ‘high nutrient concentration’, and ‘low nutrient concentration’, respectively. Different letters indicate statistical differences (<span class="html-italic">p</span> &lt; 0.05) within the experiments.</p>
Full article ">
21 pages, 7500 KiB  
Article
Analysis of Agricultural Drought Evolution Characteristics and Driving Factors in Inner Mongolia Inland River Basin Based on Three-Dimensional Recognition
by Zezhong Zhang, Hengzhi Guo, Kai Feng, Fei Wang, Weijie Zhang and Jian Liu
Water 2024, 16(3), 440; https://doi.org/10.3390/w16030440 - 29 Jan 2024
Cited by 2 | Viewed by 1368
Abstract
Agricultural drought events have become more frequent in the Inner Mongolia inland river basin in recent years, and the spatio-temporal evolution characteristics and development rules can be accurately and comprehensively understood using the three-dimensional identification method. In this paper, standardized soil moisture index [...] Read more.
Agricultural drought events have become more frequent in the Inner Mongolia inland river basin in recent years, and the spatio-temporal evolution characteristics and development rules can be accurately and comprehensively understood using the three-dimensional identification method. In this paper, standardized soil moisture index (SSMI) was used to characterize agricultural drought, and modified Mann–Kendall trend test (MMK) and 3D recognition of drought events were used to analyze the spatio-temporal evolution characteristics of agricultural drought events in this basin and reveal the drought development law. The relationships between drought and temperature (T), precipitation (P), evapotranspiration (E), and humidity (H) were analyzed using a cross-wavelet method. The results are as follows: (1) When the time scale of agricultural drought was short (monthly scale), the alternations of dry and wet were frequent, but the SSMI index of all scales showed a downward trend; (2) The spatial distribution characteristics of drought change trend in four seasons were similar, but the area with a significant downward trend of drought in spring was the largest, and the area of high frequency region was also the largest, and the drought trend was the most obvious; (3) The most serious agricultural drought event occurred from October 2000 to May 2002, and reached its maximum value in September 2001 (drought area and drought severity of 2.26 × 105 km2 and 3.61 × 105 months·km2, respectively), which mainly experienced five processes—drought onset–intensification–decay–re-intensification–termination—and the migration path of the drought center showed the characteristics of southwest–northeast transmission; (4) All the four meteorological factors were correlated with SSMI, and P had a greater impact on SSMI. This article aims to reveal the spatio-temporal evolution of agricultural drought events in the Inner Mongolia inland river basin, and provide a new way to accurately evaluate the spatio-temporal evolution of drought. Full article
(This article belongs to the Special Issue Drought Monitoring and Risk Assessment)
Show Figures

Figure 1

Figure 1
<p>Basic situation in the study area. UBB, Urat Rear Banner; UMB, Urat Middle Banner; DM, Darhan Muminggan United Banner; SG, Shiguai District; WC, Wuchuan County; SZW, Siziwang Banner; CRW, Qahar youyi houqi; EL, Erenhot City; SRB, Sonid Right Banner; XB, Xianghuang Banner; SD, Shangdu County; XH, Xilinhot City; SLB, Sonid Left Banner; ZB, Zhengxiangbai Banner; ZLB, Zhenglan Banner; AB, Abaga Banner; XH, Xinghe County; KT, Hexigten Banner; WU, Wast Ujumuqin Banner; EU, East Ujumuqin Banner.</p>
Full article ">Figure 2
<p>Schematic diagram of drought patch recognition.</p>
Full article ">Figure 3
<p>Schematic diagram of time–history connection of arid patches.</p>
Full article ">Figure 4
<p>Temporal evolution characteristics of multi-scale SSMI in IMIRB from 1960 to 2021.</p>
Full article ">Figure 5
<p>Monthly, seasonal, and annual scales SSMI time trends in IMIRB from 1960 to 2021. (<b>a</b>) Month. (<b>b</b>) Spring. (<b>c</b>) Summer. (<b>d</b>) Autumn. (<b>e</b>) Winter. (<b>f</b>) Year.</p>
Full article ">Figure 6
<p>Variation characteristics of seasonal drought intensity and area proportion in IMIRB from 1960 to 2021. (<b>a</b>) Spring. (<b>b</b>) Summer. (<b>c</b>) Autumn. (<b>d</b>) Winter.</p>
Full article ">Figure 7
<p>Spatial distribution characteristics of seasonal SSMI trends in IMIRB. (<b>a</b>) Spring. (<b>b</b>) Summer. (<b>c</b>) Autumn. (<b>d</b>) Winter.</p>
Full article ">Figure 8
<p>Spatial distribution characteristics of seasonal drought intensity of IMIRB during 1960–2021. (<b>a</b>) Spring. (<b>b</b>) Summer. (<b>c</b>) Autumn. (<b>d</b>) Winter.</p>
Full article ">Figure 9
<p>Spatial distribution characteristics of seasonal drought frequency in IMIRB from 1960 to 2021. (<b>a</b>) Spring. (<b>b</b>) Summer. (<b>c</b>) Autumn. (<b>d</b>) Winter.</p>
Full article ">Figure 10
<p>Three-dimensional perspective of the 36th agricultural drought and temporal trend of characteristic variables.</p>
Full article ">Figure 11
<p>Temporal and spatial dynamic evolution of typical agricultural drought events (<b>a</b>–<b>t</b>): October 2000–May 2002.</p>
Full article ">Figure 12
<p>Migration path of typical agricultural drought center.</p>
Full article ">Figure 13
<p>Cross-wavelet energy spectrum and condensation spectrum of SSMI with T, P, E, and H. (<b>a</b>) XWT: SSMI-T; (<b>b</b>) WTC: SSMI-T; (<b>c</b>) XWT: SSMI-P; (<b>d</b>) WTC: SSMI-P; (<b>e</b>) XWT: SSMI-E; (<b>f</b>) WTC: SSMI-E; (<b>g</b>) XWT: SSMI-H; (<b>h</b>) WTC: SSMI-H.</p>
Full article ">Figure 14
<p>The affected area of crops and the area of crop failure from 1998 to 2021.</p>
Full article ">
Back to TopTop