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

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

Years

Between: -

Subjects

remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline

Journals

Article Types

Countries / Regions

Search Results (142)

Search Parameters:
Keywords = FY3D

Order results
Result details
Results per page
Select all
Export citation of selected articles as:
24 pages, 6494 KiB  
Article
Reconstruction of Fine-Spatial-Resolution FY-3D-Based Vegetation Indices to Achieve Farmland-Scale Winter Wheat Yield Estimation via Fusion with Sentinel-2 Data
by Xijia Zhou, Tao Wang, Wei Zheng, Mingwei Zhang and Yuanyuan Wang
Remote Sens. 2024, 16(22), 4143; https://doi.org/10.3390/rs16224143 - 6 Nov 2024
Viewed by 512
Abstract
The spatial resolution (250–1000 m) of the FY-3D MERSI is too coarse for agricultural monitoring at the farmland scale (20–30 m). To achieve the winter wheat yield (WWY) at the farmland scale, based on FY-3D, a method framework is developed in this work. [...] Read more.
The spatial resolution (250–1000 m) of the FY-3D MERSI is too coarse for agricultural monitoring at the farmland scale (20–30 m). To achieve the winter wheat yield (WWY) at the farmland scale, based on FY-3D, a method framework is developed in this work. The enhanced deep convolutional spatiotemporal fusion network (EDCSTFN) was used to perform a spatiotemporal fusion on the 10 day interval FY-3D and Sentinel-2 vegetation indices (VIs), which were compared with the enhanced spatial and temporal adaptive reflectance fusion model (ESTARFM). In addition, a BP neural network was built to calculate the farmland-scale WWY based on the fused VIs, and the Aqua MODIS gross primary productivity product was used as ancillary data for WWY estimation. The results reveal that both the EDCSTFN and ESTARFM achieve satisfactory precision in the fusion of the Sentinel-2 and FY-3D VIs; however, when the period of spatiotemporal data fusion is relatively long, the EDCSTFN can achieve greater precision than ESTARFM. Finally, the WWY estimation results based on the fused VIs show remarkable correlations with the WWY data at the county scale and provide abundant spatial distribution details about the WWY, displaying great potential for accurate farmland-scale WWY estimations based on reconstructed fine-spatial-temporal-resolution FY-3D data. Full article
Show Figures

Figure 1

Figure 1
<p>Overview of the study region: (<b>a</b>) location of the Weihe Plain; (<b>b</b>) FY-3D false colour composite image for 3 May 2020; and (<b>c</b>) locations of the county-scale WWY data points used in the WWY estimation.</p>
Full article ">Figure 2
<p>Flowchart of the 10 day interval VI imagery reconstruction and farmland-scale WWY estimation.</p>
Full article ">Figure 3
<p>Flowchart of farmland-scale WWY estimation: (<b>a</b>) Y estimation model based on the cumulative GPP; and (<b>b</b>) farmland-scale Y estimation model based on multiple parameters.</p>
Full article ">Figure 4
<p>Results of the consistency analysis of the Sentinel-2 and FY-3D VIs: (<b>a</b>) R<sup>2</sup> values between the aggregated Sentinel-2 VI imagery and the FY-3D VI imagery at an SR of 250 m; and (<b>b</b>) average deviations and RMSE values of the fitting results between the aggregated Sentinel-2 VI imagery and FY-3D VI imagery. The error line in (<b>b</b>) denotes the RMSE of the fitting results.</p>
Full article ">Figure 5
<p>EVI at each WW growing stage from 2020 to 2022.</p>
Full article ">Figure 6
<p>Y estimation model based on the cumulative GPP for the main WW growing period and the Y estimation precision evaluation results in 2020 and 2021. The dotted lines in the figures denote the fitted linear functions, which are close to the diagonal solid lines, indicating that the systematic deviation in the Y estimation results is small. (<b>a</b>) Linear regression model between the cumulative GPP data for the main WW growing period and the county-scale WWY from 2014 to 2018, (<b>b</b>) linear regression results between the WWY estimation results from 2020 based on the cumulative GPP and county-scale Y statistical data, and (<b>c</b>) linear regression results between the WWY estimation results in 2021 based on the cumulative GPP and county-scale Y statistical data.</p>
Full article ">Figure 7
<p>WWY estimation results for 2020 to 2022 based on the MODIS cumulative GPP data.</p>
Full article ">Figure 8
<p>Farmland-scale WWY estimation results for the Weihe Plain from 2020 to 2022 based on multiple parameters.</p>
Full article ">Figure 9
<p>Linear regression results between the farmland-scale WWY estimation results and the Y statistical data in 2020 and 2021. The dotted lines in the figures denote the fitted linear functions, which are close to the diagonal solid lines, indicating that the systematic deviation of the Y estimation results is small.</p>
Full article ">
20 pages, 5140 KiB  
Article
Distribution-Based Approach for Efficient Storage and Indexing of Massive Infrared Hyperspectral Sounding Data
by Han Li, Mingjian Gu, Guang Shi, Yong Hu and Mengzhen Xie
Remote Sens. 2024, 16(21), 4088; https://doi.org/10.3390/rs16214088 - 1 Nov 2024
Viewed by 416
Abstract
Hyperspectral infrared atmospheric sounding data, characterized by their high vertical resolution, play a crucial role in capturing three-dimensional atmospheric spatial information. The hyperspectral infrared atmospheric detectors HIRAS/HIRAS-II, mounted on the FY3D/EF satellite, have established an initial global coverage network for atmospheric sounding. The [...] Read more.
Hyperspectral infrared atmospheric sounding data, characterized by their high vertical resolution, play a crucial role in capturing three-dimensional atmospheric spatial information. The hyperspectral infrared atmospheric detectors HIRAS/HIRAS-II, mounted on the FY3D/EF satellite, have established an initial global coverage network for atmospheric sounding. The collaborative observation approach involving multiple satellites will improve both the coverage and responsiveness of data acquisition, thereby enhancing the overall quality and reliability of the data. In response to the increasing number of channels, the rapid growth of data volume, and the specific requirements of multi-satellite joint observation applications with infrared hyperspectral sounding data, this paper introduces an efficient storage and indexing method for infrared hyperspectral sounding data within a distributed architecture for the first time. The proposed approach, built on the Kubernetes cloud platform, utilizes the Google S2 discrete grid spatial indexing algorithm to establish a grid-based hierarchical model for unified metadata-embedded documents. Additionally, it optimizes the rowkey design using the BPDS model, thereby enabling the distributed storage of data in HBase. The experimental results demonstrate that the query efficiency of the Google S2 grid-based embedded document model is superior to that of the traditional flat model, achieving a query time that is only 35.6% of the latter for a dataset of 5 million records. Additionally, this method exhibits better data distribution characteristics within the global grid compared to the H3 algorithm. Leveraging the BPDS model, the HBase distributed storage system adeptly balances the node load and counteracts the detrimental effects caused by the accumulation of time-series remote sensing images. This architecture significantly enhances both storage and query efficiency, thus laying a robust foundation for forthcoming distributed computing. Full article
Show Figures

Figure 1

Figure 1
<p>Illustration of ground field-of-view arrangement for the HIRAS series detector [<a href="#B36-remotesensing-16-04088" class="html-bibr">36</a>,<a href="#B37-remotesensing-16-04088" class="html-bibr">37</a>].</p>
Full article ">Figure 2
<p>The system architecture for the efficient storage and indexing of hyperspectral infrared sounding data.</p>
Full article ">Figure 3
<p>A schematic diagram of the design structure of rowkeys.</p>
Full article ">Figure 4
<p>The table structure of infrared hyperspectral sounding data in HBase.</p>
Full article ">Figure 5
<p>Collections on MongoDB cluster.</p>
Full article ">Figure 6
<p>The relationship between coding level and grid resolution.</p>
Full article ">Figure 7
<p>Comparison of loading times by various index methods.</p>
Full article ">Figure 8
<p>Comparison of storage structure in MongoDB.</p>
Full article ">Figure 9
<p>Grid-based heatmap visualization of different methods: (<b>a</b>) Google S2 algorithm; (<b>b</b>) H3 algorithm.</p>
Full article ">Figure 10
<p>The histogram of pixels of grids in different methods: (<b>a</b>) Google S2 algorithm; (<b>b</b>) H3 algorithm.</p>
Full article ">Figure 11
<p>Query time and storage time are utilized as comparative metrics to assess performance by different methods: (<b>a</b>) storage time and (<b>b</b>) query time.</p>
Full article ">Figure 12
<p>Distribution of regions in different models: (<b>a</b>) DP model, (<b>b</b>) FP model, and (<b>c</b>) BPDS model.</p>
Full article ">
35 pages, 7235 KiB  
Article
Change in Fractional Vegetation Cover and Its Prediction during the Growing Season Based on Machine Learning in Southwest China
by Xiehui Li, Yuting Liu and Lei Wang
Remote Sens. 2024, 16(19), 3623; https://doi.org/10.3390/rs16193623 - 28 Sep 2024
Viewed by 607
Abstract
Fractional vegetation cover (FVC) is a crucial indicator for measuring the growth of surface vegetation. The changes and predictions of FVC significantly impact biodiversity conservation, ecosystem health and stability, and climate change response and prediction. Southwest China (SWC) is characterized by complex topography, [...] Read more.
Fractional vegetation cover (FVC) is a crucial indicator for measuring the growth of surface vegetation. The changes and predictions of FVC significantly impact biodiversity conservation, ecosystem health and stability, and climate change response and prediction. Southwest China (SWC) is characterized by complex topography, diverse climate types, and rich vegetation types. This study first analyzed the spatiotemporal variation of FVC at various timescales in SWC from 2000 to 2020 using FVC values derived from pixel dichotomy model. Next, we constructed four machine learning models—light gradient boosting machine (LightGBM), support vector regression (SVR), k-nearest neighbor (KNN), and ridge regression (RR)—along with a weighted average heterogeneous ensemble model (WAHEM) to predict growing-season FVC in SWC from 2000 to 2023. Finally, the performance of the different ML models was comprehensively evaluated using tenfold cross-validation and multiple performance metrics. The results indicated that the overall FVC in SWC predominantly increased from 2000 to 2020. Over the 21 years, the FVC spatial distribution in SWC generally showed a high east and low west pattern, with extremely low FVC in the western plateau of Tibet and higher FVC in parts of eastern Sichuan, Chongqing, Guizhou, and Yunnan. The determination coefficient R2 scores from tenfold cross-validation for the four ML models indicated that LightGBM had the strongest predictive ability whereas RR had the weakest. WAHEM and LightGBM models performed the best overall in the training, validation, and test sets, with RR performing the worst. The predicted spatial change trends were consistent with the MODIS-MOD13A3-FVC and FY3D-MERSI-FVC, although the predicted FVC values were slightly higher but closer to the MODIS-MOD13A3-FVC. The feature importance scores from the LightGBM model indicated that digital elevation model (DEM) had the most significant influence on FVC among the six input features. In contrast, soil surface water retention capacity (SSWRC) was the most influential climate factor. The results of this study provided valuable insights and references for monitoring and predicting the vegetation cover in regions with complex topography, diverse climate types, and rich vegetation. Additionally, they offered guidance for selecting remote sensing products for vegetation cover and optimizing different ML models. Full article
Show Figures

Figure 1

