[go: up one dir, main page]

 
 
Sign in to use this feature.

Years

Between: -

Subjects

remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline

Journals

Article Types

Countries / Regions

Search Results (65)

Search Parameters:
Keywords = MCD12Q1

Order results
Result details
Results per page
Select all
Export citation of selected articles as:
19 pages, 9008 KiB  
Article
The Carpathian Agriculture in Poland in Relation to Other EU Countries, Ukraine and the Environmental Goals of the EU CAP 2023–2027
by Marek Zieliński, Artur Łopatka, Piotr Koza and Barbara Gołębiewska
Agriculture 2024, 14(8), 1325; https://doi.org/10.3390/agriculture14081325 - 9 Aug 2024
Viewed by 560
Abstract
This study discusses the issue of determining the direction and strength of changes taking place in the structure of agricultural land in the mountain and foothill areas of the Carpathians in Poland in comparison with Slovakia, Romania and Ukraine. The most important financial [...] Read more.
This study discusses the issue of determining the direction and strength of changes taking place in the structure of agricultural land in the mountain and foothill areas of the Carpathians in Poland in comparison with Slovakia, Romania and Ukraine. The most important financial institutional measures dedicated to the protection of the natural environment in Polish agriculture in the Areas facing Natural and other specific Constraints (ANCs) mountain and foothill in the first year of the CAP 2023–2027 were also established. Satellite data from 2001 to 2022 were used. The analyses used the land use classification MCD12Q1 provided by NASA and were made on the basis of satellite imagery collections from the MODIS sensor placed on two satellites: TERRA and AQUA. In EU countries, a decreasing trend in agricultural areas has been observed in areas below 350 m above sea level. In areas above 350 m, this trend weakened or even turned into an upward trend. Only in Ukraine was a different trend observed. It was found that in Poland, the degree of involvement of farmers from mountain and foothill areas in implementing financial institutional measures dedicated to protecting the natural environment during the study period was not satisfactory. Full article
Show Figures

Figure 1

Figure 1
<p>Scheme of the analysis of agriculture within separate groups of communes due to the fact and nuisance of ANCs mountain and foothill in Poland. Source: own study.</p>
Full article ">Figure 2
<p>Distribution of communes with different shares of ANCs mountain and foothill in Poland. Source: own study ISSPC SRI; IAFE NRI.</p>
Full article ">Figure 3
<p>Land use in the Carpathians in 2001 and 2022. Source: own study based on MODIS.</p>
Full article ">Figure 4
<p>Trends in the percentage share [%] of the total agricultural area and cropland in the total area of land in the Carpathians in 2001–2022. Source: own study based on MODIS.</p>
Full article ">Figure 5
<p>Number of farms participating in practices under eco-schemes, in organic and agri–environment–climate measures in communes with different shares of ANCs mountain and foothill in Poland in 2023. Source: own study based on ARMA.</p>
Full article ">Figure 6
<p>Share of [%] farms with eco-schemes in total number of farms in communes with ANCs mountain and foothill in 2023. Source: own study based on ARMA.</p>
Full article ">Figure 7
<p>Share of [%] farms with organic and agri–environmental–climate measure in total number of farms in communes with ANCs mountain and foothill in 2023.</p>
Full article ">Figure 8
<p>Agricultural area covered by practices under eco-schemes, ecological and agri–environment–climate measures in communes with different shares of ANCs mountain and foothill in Poland in 2023. Source: own study based on ARMA.</p>
Full article ">Figure 9
<p>Share [%] of UAA in farms with eco-schemes in total UAA in communes with ANCs mountain and foothill in 2023.</p>
Full article ">Figure 10
<p>Share [%] of UAA covered by organic and agri–environmental–climate measures in total UAA in communes with ANCs mountain and foothill in 2023.</p>
Full article ">
17 pages, 5254 KiB  
Article
Gynostemma pentaphyllum Extract Alleviates NASH in Mice: Exploration of Inflammation and Gut Microbiota
by Feng-Yan Jiang, Si-Ran Yue, Yi-Yun Tan, Nan Tang, Yue-Song Xu, Bao-Jun Zhang, Yue-Jian Mao, Zheng-Sheng Xue, Ai-Ping Lu, Bao-Cheng Liu and Rui-Rui Wang
Nutrients 2024, 16(11), 1782; https://doi.org/10.3390/nu16111782 - 6 Jun 2024
Viewed by 1164
Abstract
NASH (non-alcoholic steatohepatitis) is a severe liver disease characterized by hepatic chronic inflammation that can be associated with the gut microbiota. In this study, we explored the therapeutic effect of Gynostemma pentaphyllum extract (GPE), a Chinese herbal extract, on methionine- and choline-deficient (MCD) [...] Read more.
NASH (non-alcoholic steatohepatitis) is a severe liver disease characterized by hepatic chronic inflammation that can be associated with the gut microbiota. In this study, we explored the therapeutic effect of Gynostemma pentaphyllum extract (GPE), a Chinese herbal extract, on methionine- and choline-deficient (MCD) diet-induced NASH mice. Based on the peak area, the top ten compounds in GPE were hydroxylinolenic acid, rutin, hydroxylinoleic acid, vanillic acid, methyl vanillate, quercetin, pheophorbide A, protocatechuic acid, aurantiamide acetate, and iso-rhamnetin. We found that four weeks of GPE treatment alleviated hepatic confluent zone inflammation, hepatocyte lipid accumulation, and lipid peroxidation in the mouse model. According to the 16S rRNA gene V3–V4 region sequencing of the colonic contents, the gut microbiota structure of the mice was significantly changed after GPE supplementation. Especially, GPE enriched the abundance of potentially beneficial bacteria such as Akkerrmansia and decreased the abundance of opportunistic pathogens such as Klebsiella. Moreover, RNA sequencing revealed that the GPE group showed an anti-inflammatory liver characterized by the repression of the NF-kappa B signaling pathway compared with the MCD group. Ingenuity Pathway Analysis (IPA) also showed that GPE downregulated the pathogen-induced cytokine storm pathway, which was associated with inflammation. A high dose of GPE (HGPE) significantly downregulated the expression levels of the tumor necrosis factor-α (TNF-α), myeloid differentiation factor 88 (Myd88), cluster of differentiation 14 (CD14), and Toll-like receptor 4 (TLR4) genes, as verified by real-time quantitative PCR (RT-qPCR). Our results suggested that the therapeutic potential of GPE for NASH mice may be related to improvements in the intestinal microenvironment and a reduction in liver inflammation. Full article
(This article belongs to the Special Issue Prebiotics, Probiotics, and Gut Microbiota with Chronic Disease)
Show Figures

Figure 1

Figure 1
<p>Analysis of GPE compositions. Base peak chromatogram (BPC) of GPE in the negative (<b>A</b>) and positive (<b>B</b>) ion modes.</p>
Full article ">Figure 2
<p>GPE ameliorates liver injury in MCD diet-induced NASH mice. (<b>A</b>) H&amp;E-staining (200 magnifications, 100 μM) and Oil Red (200 magnifications, 100 μM) of liver sections; (<b>B</b>) Non-alcoholic Fatty Liver Disease Activity Score (NAS); (<b>C</b>) Sirius Red staining (200 magnifications, 100 μM) of liver sections; (<b>D</b>) quantification of the percentage of the Sirius-Red-positive area; (<b>E</b>) hepatic hydroxyproline quantification (HYP); (<b>F</b>–<b>M</b>) serum total cholesterol (TC), serum triglycerides (TG), serum alanine aminotransferase (ALT), serum aspartate aminotransferase (AST), serum high-density lipoprotein cholesterol (HDL-C), serum low-density lipoprotein cholesterol (LDL-C), hepatic TC, and hepatic TG; and (<b>N</b>,<b>O</b>) hepatic superoxide dismutase and malondialdehyde quantification (SOD and MDA). Each biological replicate (n = 5–6) had three technical replicate wells for the experiment. Data are expressed as means ± SEM. # <span class="html-italic">p</span> &lt; 0.05, ## <span class="html-italic">p</span> &lt; 0.01 vs. the MCS group; * <span class="html-italic">p</span> &lt; 0.05, ** <span class="html-italic">p</span> &lt; 0.01 vs. the MCD group.</p>
Full article ">Figure 3
<p>GPE alters the structure of the intestinal microbiota. (<b>A</b>–<b>D</b>) The alpha diversity analysis of sequence reads. (<b>A</b>) Chao1 index. (<b>B</b>) Shannon index. (<b>C</b>) Simpson. (<b>D</b>) ACE index. (<b>E</b>) Principal co-ordinate analysis at the OTU level based on the Bray–Curtis distance. (<b>F</b>) Fecal microbiota at the genus level. (<b>G</b>) Co-abundance groups interaction network of genus level among all groups based on Pearson and Spearman correlation statistical analysis. The network shows correlation relationships between nine CAGs of 174 genera. Node size represents the average abundance of each genus. Lines between nodes represent correlations of each other, with the line width representing the correlation magnitude. The red ones represent positive correlations, and the blue ones represent negative correlations. Data are expressed as means ± SEM. # <span class="html-italic">p</span> &lt; 0.05, ## <span class="html-italic">p</span> &lt; 0.01 vs. the MCS group.</p>
Full article ">Figure 4
<p>GPE enriches potential probiotics and reduces opportunistic pathogens. (<b>A</b>) Linear discriminant analysis (LDA) effect size (LEfSe) analysis for differences in abundance between MCD and LGPE groups of biomarkers; (<b>B</b>) random forest classification for top 15 bacterial genera between MCD and LGPE groups; (<b>C</b>) linear discriminant analysis (LDA) effect size (LEfSe) analysis for differences in abundance between MCD group and HGPE groups of biomarkers; (<b>D</b>) random forest classification for top 15 bacterial genera between MCD group and HGPE groups; (<b>E</b>) Venn diagrams to show the common biomarker for LEfSe analysis and random forest classification between MCD and LGPE groups; (<b>F</b>) the relative abundance of <span class="html-italic">Akkermansia</span> between MCD and LGPE groups; (<b>G</b>) Venn diagrams to show the common biomarker for LEfSe analysis and random forest classification between MCD group and HGPE groups; and (<b>H</b>) the relative abundance of <span class="html-italic">Klesiella</span> between the MCD and LGPE groups. Data are expressed as means ± SEM. * <span class="html-italic">p</span> &lt; 0.05 vs. the MCD group.</p>
Full article ">Figure 5
<p>GPE-induced alterations in the liver transcriptome. (<b>A</b>) Principal component analysis (PCA) analysis of the transcript; (<b>B</b>,<b>C</b>) volcano plots show the differentially expressed genes of MCD vs. LGPE and MCD vs. HGPE, as revealed by the transcriptome; (<b>D</b>–<b>G</b>) Gene Ontology (GO) enrichment analysis of the top 20 pathways between the MCD and LGPE groups (<b>D</b>) Gene Ontology (GO) enrichment analysis of the top 20 pathways between the MCD and HGPE groups; (<b>E</b>) Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis of the top 20 pathways between the MCD and LGPE groups; (<b>F</b>) Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis of the top 20 pathways between the MCD and HGPE groups; and (<b>H</b>,<b>I</b>) heatmap of differential gene expression in the NF-kappa B signaling pathway between the MCD and LGPE groups and between the MCD and HGPE groups.</p>
Full article ">Figure 6
<p>GPE suppresses the expression of inflammatory genes in the liver. (<b>A</b>) Pathway analysis of liver tissue of NASH mice based on transcriptomics in combination with ingenuity pathway analysis (IPA) between the MCD and LGPE groups; (<b>B</b>) pathway analysis of liver tissue of NASH mice based on transcriptomics in combination with ingenuity pathway analysis (IPA) between the MCD and HGPE groups; (<b>C</b>) the level of IL-1β in blood serum was tested by enzyme-linked immunosorbent assay (ELISA) (n = 4); (<b>D</b>) the level of LBP in blood serum was tested by ELISA (n = 4); (<b>E</b>–<b>J</b>) RT-PCR of IL-1β, Ccr5, TLR4, CD14, Myd88, and TNF-α mRNA expression in liver. The relative expression of IL-1β, Ccr5, TNF-α, Myd88, CD14, and TLR4 was adjusted with GAPDH as the housekeeping gene (n = 5). Each biological replicate (n = 4–5) had three technical replicate wells for the experiment. Data are expressed as means ± SEM. # <span class="html-italic">p</span> &lt; 0.05, ## <span class="html-italic">p</span> &lt; 0.01 vs. the MCS group; * <span class="html-italic">p</span> &lt; 0.05, ** <span class="html-italic">p</span> &lt; 0.01 vs. the MCD group.</p>
Full article ">
21 pages, 19755 KiB  
Article
Pin1 Exacerbates Non-Alcoholic Fatty Liver Disease by Enhancing Its Activity through Binding to ACC1
by Yiyi Jin, Zhaoshui Shangguan, Jiao Pang, Yuwen Chen, Suijin Lin and Hekun Liu
Int. J. Mol. Sci. 2024, 25(11), 5822; https://doi.org/10.3390/ijms25115822 - 27 May 2024
Viewed by 1012
Abstract
Non-alcoholic fatty liver disease (NAFLD) is a clinicopathological syndrome characterized by diffuse hepatocellular steatosis due to fatty deposits in hepatocytes, excluding alcohol and other known liver injury factors. However, there are no specific drugs for the clinical treatment of NAFLD. Therefore, research on [...] Read more.
Non-alcoholic fatty liver disease (NAFLD) is a clinicopathological syndrome characterized by diffuse hepatocellular steatosis due to fatty deposits in hepatocytes, excluding alcohol and other known liver injury factors. However, there are no specific drugs for the clinical treatment of NAFLD. Therefore, research on the pathogenesis of NAFLD at the cellular and molecular levels is a promising approach to finding therapeutic targets and developing targeted drugs for NAFLD. Pin1 is highly expressed during adipogenesis and contributes to adipose differentiation, but its specific mechanism of action in NAFLD is unclear. In this study, we investigated the role of Pin1 in promoting the development of NAFLD and its potential mechanisms in vitro and in vivo. First, Pin1 was verified in the NAFLD model in vitro using MCD diet-fed mice by Western Blot, RT-qPCR and immunohistochemistry (IHC) assays. In the in vitro study, we used the oleic acid (OA) stimulation-induced lipid accumulation model and examined the lipid accumulation in each group of cells by oil red O staining as well as BODIPY staining. The results showed that knockdown of Pin1 inhibited lipid accumulation in hepatocytes in an in vitro lipid accumulation model and improved lipid indices and liver injury levels. Moreover, in vivo, WT and Pin1-KO mice were fed a methionine-choline deficient (MCD) diet for 4 weeks to induce the NAFLD model. The effects of Pin1 on lipid accumulation, hepatic fibrosis, and oxidative stress were evaluated by biochemical analysis, glucose and insulin tolerance tests, histological analysis, IHC, RT-qPCR and Western blot assays. The results indicate that Pin1 knockdown significantly alleviated hepatic steatosis, fibrosis and inflammation in MCD-induced NAFLD mice, improved glucose tolerance and alleviated insulin resistance in mice. Further studies showed that the AMPK/ACC1 signalling pathway might take part in the process by which Pin1 regulates NAFLD, as evidenced by the inhibition of the AMPK/ACC1 pathway. In addition, immunofluorescence (IF), coimmunoprecipitation (Co-IP) and GST pull-down experiments also showed that Pin1 interacts directly with ACC1 and inhibits ACC1 phosphorylation levels. Our study suggests that Pin1 promotes NAFLD progression by inhibiting the activation of the AMPK/ACC1 signalling pathway, and it is possible that this effect is achieved by Pin1 interacting with ACC1 and inhibiting the phosphorylation of ACC1. Full article
(This article belongs to the Section Molecular Endocrinology and Metabolism)
Show Figures