Figure 1
<p>Geographical location (<b>a</b>), elevation distribution (<b>b</b>) and the average NDVI values for 2020 (<b>c</b>) in Southwest China.</p>
Full article ">Figure 2
<p>Density scatter plots of the annual average FVC during 2000–2020 between the FVC dataset calculated in this study and the published GLASS FVC product. (<b>a</b>) 2000, (<b>b</b>) 2007, (<b>c</b>) 2014, (<b>d</b>) 2020.</p>
Full article ">Figure 3
<p>Spatial distribution of FVC across different timescales from 2000 to 2020 in southwest China. (<b>a</b>) Annual, (<b>b</b>) growing season, (<b>c</b>) spring, (<b>d</b>) summer, (<b>e</b>) autumn, and (<b>f</b>) winter.</p>
Full article ">Figure 4
<p>Spatial change trend of FVC across different timescales from 2000 to 2020 in southwest China. (<b>a</b>) Annual, (<b>b</b>) growing season, (<b>c</b>) spring, (<b>d</b>) summer, (<b>e</b>) autumn, and (<b>f</b>) winter.</p>
Full article ">Figure 5
<p>Spatial variation of FVC classification in growing season from 2000 to 2020. (<b>a</b>) 2000, (<b>b</b>) 2007, (<b>c</b>) 2014, (<b>d</b>) 2020.</p>
Full article ">Figure 6
<p>The area ratio changes in FVC classification during the growing season from 2000 to 2020.</p>
Full article ">Figure 7
<p><span class="html-italic">R</span><sup>2</sup> score from tenfold cross-validation for various models.</p>
Full article ">Figure 8
<p>Spatial distribution of FVC predicted by different ML models and retrieved calculation from satellite products during the growing seasons from 2021 to 2023. (<b>a</b>) LightGBM-2021, (<b>b</b>) LightGBM-2022, (<b>c</b>) LightGBM-2023, (<b>d</b>) SVR-2021, (<b>e</b>) SVR-2022, (<b>f</b>) SVR-2023, (<b>g</b>) KNN-2021, (<b>h</b>) KNN-2022, (<b>i</b>) KNN-2023, (<b>j</b>) RR-2021, (<b>k</b>) RR-2022, (<b>l</b>) RR-2023, (<b>m</b>) WAHEM-2021, (<b>n</b>) WAHEM-2022, (<b>o</b>) WAHEM-2023, (<b>p</b>) MOD13A3-2021, (<b>q</b>) MOD13A3-2022, (<b>r</b>) MOD13A3-2023, (<b>s</b>) FY3D-MERSI-2021, (<b>t</b>) FY3D-MERSI-2022, (<b>u</b>) FY3D-MERSI-2023.</p>
Full article ">Figure 9
<p>Comparison of FVC changes derived from different ML model predictions and satellite product calculations across various regions during the growing seasons from 2021 to 2023. (<b>a</b>) SWC, (<b>b</b>) Sichuan, (<b>c</b>) Chongqing, (<b>d</b>) Guizhou, (<b>e</b>) Yunnan, and (<b>f</b>) Tibet.</p>
Full article ">Figure 10
<p>Feature importance scores in the LightGBM model.</p>
Full article ">
21 pages, 9876 KiB  
Article
Estimation of Leaf Area Index across Biomes and Growth Stages Combining Multiple Vegetation Indices
by Fangyi Lv, Kaimin Sun, Wenzhuo Li, Shunxia Miao and Xiuqing Hu
Sensors 2024, 24(18), 6106; https://doi.org/10.3390/s24186106 - 21 Sep 2024
Viewed by 781
Abstract
The leaf area index (LAI) is a key indicator of vegetation canopy structure and growth status, crucial for global ecological environment research. The Moderate Resolution Spectral Imager-II (MERSI-II) aboard Fengyun-3D (FY-3D) covers the globe twice daily, providing a reliable data source for large-scale [...] Read more.
The leaf area index (LAI) is a key indicator of vegetation canopy structure and growth status, crucial for global ecological environment research. The Moderate Resolution Spectral Imager-II (MERSI-II) aboard Fengyun-3D (FY-3D) covers the globe twice daily, providing a reliable data source for large-scale and high-frequency LAI estimation. VI-based LAI estimation is effective, but species and growth status impacts on the sensitivity of the VI–LAI relationship are rarely considered, especially for MERSI-II. This study analyzed the VI–LAI relationship for eight biomes in China with contrasting leaf structures and canopy architectures. The LAI was estimated by adaptively combining multiple VIs and validated using MODIS, GLASS, and ground measurements. Results show that (1) species and growth stages significantly affect VI–LAI sensitivity. For example, the EVI is optimal for broadleaf crops in winter, while the RDVI is best for evergreen needleleaf forests in summer. (2) Combining vegetation indices can significantly optimize sensitivity. The accuracy of multi-VI-based LAI retrieval is notably higher than using a single VI for the entire year. (3) MERSI-II shows good spatial–temporal consistency with MODIS and GLASS and is more sensitive to vegetation growth fluctuation. Direct validation with ground-truth data also demonstrates that the uncertainty of retrievals is acceptable (R2 = 0.808, RMSE = 0.642). Full article
(This article belongs to the Section Remote Sensors)
Show Figures

Figure 1

Figure 1
<p>Land cover map of the study area using MODIS MCD12Q1 product in 2020.</p>
Full article ">Figure 2
<p>Flowchart for MERSI-II LAI estimation.</p>
Full article ">Figure 3
<p>Time-series curves of VIs and LAI for grasses/cereal crops.</p>
Full article ">Figure 4
<p>Time-series curves of VIs and LAI for broadleaf crops.</p>
Full article ">Figure 5
<p>Time-series curves of VIs and LAI for savannahs.</p>
Full article ">Figure 6
<p>Time-series curves of VIs and LAI for deciduous broadleaf forests.</p>
Full article ">Figure 7
<p>VI–LAI density scatter plot of grasses/cereal crops: (<b>a</b>) DVI–LAI; (<b>b</b>) EVI–LAI; (<b>c</b>) MSR-LAI; (<b>d</b>) NDVI–LAI; (<b>e</b>) OSAVI–LAI; (<b>f</b>) RDVI–LAI; (<b>g</b>) RVI–LAI; (<b>h</b>) SAVI–LAI.</p>
Full article ">Figure 7 Cont.
<p>VI–LAI density scatter plot of grasses/cereal crops: (<b>a</b>) DVI–LAI; (<b>b</b>) EVI–LAI; (<b>c</b>) MSR-LAI; (<b>d</b>) NDVI–LAI; (<b>e</b>) OSAVI–LAI; (<b>f</b>) RDVI–LAI; (<b>g</b>) RVI–LAI; (<b>h</b>) SAVI–LAI.</p>
Full article ">Figure 8
<p>Optimal VIs for MERSI-II LAI estimation across different biomes and growth stages. Biomes 1–8 correspond to grasses/cereal crops, shrubs, broadleaf crops, savanna, EBF, DBF, ENF, and ENF, respectively. The seven colors represent seven different input parameters, respectively.</p>
Full article ">Figure 9
<p>Estimation results of different methods.</p>
Full article ">Figure 10
<p>Comparison of spatial distributions of LAI differences between MERSI-II, Aqua MODIS and GLASS in mainland China in 2020: (<b>a</b>) MERSI-II LAI minus MODIS LAI; (<b>b</b>) STD of the LAI differences between MERSI-II and MODIS; (<b>c</b>) MERSI-II LAI minus GLASS LAI; (<b>d</b>) STD of the LAI differences between MERSI-II and GLASS. No-data pixels in white color are observations contaminated by cloud, shadow, aerosol, etc. Gray pixels are non-vegetation areas.</p>
Full article ">Figure 11
<p>Bar chart for the proportion of different categories of LAI differences under each biome: (<b>a</b>) MERSI-II LAI minus MODIS LAI; (<b>b</b>) MERSI-II LAI minus GLASS LAI. Biomes 1–8 are grasses/cereal crops, shrubs, broadleaf crops, savanna, EBF, DBF, ENF, and ENF.</p>
Full article ">Figure 12
<p>Bar chart for the proportion of different biomes under each category of LAI difference: (<b>a</b>) MERSI-II LAI minus MODIS LAI; (<b>b</b>) MERSI-II LAI minus GLASS LAI. Biomes 1–8 are grasses/cereal crops, shrubs, broadleaf crops, savanna, EBF, DBF, ENF and ENF.</p>
Full article ">Figure 13
<p>Time series of MERSI, GLASS and MODIS (2020). The red circle indicates large fluctuations in MODIS LAI.</p>
Full article ">Figure 14
<p>Comparison between ground-truth data and FY-3D MERSI-II LAI. Biome 1 is grasses/cereal crops, Biome 2 is shrubs, Biome 4 is savanna, Biome 5 is evergreen broadleaf forests, and Biome 8 is deciduous needleleaf forests.</p>
Full article ">
17 pages, 4008 KiB  
Article
Phenotypic and Genetic Analyses of Mastitis, Endometritis, and Ketosis on Milk Production and Reproduction Traits in Chinese Holstein Cattle
by Xiaoli Ren, Haibo Lu, Yachun Wang, Lei Yan, Changlei Liu, Chu Chu, Zhuo Yang, Xiangnan Bao, Mei Yu, Zhen Zhang and Shujun Zhang
Animals 2024, 14(16), 2372; https://doi.org/10.3390/ani14162372 - 15 Aug 2024
Viewed by 735
Abstract
Mastitis (MAS), endometritis (MET), and ketosis (KET) are prevalent diseases in dairy cows that result in substantial economic losses for the dairy farming industry. This study gathered 26,014 records of the health and sickness of dairy cows and 99,102 data of reproduction from [...] Read more.
Mastitis (MAS), endometritis (MET), and ketosis (KET) are prevalent diseases in dairy cows that result in substantial economic losses for the dairy farming industry. This study gathered 26,014 records of the health and sickness of dairy cows and 99,102 data of reproduction from 13 Holstein dairy farms in Central China; the milk protein and milk fat content from 56,640 milk samples, as well as the pedigree data of 37,836 dairy cows were obtained. The logistic regression method was used to analyze the variations in the prevalence rates of MAS, MET, and KET among various parities; the mixed linear model was used to examine the effects of the three diseases on milk production, milk quality, and reproductive traits. DMU software (version 5.2) utilized the DMUAI module in conjunction with the single-trait and two-trait animal model, as well as best linear unbiased prediction (BLUP), to estimate the genetic parameters for the three diseases, milk production, milk quality, and reproductive traits in dairy cows. The primary findings of the investigation comprised the following: (1) The prevalence rates of MAS, MET, and KET in dairy farms were 20.04%, 10.68%, and 7.33%, respectively. (2) MAS and MET had a substantial impact (p < 0.01) on milk production, resulting in significant decreases of 112 kg and 372 kg in 305-d Milk Yield (305-d MY), 4 kg and 12 kg in 305-d Protein Yield (305-d PY), and 6 kg and 16 kg in 305-d Fat Yield (305-d FY). As a result of their excessive 305-d MY, some cows were diagnosed with KET due to glucose metabolism disorder. The 305-d MY of cows with KET was significantly higher than that of healthy cows (205 kg, p < 0.01). (3) All three diseases resulted in an increase in the Interval from Calving to First Service (CTFS, 0.60–1.50 d), Interval from First Service to Conception (FSTC, 0.20–16.20 d), Calving Interval (CI, 4.00–7.00 d), and Number of Services (NUMS, 0.07–0.35). (4) The heritabilities of cows with MAS, MET, and KET were found to be low, with values of 0.09, 0.01, and 0.02, respectively. The genetic correlation between these traits ranged from 0.14 to 0.44. This study offers valuable insights on the prevention and control of the three diseases, as well as feeding management and genetic breeding. Full article
(This article belongs to the Collection Advances in Cattle Breeding, Genetics and Genomics)
Show Figures

Figure 1

Figure 1
<p>The prevalence rates and odds ratios (ORs) of mastitis in different parities, compared to the first parity. The MAS Rate represents the mastitis prevalence rate. “**” represents no and highly significant differences (<span class="html-italic">p</span> &lt; 0.01).</p>
Full article ">Figure 2
<p>The prevalence rates and odds ratios (ORs) of endometritis in different parities, compared to the first parity. The MET Rate represents the endometritis prevalence rate. NS and “**”, respectively, represent no and highly significant differences (<span class="html-italic">p &gt;</span> 0.05, <span class="html-italic">p &lt;</span> 0.01).</p>
Full article ">Figure 3
<p>The prevalence rates and odds ratios (ORs) of ketosis in different parities, compared to the first parity. The KET Rate represents the ketosis prevalence rate. “**” represents no and highly significant differences (<span class="html-italic">p</span> &lt; 0.01).</p>
Full article ">Figure 4
<p>The effects of mastitis, endometritis, and ketosis on milk production in dairy cows. (<b>a</b>) 305-d MY represents 305-d Milk Yield, (<b>b</b>) 305-d PY represents 305-d Protein Yield, and (<b>c</b>) 305-d FY represents 305-d Fat Yield. MAS represents mastitis, MET represents endometritis, and KET represents ketosis. LSM represents the least squares means. NS and “**”, respectively, represent no and highly significant differences (<span class="html-italic">p &gt;</span> 0.05, <span class="html-italic">p</span> &lt; 0.01).</p>
Full article ">Figure 4 Cont.
<p>The effects of mastitis, endometritis, and ketosis on milk production in dairy cows. (<b>a</b>) 305-d MY represents 305-d Milk Yield, (<b>b</b>) 305-d PY represents 305-d Protein Yield, and (<b>c</b>) 305-d FY represents 305-d Fat Yield. MAS represents mastitis, MET represents endometritis, and KET represents ketosis. LSM represents the least squares means. NS and “**”, respectively, represent no and highly significant differences (<span class="html-italic">p &gt;</span> 0.05, <span class="html-italic">p</span> &lt; 0.01).</p>
Full article ">Figure 5
<p>The effects of mastitis, endometritis, and ketosis on milk composition in dairy cows. (<b>a</b>) ASCS represents the Average Somatic Cell Score (SCS) in lactation, (<b>b</b>) ADSCC represents the Average Differential Somatic Cell Count (DSCC) in lactation, (<b>c</b>) APROPER represents the Average Protein Percentage in lactation, and (<b>d</b>) AFATPER represents the Average Fat Percentage in Lactation. MAS represents mastitis, MET represents endometritis, and KET represents ketosis. LSM represents the least squares means. NS and “**”, respectively, represent no and highly significant differences (<span class="html-italic">p &gt;</span> 0.05, <span class="html-italic">p</span> &lt; 0.01).</p>
Full article ">Figure 5 Cont.
<p>The effects of mastitis, endometritis, and ketosis on milk composition in dairy cows. (<b>a</b>) ASCS represents the Average Somatic Cell Score (SCS) in lactation, (<b>b</b>) ADSCC represents the Average Differential Somatic Cell Count (DSCC) in lactation, (<b>c</b>) APROPER represents the Average Protein Percentage in lactation, and (<b>d</b>) AFATPER represents the Average Fat Percentage in Lactation. MAS represents mastitis, MET represents endometritis, and KET represents ketosis. LSM represents the least squares means. NS and “**”, respectively, represent no and highly significant differences (<span class="html-italic">p &gt;</span> 0.05, <span class="html-italic">p</span> &lt; 0.01).</p>
Full article ">Figure 6
<p>The effects of mastitis, endometritis, and ketosis on reproductive performance in dairy cows. (<b>a</b>) CTFS represents the Interval from Calving to First Service, (<b>b</b>) FSTC represents the Interval from First Service to Conception, (<b>c</b>) NUMS represents the Number of Services, and (<b>d</b>) CI represents Calving Interval. MAS represents mastitis, MET represents endometritis, and KET represents ketosis. LSM represents the least squares means. NS and “**”, respectively, represent no and highly significant differences (<span class="html-italic">p &gt;</span> 0.05, <span class="html-italic">p</span> &lt; 0.01).</p>
Full article ">Figure 7
<p>The heritability and genetic correlation of mastitis, endometritis, ketosis, milk yield, milk composition, and reproduction traits in dairy cattle. MAS represents mastitis, MET represents endometritis, KET represents ketosis, 305-d MY represents 305-d Milk Yield, 305-d PY represents 305-d Protein Yield, 305-d FY represents 305-d Fat Yield, ASCS represents the Average Somatic Cell Score (SCS) in lactation, ADSCC represents the Average Differential Somatic Cell Count (DSCC) in lactation, APROPER represents the Average Protein Percentage in lactation, AFATPER represents the Average Fat Percentage in lactation, CTFS represents the Interval from Calving to First Service, FSTC represents the Interval from First Service to Conception, NUMS represents the Number of Services, and CI represents the Calving Interval. Corrg, in the lower triangles, represents genetic correlation. H, in the diagonal, represents heritability. NS, “*”, and “**”, respectively, represent no, significant, and highly significant differences (<span class="html-italic">p &gt;</span> 0.05, <span class="html-italic">p &lt;</span> 0.05, <span class="html-italic">p</span> &lt; 0.01).</p>
Full article ">Figure 8
<p>The phenotypic correlation of mastitis, endometritis, ketosis, milk yield, milk composition, and reproduction traits in dairy cattle. Corrp represents phenotypic correlation. NS, “*”, and “**”, respectively, represent no, significant, and highly significant differences (<span class="html-italic">p &gt;</span> 0.05, <span class="html-italic">p &lt;</span> 0.05, <span class="html-italic">p</span> &lt; 0.01). The trait meanings are the same as in <a href="#animals-14-02372-f007" class="html-fig">Figure 7</a>.</p>
Full article ">
24 pages, 23818 KiB  
Article
Effects of Assimilating Ground-Based Microwave Radiometer and FY-3D MWTS-2/MWHS-2 Data in Precipitation Forecasting
by Bingli Wang, Wei Cheng, Yansong Bao, Shudong Wang, George P. Petropoulos, Shuiyong Fan, Jiajia Mao, Ziqi Jin and Zihui Yang
Remote Sens. 2024, 16(14), 2682; https://doi.org/10.3390/rs16142682 - 22 Jul 2024
Viewed by 685
Abstract
This study investigates the impacts of the joint assimilation of ground-based microwave radiometer (MWR) and FY-3D microwave sounder (MWTS-2/MWHS-2) observations on the analyses and forecasts for precipitation forecast. Based on the weather research and forecasting data assimilation (WRFDA) system, four experiments are conducted [...] Read more.
This study investigates the impacts of the joint assimilation of ground-based microwave radiometer (MWR) and FY-3D microwave sounder (MWTS-2/MWHS-2) observations on the analyses and forecasts for precipitation forecast. Based on the weather research and forecasting data assimilation (WRFDA) system, four experiments are conducted in this study, concerning a heavy precipitation event in Beijing on 2 July 2021, and 10-day batch experiments were also conducted. The key study findings include the following: (1) Both ground-based microwave radiometer and MWTS-2/MWHS-2 data contribute to improvements in the initial fields of the model, leading to appropriate adjustments in the thermal structure of the model. (2) The forecast fields of the experiments assimilating ground-based microwave radiometer and MWTS-2/MWHS-2 data show temperature and humidity performances closer to the true fields compared with the control experiment. (3) Separate assimilation of two types of microwave radiometer data can improve precipitation forecasts, while joint assimilation provides the most accurate forecasts among all the experiments. In the single-case, compared with the control experiment, the individual and combined assimilation of MWR and MWTS-2/MWHS-2 improves the six-hour cumulative precipitation threat score (TS) at the 25 mm level by 57.1%, 28.9%, and 38.2%, respectively. The combined assimilation also improves the scores at the 50 mm level by 54.4%, whereas individual assimilations show a decrease in performance. In the batch experiments, the MWR_FY experiment’s TS of 24 h precipitation forecast improves 28.5% at 10 mm and 330% at 25 mm based on the CTRL. Full article
Show Figures

Figure 1

Figure 1
<p>Distribution of ground microwave radiometer stations, located in Beijing, China: Haidian (Station ID 54399), Yanqing (Station ID 54406), Huairou (Station ID 54419), Pinggu (Station ID 54424), Nanjiao (Station ID 54511), and Xiayunling (Station ID 54597).</p>
Full article ">Figure 2
<p>Probability distribution function of deviations in ground-based microwave radiometer temperature data before and after preprocessing for stations (<b>a</b>) 54511, (<b>b</b>) 54597, (<b>c</b>) 54419, (<b>d</b>) 54424, (<b>e</b>) 54511, and (<b>f</b>) 54597. The red “u” and “o” represent the error mean and error standard deviation of all samples before QC and BC. The blue “u” and “o” represent the error mean and error standard deviation of all samples after QC and BC.</p>
Full article ">Figure 3
<p>Probability distribution function of observation-minus-background (OMB) for (<b>a</b>,<b>d</b>) without bias correction, (<b>b</b>,<b>e</b>) OMB with bias correction, and (<b>c</b>,<b>f</b>) OMA with bias correction of (<b>a</b>–<b>c</b>) MWTS-2 ch7, (<b>d</b>–<b>e</b>) MWHS-2 ch11, respectively. “ME” represents the mean, “STD” represents the standard deviation, and “RMSE” represents the root-mean-square error.</p>
Full article ">Figure 4
<p>Cloud distribution after cloud matching (<b>a</b>,<b>b</b>) and true cloud distribution from the moderate resolution spectral imager (MERSI) (<b>c</b>,<b>d</b>) at 1800 UTC on 1 July (<b>a</b>,<b>c</b>) and 0600 UTC on 2 July (<b>b</b>,<b>d</b>).</p>
Full article ">Figure 5
<p>Brightness temperature distribution in MWTS-2 channel 11 under clear-sky conditions (<b>a</b>,<b>c</b>) and retrieval temperature distribution at 500 hPa latitude (<b>b</b>,<b>d</b>) at 1800 UTC on 1 July (<b>a</b>,<b>b</b>) and 0600 UTC on 2 July (<b>c</b>,<b>d</b>).</p>
Full article ">Figure 6
<p>Distribution of cumulative precipitation in Beijing on 2 July 2021 (<b>a</b>) from 1800 UTC to 0000 UTC; (<b>b</b>) from 1800 UTC to 2100 UTC; (<b>c</b>) from 2100 UTC to 0000 UTC.</p>
Full article ">Figure 7
<p>Simulated area in ARW-WRF.</p>
Full article ">Figure 8
<p>Assimilation flowchart.</p>
Full article ">Figure 9
<p>MWR, retrieval, and clear-sky observations at (<b>a</b>) 850 hPa and (<b>b</b>) 500 hPa in Beijing at 0600 UTC on 2 July 2021.</p>
Full article ">Figure 10
<p>Temperature and humidity profiles for the MWR, along with its adjacent MWTS-2/WMHS-2 cloud area retrievals from sites 54399, 54419, and 54597 at 0600 UTC on 2 July 2021.</p>
Full article ">Figure 11
<p>Vertical cross-sections of the specific humidity increments after the first data assimilation along 40°N from (<b>a</b>) CTRL, (<b>b</b>) MWR, (<b>c</b>) FY, and (<b>d</b>) MWR_FY.</p>
Full article ">Figure 12
<p>Specific humidity increments after the first data assimilation at 500 hPa and 850 hPa altitude from (<b>a</b>) CTRL, (<b>b</b>) MWR, (<b>c</b>) FY, and (<b>d</b>) MWR_FY, with the black mark indicating the location of the ground-based microwave radiometer station.</p>
Full article ">Figure 13
<p>Temperature and wind superposition at 500 hPa forecast field from (<b>a</b>) GFS, (<b>b</b>) CTRL, (<b>c</b>) MWR, (<b>d</b>) FY, and (<b>e</b>) MWR_FY at 0000 UTC on 3 July 2021. The shadow represents the temperature distribution and the dashed line represents the temperature isotherm.</p>
Full article ">Figure 14
<p>Vertical cross-sections of relative humidity (shadow), temperature (red lines), and wind field (vectors) along 40°N from (<b>a</b>) GFS, (<b>b</b>) CTRL, (<b>c</b>) MWR, (<b>d</b>) FY, and (<b>e</b>) MWR_FY at 0000 UTC on 3 July 2021.</p>
Full article ">Figure 15
<p>Distribution of precipitation forecasts for the 6 h (from 1800 UTC on 2 July to 0000 UTC on 3 July), the first 3 h (from 1800 UTC to 2100 UTC on 2 July), and the last 3 h (from 2100 UTC on 2 July to 0000 UTC on 3 July) from OBS (actual observation), the CTRL, the MWR, the FY, and the MWR_FY experiments.</p>
Full article ">Figure 16
<p>Threat score (TS) for (<b>a</b>) 6 h rainfall forecast, (<b>b</b>) the first 3 h rainfall forecast, and (<b>c</b>) the last 3 h rainfall forecast from all four experiments at different rainfall thresholds.</p>
Full article ">Figure 17
<p>RMSE/bias profiles for 12 h and 24 h temperature (<b>a</b>–<b>d</b>) and relative humidity (<b>e</b>–<b>h</b>) forecast fields.</p>
Full article ">Figure 18
<p>Threat score (TS) for (<b>a</b>) 24 h rainfall forecast, (<b>b</b>) the first 12 h rainfall forecast, and (<b>c</b>) the last 12 h rainfall forecast from all four experiments of 10-days at different rainfall thresholds.</p>
Full article ">
25 pages, 11697 KiB  
Article
Improving Typhoon Muifa (2022) Forecasts with FY-3D and FY-3E MWHS-2 Satellite Data Assimilation under Clear Sky Conditions
by Feifei Shen, Xiaolin Yuan, Hong Li, Dongmei Xu, Jingyao Luo, Aiqing Shu and Lizhen Huang
Remote Sens. 2024, 16(14), 2614; https://doi.org/10.3390/rs16142614 - 17 Jul 2024
Viewed by 851
Abstract
This study investigates the impacts of assimilating the Microwave Humidity Sounder II (MWHS-2) radiance data carried on the FY-3D and FY-3E satellites on the analyses and forecasts of Typhoon Muifa in 2022 under clear-sky conditions. Data assimilation experiments are conducted using the Weather [...] Read more.
This study investigates the impacts of assimilating the Microwave Humidity Sounder II (MWHS-2) radiance data carried on the FY-3D and FY-3E satellites on the analyses and forecasts of Typhoon Muifa in 2022 under clear-sky conditions. Data assimilation experiments are conducted using the Weather Research and Forecasting (WRF) model coupled with the Three-Dimensional Variational (3D-Var) Data Assimilation method to compare the different behaviors of FY-3D and FY-3E radiances. Additionally, the data assimilation strategies are assessed in terms of the sequence of applying the conventional and MWHS-2 radiance data. The results show that assimilating MWHS-2 data is able to enhance the dynamic and thermal structures of the typhoon system. The experiment with FY-3E MWHS-2 assimilated demonstrated superior performance in terms of simulating the typhoon’s structure and providing a prediction of the typhoon’s intensity and track than the experiment with FY-3D MWHS-2 did. The two-step assimilation strategy that assimilates conventional observations before the radiance data has improved the track and intensity forecasts at certain times, particularly with the FY-3E MWHS-2 radiance. It appears that large-scale atmospheric conditions are more refined by initially assimilating the Global Telecommunication System (GTS) data, with subsequent satellite data assimilation further adjusting the model state. This strategy has also confirmed improvements in precipitation prediction as it enhances the dynamic and thermal structures of the typhoon system. Full article
Show Figures