Figure 1

Figure 1
<p>High expression of Pin1 in an in vitro lipid accumulation model and in vivo MCD-induced NAFLD model. (<b>A</b>) Protein and mRNA expression of Pin1 in the HepG2 cell model. (<b>B</b>) Protein and mRNA expression of Pin1 in the L-02 cell model. (<b>C</b>) Protein and mRNA expression of Pin1 in MCD-induced mouse liver tissue. (<b>D</b>) Immunohistochemical validation of Pin1 expression in MCD-induced mouse liver tissues (200×) (* <span class="html-italic">p</span> &lt; 0.05, ** <span class="html-italic">p</span> &lt; 0.01).</p>
Full article ">Figure 2
<p>Knockdown of Pin1 inhibited lipid accumulation in an in vitro NAFLD model. (<b>A</b>) Western blot to verify the knockdown effect of Pin1 in HepG2 cells. (<b>B</b>) Western blot to verify the knockdown effect of Pin1 in L-02 cells (* <span class="html-italic">p</span> &lt; 0.05). (<b>C</b>) Oil Red O staining in OA-induced HepG2 cells after Pin1 knockdown (100×). (<b>D</b>) Oil red O staining of OA-induced L-02 cells after Pin1 knockdown (100×). (<b>E</b>) Biochemical detection of changes in four levels of lipids and two indicators of liver injury in HepG2 cells after Pin1 knockdown. (<b>F</b>) Biochemical detection of changes in four levels of lipids and two indicators of liver function in L-02 cells after Pin1 knockdown (* indicates <span class="html-italic">p</span> &lt; 0.05 compared with the pLKO.1 group; <sup>#</sup> indicates <span class="html-italic">p</span> &lt; 0.05 compared with the pLKO.1-OA group, <span class="html-italic">p</span> &lt; 0.05).</p>
Full article ">Figure 3
<p>Pin1 overexpression exacerbates lipid accumulation in an in vitro cellular NAFLD model. (<b>A</b>) Western blot to verify the effect of Pin1 overexpression in HepG2 cells. (<b>B</b>) Western blot to verify the effect of Pin1 overexpression in L-02 cells (* <span class="html-italic">p</span> &lt; 0.05). (<b>C</b>) Oil Red O staining of OA-induced HepG2 cells after Pin1 overexpression (100×). (<b>D</b>) Oil red O staining of OA-induced L-02 cells after Pin1 overexpression (100×). (<b>E</b>) Biochemical detection of changes in four levels of lipids and liver function in HepG2 cells after Pin1 overexpression. (<b>F</b>) Biochemical detection of changes in four levels of lipids and liver function in L-02 cells after Pin1 overexpression (* indicates <span class="html-italic">p</span> &lt; 0.05 compared with the pLKO.1 group; <sup>#</sup> indicates <span class="html-italic">p</span> &lt; 0.05 compared with the pLKO.1-OA group).</p>
Full article ">Figure 4
<p>Pin1 inhibits the AMPK/ACC1 pathway in hepatocytes under metabolic stimulation. (<b>A</b>) BODIPY 493/503 staining of OA-induced HepG2 cells (200×). (<b>B</b>) BODIPY staining of OA-induced L-02 cells (200×). (<b>C</b>) Western blot to verify OA-induced changes in major molecules in the AMPK/ACC1 signalling pathway in HepG2 cells. (<b>D</b>) Western blot to verify OA-induced changes in major molecules in the AMPK/ACC1 signalling pathway in L-02 cells. (<b>E</b>) Quantification of p-ACC1/ACC1 and p-AMPK/AMPK protein levels in HepG2 cells normalized to GAPDH. (<b>F</b>) Quantification of p-ACC1/ACC1 and p-AMPK/AMPK protein levels in L-02 cells normalized to GAPDH (* <span class="html-italic">p</span> &lt; 0.05, ** <span class="html-italic">p</span> &lt; 0.01). (<b>G</b>) mRNA levels of genes in the AMPK/ACC1 pathway in the OA-induced cell model (* <span class="html-italic">p</span> &lt; 0.05, ** <span class="html-italic">p</span> &lt; 0.01).</p>
Full article ">Figure 5
<p>Pin1 inhibits the AMPK/ACC1 pathway in hepatocytes under metabolic stimulation. (<b>A</b>–<b>D</b>) Inhibition and activation of AMPK/ACC1 signalling pathway protein expression in HepG2 and L-02 cells after overexpression and knockdown of Pin1, respectively (* indicates <span class="html-italic">p</span> &lt; 0.05, ** indicates <span class="html-italic">p</span> &lt; 0.01; <sup>#</sup> indicates <span class="html-italic">p</span> &lt; 0.05, <sup>##</sup> indicates <span class="html-italic">p</span> &lt; 0.01). (<b>E</b>,<b>F</b>) mRNA levels of genes involved in the AMPK/ACC1 pathway in HepG2 and L-02 cells in OA-induced cell models after overexpression and knockdown of Pin1 (* indicates <span class="html-italic">p</span> &lt; 0.05, ** indicates <span class="html-italic">p</span> &lt; 0.01; <sup>#</sup> indicates <span class="html-italic">p</span> &lt; 0.05, <sup>##</sup> indicates <span class="html-italic">p</span> &lt; 0.01).</p>
Full article ">Figure 6
<p>Pin1 inhibits the AMPK/ACC1 pathway in hepatocytes under metabolic stimulation. (<b>A</b>,<b>B</b>) Overexpression and knockdown of Pin1 in HepG2 and L-02 cells after the addition of the AMPK pathway inhibitor Compound C to validate AMPK/ACC1 signalling pathway protein expression and quantitative statistical plots, respectively. (* indicates <span class="html-italic">p</span> &lt; 0.05, ** indicates <span class="html-italic">p</span> &lt; 0.01; <sup>##</sup> indicates <span class="html-italic">p &lt;</span> 0.01). Note: In the statistical graph of overexpressed cells, A indicates the pBybe-OA group, B indicates the Pin1-OA group, C indicates the pBybe-OA + CC group, and D indicates the Pin1-OA + CC group. In the statistical graph of knockdown cells, A indicates the pLKO.1-OA group, B indicates the shPin1-OA group, C indicates the pLKO.1-OA + CC group, and D indicates the shPin1-OA + CC group.</p>
Full article ">Figure 7
<p>Pin1 deficiency alleviates MCD diet-induced NAFLD progression. (<b>A</b>) Identification of mouse genotypes. Pin1 wild-type (Pin1 +/+) mice had a band at 1000 bp, Pin1 heterozygous (Pin1 +/−) mice had two bands at 1000 bp and 1400 bp, and Pin1 knockout (Pin1 −/−) mice had a band at 1400 bp. (<b>B</b>) Pin1 knockout reduces the MCD diet-induced increase in liver index (* <span class="html-italic">p &lt;</span> 0.05, <sup>#</sup><span class="html-italic">p &lt;</span> 0.05). (<b>C</b>) Knockout of Pin1 improves glucose tolerance in mice after MCD diet feeding (** indicates <span class="html-italic">p &lt;</span> 0.01 compared with the WT-ND group; <sup>#</sup> indicates <span class="html-italic">p &lt;</span> 0.05 compared with the WT-MCD group; <sup>##</sup> indicates <span class="html-italic">p &lt;</span> 0.01 compared with the WT-MCD group). (<b>D</b>) Knockout of Pin1 improves insulin sensitivity in mice after MCD diet feeding (* indicates compared with the WT-ND group; <span class="html-italic">p &lt;</span> 0.05; <sup>#</sup> indicates <span class="html-italic">p &lt;</span> 0.05 compared with the WT-MCD group). (<b>E</b>) Upper layer: HE staining to observe the effect of Pin1 knockout on histopathological changes in mouse liver tissues (400×); middle layer: Oil Red O staining to observe the effect of Pin1 knockout on lipid accumulation in mouse liver tissues (400×); lower layer: Sirius Scarlet staining to observe the effect of Pin1 knockout on the degree of fibrosis in mouse liver tissues (400×).</p>
Full article ">Figure 8
<p>Pin1 deficiency alleviates MCD diet-induced NAFLD progression. (<b>A</b>) Knockout of Pin1 improved MCD-induced lipid levels and liver function impairment in mouse serum (<span class="html-italic">* p &lt;</span> 0.05, <sup>#</sup><span class="html-italic">p &lt;</span> 0.05). (<b>B</b>) Western blot detection and quantification of AMPK/ACC1 pathway protein expression in liver tissues of mice induced by MCD diet. (<b>C</b>) RT-qPCR was used to detect mRNA levels of related genes in AMPK/ACC1 pathway in liver tissue of mice induced by MCD diet (* indicates compared with WT-ND group, <span class="html-italic">p &lt;</span> 0.05; ** indicates <span class="html-italic">p &lt;</span> 0.01 compared with WT-ND group; <sup>#</sup> indicates <span class="html-italic">p &lt;</span> 0.05 compared with WT-MCD group; <sup>##</sup> indicates <span class="html-italic">p &lt;</span> 0.01 compared with WT-MCD group). (<b>D</b>) Immunohistochemical detection of p-ACC1, ACC1 and p-AMPK, AMPK protein expression levels in mouse liver tissues (200×). (<b>E</b>) Immunohistochemical protein expression quantification plots (** indicates <span class="html-italic">p &lt;</span> 0.01 compared with the WT-ND group; <sup>#</sup> indicates <span class="html-italic">p &lt;</span> 0.05 compared with the WT-MCD group; <sup>##</sup> indicates <span class="html-italic">p &lt;</span> 0.01 compared with the WT-MCD group).</p>
Full article ">Figure 9
<p>Pin1 interacts with ACC1 and blocks ACC1 phosphorylation. (<b>A</b>) Immunofluorescence assay results show the colocalization of Pin1 with ACC1 in cells (200×). (<b>B</b>) Co-IP detection of the interaction between Pin1 and ACC1 in cells.</p>
Full article ">Figure 10
<p>Pin1 interacts with ACC1 and blocks ACC1 phosphorylation. (<b>A</b>) GST pull-down assay detected the direct interaction between Pin1 and ACC1. (<b>B</b>) Schematic representation of the Pin1 WW and PPI structural domains. (<b>C</b>) GST-pull down detection of full-length ACC1 and Pin1 and Pin1 WW and PPI structural domain interactions. (<b>D</b>) Schematic representation of the 5 truncated bodies of ACC1. (<b>E</b>) GST pull-down detection of Pin1 and ACC1 full-length and 5 truncated body interactions. (<b>F</b>) The 1705-2380 domain containing T1791A and T2229A mutant ACC1 cannot bind to GST-Pin1 at all.</p>
Full article ">Figure 11
<p>Diagram of the mechanism by which Pin1 regulates NAFLD. On the one hand, Pin1 promotes liver lipid aggregation, inflammation, and fibrosis by inhibiting AMPK mediated phosphorylation of ACC1 and downstream signaling pathway proteins, exacerbating non-alcoholic fatty liver disease; On the other hand, Pin1 regulates the activation of SREBP1, promotes the expression of fat synthase, aggravates liver lipid deposition, and promotes the process of NAFLD. During this process, Pin1 can enhance the activity of ACC1 by directly binding to block its phosphorylation, thereby accelerating the synthesis of fatty acids and promoting liver steatosis by acting on the AMPK/ACC1 pathway.</p>
Full article ">
19 pages, 34859 KiB  
Article
Role of Probiotics in Gut Microbiome and Metabolome in Non-Alcoholic Fatty Liver Disease Mouse Model: A Comparative Study
by Tian Wu, Zheng Zeng and Yanyan Yu
Microorganisms 2024, 12(5), 1020; https://doi.org/10.3390/microorganisms12051020 - 17 May 2024
Cited by 1 | Viewed by 1280
Abstract
Non-alcoholic fatty liver disease (NAFLD) is the most prevalent chronic liver condition worldwide. Numerous studies conducted recently have demonstrated a connection between the dysbiosis of the development of NAFLD and gut microbiota. Rebuilding a healthy gut ecology has been proposed as a strategy [...] Read more.
Non-alcoholic fatty liver disease (NAFLD) is the most prevalent chronic liver condition worldwide. Numerous studies conducted recently have demonstrated a connection between the dysbiosis of the development of NAFLD and gut microbiota. Rebuilding a healthy gut ecology has been proposed as a strategy involving the use of probiotics. The purpose of this work is to investigate and compare the function of probiotics Akkermansia muciniphila (A. muciniphila) and VSL#3 in NAFLD mice. Rodent NAFLD was modeled using a methionine choline-deficient diet (MCD) with/without oral probiotic delivery. Subsequently, qPCR, histological staining, and liver function tests were conducted. Mass spectrometry-based analysis and 16S rDNA gene sequencing were used to investigate the liver metabolome and gut microbiota. We found that while both A. muciniphila and VSL#3 reduced hepatic fat content, A. muciniphila outperformed VSL#3. Furthermore, probiotic treatment restored the β diversity of the gut flora and A. muciniphila decreased the abundance of pathogenic bacteria such as Ileibacterium valens. These probiotics altered the metabolism in MCD mice, especially the glycerophospholipid metabolism. In conclusion, our findings distinguished the role of A. muciniphila and VSL#3 in NAFLD and indicated that oral-gavage probiotics remodel gut microbiota and improve metabolism, raising the possibility of using probiotics in the cure of NAFLD. Full article
(This article belongs to the Special Issue Probiotics, Prebiotics, and Gut Microbes)
Show Figures

Figure 1

Figure 1
<p>Probiotics improve liver functions and reduce lipid accumulation. (<b>A</b>) Alterations in body weight over 8 weeks; (<b>B</b>) Liver weight of mice in four groups; (<b>C</b>) The ratio of liver weight/body weight; (<b>D</b>) Liver function ALT for mice in each group; (<b>E</b>) Serum TG levels, (<b>F</b>) serum LDL levels, and (<b>G</b>) serum HDL levels for mice in each group. The above values are expressed as mean ± SD, <span class="html-italic">n</span> = 10; ns, no significant difference; * <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; **** <span class="html-italic">p</span> &lt; 0.0001. ALT, alanine aminotransferase; TG, triglyceride.</p>
Full article ">Figure 2
<p><span class="html-italic">A. muciniphila</span> and VSL#3 reduce lipid deposition and inflammation in the liver. (<b>A</b>,<b>B</b>) Hepatic histological presentation by Oil Red O and HE staining. White arrows indicate hepatocyte ballooning and hepatocyte steatosis. (<b>C</b>) Quantitative analysis of Oil Red O staining. (<b>D</b>–<b>F</b>) The expression of IL-1b, IL-6, and IL-10 in livers of mice. The above values are expressed as mean ± SD, <span class="html-italic">n</span> = 10; ns, no significant difference; * <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; **** <span class="html-italic">p</span> &lt; 0.0001.</p>
Full article ">Figure 3
<p>Probiotics improve gut barrier function. (<b>A</b>) Morphological analysis of intestine tissues by HE staining; (<b>B</b>). The expression of Cldn-1 in livers of mice; (<b>C</b>) The expression of ZO-2 in livers of mice; (<b>D</b>–<b>F</b>) Expression levels of hepatic cytokines TNF-α, IL-1β, and IL-6, as determined by RT-PCR; (<b>G</b>) The expression of IL-10 in livers of mice. Values are expressed as mean ± SD, <span class="html-italic">n</span> = 10; * <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; **** <span class="html-italic">p</span> &lt; 0.0001.</p>
Full article ">Figure 4
<p>Intestinal microbial diversity and composition among the four groups. (<b>A</b>) Rank abundance curves of the four groups; (<b>B</b>) Venn diagram of the makeup of ASVs in intestinal flora; (<b>C</b>) Average relative abundances of dominant bacterial phyla; (<b>D</b>) Relative abundance of Proteobacteria; (<b>E</b>) α diversity comparison: Shannon index; (<b>F</b>,<b>G</b>) β diversity of clustering analysis: PCoA and NMDS analysis of 4 groups at the amplicon sequence variant (ASV) level. (<b>H</b>) Statistical analysis based on unweighted UniFrac metrics. PCoA, principal coordinates analysis; NMDS, nonmetric multidimensional scaling analysis. * <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 5
<p>LEfSe analysis among the four groups: Cladogram representation of the significantly different taxa features, from phylum (inner circle) to genus (outer circle). LEfSe: linear discriminant analysis (LDA) effect size analysis.</p>
Full article ">Figure 6
<p>Relative abundance functional comparisons by PICRUSt2 among the four groups, based on the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. (<b>A</b>) Level 2 function categories. (<b>B</b>) Level 3 function categories.</p>
Full article ">Figure 7
<p>Metabolic characteristics in the four groups. (<b>A</b>) Venn diagram of the changed metabolites corresponding to groups MCD_A versus MCD, MCD_V versus MCD, and MCD versus SD; (<b>B</b>) OPLS-DA score plots of metabolites. Dynamic distribution of metabolite content differences: (<b>C</b>) MCD_A versus MCD and (<b>D</b>–<b>F</b>) MCD_V versus MCD. Metabolic analysis is based on the KEGG database, and enriched pathways are displayed by bubble plots (MCD_A versus MCD, MCD_V versus MCD).</p>
Full article ">Figure 8
<p>The close correlations between the relative abundance of gut microbiota and altered metabolites. (<b>A</b>) The close correlations between the relative abundance of the main family in the gut microbiota and the altered metabolites between the MCD_V and MCD groups; (<b>B</b>) The close correlations between the relative abundance of the main family in the gut microbiota and the differential metabolites between the MCD_V and MCD groups. * <span class="html-italic">p</span> &lt; 0.05; ** <span class="html-italic">p</span> &lt; 0.01.</p>
Full article ">
24 pages, 5507 KiB  
Article
Comparison and Evaluation of Five Global Land Cover Products on the Tibetan Plateau
by Yongjie Pan, Danyun Wang, Xia Li, Yong Liu and He Huang
Land 2024, 13(4), 522; https://doi.org/10.3390/land13040522 - 14 Apr 2024
Viewed by 1033
Abstract
The Tibetan Plateau (TP) region contains maximal alpine grassland ecology at the mid-latitudes. This region is also recognized as an ecologically fragile and sensitive area under the effects of global warming. Regional climate modeling and ecosystem research depend on accurate land cover (LC) [...] Read more.
The Tibetan Plateau (TP) region contains maximal alpine grassland ecology at the mid-latitudes. This region is also recognized as an ecologically fragile and sensitive area under the effects of global warming. Regional climate modeling and ecosystem research depend on accurate land cover (LC) information. In order to obtain accurate LC information over the TP, the reliability and precision of five moderate/high-resolution LC products (MCD12Q1, C3S-LC, GlobeLand30, GLC_FCS30, and ESA2020 in 2020) were analyzed and evaluated in this study. The different LC products were compared with each other in terms of areal/spatial consistency and assessed with four reference sample datasets (Geo-Wiki, GLCVSS, GOFC-GOLD, and USGS) using the confusion matrix method for accuracy evaluation over the TP. Based on the paired comparison of these five LC datasets, all five LC products show that grass is the major land cover type on the TP, but the range of grass coverage identified by the different products varies noticeably, from 43.35% to 65.49%. The fully consistent spatial regions account for 43.72% of the entire region of the TP, while, in the transition area between grass and bare soil, there is still a large area of medium-to-low consistency. In addition, a comparison of LC datasets using integrated reference datasets shows that the overall accuracies of MCD12Q1, C3S-LC, GlobeLand30, GLC_FCS30, and ESA2020 are 54.29%, 49.32%, 53.03%, 53.73%, and 60.11%, respectively. The producer accuracy of the five products is highest for grass, while glaciers have the most reliable and accurate characteristics among all LC products for users. These findings provide valuable insights for the selection of rational and appropriate LC datasets for studying land-atmosphere interactions and promoting ecological preservation in the TP. Full article
Show Figures

Figure 1

Figure 1
<p>Digital elevation model (DEM) map and the location of Tibetan Plateau.</p>
Full article ">Figure 2
<p>(<b>a</b>) Spatial distribution of the four reference validation samples, (<b>b</b>) and their land cover classification.</p>
Full article ">Figure 3
<p>Spatial patterns for the five LC products on TP: (<b>a</b>) MCD12Q1, (<b>b</b>) C3S-LC, (<b>c</b>) GlobeLand30, (<b>d</b>) GLC_FCS30, and (<b>e</b>) ESA2020.</p>
Full article ">Figure 4
<p>Percentage of separate classes area for five LC products in 2020.</p>
Full article ">Figure 5
<p>Spatial confusion of LC types for each of two different products (ft stands for confounding relationships between two LC datasets).</p>
Full article ">Figure 6
<p>(<b>a</b>) Distribution of consistency of five LC products on the TP, and (<b>b</b>) percentage of spatial consistency at different elevations.</p>
Full article ">Figure 7
<p>Spatial consistency distribution with individual LC types in the TP.</p>
Full article ">Figure 8
<p>Spatial distribution of low and moderate consistency in grass identified by five land cover products (<b>a</b>), and other land cover types besides grass (<b>b</b>).</p>
Full article ">Figure 9
<p>Spatial consistency distribution between different LC datasets and reference validation samples in the TP.</p>
Full article ">
30 pages, 14594 KiB  
Article
Analysis of Factors Driving Subtropical Forest Phenology Differentiation, Considering Temperature and Precipitation Time-Lag Effects: A Case Study of Fujian Province
by Menglu Ma, Hao Zhang, Jushuang Qin, Yutian Liu, Baoguo Wu and Xiaohui Su
Forests 2024, 15(2), 334; https://doi.org/10.3390/f15020334 - 8 Feb 2024
Cited by 1 | Viewed by 1075
Abstract
Subtropical forest phenology differentiation is affected by temperature, precipitation, and topography. Understanding the primary contributing elements and their interactions with forest phenology can help people better comprehend the subtropical forest growth process and its response to climate. Meanwhile, the temporal and spatial variations [...] Read more.
Subtropical forest phenology differentiation is affected by temperature, precipitation, and topography. Understanding the primary contributing elements and their interactions with forest phenology can help people better comprehend the subtropical forest growth process and its response to climate. Meanwhile, the temporal and spatial variations of phenological rhythms are important indicators of climatic impacts on forests. Therefore, this study aimed to analyze both a total area and different forest growth environments within the whole (i.e., coastal site areas (II, IV) and inland site areas (I, III)) as to spatiotemporal patterns associated with subtropical forests in Fujian Province, which is located at the boundary between the middle and south subtropical zones. Considering the asymmetric effects of climate and forest growth, this study chose pre-seasonal and cumulative temperature and precipitation factors and utilized the GeoDetector model to analyze the dominant drivers and interactions within phenology differentiation in Fujian Province. The results show the following: (1) All of the phenological parameters were advanced or shortened over the 19-year observation period; those of shrubland and deciduous broadleaf forests fluctuated greatly, and their stability was poor. (2) The phenological parameters were more distinct at the borders of the site areas. Additionally, the dates associated with the end of the growth season (EOS) and the date-position of peak value (POP) in coastal areas (i.e., II and IV) were later than those in inland areas (i.e., I and III). Among the parameters, the length of the growth season (LOS) was most sensitive to altitude. (3) Precipitation was the main driving factor affecting the spatial heterogeneity of the start of the growth season (SOS) and the EOS. The relatively strong effects of preseason and current-month temperatures on the SOS may be influenced by the temperature threshold required to break bud dormancy, and the relationship between the SOS and temperature was related to the lag time and the length of accumulation. The EOS was susceptible to the hydrothermal conditions of the preseason accumulation, and the variation trend was negatively correlated with temperature and precipitation. Spatial attribution was used to analyze the attribution of phenology differentiation from the perspectives of different regions, thus revealing the relationships between forest phenology and meteorological time-lag effects, the result which can contribute to targeted guidance and support for scientific forest management. Full article
(This article belongs to the Special Issue Application of Remote Sensing in Vegetation Dynamic and Ecology)
Show Figures

Figure 1

Figure 1
<p>The geographical location of the study area: (<b>a</b>) the geographical location of the study area in China; (<b>b</b>) the distribution of forest type, site area, demarcation line, and main rivers in Fujian Province. I, II, III, and IV, respectively, represent the mountainous region of Wuyi Mountain; the coastal low mountain and hill region of Zhejiang and Fujian; the mountainous and hill region in Jiangxi, Fujian, and Guangdong; and the coastal hill and plain region of Fujian and Guangdong. The demarcation line is that between the middle and southern subtropical zones.</p>
Full article ">Figure 2
<p>(<b>a</b>) The proportions of individual regions in the area as to distinguished forest types. (<b>b</b>) The elevation distribution of different vegetation types in Fujian Province. In the figure, ENF, EBF, DBF1, DBF2, and SHL, respectively, represent evergreen needleleaf forests, evergreen broadleaf forests, deciduous broadleaf forests, closed (&gt;40%) deciduous broadleaf forests, and shrubland.</p>
Full article ">Figure 3
<p>The division of the (<b>a</b>) elevation data, (<b>b</b>) slope data, and (<b>c</b>) aspect data in Fujian Province.</p>
Full article ">Figure 4
<p>The research framework of this study. SRTM is the elevation data that could produce slope, as well as the aspect data. ESA CCI represents the land-cover data from the European Space Agency’s Climate Change Initiative product. SOS, EOS, LOS, and POP represent the start of the growth season, the end of the growth season, the length of the growth season, and the date-position of peak value, respectively. Sen + MK represents the methods of Theil–Sen median analysis and the Mann–Kendall test.</p>
Full article ">Figure 5
<p>The spatial distributions of the (<b>a</b>) SOS, (<b>b</b>) EOS, (<b>c</b>) LOS, and (<b>d</b>) POP in the Fujian province. The histograms present the proportion of pixels in the phenological interval. DOY indicates the day of the year.</p>
Full article ">Figure 6
<p>The fluctuations of the (<b>a</b>) SOS, (<b>b</b>) EOS, (<b>c</b>) LOS, and (<b>d</b>) POP among different forest types and site areas.</p>
Full article ">Figure 7
<p>The phenological time trends of (<b>a</b>) evergreen needleleaf forests, (<b>b</b>) evergreen broadleaf forests, (<b>c</b>) deciduous broadleaf forests, (<b>d</b>) deciduous broadleaf forests mixed with other vegetation, and (<b>e</b>) shrubland, in the four site areas, as well as that of (<b>f</b>) all forest pixels.</p>
Full article ">Figure 8
<p>The trend figures of the (<b>a</b>) SOS, (<b>b</b>) EOS, (<b>c</b>) LOS, and (<b>d</b>) POP. In this figure, ESA (ESD), SA (SD), SSA (SSD), NA (ND), and NF, respectively, represent the trends of extremely significantly advanced (or delayed), significantly advanced (or delayed), slightly significantly advanced (or delayed), not significantly advanced (or delayed), and no fluctuation.</p>
Full article ">Figure 9
<p>The factor detection <span class="html-italic">q</span>-values of the SOS and EOS in regions (<b>a</b>) I, (<b>b</b>) II, (<b>c</b>) III, and (<b>d</b>) IV. In the figure, “*” means that the factor has passed the significance test. The lagged impacts of precipitation and temperature are, respectively, presented in the formats “#month_pre” and “#month_tmp,” where “#” represents the number of preseason months. Moreover, the cumulative impacts are, respectively, presented in the formats “acc#_pre” and “acc#_tmp,” where “#” denotes the cumulative length. Finally, DEM is the elevation, Slope and Aspect are terrain factors, and FC_type is the forest type.</p>
Full article ">Figure 9 Cont.
<p>The factor detection <span class="html-italic">q</span>-values of the SOS and EOS in regions (<b>a</b>) I, (<b>b</b>) II, (<b>c</b>) III, and (<b>d</b>) IV. In the figure, “*” means that the factor has passed the significance test. The lagged impacts of precipitation and temperature are, respectively, presented in the formats “#month_pre” and “#month_tmp,” where “#” represents the number of preseason months. Moreover, the cumulative impacts are, respectively, presented in the formats “acc#_pre” and “acc#_tmp,” where “#” denotes the cumulative length. Finally, DEM is the elevation, Slope and Aspect are terrain factors, and FC_type is the forest type.</p>
Full article ">Figure 10
<p>The mean SOS and EOS for (<b>a</b>) forest type, (<b>b</b>) elevation, (<b>c</b>) slope, and (<b>d</b>) aspect factors.</p>
Full article ">Figure 11
<p>The ecological detection results of the (<b>a</b>) SOS and (<b>b</b>) EOS in site areas. Variables <span class="html-italic">X</span><sub>0</sub> to <span class="html-italic">X</span><sub>15</sub> represent the driving factors, Y indicates that the factor had a significant difference in the 95% confidence interval for the SOS or EOS, and N represents “no significant difference.” The factors with no significant difference were labeled according to different regions, and all the other factors except the labeled ones had significant differences.</p>
Full article ">Figure 12
<p>The results of the relationships of lagged and accumulated (<b>a</b>) temperature and (<b>b</b>) precipitation with phenology. TMP indicates temperature and PRE represents precipitation.</p>
Full article ">Figure A1
<p>The relationships between altitude and the average (<b>a</b>) SOS, (<b>b</b>) EOS, (<b>c</b>) LOS, and (<b>d</b>) POP.</p>
Full article ">Figure A2
<p><math display="inline"><semantics> <mrow> <mi>β</mi> </mrow> </semantics></math> values of each phenological parameter in the trend analysis. The horizontal axis is arranged from top to bottom to represent the <math display="inline"><semantics> <mrow> <mi>β</mi> </mrow> </semantics></math> values of SOS (<b>a</b>–<b>c</b>), EOS (<b>d</b>–<b>f</b>), LOS (<b>g</b>–<b>i</b>), and POP (<b>j</b>–<b>l</b>), respectively. The vertical axis, from left to right, represents the spatial distribution of the <math display="inline"><semantics> <mrow> <mi>β</mi> </mrow> </semantics></math> values as follows: less than 0 (<b>a</b>,<b>d</b>,<b>g</b>,<b>j</b>), equal to 0 (<b>b</b>,<b>e</b>,<b>h</b>,<b>k</b>), and greater than 0 (<b>c</b>,<b>f</b>,<b>i</b>,<b>l</b>).</p>
Full article ">Figure A3
<p>The relationships between altitude and the (<b>a</b>) average annual temperature and (<b>b</b>) precipitation. Subfigures (<b>c</b>–<b>f</b>) show the relationship between elevation and precipitation in different elevation gradients.</p>
Full article ">
22 pages, 14409 KiB  
Article
Study on the Susceptibility of Drifting Snow in Ya’an–Qamdo Section of the Railway in Southwest China
by Xue Zhou, Zhen Zhang, Weidong Yang and Qingkuan Liu
Appl. Sci. 2024, 14(2), 475; https://doi.org/10.3390/app14020475 - 5 Jan 2024
Viewed by 760
Abstract
To investigate the susceptibility of drifting snow along the Ya’an–Qamdo section of the railway, which is located in a high-altitude and cold plateau in Southwest China with scarce meteorological information, the Weather Research and Forecasting Model (WRF) is used in this paper to [...] Read more.
To investigate the susceptibility of drifting snow along the Ya’an–Qamdo section of the railway, which is located in a high-altitude and cold plateau in Southwest China with scarce meteorological information, the Weather Research and Forecasting Model (WRF) is used in this paper to simulate the spatio-temporal distribution of meteorological data. According to the varying terrain, the railway section from Ya’an to Qamdo is divided into two regions along 100.8° E for double-layer nested simulation. The original land use data of the WRF model are used in region 1. Due to the increased number of mountains in region 2, the original data are replaced by the MCD12Q1v006 land use data, and the vertical direction layers are densified near the ground to increase simulation accuracy. The simulated results are compared with the observation data. It is found that after densification, the results have been significantly improved. The results obtained by the WRF model can accurately simulate the change trends of temperature, rainfall, and wind speed, and the correlation coefficients are relatively high, which verifies the accuracy of WRF for simulating complex terrain regions. The simulation results further indicate that approximately 300 km of the Ya’an–Qamdo railway may experience drifting snow. Among them, no drifting snow events occur in Ya’an County, and the areas with higher probability are located at the border between Luding County and Tianquan County, followed by Kangding area. The remaining areas have a probability of less than 10%. The WRF model demonstrates its capability in the drifting snow protection of railways with limited meteorological data. Full article
Show Figures

Figure 1

Figure 1
<p>The railway of Ya’an–Qamdo section.</p>
Full article ">Figure 2
<p>Elevation of Ya’an–Qamdo section.</p>
Full article ">Figure 3
<p>Simulation range for both regions. (<b>a</b>) The simulation range of region 1; (<b>b</b>) the simulation range of region 2.</p>
Full article ">Figure 4
<p>Vertical stratification of region 1 (<b>left</b>) and region 2 (<b>right</b>).</p>
Full article ">Figure 5
<p>Wind speed variation trends at different ETA levels.</p>
Full article ">Figure 6
<p>Comparison of simulation results in Kangding. (<b>a</b>) Comparison of precipitation result; (<b>b</b>) comparison of wind speed result.</p>
Full article ">Figure 7
<p>The results of the WRF fitting process and a comparison between the simulated and observed values. (<b>a</b>) The fitted line of temperature in Kangding; (<b>b</b>) the variation in temperature in Kangding; (<b>c</b>) the fitted line of temperature in Litang; (<b>d</b>) the variation in temperature in Litang; (<b>e</b>) the fitted line of wind speed in Kangding; (<b>f</b>) the variation in wind speed in Kangding; (<b>g</b>) the fitted line of wind speed in Litang; (<b>h</b>) the variation in wind speed in Litang.</p>
Full article ">Figure 7 Cont.
<p>The results of the WRF fitting process and a comparison between the simulated and observed values. (<b>a</b>) The fitted line of temperature in Kangding; (<b>b</b>) the variation in temperature in Kangding; (<b>c</b>) the fitted line of temperature in Litang; (<b>d</b>) the variation in temperature in Litang; (<b>e</b>) the fitted line of wind speed in Kangding; (<b>f</b>) the variation in wind speed in Kangding; (<b>g</b>) the fitted line of wind speed in Litang; (<b>h</b>) the variation in wind speed in Litang.</p>
Full article ">Figure 8
<p>Comparison of daily rainfall in Litang.</p>
Full article ">Figure 9
<p>The distribution of average temperature in the simulation period. (<b>a</b>) Temperature of region 1 in December; (<b>b</b>) temperature of region 2 in December; (<b>c</b>) temperature of region 1 in January; (<b>d</b>) temperature of region 2 in January; (<b>e</b>) temperature of region 1 in February; (<b>f</b>) temperature of region 2 in February.</p>
Full article ">Figure 10
<p>The distribution of cumulative precipitation in the simulation period. (<b>a</b>) Precipitation of region 1 in December; (<b>b</b>) precipitation of region 2 in December; (<b>c</b>) precipitation of region 1 in January; (<b>d</b>) precipitation of region 2 in January; (<b>e</b>) precipitation of region 1 in February; (<b>f</b>) precipitation of region 2 in February.</p>
Full article ">Figure 11
<p>The distribution of average wind field in the simulation period. (<b>a</b>) Wind field of region 1 in December; (<b>b</b>) wind field of region 2 in December; (<b>c</b>) wind field of region 1 in January; (<b>d</b>) wind field of region 2 in January; (<b>e</b>) wind field of region 1 in February; (<b>f</b>) wind field of region 2 in February.</p>
Full article ">Figure 12
<p>Wind rose diagram and the wind speed probability density graph. (<b>a</b>) the wind rose diagram of P1; (<b>b</b>) the wind speed probability density graph of P1; (<b>c</b>) the wind rose diagram of P2; (<b>d</b>) the wind speed probability density graph of P2; (<b>e</b>) the wind rose diagram of P3; (<b>f</b>) the wind speed probability density graph of P3; (<b>g</b>) the wind rose diagram of P4; (<b>h</b>) the wind speed probability density graph of P4; (<b>i</b>) the wind rose diagram of P5; (<b>j</b>) the wind speed probability density graph of P5.</p>
Full article ">Figure 12 Cont.
<p>Wind rose diagram and the wind speed probability density graph. (<b>a</b>) the wind rose diagram of P1; (<b>b</b>) the wind speed probability density graph of P1; (<b>c</b>) the wind rose diagram of P2; (<b>d</b>) the wind speed probability density graph of P2; (<b>e</b>) the wind rose diagram of P3; (<b>f</b>) the wind speed probability density graph of P3; (<b>g</b>) the wind rose diagram of P4; (<b>h</b>) the wind speed probability density graph of P4; (<b>i</b>) the wind rose diagram of P5; (<b>j</b>) the wind speed probability density graph of P5.</p>
Full article ">Figure 13
<p>Probability of drifting snow along the railway. (<b>a</b>) The susceptibility of drifting snow in region 1; (<b>b</b>) the susceptibility of drifting snow in region 2.</p>
Full article ">
27 pages, 9205 KiB  
Article
Seasonal Effect of the Vegetation Clumping Index on Gross Primary Productivity Estimated by a Two-Leaf Light Use Efficiency Model
by Zhilong Li, Ziti Jiao, Chenxia Wang, Siyang Yin, Jing Guo, Yidong Tong, Ge Gao, Zheyou Tan and Sizhe Chen
Remote Sens. 2023, 15(23), 5537; https://doi.org/10.3390/rs15235537 - 28 Nov 2023
Cited by 2 | Viewed by 1100
Abstract
Recently, light use efficiency (LUE) models driven by remote sensing data have been widely employed to estimate the gross primary productivity (GPP) of different terrestrial ecosystems at global or regional scales. Furthermore, the two-leaf light use efficiency (TL-LUE) model has been reported to [...] Read more.
Recently, light use efficiency (LUE) models driven by remote sensing data have been widely employed to estimate the gross primary productivity (GPP) of different terrestrial ecosystems at global or regional scales. Furthermore, the two-leaf light use efficiency (TL-LUE) model has been reported to improve the accuracy of GPP estimation, relative to the big-leaf MOD17 model, by separating the entire canopy into sunlit and shaded leaves through the use of constant clumping index estimation (Ω). However, ignoring obvious seasonal changes in the vegetation clumping index (CI) most likely results in GPP estimation errors since the CI tends to present seasonal changes, especially with respect to the obvious presence or absence of leaves within the canopy of deciduous vegetation. Here, we propose a TL-CLUE model that considers the seasonal difference in the CI based on the TL-LUE model to characterize general changes in canopy seasonality. This method composites monthly CI values into two or three Ω values to capture the general seasonal changes in CI while attempting to reduce the potential uncertainty caused during CI inversion. In theory, CI seasonality plays an essential role in the distribution of photosynthetically active radiation absorbed by the canopy (APAR). Specifically, the seasonal difference in CI values mainly considers the state of leaf growth, which is determined by the MODIS land surface phenology (LSP) product (MCD12Q2). Therefore, the one-year cycle (OYC) of leaf life is divided into two (leaf-off and leaf-on) or three seasons (leaf-off, leaf-scattering, and leaf-gathering) according to this MODIS LSP product, and the mean CI of each corresponding season for each vegetation class is computed to smoothen the uncertainties within each seasonal section. With these two or three seasonal Ω values as inputs, the TL-CLUE model by which the seasonal differences in CI are incorporated into the TL-LUE model is run and evaluated based on observations from 84 eddy covariance (EC) tower sites across North America. The results of the analysis reveal that the TL-LUE model widely overestimates GPP for most vegetation types during the leaf-on season, particularly during the growth peak. Although the TL-LUE model shows that the temporal characteristics of GPP agree with the EC observations in terms of general trends, the TL-CLUE model further improves the accuracy of GPP estimation by considering the seasonal changes in the CI. The result of GPP estimation from the TL-CLUE model shows a lower error (RMSE = 2.46 g C m−2 d−1) than the TL-LUE model (RMSE = 2.75 g C m−2 d−1) and somewhat decreases the eight-day GPP overestimation in the TL-LUE model with a constant Ω by approximately 9.76 and 8.970% when adapting three and two Ωs from different seasons, respectively. The study demonstrates that the uncertainty of seasonal disturbance in the CI, quantified by a standard deviation of approximately 0.071 relative to the mean CI of 0.746, is diminished through simple averaging. The seasonal difference in CI should be considered in GPP estimation of terrestrial ecosystems, particularly for vegetation with obvious canopy changes, where leaves go through the complete physiological processes of germination, stretching, maturity, and falling within a year. This study demonstrates the potential of the MODIS CI application in developing ecosystem and hydrological models. Full article
(This article belongs to the Section Ecological Remote Sensing)
Show Figures

Figure 1

Figure 1
<p>Site distribution of GPP measurement eddy covariance towers across North America in different vegetation types obtained from the MODIS Year Land Cover Version 061 (MCD12Q1.061) product. There are 35 sites from the AmeriFlux dataset indicated by orange rectangles, and 49− sites from the FLUXNET2015 dataset represented by red circles.</p>
Full article ">Figure 2
<p>Flow chart of the improvement of GPP estimation by considering seasonal differences in the CI (TL-CLUE). The total sensitivity of GPP to the CI is explored by the EFAST method. The OYC for leaf growth is divided into LFS and LSS, consisting of LOS and LGS, labeled according to MODIS LC (land cover) data (MCD12Q1.061) and MODIS LSP (land surface phenology) data (MCD12Q2 V061). Then, the corresponding average CI (Ω) of each season for each vegetation type is estimated based on MODIS CI data and published literature.</p>
Full article ">Figure 3
<p>Comparison of the eight-day change in GPP (red point line) observed at the EC site with GPP (deep green point line) estimated by the TL-LUE model for DBF (<b>a</b>): deciduous broadleaf forest, ENF (<b>b</b>): evergreen needleleaf forest, MF (<b>c</b>): mixed forest, OSH (<b>d</b>): open shrublands, CRO (<b>e</b>): croplands, GRA (<b>f</b>): grasslands, WSA (<b>g</b>): woody savannas, SAV (<b>h</b>): savannas during 2002~2020, respectively. There is an evident overestimation of GPP in the TL-LUE model, particularly in the growing peak, by approximately 5.70 for MF, 3.78 for OSH, and 2.37 g C m<sup>−2</sup> d<sup>−1</sup> for CRO.</p>
Full article ">Figure 4
<p>Eight-day changes in the mean CI (black point line) for different vegetation types during 2002~2020. The gray solid vertical line represents the uncertainty of the CI, as measured by the standard deviation (SD), the black dashed line is the default of Ω currently used in the TL-LUE model, and the black dots are eight-day CI mean in each vegetation type by averaging CI values from 2002 to 2020 in all sites extracted from the MODIS CI. This figure shows that although the eight-day average CIs capture general seasonal trends for some vegetation classes, some large uncertainties exist due to the changes in the background and ambient reflectivity, abnormal land surface shadows, insufficient BRDF information, and even some systematic vibrations that cannot be well understood in some classes, e.g., approximately from 0.016 to 0.161 for MF, from 0.032 to 0.132 for WSA, and from 0.048 to 0.093 for GRA.</p>
Full article ">Figure 5
<p>The total sensitivity of the eight-day GPP to the CI in different vegetation types. This total sensitivity is more than 0.72 for various vegetation types in the figure, which indicates that changes in the CI would affect the magnitude of GPP to some extent.</p>
Full article ">Figure 6
<p>Comparison of the consistency of the observed GPP with GPP estimated by the TL-CLUE (<b>a1</b>–<b>a5</b>) and TL-LUE (<b>b1</b>–<b>b5</b>) models in the LOS (<b>a1</b>,<b>b1</b>), LGS (<b>a2</b>,<b>b2</b>), LSS (<b>a3</b>,<b>b3</b>), LFS (<b>a4</b>,<b>b4</b>), and OYC (<b>a5</b>,<b>b5</b>), respectively. The 1:1 theory line is shown by the black dashed line, and the regression line is displayed by the blue dashed line.</p>
Full article ">Figure 7
<p>R<sub>bias</sub> difference of eight-day GPP estimation between the TL-CLUE model and the TL-LUE model. The TL-CLUE model adopts 2 (yellow par) and 3 (light blue par) Ωs from different seasons.</p>
Full article ">Figure 8
<p>The spatial difference in GPP simulation between the TL-CLUE and TL-LUE models in the leaf-on season (<b>a</b>), leaf-gathering season (<b>b</b>), leaf-scattering season (<b>c</b>), leaf-off season (<b>d</b>), and one-year cycle (<b>e</b>) in different vegetation types (<b>f</b>) obtained from the MODIS Year Land Cover (MCD12Q1.061) across North America during 2002~2020. The magnitude and spatial characteristics of the difference in GPP simulation between the TL-CLUE and TL-LUE models in the leaf-on season are similar to those in the OYC for most areas of North America. The difference is less than 0, approximately from −100 to 0 g C m<sup>−2</sup> yr<sup>−1</sup>, in the leaf-on season and OYC across North America.</p>
Full article ">Figure 9
<p>The spatial characteristics of GPP simulated by the TL-CLUE model in the leaf-on season (<b>a</b>), leaf-gathering season (<b>b</b>), leaf-scattering season (<b>c</b>), and leaf-off season (<b>d</b>) across North America during 2002~2020. GPP in the leaf-on season shows obvious spatial heterogeneity; the lowest value, ranging from approximately 0 to 400 g C m<sup>−2</sup> yr<sup>−1</sup>, occurs in the high elevation areas of the Rocky Mountains and high latitude areas north of 52°N, and the highest value, ranging from approximately 1800 to 2200 g C m<sup>−2</sup> yr<sup>−1</sup>, occurs on the east coast south of 48°N. GPP is near 0 in the leaf-off season for most areas of North America.</p>
Full article ">Figure 10
<p>The annual spatial pattern of GPP simulated by the TL-CLUE model across North America during 2002~2020. The magnitude and spatial patterns of the annual GPP are close to those of the leaf-on season shown in <a href="#remotesensing-15-05537-f009" class="html-fig">Figure 9</a>a, which could be attributed to the fact that GPP is mainly generated by vegetation in the leaf-on season.</p>
Full article ">
20 pages, 10522 KiB  
Article
Phenology of Vegetation in Arid Northwest China Based on Sun-Induced Chlorophyll Fluorescence
by Zhizhong Chen, Mei Zan, Jingjing Kong, Shunfa Yang and Cong Xue
Forests 2023, 14(12), 2310; https://doi.org/10.3390/f14122310 - 24 Nov 2023
Viewed by 1074
Abstract
The accurate monitoring of vegetation phenology is critical for carbon sequestration and sink enhancement. Vegetation phenology in arid zones is more sensitive to climate responses; therefore, it is important to conduct research on phenology in arid zones in response to global climate change. [...] Read more.
The accurate monitoring of vegetation phenology is critical for carbon sequestration and sink enhancement. Vegetation phenology in arid zones is more sensitive to climate responses; therefore, it is important to conduct research on phenology in arid zones in response to global climate change. This study compared the applicability of the enhanced vegetation index (EVI), which is superior in arid zones, and global solar-induced chlorophyll fluorescence (GOSIF), which has a high spatial resolution, in extracting vegetation phenology in arid zones, and explored the mechanism of the differences in the effects of environmental factors on the phenology of different vegetation types. Therefore, this study employed a global solar-induced chlorophyll fluorescence (GOSIF) dataset to determine the start and end of the vegetation growth season (SOSSIF and EOSSIF, respectively) in the arid zone of Northwest China from 2001 to 2019. The results were compared with those from the EVI-based MODIS climate product MCD12Q2 (SOSEVI and EOSEVI). Variations in the sensitivity of these climatic datasets concerning temperature, precipitation, and standardised precipitation evapotranspiration index (SPEI) were assessed through partial correlation analysis. Results: Compared to the MCD12Q2 climatic products, SOSSIF and EOSSIF closely matched the observed climate data in the study area. Spring onset was delayed at higher altitudes and latitudes, and the end of the growing season occurred earlier in these areas. Both SOSSIF and EOSSIF significantly advanced from 2001 to 2019 (trend degrees −0.22 and −0.48, respectively). Spring vegetation phenology was chiefly influenced by precipitation while autumn vegetation phenology was driven by both precipitation and SPEI. GOSIF-based climate data provides a more accurate representation of vegetation phenology compared to traditional vegetation indices. The findings of this study contribute to a deeper understanding of the potential ability of EVI and SIF to reveal the influence of vegetation phenology on the carbon cycle. Full article
(This article belongs to the Special Issue Woody Plant Phenology in a Changing Climate)
Show Figures

Figure 1

Figure 1
<p>Overview map of the study area.</p>
Full article ">Figure 2
<p>Spatial distribution map of vegetation phenology for SIF and MODIS.</p>
Full article ">Figure 3
<p>Time series curves of month-by-month SIF, EVI, and GPP averages from 2001 to 2019.</p>
Full article ">Figure 4
<p>Time series of SEN trend of vegetation phenology in the arid zone.</p>
Full article ">Figure 5
<p>Spatial characteristics of Mann–Kendall significance test for SEN trend of vegetation phenology in the arid zone.</p>
Full article ">Figure 6
<p>(<b>a</b>) Spatial representation of the bias correlation between vegetation phenology and temperature. (<b>b</b>) Spatial depiction of the bias correlation between vegetation phenology and precipitation. (<b>c</b>) Frequency plot of the bias correlation coefficient among vegetation phenology, temperature, and precipitation.</p>
Full article ">Figure 6 Cont.
<p>(<b>a</b>) Spatial representation of the bias correlation between vegetation phenology and temperature. (<b>b</b>) Spatial depiction of the bias correlation between vegetation phenology and precipitation. (<b>c</b>) Frequency plot of the bias correlation coefficient among vegetation phenology, temperature, and precipitation.</p>
Full article ">Figure 7
<p>Frequency plot of significantly biased correlation of vegetation phenology with temperature and precipitation.</p>
Full article ">Figure 8
<p>Bias correlation analysis between SIF- and MODIS-based phenology with SPEI changes in the arid zone. (<b>a</b>) Spatial pattern of bias correlation coefficients between vegetation phenology and SPEI. (<b>b</b>) Frequency plot displaying bias correlation coefficients between vegetation phenology and SPEI. (<b>c</b>) Frequency distribution of image elements with a significant bias correlation between vegetation phenology and SPEI.</p>
Full article ">Figure 9
<p>Partial correlation coefficients of phenological response to climate for different vegetation types. The orange, green and mauve colours in the figure indicate the correlation of phenology with temperature, precipitation and SPEI, respectively, for different vegetation types.</p>
Full article ">Figure 10
<p>Validation of physical weather results based on site data.</p>
Full article ">
23 pages, 8208 KiB  
Article
Spatiotemporal Vegetation Variability and Linkage with Snow-Hydroclimatic Factors in Western Himalaya Using Remote Sensing and Google Earth Engine (GEE)
by Dhiraj Kumar Singh, Kamal Kant Singh, George P. Petropoulos, Priestly Shan Boaz, Prince Jain, Sartajvir Singh, Dileep Kumar Gupta and Vishakha Sood
Remote Sens. 2023, 15(21), 5239; https://doi.org/10.3390/rs15215239 - 4 Nov 2023
Cited by 1 | Viewed by 1700
Abstract
The mountain systems of the Himalayan regions are changing rapidly due to climatic change at a local and global scale. The Indian Western Himalaya ecosystem (between the tree line and the snow line) is an underappreciated component. Yet, knowledge of vegetation distribution, rates [...] Read more.
The mountain systems of the Himalayan regions are changing rapidly due to climatic change at a local and global scale. The Indian Western Himalaya ecosystem (between the tree line and the snow line) is an underappreciated component. Yet, knowledge of vegetation distribution, rates of change, and vegetation interactions with snow-hydroclimatic elements is lacking. The purpose of this study is to investigate the linkage between the spatiotemporal variability of vegetation (i.e., greenness and forest) and related snow-hydroclimatic parameters (i.e., snow cover, land surface temperature, Tropical Rainfall Measuring Mission (TRMM), and Evapotranspiration (ET)) in Himachal Pradesh (HP) Basins (i.e., Beas, Chandra, and Bhaga). Spatiotemporal variability in forest and grassland has been estimated from MODIS land cover product (MCD12Q1) using Google Earth Engine (GEE) for the last 19 years (2001–2019). A significant inter- and intra-annual variation in the forest, grassland, and snow-hydroclimatic factors have been observed during the data period in HP basins (i.e., Beas, Chandra, and Bhaga basin). The analysis demonstrates a significant decrease in the forest cover (214 ha/yr.) at the Beas basin; however, a significant increase in grassland cover is noted at the Beas basin (459 ha/yr.), Chandra (176.9 ha/yr.), and Bhaga basin (9.1 ha/yr.) during the data period. Spatiotemporal forest cover loss and gain in the Beas basin have been observed at ~7504 ha (6.6%) and 1819 ha (1.6%), respectively, from 2001 to 2019. However, loss and gain in grassland cover were observed in 3297 ha (2.9%) and 10,688 ha (9.4%) in the Beas basin, 1453 ha (0.59%) and 3941 ha (1.6%) in the Chandra basin, and 1185 ha (0.92%) and 773 ha (0.60%) in the Bhaga basin, respectively. Further, a strong negative correlation (r = −0.65) has been observed between forest cover and evapotranspiration (ET). However, a strong positive correlation (r = 0.99) has been recorded between grassland cover and ET as compared to other factors. The main outcome of this study in terms of spatiotemporal loss and gain in forest and grassland shows that in the Bhaga basin, very little gain and loss have been observed as compared to the Chandra and Beas basins. The present study findings may provide important aid in the protection and advancement of the knowledge gap of the natural environment and the management of water resources in the HP Basin and other high-mountain regions of the Himalayas. For the first time, this study provides a thorough examination of the spatiotemporal variability of forest and grassland and their interactions with snow-hydroclimatic factors using GEE for Western Himalaya. Full article
Show Figures

Figure 1

Figure 1
<p>(<b>a</b>) Geographical location of the study area, (<b>b</b>) spatial elevation map of three HP basins (i.e., Beas, Chandra, and Bhaga), and (<b>c</b>) hypsometric curve of Beas, Chandra, and Bhaga basins.</p>
Full article ">Figure 2
<p>Methodology adopted in this study.</p>
Full article ">Figure 3
<p>Annual spatiotemporal variation of forest and grassland in the Beas basin from 2001 to 2019.</p>
Full article ">Figure 4
<p>Annual spatiotemporal variation of grassland in the Chandra basin from 2001 to 2019.</p>
Full article ">Figure 5
<p>Annual spatiotemporal variation of grassland in the Bhaga basin from 2001 to 2019.</p>
Full article ">Figure 6
<p>Inter- and intra-annual variation of (<b>a</b>) PPT, (<b>b</b>) ET, (<b>c</b>) LST, and (<b>d</b>) SCA in the Beas basin from 2001 to 2019.</p>
Full article ">Figure 7
<p>Intra- and inter-annual variation of (<b>a</b>) PPT, (<b>b</b>) ET, (<b>c</b>) LST, and (<b>d</b>) SCA in the Chandra basin from 2001 to 2019.</p>
Full article ">Figure 8
<p>Intra- and inter-annual variation of (<b>a</b>) PPT, (<b>b</b>) ET, (<b>c</b>) LST, and (<b>d</b>) SCA in the Bhaga basin during 2001–2019.</p>
Full article ">Figure 9
<p>Temporal variation of annual mean forest cover, grassland, ET, PPT, LST, and SCA and statistical significance (α) of the trend in the Beas basin during 2001–2019.</p>
Full article ">Figure 10
<p>Temporal variation of annual mean grassland, ET, PPT, LST, and SCA and statistical significance (α) of the trend in the Chandra basin during 2001–2019.</p>
Full article ">Figure 11
<p>Temporal variation of annual mean grassland, ET, PPT, LST, and SCA, and statistical significance (α) of the trend in the Bhaga basin during 2001–2019.</p>
Full article ">Figure 12
<p>Spatiotemporal gain and loss of forest and grassland in the Beas basin (<b>a</b>) during 2001–2019, (<b>b</b>) during 2001–2010, and (<b>c</b>) during 2010–2019. Grassland gain and loss in the Chandra basin (<b>d</b>) during 2001–2019, (<b>e</b>) during 2001–2010, and (<b>f</b>) during 2010–2019. Grassland gain and loss in the Bhaga basin (<b>g</b>) during 2001–2019, (<b>h</b>) during 2001–2010, and (<b>i</b>) during 2010–2019.</p>
Full article ">Figure 13
<p>Heat map of annual Pearson correlation coefficient, <span class="html-italic">p</span>-value, and * represents statistical significance at the 0.05 level for forest and grassland with ET, PPT, LST, and SCA in the study area.</p>
Full article ">
21 pages, 12947 KiB  
Article
Spatio-Temporal Variability Analysis of Vegetation Dynamics in China from 2000 to 2022 Based on Leaf Area Index: A Multi-Temporal Image Classification Perspective
by Tianchi Xu, Kai Yan, Yuanpeng He, Si Gao, Kai Yang, Jingrui Wang, Jinxiu Liu and Zhao Liu
Remote Sens. 2023, 15(12), 2975; https://doi.org/10.3390/rs15122975 - 7 Jun 2023
Cited by 1 | Viewed by 1980
Abstract
Leaf Area Index (LAI) is one of the most important biophysical parameters of vegetation, and its dynamic changes can be used as a reflective indicator and differentiation basis of vegetation function. In this study, a VCA–MLC (Vertex Component Analysis–Maximum Likelihood Classification) algorithm is [...] Read more.
Leaf Area Index (LAI) is one of the most important biophysical parameters of vegetation, and its dynamic changes can be used as a reflective indicator and differentiation basis of vegetation function. In this study, a VCA–MLC (Vertex Component Analysis–Maximum Likelihood Classification) algorithm is proposed from the perspective of multi-temporal satellite LAI image classification to monitor and quantify the spatial and temporal variability of vegetation dynamics in China since 2000. The algorithm extracts the vegetation endmembers from 46 multi-temporal images of MODIS LAI in 2011 without the aid of other a priori knowledge and uses the maximum likelihood classification method to select the categories that satisfy the requirements of the number of missing periods, absolute distance, and relative distance for the rest pixels to be classified, ultimately dividing the vegetation area of China into 10 vegetation zones called China Vegetation Functional Zones (CVFZ). CVFZ outperforms MCD12Q1 and CLCD land cover datasets in the overall differentiation of vegetation functions and can be used synergistically with other land cover datasets. In this study, CVFZ is used to cut the constant vegetation-type pixels of MCD12Q1 during 2001–2022. The results of the LAI mean time series decomposition of each subregion using the STL (Seasonal-Trend Decomposition based on Loess) method show that the rate of vegetation greening ranges from 9.02 × 10−4 m2m−2yr−1 in shrubland subregions to 2.34 × 10−2 m2m−2yr−1 in savanna subregions. In relative terms, the average greening speed of forests is moderate, and savannas tend to have the fastest average greening speed. The greening speed of grasslands and croplands in different zones varies widely. In contrast, the average greening speed of shrublands is the slowest. In addition, CVFZ detected grasslands with one or two phenological cycles, broadleaf croplands with one or two phenological cycles, and shrublands with no apparent or one phenological cycle. Full article
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Temporal average of MODIS LAI in China, 25 June 2000 to 25 June 2022.</p>
Full article ">Figure 2
<p>VCA–MLC algorithm flowchart. The process of extracting vegetation endmembers was conducted 1000 times for the sample of MOD-2011, which was also used as the sample for subsequent MLC.</p>
Full article ">Figure 3
<p>China Vegetation Functional Zones (CVFZ) and location of the selected vegetation endmembers. The curve in the lower left corner corresponds to the color of the classification chart.</p>
Full article ">Figure 4
<p>Validation results of CVFZ and several land cover products. The proportion of each vegetation type was calculated to provide more information. (<b>a</b>) MCD-2011, (<b>b</b>) CLCD-2011, (<b>c</b>) MCD-unchanged, and (<b>d</b>) MCD-CVFZ (the legend of <a href="#remotesensing-15-02975-f004" class="html-fig">Figure 4</a>d only shows areas with more than 10,000 pixels).</p>
Full article ">Figure 5
<p>The performance of CVFZ in distinguishing the LAI mean and LAI standard deviation of MOD-2011. The LAI mean for each vegetation type is represented by black dots, and the legend’s color scheme is consistent with <a href="#remotesensing-15-02975-f003" class="html-fig">Figure 3</a> (some zones are not distributed in all vegetation types).</p>
Full article ">Figure 6
<p>The performance of CVFZ in distinguishing LAI time series, 25 June 2000 to 25 June 2022. The legend’s color scheme is consistent with <a href="#remotesensing-15-02975-f003" class="html-fig">Figure 3</a>.</p>
Full article ">Figure 7
<p>LAI time series and linear regression results of annual fluctuations, 25 June 2000 to 25 June 2022. The legend’s color scheme is consistent with <a href="#remotesensing-15-02975-f003" class="html-fig">Figure 3</a>.</p>
Full article ">Figure 8
<p>LAI time series of intra-annual fluctuations, 25 June 2000 to 25 June 2022. The legend’s color scheme is consistent with <a href="#remotesensing-15-02975-f003" class="html-fig">Figure 3</a>.</p>
Full article ">Figure 9
<p>The LAI mean and <span class="html-italic">D</span><sub>2</sub> (indicated by the error bar) in MCD-2011 and MCD-CVFZ in the same vegetation. The red curve represents the proportion of pixels remaining after MCD is cut.</p>
Full article ">Figure 10
<p>Spatial variability in vegetation dynamics revealed by CVFZ: a case study of the Hengduan Mountains, the Chengdu Plain, and the Sichuan Basin.</p>
Full article ">
22 pages, 9765 KiB  
Article
Comparison of Various Annual Land Cover Datasets in the Yellow River Basin
by Bo Liu, Zemin Zhang, Libo Pan, Yibo Sun, Shengnan Ji, Xiao Guan, Junsheng Li and Mingzhu Xu
Remote Sens. 2023, 15(10), 2539; https://doi.org/10.3390/rs15102539 - 12 May 2023
Cited by 4 | Viewed by 1579
Abstract
Accurate land cover (LC) datasets are the basis for global environmental and climate change studies. Recently, numerous open-source annual LC datasets have been created due to advances in remote sensing technology. However, the agreements and sources of error that affect the accuracy of [...] Read more.
Accurate land cover (LC) datasets are the basis for global environmental and climate change studies. Recently, numerous open-source annual LC datasets have been created due to advances in remote sensing technology. However, the agreements and sources of error that affect the accuracy of current annual LC datasets are not well understood, which limits the widespread use of these datasets. We compared four annual LC datasets, namely the CLCD, MCD12Q1, CCI-LC, and GLASS-LC, in the Yellow River Basin (YRB) to identify their spatial and temporal agreement for nine LC classes and to analyze their sources of error. The Mann–Kendall test, Sen’s slope analysis, Taylor diagram, and error decomposition analysis were used in this study. Our results showed that the main LC classes in the four datasets were grassland and cropland (total area percentage > 80%), but their trends in area of change were different. For the main LC classes, the temporal agreement was the highest between the CCI-LC and CLCD (0.85), followed by the MCD12Q1 (0.21), while the lowest was between the GLASS-LC and CLCD (−0.11). The spatial distribution of area for the main LC classes was largely similar between the four datasets, but the spatial agreement in their trends in area of change varied considerably. The spatial variation in the trends in area of change for the cropland, forest, grassland, barren, and impervious LC classes were mainly located in the upstream area region (UA) and the midstream area region (MA) of the YRB, where the percentage of systematic error was high (>68.55%). This indicated that the spatial variation between the four datasets was mainly caused by systematic errors. Between the four datasets, the total error increased along with landscape heterogeneity. These results not only improve our understanding of the spatial and temporal agreement and sources of error between the various current annual LC datasets, but also provide support for land policy making in the YRB. Full article
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Location and geomorphology of the Yellow River Basin (YRB).</p>
Full article ">Figure 2
<p>Flowchart of the comparison between the four annual LC datasets.</p>
Full article ">Figure 3
<p>Comparison of the annual area of change for the nine LC classes between the four datasets in the YRB from 2001 to 2015. (<b>a</b>) Cropland, (<b>b</b>) forest, (<b>c</b>) shrubland, (<b>d</b>) grassland, (<b>e</b>) water, (<b>f</b>) snow/ice, (<b>g</b>) barren, (<b>h</b>) impervious, (<b>i</b>) wetland. The inset graph is the Sen’s slope estimate (β) of the area of change in the four datasets and * indicates that the result passed the Mann–Kendall test with α = 0.05.</p>
Full article ">Figure 4
<p>Taylor diagrams of the area of change for the nine LC classes between the four datasets in the YRB from 2001 to 2015. (<b>a</b>) Cropland, (<b>b</b>) forest, (<b>c</b>) shrubland, (<b>d</b>) grassland, (<b>e</b>) water, (<b>f</b>) snow/ice, (<b>g</b>) barren, (<b>h</b>) impervious, (<b>i</b>) wetland.</p>
Full article ">Figure 5
<p>Comparison of the spatial variation of the Sen’s slope estimate (<span class="html-italic">β</span>, M ha/a) for nine LC classes between the four datasets in the YRB from 2001 to 2015. Values are summary results in a 0.5° × 0.5° grid, blank areas have no values, and the inset graph shows the differences in <span class="html-italic">β</span> for each LC class in the SA, UA, MA, and DA watersheds, respectively.</p>
Full article ">Figure 6
<p>Comparison of the spatial distribution of the correlation (<span class="html-italic">R</span> &gt; 0, <span class="html-italic">p</span> &lt; 0.05) for nine LC classes between the four other datasets in the YRB from 2001 to 2015. Values are summary results in a 0.5° × 0.5° grid, blank areas have no values, and the inset graph shows the percentage of correlated grids for each LC class in the SA, UA, MA, and DA watersheds, respectively.</p>
Full article ">Figure 7
<p>Comparison of the spatial variation of the percentage of systematic and random error for nine LC classes between the four datasets in the YRB from 2001 to 2015. Values are summary results in a 0.5° × 0.5° grid, blank areas have no values, and the inset graph shows the percentage of systematic and random error for each LC class in the SA, UA, MA, and DA watersheds, respectively.</p>
Full article ">Figure 8
<p>Comparison of the variations in spatial agreement (<span class="html-italic">R</span> &gt; 0, <span class="html-italic">p</span> &lt; 0.05) (<b>left</b> panel) and total errors (<b>right</b> panel) with landscape diversity (<span class="html-italic">H</span>) for nine LC classes between the four datasets in the YRB from 2001 to 2015.</p>
Full article ">
16 pages, 5294 KiB  
Article
Improving the Spatial Prediction of Soil Organic Carbon Content Using Phenological Factors: A Case Study in the Middle and Upper Reaches of Heihe River Basin, China
by Xinyu Liu, Jian Wang and Xiaodong Song
Remote Sens. 2023, 15(7), 1847; https://doi.org/10.3390/rs15071847 - 30 Mar 2023
Cited by 1 | Viewed by 1635
Abstract
The accurate mapping of soil organic carbon (SOC) distribution is important for carbon sequestration and land management strategies, contributing to mitigating climate change and ensuring agricultural productivity. The Heihe River Basin in China is an important region that has immense potential for SOC [...] Read more.
The accurate mapping of soil organic carbon (SOC) distribution is important for carbon sequestration and land management strategies, contributing to mitigating climate change and ensuring agricultural productivity. The Heihe River Basin in China is an important region that has immense potential for SOC storage. Phenological variables are effective indicators of vegetation growth, and hence are closely related to SOC. However, few studies have incorporated phenological variables in SOC prediction, especially in alpine areas such as the Heihe River Basin. This study used random forest (RF) and extreme gradient boosting (XGBoost) to study the effects of phenological variables (e.g., Greenup, Dormancy, etc.) obtained from MODIS (i.e., Moderate Resolution Imaging Spectroradiometer) product (MCD12Q2) on SOC content prediction in the middle and upper reaches of Heihe River Basin. The current study also identified the dominating variables in SOC prediction and compared model performance using a cross validation procedure. The results indicate that: (1) when phenological variables were considered, the R2 (coefficient of determination) of RF and XGBoost were 0.68 and 0.56, respectively, and RF consistently outperforms XGBoost in various cross validation experiments; (2) the environmental variables MAT, MAP, DEM and NDVI play the most important roles in SOC prediction; (3) the phenological variables can account for 32–39% of the spatial variability of SOC in both the RF and XGBoost models, and hence were the most important factor among the five categories of predictive variables. This study proved that the introduction of phenological variables can significantly improve the performance of SOC prediction. They should be used as indispensable variables for accurately modeling SOC in related studies. Full article
Show Figures

Figure 1

Figure 1
<p>Location of study area and distribution of sample points.</p>
Full article ">Figure 2
<p>The spatial distribution of phenological variables Greenup (<b>a</b>) and Dormancy (<b>b</b>) in 2010. Note that values of the phenological variables represent accumulated days since 1 January 1970.</p>
Full article ">Figure 3
<p>Frequency distribution of SOC (<b>a</b>) and lnSOC (<b>b</b>).</p>
Full article ">Figure 4
<p>Model accuracy under <span class="html-italic">k</span>-fold cross validation with commonly used variables (<b>a</b>) and commonly used variables and phenological variables (<b>b</b>).</p>
Full article ">Figure 5
<p>The histogram of model accuracy of <math display="inline"><semantics> <mrow> <msup> <mi mathvariant="normal">R</mi> <mn>2</mn> </msup> </mrow> </semantics></math> (<b>a</b>), RMSE (<b>b</b>), and MAE (<b>c</b>), when different environmental variables were used.</p>
Full article ">Figure 6
<p>Model fitting result by RF (<b>a</b>) and XGBoost (<b>b</b>) when phenological variables and normal variables were used.</p>
Full article ">Figure 7
<p>The spatial distribution of SOC content predicted by RF using commonly used variables (<b>a</b>) and commonly-used variables and phenological variables (<b>b</b>).</p>
Full article ">Figure 8
<p>Importance ranking of environmental variables by RF using commonly used variables (<b>a</b>) and commonly used variables and phenological variables (<b>b</b>).</p>
Full article ">Figure 9
<p>The mean SOC content and relative total SOC content of different land use types (<b>a</b>) and the spatial distribution of different land use types (<b>b</b>).</p>
Full article ">
22 pages, 6429 KiB  
Article
Monitoring of Land Degradation in Greece and Tunisia Using Trends.Earth with a Focus on Cereal Croplands
by Ines Cherif, Eleni Kolintziki and Thomas K. Alexandridis
Remote Sens. 2023, 15(7), 1766; https://doi.org/10.3390/rs15071766 - 25 Mar 2023
Cited by 4 | Viewed by 2626
Abstract
Land degradation (LD) processes are widespread in drylands worldwide and are accelerated by climate change. As a result, food security and livelihoods are at risk. Thus, there is a need to monitor LD trends, especially in agricultural areas. Mediterranean countries, including Tunisia and [...] Read more.
Land degradation (LD) processes are widespread in drylands worldwide and are accelerated by climate change. As a result, food security and livelihoods are at risk. Thus, there is a need to monitor LD trends, especially in agricultural areas. Mediterranean countries, including Tunisia and Greece, are concerned due to the presence of drivers and pressures causing land degradation. Through the Trends.Earth plugin, the SDG 15.3.1 indicator can be implemented to map LD status. In this study, we mapped LD in Greece and Tunisia for the recommended baseline period of 2001–2015 and the selected reporting period of 2016–2020. The land productivity was assessed within Trends.Earth using the MODIS MOD13Q1 product, while the default datasets were used for the other sub-indicators. The main findings are: (i) the percentage of degraded land decreased from the baseline to the reporting period from 4.83% to 2.62% of total area in Greece and 9.97% to 6.26% in Tunisia—degradation rates that differ from those reported to the UNCCD (United Nations Convention to Combat Desertification) by the respective national authorities; (ii) the dominant land condition in Greece was improved, while in Tunisia, it was stable; (iii) land productivity presented a similar trend through the SDG 15.3.1 indicator over both countries, including the net land productivity dynamics over croplands; (iv) based on analysis using plant functional types performed with MODIS MCD12Q1, the highest portion of degraded land in Greece was located in grasslands and in Tunisia in cereal croplands (after desert areas); and (v) with a focus on LD over cereal croplands, the portion of degraded areas appeared to decrease in both Greece and Tunisia. The percentage was higher in Tunisia, representing 16.52% of the total degraded land during the reporting period compared to 10.83% in Greece. All the above stress the need to foster the adoption of sustainable land management practices, especially in Tunisia, and speed up the implementation of measures to achieve LD neutrality. Full article
Show Figures

Figure 1

Figure 1
<p>Location of the study areas.</p>
Full article ">Figure 2
<p>Spatial distribution of SDG 15.3.1 indicator in Greece during the baseline period (<b>left</b>) and the reporting period (<b>right</b>).</p>
Full article ">Figure 3
<p>Spatial distribution of SDG 15.3.1 indicator in Tunisia during the baseline period (<b>left</b>) and the reporting period (<b>right</b>).</p>
Full article ">Figure 4
<p>Spatial distribution of land productivity in Greece during the baseline period (<b>left</b>) and the reporting period (<b>right</b>).</p>
Full article ">Figure 5
<p>Spatial distribution of land productivity in Tunisia during the baseline period (<b>left</b>) and the reporting period (<b>right</b>).</p>
Full article ">Figure 6
<p>Distribution map of functional types of plants that have remained unchanged during the baseline period (<b>left</b>: 2001–2015) and during the reporting period (<b>right</b>: 2016–2020). In Greece (<b>top</b>) and in Tunisia (<b>bottom</b>).</p>
Full article ">Figure 7
<p>Spatial distribution of SDG 15.3.1 indicator in cereal croplands during the reporting period in Greece (<b>right</b>) and in Tunisia (<b>left</b>).</p>
Full article ">Figure 8
<p>Spatial distribution of JRC land productivity in 2005–2019 in Greece (<b>left</b>) and Tunisia (<b>right</b>).</p>
Full article ">
18 pages, 8558 KiB  
Article
Microwave Dual-Crack Sensor with a High Q-Factor Using the TE20 Resonance of a Complementary Split-Ring Resonator on a Substrate-Integrated Waveguide
by Yelim Kim, Eiyong Park, Ahmed Salim, Junghyeon Kim and Sungjoon Lim
Micromachines 2023, 14(3), 578; https://doi.org/10.3390/mi14030578 - 28 Feb 2023
Cited by 4 | Viewed by 1943
Abstract
Microwave sensors have attracted interest as non-destructive metal crack detection (MCD) devices due to their low cost, simple fabrication, potential miniaturization, noncontact nature, and capability for remote detection. However, the development of multi-crack sensors of a suitable size and quality factor (Q-factor) remains [...] Read more.
Microwave sensors have attracted interest as non-destructive metal crack detection (MCD) devices due to their low cost, simple fabrication, potential miniaturization, noncontact nature, and capability for remote detection. However, the development of multi-crack sensors of a suitable size and quality factor (Q-factor) remains a challenge. In the present study, we propose a multi-MCD sensor that combines a higher-mode substrate-integrated waveguide (SIW) and complementary split-ring resonators (CSRRs). In order to increase the Q-factor, the device is miniaturized; the MCD is facilitated; and two independent CSRRs are loaded onto the SIW, where the electromagnetic field is concentrated. The concentrated electromagnetic field of the SIW improves the Q-factor of the CSRRs, and each CSRR creates its own resonance and produces a miniaturizing effect by activating the sensor below the cut-off frequency of the SIW. The proposed multi-MCD sensor is numerically and experimentally demonstrated for cracks with different widths and depths. The fabricated sensor with a TE20-mode SIW and CSRRs is able to efficiently detect two sub-millimeter metal cracks simultaneously with a high Q-factor of 281. Full article
(This article belongs to the Special Issue Design and Fabrication of Micro/Nano Sensors and Actuators, Volume II)
Show Figures

Figure 1

Figure 1
<p>Proposed multi-crack detecting sensor with a high Q-factor: (<b>a</b>) TE<sub>20</sub> mode SIW only E-field vector distribution. E-field vector distribution after the integration of the CSRRs (<b>b</b>) at 4.69 GHz and (<b>c</b>) 5.49 GHz. The E-field intensity range of (<b>a</b>–<b>c</b>) is the same. Concept illustration of (<b>d</b>) independent multi-crack detection and (<b>e</b>) the reflective coefficients before and after CSRR loading.</p>
Full article ">Figure 2
<p>Proposed dual-crack detection sensor: (<b>a</b>) top view, (<b>b</b>) CSRR structure and equivalent circuit, (<b>c</b>) bottom view, and (<b>d</b>) side view.</p>
Full article ">Figure 3
<p>Simulation results for the proposed crack sensor: (<b>a</b>) S11 without an Al sheet and with an uncracked Al sheet, (<b>b</b>) S11 for different crack numbers and positions (two cracks, single crack at CSRR1, single crack at CSRR2, no crack), S11 results (<b>c</b>) in the CSRR2 resonance frequency band, and (<b>d</b>) in the CSRR1 resonance frequency band (crack width, depth, and length of 0.5 <math display="inline"><semantics> <mo>×</mo> </semantics></math> 0.5 <math display="inline"><semantics> <mo>×</mo> </semantics></math> 6.1 mm).</p>
Full article ">Figure 4
<p>Simulated S11 results for the proposed dual-crack sensor for different crack widths above (<b>a</b>) CSRR1 and (<b>b</b>) CSRR2. (<b>c</b>) Fundamental metal sheet structure with two cracks with a width, depth, and length of 0.5 <math display="inline"><semantics> <mo>×</mo> </semantics></math> 0.9 <math display="inline"><semantics> <mo>×</mo> </semantics></math> 6.1 mm. S11 results for different crack depths above (<b>d</b>) CSRR1 and (<b>e</b>) CSRR2. (<b>f</b>) S11 results for the rotation of a straight crack above CSRR1. (<b>g</b>) S11 results for five crack shapes above CSRR2.</p>
Full article ">Figure 5
<p>Simulation results for different crack locations (<span class="html-italic">x</span>, <span class="html-italic">y</span>) on the Al sheet: (<b>a</b>) coordinate of the crack; (<b>b</b>–<b>e</b>) S11 when the crack moves <span class="html-italic">x</span>-axis; (<b>b</b>) from (−9, 0) to (-2, 0), (<b>c</b>) from (−1, 0) to (8, 0), (<b>d</b>) from (−9, 3) to (8, 3), and (<b>e</b>) from (−9, 4) to (8, 4); (<b>f</b>) S11 when the crack move <span class="html-italic">y</span>-axis from (4, −6) to (4, 6).</p>
Full article ">Figure 6
<p>Simulation result of the proposed sensor for different depths from a metal surface: (<b>a</b>) side view and top view of structure and crack; (<b>b</b>,<b>c</b>) the simulated S<sub>11</sub> for different crack positions: depth from a metal surface D<sub>ud</sub> = 0–2.5 mm.</p>
Full article ">Figure 7
<p>Simulation results of the proposed sensor for (<b>a</b>) different Al sheet thicknesses and (<b>b</b>,<b>c</b>) crack lengths.</p>
Full article ">Figure 8
<p>Fabricated prototype of the proposed dual-crack sensor: (<b>a</b>) top view, (<b>b</b>) bottom view, and (<b>c</b>) side view with the Al sheet. Simulation and measurement results for (<b>d</b>) the sensor with no Al sheet and (<b>e</b>) the sensor with an uncracked Al sheet.</p>
Full article ">Figure 9
<p>Measured S-parameters for the proposed dual-crack sensor (<b>a</b>) with various slot widths and depths at (<b>b</b>) CSRR 2 and (<b>c</b>) CSRR 1. Measured relationship between the frequency and various crack sizes: (<b>d</b>) depth and (<b>e</b>) width.</p>
Full article ">
Back to TopTop