Figure 1

Figure 1
<p>The evolution of Typhoon Muifa’s intensity levels throughout its track, recorded from 0600 on 6 September to 2100 on 16 September UTC. Time (day and hour) and central pressure of typhoon are noted at nodes where the typhoon intensity shifts.</p>
Full article ">Figure 2
<p>Terrain height (filled colors, unit: m) and the best track of Typhoon Muifa 24 h interval (black dot), from 0000 UTC on 8 September to 0000 UTC 16 September within the model domain. Blue points represent FY-3D MWHS-2 data, while red points represent FY-3E MWHS-2 data.</p>
Full article ">Figure 3
<p>The flow chart for the data assimilation experiments. One-step experiments including GTS_DA (yellow solid line), 3D_DA (blue solid line), and 3E_DA (red solid line) are depicted in (<b>a</b>), two-step experiments including 3D_R_DA (blue dashed line) and 3E_R_DA (red dashed line) are depicted in (<b>b</b>).</p>
Full article ">Figure 4
<p>OMB of (<b>a</b>) 3D_DA, (<b>c</b>) 3D_R_DA, (<b>e</b>) 3E_DA, and (<b>g</b>) 3E_R_DA; OMA of (<b>b</b>) 3D_DA, (<b>d</b>) 3D_R_DA, (<b>f</b>) 3E_DA, and (<b>h</b>) 3E_R_DA for the brightness temperature (units: K) of channel 11. The blue dot and red dot represent the location of the typhoon center at 0600 UTC on 14 September 2022 and 0900 UTC on 14 September 2022, respectively.</p>
Full article ">Figure 5
<p>Scatter distribution of the brightness temperature (unit: K) of channel 11 simulated from (<b>a</b>,<b>d</b>,<b>g</b>,<b>j</b>) the background before the bias correction, (<b>b</b>,<b>e</b>,<b>h</b>,<b>k</b>) the background after bias correction, and (<b>c</b>,<b>f</b>,<b>i</b>,<b>l</b>) the analysis after the bias correction in the y-axis versus the observed radiances in the x-axis. (<b>a</b>–<b>c</b>) 3D_DA experiment, (<b>d</b>–<b>f</b>) 3D_R_DA experiment, (<b>g</b>–<b>i</b>) 3E_DA experiment, and (<b>j</b>–<b>l</b>) 3E_R_DA experiment.</p>
Full article ">Figure 6
<p>(<b>a</b>) The number of assimilated satellite data for the experiment; (<b>b</b>,<b>c</b>) Mean for the experiment and (<b>d</b>,<b>e</b>) stdv for the experiment.</p>
Full article ">Figure 7
<p>The 500 hPa water vapor mixing ratio differences (analysis minus background, shading, unit: g/kg) and wind speed (vector, unit: m/s) for (<b>a</b>) GTS_DA, (<b>b</b>) 3D_DA, and (<b>c</b>) 3E_DA. The differences in the (<b>d</b>) first step, (<b>e</b>) second step, and (<b>f</b>) two-steps in 3D_R_DA at 0600 UTC 14 September. The differences in the (<b>g</b>) first step, (<b>h</b>) second step, and (<b>i</b>) two-steps in 3E_R_DA at 0900 UTC 14 September. The red symbol represents the position of the typhoon sourced from the CMA at 0600 UTC 14 September (<b>a</b>,<b>b</b>,<b>d</b>–<b>f</b>) and 0900 UTC 14 September (<b>c</b>,<b>g</b>–<b>i</b>).</p>
Full article ">Figure 8
<p>Same as <a href="#remotesensing-16-02614-f007" class="html-fig">Figure 7</a>, but for the temperature (unit: K).</p>
Full article ">Figure 9
<p>Same as <a href="#remotesensing-16-02614-f007" class="html-fig">Figure 7</a>, but for geopotential height (contours, units: m<sup>2</sup> s<sup>−2</sup>; contours every 2 m<sup>2</sup> s<sup>−2</sup>) and geopotential height differences (analysis minus background, shading, unit: m).</p>
Full article ">Figure 10
<p>RMSE vertical profiles of (<b>a</b>) u-wind (m/s), (<b>b</b>) v-wind (m/s), (<b>c</b>) temperature (K), and (<b>d</b>) specific humidity (g/kg) forecasts versus sounding and surface synoptic (SYNOP) observations at 1200 UTC on 15 September 2022. Error dots indicate statistical significance at the 95% confidence level.</p>
Full article ">Figure 11
<p>(<b>a</b>) Typhoon Muifa’s track from 0900 UTC on 14 September to 0000 UTC on 16 September 2022 and (<b>b</b>) the corresponding track error, the numbers on the x-axis represent day and hour, respectively.</p>
Full article ">Figure 12
<p>The 39 h forecast errors of (<b>a</b>) MSLP (unit: hPa) and (<b>b</b>) MSWS (unit: m/s), initialized at 0900 UTC on 14 September 2022.</p>
Full article ">Figure 13
<p>The vertical cross section of the wind speed (shading, unit: m/s) and potential temperature (contours, unit: K, intervals of 4 K) for (<b>a</b>) CTRL, (<b>b</b>) GTS_DA, (<b>c</b>) 3D_DA, and (<b>d</b>) 3D_R_DA at 0600 UTC 14 September, (<b>e</b>) 3E_DA, and (<b>f</b>) 3E_R_DA at 0900 UTC 14 September. (<b>g</b>) is the cross section with the wind speed (blue wind barbs, unit: m/s) and potential temperature (shading, unit: K). The red symbol is the position of the typhoon at 0600 UTC 14 September.</p>
Full article ">Figure 14
<p>The 24 h accumulated precipitation distribution (unit: mm) from 0900 UTC on 14 September to 0900 UTC on 15 September in the experiments (<b>a</b>) CTRL, (<b>b</b>) GTS_DA, (<b>c</b>) 3D_DA, (<b>d</b>) 3D_R_DA, (<b>e</b>) 3E_DA, (<b>f</b>) 3E_R_DA, and in (<b>g</b>) the observation.</p>
Full article ">Figure 15
<p>ETS of 24 h accumulated precipitation at different thresholds.</p>
Full article ">
14 pages, 2215 KiB  
Article
DNA Methylation of the Autonomous Pathway Is Associated with Flowering Time Variations in Arabidopsis thaliana
by Hongjie Xie, Xinchen Li, Yuli Sun, Lei Lin, Keke Xu, Huan Lu, Biao Cheng, Siming Xue, Dan Cheng and Sheng Qiang
Int. J. Mol. Sci. 2024, 25(13), 7478; https://doi.org/10.3390/ijms25137478 - 8 Jul 2024
Viewed by 1024
Abstract
Plant flowering time is affected by endogenous and exogenous factors, but its variation patterns among different populations of a species has not been fully established. In this study, 27 Arabidopsis thaliana accessions were used to investigate the relationship between autonomous pathway gene methylation, [...] Read more.
Plant flowering time is affected by endogenous and exogenous factors, but its variation patterns among different populations of a species has not been fully established. In this study, 27 Arabidopsis thaliana accessions were used to investigate the relationship between autonomous pathway gene methylation, gene expression and flowering time variation. DNA methylation analysis, RT-qPCR and transgenic verification showed that variation in the flowering time among the Arabidopsis populations ranged from 19 to 55 days and was significantly correlated with methylation of the coding regions of six upstream genes in the autonomous pathway, FLOWERING LOCUS VE (FVE), FLOWERING LOCUS Y (FY), FLOWERING LOCUS D (FLD), PEPPER (PEP), HISTONE DEACETYLASE 5 (HAD5) and Pre-mRNA Processing Protein 39-1 (PRP39-1), as well as their relative expression levels. The expression of FVE and FVE(CS) was modified separately through degenerate codon substitution of cytosine and led to earlier flowering of transgenic plants by 8 days and 25 days, respectively. An accurate determination of methylated sites in FVE and FVE(CS) among those transgenic plants and the recipient Col-0 verified the close relationship between the number of methylation sites, expression and flowering time. Our findings suggest that the methylation variation of these six key upstream transcription factors was associated with the gene expression level of the autonomous pathway and flowering time in Arabidopsis. The FVE(CS) and FVE genes in transgenic plants tended to be hypermethylated, which could be a protective mechanism for plants. However, modification of gene sequences through degenerate codon substitution to reduce cytosine can avoid hypermethylated transferred genes in transgenic plants. It may be possible to partially regulate the flowering of plants by modified trans-epigenetic technology. Full article
(This article belongs to the Special Issue Molecular and Structural Research Advances in Model Plants)
Show Figures

Figure 1

Figure 1
<p>Methylation levels of gene coding regions in <span class="html-italic">A. thaliana</span> autonomous pathway. Ellipses: genes with methylation related to flowering time and leaf number; Rectangles: gene coding regions universally methylated in the 27 <span class="html-italic">A. thaliana</span> accessions but not related to flowering time or leaf number; and Hexagons: gene coding regions generally not methylated in 27 <span class="html-italic">A. thaliana</span> accessions and not related to flowering time or leaf number.</p>
Full article ">Figure 2
<p>Correlation plot of methylation levels of the coding regions of <span class="html-italic">A. thaliana</span> autonomous pathway genes, flowering time and leaf number. mC: total number of methylated sites of coding regions.</p>
Full article ">Figure 3
<p>Relative expression levels of the <span class="html-italic">FVE</span> (<b>A</b>), <span class="html-italic">FY</span> (<b>B</b>), <span class="html-italic">FLD</span> (<b>C</b>), <span class="html-italic">PEP</span> (<b>D</b>), <span class="html-italic">HDA5</span> (<b>E</b>) and <span class="html-italic">PRP39-1</span> (<b>F</b>) genes and correlations with coding region methylation levels among Col-0, Br-0 and Tscha-1 (<b>G</b>–<b>L</b>). The different letters indicate significant differences among the different accessions, <span class="html-italic">p</span> &lt; 0.05.</p>
Full article ">Figure 4
<p>Expression analysis of <span class="html-italic">FVE</span>(<span class="html-italic">CS</span>) and <span class="html-italic">FVE</span> in <span class="html-italic">fve-3</span> and their flowering phenotype. (<b>A</b>), Construction of <span class="html-italic">FVE</span> and <span class="html-italic">FVE</span>(<span class="html-italic">CS</span>) plant expression vectors; (<b>B</b>), Phenotypes of six-week-old seedlings of Col-0, <span class="html-italic">fve-3</span>+p<span class="html-italic">FVE</span>::<span class="html-italic">FVE</span>(<span class="html-italic">CS</span>), <span class="html-italic">fve-3</span>+p<span class="html-italic">FVE</span>::<span class="html-italic">FVE</span> and <span class="html-italic">fve-3</span>; (<b>C</b>), Flowering time differences among these plants; and (<b>D</b>), Leaf number differences among these plants. The different letters indicate significant differences among the different plants, <span class="html-italic">p</span> &lt; 0.05.</p>
Full article ">Figure 5
<p>RT-qPCR analysis of the autonomous pathway genes <span class="html-italic">FVE</span>, <span class="html-italic">FLC</span>, <span class="html-italic">SOC1</span>, <span class="html-italic">LFY</span> and <span class="html-italic">AP1</span> of Col-0, <span class="html-italic">fve-3</span>+p<span class="html-italic">FVE</span>::<span class="html-italic">FVE</span>(<span class="html-italic">CS</span>), <span class="html-italic">fve-3</span>+p<span class="html-italic">FVE</span>::<span class="html-italic">FVE</span> and <span class="html-italic">fve-3</span>. The letters indicate significant differences among the different transgenic plants, <span class="html-italic">p</span> &lt; 0.05.</p>
Full article ">Figure 6
<p>Methylated sites of <span class="html-italic">FVE</span> and <span class="html-italic">FVE</span>(<span class="html-italic">CS</span>) genes in Col-0, <span class="html-italic">fve-3</span>+p<span class="html-italic">FVE</span>::<span class="html-italic">FVE</span> and <span class="html-italic">fve-3</span>+p<span class="html-italic">FVE</span>::<span class="html-italic">FVE</span>(<span class="html-italic">CS</span>) plants.</p>
Full article ">
20 pages, 7017 KiB  
Article
Inter-Comparison of SST Products from iQuam, AMSR2/GCOM-W1, and MWRI/FY-3D
by Yili Zhao, Ping Liu and Wu Zhou
Remote Sens. 2024, 16(11), 2034; https://doi.org/10.3390/rs16112034 - 6 Jun 2024
Cited by 1 | Viewed by 1030
Abstract
Evaluating sea surface temperature (SST) products is essential before their application in marine environmental monitoring and related studies. SSTs from the in situ SST Quality Monitor (iQuam) system, Advanced Microwave Scanning Radiometer 2 (AMSR2) aboard the Global Change Observation Mission 1st-Water, and the [...] Read more.
Evaluating sea surface temperature (SST) products is essential before their application in marine environmental monitoring and related studies. SSTs from the in situ SST Quality Monitor (iQuam) system, Advanced Microwave Scanning Radiometer 2 (AMSR2) aboard the Global Change Observation Mission 1st-Water, and the Microwave Radiation Imager (MWRI) aboard the Chinese Fengyun-3D satellite are intercompared utilizing extended triple collocation (ETC) and direct comparison methods. Additionally, error characteristic variations with respect to time, latitude, SST, sea surface wind speed, columnar water vapor, and columnar cloud liquid water are analyzed comprehensively. In contrast to the prevailing focus on SST validation accuracy, the random errors and the capability to detect SST variations are also evaluated in this study. The result of ETC analysis indicates that iQuam SST from ships exhibits the highest random error, above 0.83 °C, whereas tropical mooring SST displays the lowest random error, below 0.28 °C. SST measurements from drifters, tropical moorings, Argo floats, and high-resolution drifters, which possess random errors of less than 0.35 °C, are recommended for validating remotely sensed SST. The ability of iQuam, AMSR2, and MWRI to detect SST variations diminishes significantly in ocean areas between 0°N and 20°N latitude and latitudes greater than 50°N and 50°S. AMSR2 and iQuam demonstrate similar random errors and capabilities for detecting SST variations, whereas MWRI shows a high random error and weak capability. In comparison to iQuam SST, AMSR2 exhibits a root-mean-square error (RMSE) of about 0.51 °C with a bias of −0.05 °C, while MWRI shows an RMSE of about 1.26 °C with a bias of −0.14 °C. Full article
Show Figures

Figure 1

Figure 1
<p>Spatial distribution of triple collocations.</p>
Full article ">Figure 2
<p>Comparison between AMSR2 SST and iQuam SST in the daytime and nighttime. (<b>a</b>) Daytime. (<b>b</b>) Nighttime.</p>
Full article ">Figure 3
<p>Comparison between MWRI SST and iQuam SST in the daytime and nighttime. (<b>a</b>) Daytime. (<b>b</b>) Nighttime.</p>
Full article ">Figure 4
<p>Spatial distribution of AMSR2 SST bias referring to iQuam SST. (<b>a</b>) Daytime. (<b>b</b>) Nighttime.</p>
Full article ">Figure 5
<p>Spatial distribution of SST difference between MWRI and iQuam. (<b>a</b>) Daytime. (<b>b</b>) Nighttime.</p>
Full article ">Figure 6
<p>Temporal variation in error characteristics. (<b>a</b>) ESD. (<b>b</b>) <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>SNR</mi> </mrow> <mrow> <mi>SUb</mi> </mrow> </msub> </mrow> </semantics></math>. (<b>c</b>) Bias. (<b>d</b>) RMSE.</p>
Full article ">Figure 7
<p>Latitudinal variation in error characteristics. (<b>a</b>) ESD. (<b>b</b>) <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>SNR</mi> </mrow> <mrow> <mi>SUb</mi> </mrow> </msub> </mrow> </semantics></math>. (<b>c</b>) Bias. (<b>d</b>) RMSE.</p>
Full article ">Figure 8
<p>Variation in error characteristics along SST. (<b>a</b>) ESD. (<b>b</b>) <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>SNR</mi> </mrow> <mrow> <mi>SUb</mi> </mrow> </msub> </mrow> </semantics></math>. (<b>c</b>) Bias. (<b>d</b>) RMSE.</p>
Full article ">Figure 9
<p>Variation in error characteristics along sea surface wind speed. (<b>a</b>) ESD. (<b>b</b>) <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>SNR</mi> </mrow> <mrow> <mi>SUb</mi> </mrow> </msub> </mrow> </semantics></math>. (<b>c</b>) Bias. (<b>d</b>) RMSE.</p>
Full article ">Figure 10
<p>Variation in error characteristics along columnar water vapor. (<b>a</b>) ESD. (<b>b</b>) <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>SNR</mi> </mrow> <mrow> <mi>SUb</mi> </mrow> </msub> </mrow> </semantics></math>. (<b>c</b>) Bias. (<b>d</b>) RMSE.</p>
Full article ">Figure 11
<p>Variation in error characteristics along columnar cloud liquid water (<b>a</b>) ESD. (<b>b</b>) <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>SNR</mi> </mrow> <mrow> <mi>SUb</mi> </mrow> </msub> </mrow> </semantics></math>. (<b>c</b>) Bias. (<b>d</b>) RMSE.</p>
Full article ">
19 pages, 7359 KiB  
Article
Evolutionary and Phylogenetic Dynamics of SARS-CoV-2 Variants: A Genetic Comparative Study of Taiyuan and Wuhan Cities of China
by Behzad Hussain and Changxin Wu
Viruses 2024, 16(6), 907; https://doi.org/10.3390/v16060907 - 3 Jun 2024
Viewed by 974
Abstract
Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a positive-sense, single-stranded RNA genome-containing virus which has infected millions of people all over the world. The virus has been mutating rapidly enough, resulting in the emergence of new variants and sub-variants which have reportedly [...] Read more.
Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a positive-sense, single-stranded RNA genome-containing virus which has infected millions of people all over the world. The virus has been mutating rapidly enough, resulting in the emergence of new variants and sub-variants which have reportedly been spread from Wuhan city in China, the epicenter of the virus, to the rest of China and all over the world. The occurrence of mutations in the viral genome, especially in the viral spike protein region, has resulted in the evolution of multiple variants and sub-variants which gives the virus the benefit of host immune evasion and thus renders modern-day vaccines and therapeutics ineffective. Therefore, there is a continuous need to study the genetic characteristics and evolutionary dynamics of the SARS-CoV-2 variants. Hence, in this study, a total of 832 complete genomes of SARS-CoV-2 variants from the cities of Taiyuan and Wuhan in China was genetically characterized and their phylogenetic and evolutionary dynamics studied using phylogenetics, genetic similarity, and phylogenetic network analyses. This study shows that the four most prevalent lineages in Taiyuan and Wuhan are as follows: the Omicron lineages EG.5.1.1, followed by HK.3, FY.3, and XBB.1.16 (Pangolin classification), and clades 23F (EG.5.1), followed by 23H (HK.3), 22F (XBB), and 23D (XBB.1.9) (Nextclade classification), and lineage B followed by the Omicron FY.3, lineage A, and Omicron FL.2.3 (Pangolin classification), and the clades 19A, followed by 22F (XBB), 23F (EG.5.1), and 23H (HK.3) (Nextclade classification), respectively. Furthermore, our genetic similarity analysis show that the SARS-CoV-2 clade 19A-B.4 from Wuhan (name starting with 412981) has the least genetic similarity of about 95.5% in the spike region of the genome as compared to the query sequence of Omicron XBB.2.3.2 from Taiyuan (name starting with 18495234), followed by the Omicron FR.1.4 from Taiyuan (name starting with 18495199) with ~97.2% similarity and Omicron DY.3 (name starting with 17485740) with ~97.9% similarity. The rest of the variants showed ≥98% similarity with the query sequence of Omicron XBB.2.3.2 from Taiyuan (name starting with 18495234). In addition, our recombination analysis results show that the SARS-CoV-2 variants have three statistically significant recombinant events which could have possibly resulted in the emergence of Omicron XBB.1.16 (recombination event 3), FY.3 (recombination event 5), and FL.2.4 (recombination event 7), suggesting some very important information regarding viral evolution. Also, our phylogenetic tree and network analyses show that there are a total of 14 clusters and more than 10,000 mutations which may have probably resulted in the emergence of cluster-I, followed by 47 mutations resulting in the emergence of cluster-II and so on. The clustering of the viral variants of both cities reveals significant information regarding the phylodynamics of the virus among them. The results of our temporal phylogenetic analysis suggest that the variants of Taiyuan have likely emerged as independent variants separate from the variants of Wuhan. This study, to the best of our knowledge, is the first ever genetic comparative study between Taiyuan and Wuhan cities in China. This study will help us better understand the virus and cope with the emergence and spread of new variants at a local as well as an international level, and keep the public health authorities informed for them to make better decisions in designing new viral vaccines and therapeutics. It will also help the outbreak investigators to better examine any future outbreak. Full article
(This article belongs to the Section Coronaviruses)
Show Figures

Figure 1

Figure 1
<p>Genome annotation of SARS-CoV-2 Wuhan-Hu-1. The figure shows the circular form of the annotated genome of the SARS-CoV-2 original Wuhan strain (Wuhan-Hu-1; accession number: NC_045512). The complete genome is 29903 nucleotide bases in size. It shows different open reading frames (ORFs), including the spike region (from about 22 k to 24.5 k) which is one of the most important viral genes. The figure was created with the Unipro UGENE v48.0 software.</p>
Full article ">Figure 2
<p>Maximum Likelihood (ML) Phylogenetic tree. The figure shows the maximum likelihood (ML) phylogenetic tree constructed using the IQ-Tree2 multicore software with the best-fit substitution model GTR+F+I+G4. The tree is rooted on the mid-point and shows 14 clusters containing the different variants of SARS-CoV-2. Each cluster is given a different color randomly. Some sequences have been collapsed to better visualize the tree. The tree was edited and visualized by the iTOL webserver v6. Bootstrap values have been shown on each node.</p>
Full article ">Figure 3
<p>Time-based phylogenetic tree. The figure shows the time-based phylogenetic tree constructed using BEAST software with the substitution model GTR+F+I+G4. Some sequences have been collapsed to better visualize the tree. The tree was edited and visualized by iTOL webserver v6. The median height values are shown on each node and branch of the tree to indicate the time of distribution of each variant among the two cities.</p>
Full article ">Figure 4
<p>Overall temporal trend of the variants. The figure shows the overall trend of the SARS-CoV-2 variants based on time. The graph shows the designation date on the X-axis and clade name on the Y-axis. From left to right, the X-axis shows the number of years from 2020 to 2024, and each colored dot represents an individual genome sequence with a particular year. The dark blue and purple dots show the genome sequences of year 2020–2021 followed by the light blue and aqua dots (year 2021), light green dots (year 2022), yellow and orange dots (year 2023), red dots (later time in the year 2023) and the gray dots (year 2024). As the analysis was conducted based on the spike region, some sequences might have been misidentified because of the inherent errors of the classification system. The colors have been assigned randomly by the webserver.</p>
Full article ">Figure 5
<p>Genetic similarity plot of representative SARS-CoV-2 variants. The figure shows the genetic similarity plots of the selected variants of SARS-CoV-2 in comparison with the Omicron XBB.2.3.2 from Taiyuan (18495234) as query sequence. Genetic similarity (%) has been shown on the Y-axis and the nucleotide position has been shown on the X-axis. Clade 19A-B.4 from Wuhan (412981) has the lowest genetic similarity of about 95.5% between the nucleotide position 22000 and 24000, which is the spike protein’s nucleotide position. Most of the variants show a genetic similarity of 98% or more.</p>
Full article ">Figure 6
<p>Phylogenetic network of SARS-CoV-2 variants from Taiyuan and Wuhan. The figure shows the phylogenetic network of all the 832 complete genomes of SARS-CoV-2 from Taiyuan and Wuhan, plus the reference and bat-CoV genomes, constructed by the PopArt v1.7 software. The network was constructed using the TCS method and shows that the SARS-CoV-2 genomes probably originated from the bat-CoV with a result of more than 10,000 mutations. The network was divided into 14 clusters, consistent with the phylogenetic tree. The numbers indicate the number of mutations, and each color represents a different main variant. Some interesting variants have been labeled on the network.</p>
Full article ">
17 pages, 32228 KiB  
Article
Precipitation Characteristics at Different Developmental Stages of the Tibetan Plateau Vortex in July 2021 Based on GPM-DPR Data
by Bingyun Yang, Suling Ren, Xi Wang and Ning Niu
Remote Sens. 2024, 16(11), 1947; https://doi.org/10.3390/rs16111947 - 28 May 2024
Viewed by 794
Abstract
The Tibetan Plateau vortex (TPV), as an α-scale mesoscale weather system, often brings severe weather conditions like torrential rain and severe convective storms. Based on the detections from the Global Precipitation Measurement (GPM) Core Observatory’s Dual-frequency Precipitation Radar (DPR) and the FY-4A satellite’s [...] Read more.
The Tibetan Plateau vortex (TPV), as an α-scale mesoscale weather system, often brings severe weather conditions like torrential rain and severe convective storms. Based on the detections from the Global Precipitation Measurement (GPM) Core Observatory’s Dual-frequency Precipitation Radar (DPR) and the FY-4A satellite’s Advanced Geostationary Radiation Imager (AGRI), combined with ERA5 reanalysis data, the precipitation characteristics of a TPV moving eastward during 8–13 July 2021 at different developmental stages are explored in this study. It was clear that the near-surface precipitation rate of the TPV during the initial stage at the eastern Tibetan Plateau (TP) was below 1 mm·h−1, implying overall weak precipitation dominated by stratiform clouds. After moving out of the TP, the radar reflectivity factor (Ze), precipitation rate, and normalized intercept parameter (dBNw) significantly increased, while the proportion of convective clouds gradually rose. Following the TPV movement, the distribution range and vertical thickness of Ze, mass-weighted mean diameter (Dm), and dBNw tended to increase. The high-frequency region of Ze appeared at 15–20 dBZ, while Dm and dBNw occurred at around 1 mm and 33 mm−1·m−3, respectively. Near the melting layer, Ze was characterized by a significant increase due to the aggregation and melting of ice crystals. The precipitation rate of convective clouds was generally greater than that of stratiform clouds, whilst both of them increased during the movement of the TPV. Particularly, at 01:00 on 12 July, there was a significant increase in the precipitation rate and Dm of convective clouds, while dBNw noticeably decreased. These findings could provide valuable insights into the three-dimensional structure and microphysical characteristics of the precipitation during the movement of the TPV, contributing to a better understanding of cloud precipitation mechanisms. Full article
Show Figures

Figure 1

Figure 1
<p>Path of the TPV during 8–13 July 2021 overlaying topography. The blue dots indicated the locations of the TPV centers from 15:00 on 8 July to 03:00 on 13 July. The red dots represent the locations of the TPV centers from 00:00 on 9 July to 00:00 on 13 July.</p>
Full article ">Figure 2
<p>Daily gauge precipitation at (<b>a</b>) 00:00 on 8 July to 00:00 on 9 July, (<b>b</b>) 00:00 on 9 July to 00:00 on 10 July, (<b>c</b>) 00:00 on 10 July to 00:00 on 11 July and (<b>d</b>) 00:00 on 11 July to 12:00 on 12 July.</p>
Full article ">Figure 3
<p>The brightness temperature of FY-4A/AGRI 10.8 µm channel (shading, K) with 500 hPa geopotential height (contours, 10 gpm) of ERA5 data at (<b>a</b>) 04:00 on 9 July, (<b>b</b>) 17:00 on 9 July, (<b>c</b>) 02:00 on 11 July, and (<b>d</b>) 01:00 on 12 July. The paralleled black lines represent the boundaries of GPM/DPR scanning path. The black box areas represent the study regions.</p>
Full article ">Figure 4
<p>The total column vertically-integrated water vapor (shading, kg·m<sup>−2</sup>) of ERA5 data at (<b>a</b>) 04:00 on 9 July, (<b>b</b>) 17:00 on 9 July, (<b>c</b>) 02:00 on 11 July and (<b>d</b>) 01:00 on 12 July. The paralleled black lines represent the boundaries of GPM/DPR scanning path.</p>
Full article ">Figure 4 Cont.
<p>The total column vertically-integrated water vapor (shading, kg·m<sup>−2</sup>) of ERA5 data at (<b>a</b>) 04:00 on 9 July, (<b>b</b>) 17:00 on 9 July, (<b>c</b>) 02:00 on 11 July and (<b>d</b>) 01:00 on 12 July. The paralleled black lines represent the boundaries of GPM/DPR scanning path.</p>
Full article ">Figure 5
<p>Average (<b>a</b>) divergence(s<sup>−1</sup>) and (<b>b</b>) vertical velocity(Pa·s<sup>−1</sup>) profiles of ERA5 data in the study regions in <a href="#remotesensing-16-01947-f003" class="html-fig">Figure 3</a>. The green, blue, cyan, and red solid lines represent 04:00 on 9 July, 17:00 on 9 July, 02:00 on 11 July, and 01:00 on 12 July, respectively.</p>
Full article ">Figure 6
<p>Horizontal distributions of (<b>a</b>,<b>d</b>,<b>g</b>,<b>j</b>) the near-surface precipitation rate (mm·h<sup>−1</sup>), (<b>b</b>,<b>e</b>,<b>h</b>,<b>k</b>) precipitation type, and (<b>c</b>,<b>f</b>,<b>i</b>,<b>l</b>) storm top height (km). The gray colors represent the brightness temperature of FY-4A/AGRI 10.8 µm channel. The paralleled black lines represent the boundaries of GPM/DPR scanning path. (<b>a</b>–<b>c</b>) 04:00 on 9 July, (<b>d</b>–<b>f</b>) 17:00 on 9 July, (<b>g</b>–<b>i</b>) 02:00 on 11 July, and (<b>j</b>–<b>l</b>) 01:00 on 12 July.</p>
Full article ">Figure 6 Cont.
<p>Horizontal distributions of (<b>a</b>,<b>d</b>,<b>g</b>,<b>j</b>) the near-surface precipitation rate (mm·h<sup>−1</sup>), (<b>b</b>,<b>e</b>,<b>h</b>,<b>k</b>) precipitation type, and (<b>c</b>,<b>f</b>,<b>i</b>,<b>l</b>) storm top height (km). The gray colors represent the brightness temperature of FY-4A/AGRI 10.8 µm channel. The paralleled black lines represent the boundaries of GPM/DPR scanning path. (<b>a</b>–<b>c</b>) 04:00 on 9 July, (<b>d</b>–<b>f</b>) 17:00 on 9 July, (<b>g</b>–<b>i</b>) 02:00 on 11 July, and (<b>j</b>–<b>l</b>) 01:00 on 12 July.</p>
Full article ">Figure 7
<p>Vertical cross sections of (<b>a</b>,<b>d</b>,<b>g</b>,<b>j</b>) <span class="html-italic">Ze</span> (dBZ), (<b>b</b>,<b>e</b>,<b>h</b>,<b>k</b>) <span class="html-italic">Dm</span> (mm), and (<b>c</b>,<b>f</b>,<b>i</b>,<b>l</b>) dB<span class="html-italic">Nw</span> (mm<sup>−1</sup>·m<sup>−3</sup>) along the AB lines in <a href="#remotesensing-16-01947-f006" class="html-fig">Figure 6</a>. The black solid lines represent surface heights. The dotted lines represent melting layer heights. The gray solid lines represent bright band heights. (<b>a</b>–<b>c</b>) 04:00 on 9 July, (<b>d</b>–<b>f</b>) 17:00 on 9 July, (<b>g</b>–<b>i</b>) 02:00 on 11 July, and (<b>j</b>–<b>l</b>) 01:00 on 12 July.</p>
Full article ">Figure 7 Cont.
<p>Vertical cross sections of (<b>a</b>,<b>d</b>,<b>g</b>,<b>j</b>) <span class="html-italic">Ze</span> (dBZ), (<b>b</b>,<b>e</b>,<b>h</b>,<b>k</b>) <span class="html-italic">Dm</span> (mm), and (<b>c</b>,<b>f</b>,<b>i</b>,<b>l</b>) dB<span class="html-italic">Nw</span> (mm<sup>−1</sup>·m<sup>−3</sup>) along the AB lines in <a href="#remotesensing-16-01947-f006" class="html-fig">Figure 6</a>. The black solid lines represent surface heights. The dotted lines represent melting layer heights. The gray solid lines represent bright band heights. (<b>a</b>–<b>c</b>) 04:00 on 9 July, (<b>d</b>–<b>f</b>) 17:00 on 9 July, (<b>g</b>–<b>i</b>) 02:00 on 11 July, and (<b>j</b>–<b>l</b>) 01:00 on 12 July.</p>
Full article ">Figure 8
<p>CFAD distributions of (<b>a</b>,<b>d</b>,<b>g</b>,<b>j</b>) <span class="html-italic">Ze</span> (dBZ), (<b>b</b>,<b>e</b>,<b>h</b>,<b>k</b>) <span class="html-italic">Dm</span> (mm), and (<b>c</b>,<b>f</b>,<b>i</b>,<b>l</b>) dB<span class="html-italic">Nw</span> (mm<sup>−1</sup>·m<sup>−3</sup>) in the study regions. The solid lines represent average lines. The long-dashed lines represent median lines. The dotted lines represent melting layer heights. (<b>a</b>–<b>c</b>) 04:00 on 9 July, (<b>d</b>–<b>f</b>) 17:00 on 9 July, (<b>g</b>–<b>i</b>) 02:00 on 11 July, and (<b>j</b>–<b>l</b>) 01:00 on 12 July.</p>
Full article ">Figure 8 Cont.
<p>CFAD distributions of (<b>a</b>,<b>d</b>,<b>g</b>,<b>j</b>) <span class="html-italic">Ze</span> (dBZ), (<b>b</b>,<b>e</b>,<b>h</b>,<b>k</b>) <span class="html-italic">Dm</span> (mm), and (<b>c</b>,<b>f</b>,<b>i</b>,<b>l</b>) dB<span class="html-italic">Nw</span> (mm<sup>−1</sup>·m<sup>−3</sup>) in the study regions. The solid lines represent average lines. The long-dashed lines represent median lines. The dotted lines represent melting layer heights. (<b>a</b>–<b>c</b>) 04:00 on 9 July, (<b>d</b>–<b>f</b>) 17:00 on 9 July, (<b>g</b>–<b>i</b>) 02:00 on 11 July, and (<b>j</b>–<b>l</b>) 01:00 on 12 July.</p>
Full article ">Figure 9
<p>Vertical average profiles of (<b>a</b>,<b>d</b>,<b>g</b>,<b>j</b>) precipitation rate (mm·h<sup>−1</sup>), (<b>b</b>,<b>e</b>,<b>h</b>,<b>k</b>) <span class="html-italic">Dm</span> (mm), and (<b>c</b>,<b>f</b>,<b>i</b>,<b>l</b>) dB<span class="html-italic">Nw</span> (mm<sup>−1</sup>·m<sup>−3</sup>) in the study regions. The blue lines represent stratiform clouds, and the green lines represent convective clouds. (<b>a</b>–<b>c</b>) 04:00 on 9 July, (<b>d</b>–<b>f</b>) 17:00 on 9 July, (<b>g</b>–<b>i</b>) 02:00 on 11 July and (<b>j</b>–<b>l</b>) 01:00 on 12 July.</p>
Full article ">Figure 10
<p>Horizontal distributions of (<b>a</b>,<b>d</b>,<b>g</b>,<b>j</b>) GPM-DPR near-surface precipitation rate (mm·h<sup>−1</sup>), (<b>b</b>,<b>e</b>,<b>h</b>,<b>k</b>) hourly gauge precipitation rate (mm·h<sup>−1</sup>), and (<b>c</b>,<b>f</b>,<b>i</b>,<b>l</b>) ERA5 hourly precipitation rate (mm·h<sup>−1</sup>). The paralleled black lines represent the boundaries of GPM/DPR scanning path at 04:00 on 9 July. (<b>a</b>–<b>c</b>) 04:00 on 9 July, (<b>d</b>–<b>f</b>) 17:00 on 9 July, (<b>g</b>–<b>i</b>) 02:00 on 11 July, and (<b>j</b>–<b>l</b>) 01:00 on 12 July.</p>
Full article ">Figure 10 Cont.
<p>Horizontal distributions of (<b>a</b>,<b>d</b>,<b>g</b>,<b>j</b>) GPM-DPR near-surface precipitation rate (mm·h<sup>−1</sup>), (<b>b</b>,<b>e</b>,<b>h</b>,<b>k</b>) hourly gauge precipitation rate (mm·h<sup>−1</sup>), and (<b>c</b>,<b>f</b>,<b>i</b>,<b>l</b>) ERA5 hourly precipitation rate (mm·h<sup>−1</sup>). The paralleled black lines represent the boundaries of GPM/DPR scanning path at 04:00 on 9 July. (<b>a</b>–<b>c</b>) 04:00 on 9 July, (<b>d</b>–<b>f</b>) 17:00 on 9 July, (<b>g</b>–<b>i</b>) 02:00 on 11 July, and (<b>j</b>–<b>l</b>) 01:00 on 12 July.</p>
Full article ">
11 pages, 21400 KiB  
Communication
Arctic Winds Retrieved from FY-3D Microwave Humidity Sounder-II 183.31 GHz Brightness Temperature Using Atmospheric Motion Vector Method
by Bingxu Li, Xi Guo, Hao Liu, Donghao Han, Gang Li and Ji Wu
Remote Sens. 2024, 16(10), 1715; https://doi.org/10.3390/rs16101715 - 12 May 2024
Viewed by 1168
Abstract
In this study, we develop an Atmospheric Motion Vector (AMV)-based method for retrieving wind vectors using 183.31 GHz water-vapor absorption channels. The method involves tracking water-vapor features from image triplets and subsequently deriving wind fields from motion vectors. The height of the derived [...] Read more.
In this study, we develop an Atmospheric Motion Vector (AMV)-based method for retrieving wind vectors using 183.31 GHz water-vapor absorption channels. The method involves tracking water-vapor features from image triplets and subsequently deriving wind fields from motion vectors. The height of the derived wind for each channel is determined by calculating the weighing function peak using monthly averaged ERA5 reanalysis data. By utilizing Microwave Humidity Sounder-II (MWHS-II) brightness temperatures from the five channels centered around 183.31 GHz, wind vectors are retrieved within the Arctic region for the entire year of 2022. The retrieval quality is evaluated through comparative analysis with ERA5 reanalysis data and the Visible Infrared Imaging Radiometer Suite (VIIRS) wind product. The resultant vector root mean square errors (RMSEs) are approximately 4.5 m/s for the three lower-height channels and 5.5 m/s for the two upper-height channels. These findings demonstrate a wind retrieval performance comparable to the existing methods, highlighting its potential for augmenting wind availability at lower height levels. Full article
(This article belongs to the Special Issue Advancements in Microwave Radiometry for Atmospheric Remote Sensing)
Show Figures

Figure 1

Figure 1
<p>Sample overlapping areas from which MWHS-II AMV-based winds are extracted in image triplets.</p>
Full article ">Figure 2
<p>Illustration of AMV-based wind tracking process. (<b>a</b>) The target box selected following the process described in <a href="#sec2dot4-remotesensing-16-01715" class="html-sec">Section 2.4</a>. (<b>b</b>) The target box tracked in the next image following the process described in <a href="#sec2dot5-remotesensing-16-01715" class="html-sec">Section 2.5</a>.</p>
Full article ">Figure 3
<p>The weighting functions of the five water-vapor channels: (<b>a</b>) June 2022, (<b>b</b>) December 2022.</p>
Full article ">Figure 4
<p>Sample retrieved wind fields and brightness temperatures in the region 70°N to 90°N, 135°W to 150°E at 19:28:00 UTC, 3 July 2022.</p>
Full article ">Figure 5
<p>Scatterplots of the retrieved MWHS-II winds in U (<b>left column</b>) and V (<b>right column</b>) components compared with ERA5 data in the Arctic region for June 2022: (<b>a</b>,<b>b</b>) 183 ± 1 GHz (400 hPa), (<b>c</b>,<b>d</b>) 183 ± 1.8 GHz (450 hPa), (<b>e</b>,<b>f</b>) 183 ± 3 GHz (550 hPa), (<b>g</b>,<b>h</b>) 183 ± 4.5 GHz (600 hPa), (<b>i</b>,<b>j</b>) represents 183 ± 7 GHz (750 hPa). Regions with a more yellow hue represent areas of higher data density, while regions with a more blue hue indicate areas of lower data density.</p>
Full article ">Figure 6
<p>Spatial distribution and data count of (<b>a</b>) VIIRS and (<b>b</b>) MWHS-II AMVs in June 2022.</p>
Full article ">
13 pages, 678 KiB  
Article
CSN1S1, CSN3 and LPL: Three Validated Gene Polymorphisms Useful for More Sustainable Dairy Production in the Mediterranean River Buffalo
by Alfredo Pauciullo, Giustino Gaspa, Yi Zhang, Qingyou Liu and Gianfranco Cosenza
Animals 2024, 14(10), 1414; https://doi.org/10.3390/ani14101414 - 9 May 2024
Cited by 1 | Viewed by 1296
Abstract
The search for DNA polymorphisms useful for the genetic improvement of dairy farm animals has spanned more than 40 years, yielding relevant findings in cattle for milk traits, where the best combination of alleles for dairy processing has been found in casein genes [...] Read more.
The search for DNA polymorphisms useful for the genetic improvement of dairy farm animals has spanned more than 40 years, yielding relevant findings in cattle for milk traits, where the best combination of alleles for dairy processing has been found in casein genes and in DGAT1. Nowadays, similar results have not yet been reached in river buffaloes, despite the availability of advanced genomic technologies and accurate phenotype records. The aim of the present study was to investigate and validate the effect of four single nucleotide polymorphisms (SNP) in the CSN1S1, CSN3, SCD and LPL genes on seven milk traits in a larger buffalo population. These SNPs have previously been reported to be associated with, or affect, dairy traits in smaller populations often belonging to one farm. A total of 800 buffaloes were genotyped. The following traits were individually recorded, monthly, throughout each whole lactation period from 2010 to 2021: daily milk yield (dMY, kg), protein yield (dPY, kg) and fat yield (dFY, kg), fat and protein contents (dFP, % and dPP, %), somatic cell count (SCC, 103 cell/mL) and urea (mg/dL). A total of 15,742 individual milk test day records (2496 lactations) were available for 680 buffalo cows, with 3.6 ± 1.7 parities (from 1 to 13) and an average of 6.1 ± 1.2 test day records per lactation. Three out four SNPs in the CSN1S1, CSN3 and LPL genes were associated with at least one of analyzed traits. In particular, the CSN1S1 (AJ005430:c.578C>T) gave favorable associations with all yield traits (dMY, p = 0.022; dPY, p = 0.014; dFY, p = 0.029) and somatic cell score (SCS, p = 0.032). The CSN3 (HQ677596: c.536C>T) was positively associated with SCS (p = 0.005) and milk urea (p = 0.04). Favorable effects on daily milk yield (dMY, p = 0.028), fat (dFP, p = 0.027) and protein (dPP, p = 0.050) percentages were observed for the LPL. Conversely, the SCD did not show any association with milk traits. This is the first example of a confirmation study carried out in the Mediterranean river buffalo for genes of economic interest in the dairy field, and it represents a very important indication for the preselection of young bulls destined for breeding programs aimed at more sustainable dairy production. Full article
Show Figures

Figure 1

Figure 1
<p>Plot of LS means of FP and PP over Days in milk within genotypic classes for the polymorphisms in the <span class="html-italic">LPL</span> and <span class="html-italic">CSN3</span> genes.</p>
Full article ">
15 pages, 6871 KiB  
Article
FY-4A Measurement of Cloud-Seeding Effect and Validation of a Catalyst T&D Algorithm
by Liangrui Yan, Yuquan Zhou, Yixuan Wu, Miao Cai, Chong Peng, Can Song, Shuoyin Liu and Yubao Liu
Atmosphere 2024, 15(5), 556; https://doi.org/10.3390/atmos15050556 - 30 Apr 2024
Viewed by 1040
Abstract
The transport and dispersion (T&D) of catalyst particles seeded by weather modification aircraft is crucial for assessing their weather modification effects. This study investigates the capabilities of the Chinese geostationary weather satellite FY-4A for identifying the physical response of cloud seeding with AgI-based [...] Read more.
The transport and dispersion (T&D) of catalyst particles seeded by weather modification aircraft is crucial for assessing their weather modification effects. This study investigates the capabilities of the Chinese geostationary weather satellite FY-4A for identifying the physical response of cloud seeding with AgI-based catalysts and continuously monitoring its evolution for a weather event that occurred on 15 December 2019 in Henan Province, China. Satellite measurements are also used to verify an operational catalyst T&D algorithm. The results show that FY-4A exhibits a remarkable capability of identifying the cloud-seeding tracks and continuously tracing their evolution for a period of over 3 h. About 60 min after the cloud seeding, the cloud crystallization track became clear in the FY-4A tri-channel composite cloud image and lasted for about 218 min. During this time period, the cloud track moved with the cloud system about 153 km downstream (northeast of the operation area). An operational catalyst T&D model was run to simulate the cloud track, and the outputs were extensively compared with the satellite observations. It was found that the forecast cloud track closely agreed with the satellite observations in terms of the track widths, morphology, and movement. Finally, the FY-4A measurements show that there were significant differences in the microphysical properties across the cloud track. The effective cloud radius inside the cloud track was up to 15 μm larger than that of the surrounding clouds; the cloud optical thickness was about 30 μm smaller; and the cloud-top heights inside the cloud track were up to 1 km lower. These features indicate that the cloud-seeding catalysts led to the development of ice-phase processes within the supercooled cloud, with the formation of large ice particles and some precipitation sedimentation. Full article
(This article belongs to the Section Atmospheric Techniques, Instruments, and Modeling)
Show Figures

Figure 1

Figure 1
<p>(<b>a</b>) The aircraft flight trajectory on 15 December 2019; (<b>b</b>) the flight altitudes. The take-off and climbing phases are in green (C–A), the cloud-seeding phase is in purple (A–B), and the aircraft returning phase is in blue (B–C); The red triangle is Zhengzhou Xinzheng Airport.</p>
Full article ">Figure 2
<p>(<b>a</b>) FY-4A cloud top height and a portion of the flight path (black line) at 00:00 on 15 December 2019; (<b>b</b>) L-band radiosonde plot at Xingtai at 0800 h. The green line is the temperature, the blue line is the relative humidity curve, and the black lines are wind barbs. Green shades mark the cloud cover.</p>
Full article ">Figure 3
<p>The FY-4A cloud-top height and the aircraft flight track at different times during the cloud-seeding operation on 15 December 2019: (<b>a</b>) 03:00; (<b>b</b>) 03:15; (<b>c</b>) 03:30; (<b>d</b>) 03:49; (<b>e</b>) 04:23; (<b>f</b>) 05:30. The red line is for the period with supercooled water reported.</p>
Full article ">Figure 4
<p>FY4A tri-channel fusion cloud image on 15 December 2019: (<b>a</b>) 04:15; (<b>b</b>) 04:45; (<b>c</b>) 05:15; (<b>d</b>) 05:45; (<b>e</b>) 06:15; (<b>f</b>) 06:45; (<b>g</b>) 07:15; (<b>h</b>) 07:53. The black line shows the flight path of the aircraft.</p>
Full article ">Figure 5
<p>(<b>a</b>) Tri-channel fusion cloud image of the FY-3D polar-orbiting satellite at 05:40 on 15 December 2019; (<b>b</b>) that of the FY-4A satellite at 05:38. The black line marks the flight path of the aircraft, and the red lines represent the four cross-sections examined (see the red lines).</p>
Full article ">Figure 6
<p>The AgI particle number concentration calculated by the aircraft cloud-seeding T&amp;D model ((<b>a</b>): 03:30, (<b>b</b>): 03:49, (<b>c</b>): 04:23, (<b>d</b>): 05:30, (<b>e</b>): 05:45, (<b>f</b>): 06:15), the red line is for the second segment (cloud seeding between 03:15 and 03:49).</p>
Full article ">Figure 7
<p>The calculated catalyst dispersion region (red contours), overlapped on the FY-4A 0.47 channel cloud track: (<b>a</b>) 04:15; (<b>b</b>) 04:45; (<b>c</b>) 05:15; (<b>d</b>) 05:30; (<b>e</b>) 05:45; (<b>f</b>) 06:15. The blue line is the flight trajectory, and the yellow line marks the calculated dispersion area of the catalyst.</p>
Full article ">Figure 8
<p>FY-4A tri-channel fused cloud images at different times on 15 December 2019: (<b>a</b>) 04:35; (<b>b</b>) 04:45; (<b>c</b>) 05:15; (<b>d</b>) 05:30. The red line is for the fixed space position selected. The black line marks the aircraft’s flight path.</p>
Full article ">Figure 9
<p>Cloud effective particle radius at 163 min (05:30), and the locations of the four selected segments (red lines) across the cloud-seeding-induced cloud track.</p>
Full article ">Figure 10
<p>Microphysical properties across the cloud track for the 4 cross-sections shown in <a href="#atmosphere-15-00556-f009" class="html-fig">Figure 9</a> correspond to (<b>a</b>–<b>d</b>) of <a href="#atmosphere-15-00556-f010" class="html-fig">Figure 10</a>. The cloud effective radius is in blue; the cloud optical thickness is in purple; the cloud-top height is in orange; and the cloud-top temperature is in green.</p>
Full article ">
19 pages, 6697 KiB  
Article
Methane Retrieval from Hyperspectral Infrared Atmospheric Sounder on FY3D
by Xinxin Zhang, Ying Zhang, Fan Meng, Jinhua Tao, Hongmei Wang, Yapeng Wang and Liangfu Chen
Remote Sens. 2024, 16(8), 1414; https://doi.org/10.3390/rs16081414 - 16 Apr 2024
Viewed by 993
Abstract
This study utilized an infrared spotlight Hyperspectral infrared Atmospheric Sounder (HIRAS) and the Medium Resolution Spectral Imager (MERSI) mounted on FY3D cloud products from the National Satellite Meteorological Center of China to obtain methane profile information. Methane inversion channels near 7.7 μm were [...] Read more.
This study utilized an infrared spotlight Hyperspectral infrared Atmospheric Sounder (HIRAS) and the Medium Resolution Spectral Imager (MERSI) mounted on FY3D cloud products from the National Satellite Meteorological Center of China to obtain methane profile information. Methane inversion channels near 7.7 μm were selected based on the different distribution of methane weighting functions across different seasons and latitudes, and the selected retrieval channels had a great sensitivity to methane but not to other parameters. The optimization method was employed to retrieve methane profiles using these channels. The ozone profiles, temperature, and water vapor of the European Centre for Medium-Range Weather Forecasts (ECMWF) fifth-generation reanalysis data (ERA5) were applied to the retrieval process. After validating the methane profile concentrations retrieved by HIRAS, the following conclusions were drawn: (1) compared with Civil Aircraft for the Regular Investigation of the Atmosphere Based on an Instrument Container (CARIBIC) flight data, the average correlation coefficient, relative difference, and root mean square error were 0.73, 0.0491, and 18.9 ppbv, respectively, with lower relative differences and root mean square errors in low-latitude regions than in mid-latitude regions. (2) The methane profiles retrieved from May 2019 to September 2021 showed an average error within 60 ppbv compared with the Fourier transform infrared spectrometer (FTIR) station observations of the Infrared Working Group (IRWG) of the Network for the Detection of Atmospheric Composition Change (NDACC). The errors between the a priori and retrieved values, as well as between the retrieved and smoothed values, were larger by around 400–500 hPa. Apart from Toronto and Alzomoni, which had larger peak values in autumn and spring respectively, the mean column averaging kernels typically has a larger peak in summer. Full article
(This article belongs to the Section Atmospheric Remote Sensing)
Show Figures

Figure 1

Figure 1
<p>Data from four flights of the Civil Aircraft for the Regular Investigation of the Atmosphere Based on an Instrument Container (CARIBIC). The light blue line is flight 569 (1 May 2019), the light green line is flight 575 (14 August 2019), the red line is flight 586 (9 January 2020), and the light pink line is flight 587 (10 January 2020).</p>
Full article ">Figure 2
<p>Sensitivities of HIRAS channels located in 1200−1750 cm<sup>−1</sup> to variations of 1 K temperature (blue), 10% H<sub>2</sub>O (grey), 1% CO<sub>2</sub> (yellow),10% CH<sub>4</sub> (red), 10% N<sub>2</sub>O (pink), 10% O<sub>3</sub> (dark red) and 10% CO (light blue).</p>
Full article ">Figure 3
<p>CH<sub>4</sub> Jacobians (K/ppmv) from 1210 to 1400 cm<sup>−1</sup> with a spectral resolution of 0.625 cm<sup>−1</sup> for (<b>a</b>) high latitude between 90°S–30°S and 30–90°N in spring and summer (HS), (<b>b</b>) high latitude between 90°S–30°S and 30–90°N in spring and summer (HS), (<b>c</b>) low latitude (30°S–30°N) in spring and summer (LS) and (<b>d</b>) low latitude (30°S–30°N) in autumn and winter (LW). The different colors of the lines in the figure represent different channels (1210 to 1400 cm<sup>−1</sup>).</p>
Full article ">Figure 4
<p>CH<sub>4</sub> retrieval channels in (<b>a</b>) HS, (<b>b</b>) HW (light blue line), (<b>c</b>) LS and (<b>d</b>) LW.</p>
Full article ">Figure 5
<p>CH<sub>4</sub> comparison between HIRAS retrievals and smoothed CARIBIC with a 0.5° × 0.5° spatial resolution. (<b>a</b>) Flight no. 569, (<b>b</b>) flight no. 575, (<b>c</b>) flight no. 586, and (<b>d</b>) flight no. 587.</p>
Full article ">Figure 6
<p>CH<sub>4</sub> profile comparisons between a priori (pink line), smoothing Fourier transform infrared spectrometer (FTIR) (light green line) observation products, and HIRAS retrieval (light blue line), in FTIR sites.</p>
Full article ">Figure 7
<p>The mean column averaging kernels of CH<sub>4</sub> retrievals in FTIR stations. The blue, green, red, and black lines refer to spring, summer, autumn, and winter, respectively.</p>
Full article ">Figure 8
<p>Means of variations of CH<sub>4</sub> induced by variations of T (1–10 hPa: 2 K; 10–1000 hPa: 0.2 K) and H<sub>2</sub>O (5%) in Altzomoni and Rikubetsu.</p>
Full article ">Figure 9
<p>Comparison of Altzomoni and Rikubetsu observations with inversion results (<b>a</b>) and inversion results using changed channels. (<b>a</b>) Channels in opposite seasons in Altzomoni and (<b>b</b>) channels in high latitudes in Rikubetsu.</p>
Full article ">
Back to TopTop