Fengyun-3D Research Articles

Winter wheat identification in China through FY-3D NDVI&EVI satellite data and DNN Fusion

ABSTRACT Medium Resolution Spectral Imager-II (MERSI-II) on FengYun 3D (FY-3D) meteorological satellite is a visible and infrared spectral imaging instrument comparable with Moderate Resolution Imaging Spectroradiometer (MODIS), which provides ample opportunities for regional crop mapping. However, current research may have overlooked the potential of combining FY-3D MERSI-II data with deep neural network (DNN) algorithms for winter wheat identification. This research utilizes a DNN method to train a nonlinear model with the synthetic normalized difference vegetation index (NDVI) and enhanced vegetation index (EVI) at a 250 m resolution from FY-3D MERSI-II as the eigenvalues, enabling to accurately identify spatial distribution of winter wheat in China’s main production regions and filling the gap in the application of FY-3D MERSI-II data. The results indicate that the winter wheat identification with the combination of NDVI and EVI achieved an accuracy of 95.32%, with a Kappa coefficient of 0.87 and a determination coefficient (R 2 ) of 0.73, improving the overall accuracy by 1.27% and 1.07% compared to models trained solely using NDVI or EVI, respectively. In terms of spatial distribution, the overlap rate with the winter wheat dataset of China from MODIS at a 500 m resolution exceeded 95%, with an area bias of 2.67%. This study indicates that the medium resolution FY-3D MERSI-II remote sensing dataset combined with the DNN algorithm can accurately identify crop information over large areas.

Inter-Calibration of Passive Microwave Satellite Brightness Temperature Observations between FY-3D/MWRI and GCOM-W1/AMSR2

Microwave sensors possess the capacity to effectively penetrate through clouds and fog and are widely used in obtaining soil moisture, atmospheric water vapor, and surface temperature measurements. Long time-series datasets play a pivotal role in climate change studies. Unfortunately, the lifespan of operational satellites often falls short of the needs of these extensive datasets. Hence, comparing and cross-calibrating sensors with similar configurations is paramount. The Microwave Radiation Imager (MWRI) onboard Fengyun-3D (FY-3D) is the latest generation of satellite-based microwave remote sensing instruments in China, and its data quality and application prospects have attracted widespread attention. To comprehensively assess the data quality of MWRI, a comparison of the orbital brightness temperature (TB) data between FY-3D/MWRI and Global Change Observation Mission 1st-Water (GCOM-W1)/Advanced Microwave Scanning Radiometer 2 (AMSR2) is conducted, and then a calibration model is established. The results indicate a strong correlation between the two sensors, with a correlation coefficient exceeding 0.9 across all channels. The mean bias ranges from −1.5 K to 0.15 K. Notably, the bias of vertical polarization is more pronounced than that of horizontal polarization. The TB distribution patterns and temporal evolutions are highly consistent for both sensors, particularly under snow and ice. The small intercepts and close-to-1 slopes obtained during calibration further demonstrate the minor data differences between the two sensors. However, the calibration process effectively reduces the existing errors, and the calibrated FY-3D/MWRI TB data are closer to GCOM-W1/AMSR2, with a mean bias approximately equal to 0 K and a correlation coefficient exceeding 0.99. The excellent consistency of the TB data between the two sensors provides a vital data basis for retrieving surface parameters and establishing long time-series datasets.

A Physics-Based Method for Retrieving Land Surface Emissivities from FengYun-3D Microwave Radiation Imager Data

The passive microwave land surface emissivity (MLSE) plays a crucial role in retrieving various land surface and atmospheric parameters and in Numerical Weather Prediction models. The retrieval accuracy of MLSE depends on many factors, including the consistency of the input data acquisition time. The FengYun-3D (FY-3D) polar-orbiting meteorological satellite, equipped with passive microwave and infrared bands, offers time-consistent data crucial for MLSE retrieval. This study proposes a physics-based MLSE retrieval algorithm using all the input data from the FY-3D satellite. Based on the retrieved MLSE, the spatial distribution of the MLSE is closely correlated with the land cover types and topography. Lower emissivities prevailed over barren or sparsely vegetated regions, river basins, and coastal areas. Higher emissivities dominated densely vegetated regions and mountainous areas. Moderate emissivities dominated grasslands and croplands. Lower-frequency channels showed larger emissivity differences with different polarizations than those of higher-frequency channels in barren or sparsely vegetated regions. The MLSE across densely vegetated land areas, mountainous areas, and deserts showed small seasonal variations. However, woody savannas, grasslands, croplands, and seasonal snow-covered areas showed noticeable seasonal variations. For most land cover types, the differences between vertically and horizontally polarized emissivities remained relatively constant across seasons. However, certain grasslands in eastern Inner Mongolia and southern Mongolia showed clear seasonal variations. It is very difficult to verify the MLSE on a large scale. Consequently, the possible error sources in the retrieved MLSE were analyzed, including the brightness temperature errors (correlation coefficient ranging from 0.92 to 0.99) and the retrieved land surface temperature errors (Root Mean Square Error was 3.34 K and relation coefficient was 0.958).

Enhanced precise orbit determination for GPS and BDS-2/3 with real LEO onboard and ground observations

Due to the current geographical situation in China that it is difficult to achieve globally well distributed ground monitoring stations for Chinese BeiDou Navigation System (BDS), the BDS Medium Earth Orbit (MEO) satellite orbital accuracy can be seriously affected by the insufficient length of continuous observation arc based on regional ground network. What is more, the highly static characteristics of BDS Geosynchronous Orbiting (GEO) satellites resulting in meter-level orbit determination accuracy is far from satisfactory for precise positioning services. In order to overcome these problems, the low Earth orbit (LEO) satellites serving as fast-moving space monitoring stations can be used as a good supplement for the current BDS ground monitoring stations layout. In this study, real measured spaceborne data from two LEO satellites, namely Fengyun-3D (FY-3D) and Tianjin University No. 1 (TJU-01), are collected and studied for the augmentation of BDS precise orbit determination (POD). The multipath errors of spaceborne code observations are comprehensively analyzed. It is found that not only BDS-2 relevant satellite-induced code bias but also near-field multipath errors, especially distinct in GPS, exist in a spaceborne environment. The piece-wise linear model is therefore established and applied to correct the code observations to reduce the multipath errors in the integrated POD processing. The GPS and BDS-2/3 POD solutions are compared with orbital overlapping differences and Satellite Laser Ranging (SLR) residuals validation. It is demonstrated that including the two LEO satellites can improve the orbit determination accuracy of GPS and BDS-2/3 significantly under a regional ground network, especially in the along component. The average overlapping orbital Root Mean Square (RMS) is decreased by 18.4 %, 34.3 %, 24.9 %, 21.4 %, and 7.8 % for the GPS, GEO, Inclined Geosynchronous Orbit (IGSO), BDS-2 MEO, and BDS-3 MEO satellites, respectively, in the along component. The RMS of SLR residuals is decreased from 37.1 cm to 17.8 cm for the GEO satellites, and the overall improvement is about 8 %-12 % for the BDS MEO satellites. There has also been a slight improvement under a global network, with an average improvement of 2–3 cm in the SLR residuals RMS. However, significant systematic errors are also observed in the SLR residuals in BDS GEO, IGSO, and MEO satellites, and are expected to be improved with refined observation models based on more LEO satellites, i.e., the phase center variation (PCV).

Impact of the Detection Channels Added by Fengyun Satellite MWHS-II at 183 GHz on Global Numerical Weather Prediction

Fine spectral detection can basically solve the problem of low vertical resolution at the 183 GHz water-vapor absorption line, and it is expected to become one of the main methods for next-generation geostationary and polar-orbiting satellites. Here, using data from Microwave Humidity Sounder II (MWHS-II) onboard the Chinese Fengyun 3D (FY-3D) satellite in the Global/Regional Assimilation and Prediction System (GRAPES) Four-Dimensional Variational (4D-Var) system of the China Meteorological Administration (CMA), we explore the assimilation application of the water-vapor absorption line at 183.31 ± 1 GHz, 183.31 ± 3 GHz and 183.31 ± 7 GHz, as well as 183.31 ± 1.8 GHz and 183.31 ± 4.5 GHz, two added channels, to assess the impact of adding the 183.31 ± 1.8 GHz and 183.31 ± 4.5 GHz sampling channels on data assimilation and numerical weather prediction. Our findings reveal a significant increase in the specific-humidity increment, which in the middle–upper troposphere is numerically much larger than in the lower troposphere. Specifically, the assimilation of 183.31 ± 1.8 GHz observations, positioned near the center of the water-vapor absorption line, results in a pronounced adjustment compared with the 183.31 ± 4.5 GHz observations. And under the strong constraint of the numerical model, the Root Mean Square Error (RMSE) of the wind field diminishes more significantly (by an average of 2–4%) after assimilating the water-vapor observations at greater heights.

Change assessment of the Tonga island surface temperature before and after volcanic eruption based on FY-3D time series dataset

Taking advantage of the 250 m resolution and daily observation of dual thermal infrared channels data of the Medium Resolution Spectral Imager (MERSI) onboard the Fengyun 3D (FY-3D) satellite, a refined thermal infrared island surface temperature (ILST) inversion model is established through radiation transmission simulation. And the Tonga ILST change maps before and after the volcanic eruption are produced and analyzed based on ILST time series dataset. The results show that the ILST decreases by 2 K –3 K within the five days after the eruption due to the amount of volcanic ash deposition, but increases significantly by 2 K –6 K within the fifth to tenth days after the eruption. The trend of time-limited cooling, warming, and then gradual recovery after volcanic eruption is very obvious. Through the manual cleaning and the scouring by seawater or rainwater to the volcanic ash covered ground, more parts of the island surface restored its original state. Hence, the Tonga ILST also increased gradually within a week thereafter. However, in regions with local vegetation, the recovery time of ILST would be longer. The temporal and spatial resolution of FY-3D MERSI meets the need for monitoring the changes of the island surface states before and after the volcanic eruption.

Comparative Study of the Atmospheric Gas Composition Detection Capabilities of FY-3D/HIRAS-I and FY-3E/HIRAS-II Based on Information Capacity

Fengyun-3E (FY-3E)/Hyperspectral Infrared Atmospheric Sounder-II (HIRAS-II) is an extension Fengyun-3D (FY-3D)/HIRAS-I. It is crucial to fully explore and analyze the detection capabilities of these two instruments for atmospheric gas composition. Based on the observed spectral data from the infrared hyperspectral detection instruments FY-3D/HIRAS-I and FY-3E/HIRAS-II, simulated radiance data and Jacobian matrices are obtained using the Rapid Radiative Transfer Model RTTOV (Radiative Transfer for TOVS (TIROS Operational Vertical Sounder)). By perturbing temperature (T), surface temperature (Tsurf), water vapor (H2O), ozone (O3), carbon dioxide (CO2), methane (CH4), carbon monoxide (CO), and nitrous oxide (N2O), the brightness temperature differences before and after the perturbations are calculated to analyze the sensitivity of temperature and various atmospheric gas components. The Improved Optimal Sensitivity Profile (OSP) algorithm is used to select the channels for atmospheric gas retrieval. The observation error covariance and background error covariance matrices are calculated, and then the information capacity is calculated, specifically the degrees of freedom for signal(DFS) and the entropy reduction (ER). Based on this, a comparative analysis is conducted on the information capacity of atmospheric water vapor and ozone components contained in the hyperspectral detection data from HIRAS-I and HIRAS-II instruments, respectively, to explore the retrieval capabilities of the two instruments for atmospheric gas components. We selected clear-sky data from the African oceanic region and the Chinese Yangtze River Delta terrestrial region for quantitative analysis of the information capacity of HIRAS-I and HIRAS-II. The results show that FY-3D/HIRAS-I and FY-3E/HIRAS-II exhibit different sensitivities to atmospheric gas components. In different experimental regions, temperature and water vapor show the most dramatic sensitivity changes, followed by ozone, methane, and nitrous oxide, while carbon monoxide and carbon dioxide exhibit the lowest variability. Regarding channel selection, HIRAS-II identifies more gas channels compared to HIRAS-I. The experiments concluded that HIRAS-II has a significantly higher information capacity than HIRAS-I, and the information capacity of atmospheric gas components varies across different experimental regions. Water vapor and ozone exhibit the highest information capacity, followed by nitrous oxide and methane, while carbon monoxide and carbon dioxide demonstrate the lowest capacity. The H2O ER (DFS) contained in FY-3E/HIRAS-II is 1.51 (0.35) higher than that in FY-3D/HIRAS-I, the O3 ER (DFS) in FY-3E/HIRAS-II is 1.51 (0.36) higher than that in FY-3D/HIRAS-I, while the N2O ER (DFS) in FY-3E/HIRAS-II is 0.17 (0.19) higher and the CH4 ER (DFS) is 0.07 (0.04) higher than that in FY-3D/HIRAS-I.

The Uncertainty of SNO Cross-Calibration for Satellite Infrared Channels

The on-orbit radiometric calibration is a fundamental task in quantitative remote sensing applications. A widely used calibration method is the cross-calibration based on Simultaneous Nadir Observation (SNO), which involves using high-precision reference instruments to calibrate lower-precision onboard instruments. However, despite efforts to match the observation time, spatial location, field geometry, and instrument spectra, errors can still be introduced during the matching processes and linear regression analysis. This paper focuses on the error generated by sample matching and the error fitting method generated by the sample fitting method. An error propagation analysis is performed to develop a generic model for assessing the uncertainty of the SNO cross-calibration method itself in meteorological satellite infrared channels. The model is validated using the payload parameters of the Hyperspectral Infrared Atmospheric Sounder (HIRAS) and the Medium Resolution Spectral Imager (MERSI) instruments aboard the FengYun-3D (FY-3D). Simulation experiments are performed considering typical bright temperatures, different background fields, and varying matching threshold conditions. The results demonstrate the effectiveness of the proposed model in capturing the error propagation chain in the SNO cross-calibration process. The model provides valuable insight into error analysis in the SNO cross-calibration method and can assist in determining the optimal sample matching threshold for achieving radiometric calibration accuracy.

Research on infrared hyperspectral remote sensing cloud detection method based on deep learning

ABSTRACT Infrared hyperspectral is susceptible to clouds, and accurately identifying whether the hyperspectral infrared sounder data is polluted by clouds is of great significance for numerical weather prediction and atmospheric parameter inversion. Since the complex spectral characteristics of clouds, the existing spectral threshold methods and machine learning methods have the difficulties of undetermined threshold and clear field of view (FOV) missed and false detections. In order to improve the cloud recognition accuracy of infrared hyperspectral data, three end-to-end cloud detection models combining deep neural network (DNN) and convolutional neural network (CNN) and long short-term memory network (LSTM) are proposed. In this paper, taking the High Spectral Infrared Atmospheric Sounder (HIRAS) equipped with Fengyun-3D (FY-3D) satellite as the research object, based on the same platform Moderate Resolution Spectral Imager-II (MERSI-II) cloud mask (CLM) product, the HIRAS Cloud dataset is established, and the accuracy test and qualitative analysis are carried out by using the test datasets and Typhoon Siamba, July 3, 2022, as well as the earth observation scene under the conditions of ice and snow surface. The test datasets analysis results show that the cloud detection accuracy of CNN and CNN-LSTM model is stable at 0.96, and the false alarm rate of cloud is 0.035 and 0.036, respectively, and the detection ability of DNN model is slightly inferior to the former two in the same hidden layer, with an accuracyof 0.94. In further qualitative research, we found that the CNN-LSTM model has high accuracy and robustness in infrared hyperspectral cloud detection, and the detection results in a variety of surface scenarios are consistent with the actual situation of whether clouds occur in the FOV of the instrument. Compared with CLM products, it can better identify clear ocean scenes, and provide fast and efficient cloud detection reference for data assimilation systems.

The Capabilities of FY-3D/MERSI-II Sensor to Detect and Quantify Thermal Volcanic Activity: The 2020–2023 Mount Etna Case Study

Satellite data provide crucial information to better understand volcanic processes and mitigate associated risks. In recent years, exploiting the growing number of spaceborne polar platforms, several automated volcanic monitoring systems have been developed. These, however, rely on good geometrical and meteorological conditions, as well as on the occurrence of thermally detectable activity at the time of acquisition. A multiplatform approach can thus increase the number of volcanological-suitable scenes, minimise the temporal gap between acquisitions, and provide crucial information on the onset, evolution, and conclusion of both transient and long-lasting volcanic episodes. In this work, we assessed the capabilities of the MEdium Resolution Spectral Imager-II (MERSI-II) sensor aboard the Fengyun-3D (FY-3D) platform to detect and quantify heat flux sourced from volcanic activity. Using the Middle Infrared Observation of Volcanic Activity (MIROVA) algorithm, we processed 3117 MERSI-II scenes of Mount Etna acquired between January 2020 and February 2023. We then compared the Volcanic Radiative Power (VRP, in Watt) timeseries against those obtained by MODIS and VIIRS sensors. The remarkable agreement between the timeseries, both in trends and magnitudes, was corroborated by correlation coefficients (ρ) between 0.93 and 0.95 and coefficients of determination (R2) ranging from 0.79 to 0.84. Integrating the datasets of the three sensors, we examined the effusive eruption of Mount Etna started on 27 November 2022, and estimated a total volume of erupted lava of 8.15 ± 2.44 × 106 m3 with a Mean Output Rate (MOR) of 1.35 ± 0.40 m3 s−1. The reduced temporal gaps between acquisitions revealed that rapid variations in cloud coverage as well as geometrically unfavourable conditions play a major role in thermal volcano monitoring. Evaluating the capabilities of MERSI-II, we also highlight how a multiplatform approach is essential to enhance the efficiency of satellite-based systems for volcanic surveillance.

Real‐time precise orbit determination for FY‐3C and FY‐3D based on BDS and GPS onboard observation

AbstractThe code multipath error associated with elevation in BDS‐2 is one of the main factors whose impact on the precision of real‐time reduced‐dynamic precise orbit determination (POD) using spaceborne multi‐global navigation satellite systems (GNSS) observation data for FengYun‐3C (FY‐3C) and FengYun‐3D (FY‐3D) satellites. The aim is to construct a code multipath piece‐wise linear error model and analyse the contribution of BDS‐2 to the performance of low earth orbit POD. Conclusions are drawn out from the experimental results: (1) the real‐time POD average precision of FY‐3C and FY‐3D is improved by about 11% and 12% after the correction of BDS‐2 code multipath error; (2) Due to the small number of available BDS‐2 satellites and the limited accuracy of geostationary earth orbit (GEO) satellite real‐time orbit products, the precision of the real‐time overlapping comparison for FY‐3C/FY‐3D by using the BDS‐2 onboard observation data can only reach the decimetre level. However, with the BDS‐2 signal capture capability of the FY‐3D onboard GNSS Occultation Sounder upgraded, the average precision of the BDS‐2‐based POD for FY‐3D is better than that of FY‐3C in each direction and with 10% improvement in the average 3 dimensional‐root mean square (3D‐RMS); (3) Cases of onboard global positioning system (GPS), GPS + BDS‐2 (with GEO), and GPS + BDS‐2 (without GEO) for real‐time POD can meet 5–8 cm precision requirement. The average 3D‐RMS of the real‐time POD for FY‐3C/FY‐3D, which utilises onboard GPS + BDS‐2 data is inferior only to GPS data, owing to the influence of the orbit accuracy of GEO satellites. However, after excluding the influence of GEO satellites, the average result (3D‐RMS) of the real‐time POD for FY‐3C (6.11 cm) and FY‐3D (5.95 cm) is better than those of other cases.

Automatic Detection of Daytime Sea Fog Based on Supervised Classification Techniques for FY-3D Satellite

Polar-orbiting satellites have been widely used for detecting sea fog because of their wide coverage and high spatial and spectral resolution. FengYun-3D (FY-3D) is a Chinese satellite that provides global sea fog observation. From January 2021 to October 2022, the backscatter and virtual file manager products from CALIPSO (Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observation) were used to label samples of different atmospheric conditions in FY-3D images, including clear sky, sea fog, low stratus, fog below low stratus, mid–high-level clouds, and fog below the mid–high-level clouds. A 13-dimensional feature matrix was constructed after extracting and analyzing the spectral and texture features of these samples. In order to detect daytime sea fog using a 13-dimensional feature matrix and CALIPSO sample labels, four supervised classification models were developed, including Decision Tree (DT), Support Vector Machine (SVM), K-Nearest Neighbor (KNN), and Neural Network. The accuracy of each model was evaluated and compared using a 10-fold cross-validation procedure. The study found that the SVM, KNN, and Neural Network performed equally well in identifying low stratus, with 85% to 86% probability of detection (POD). As well as identifying the basic components of sea fog, the SVM model demonstrated the highest POD (93.8%), while the KNN had the lowest POD (92.4%). The study concludes that the SVM, KNN, and Neural Network can effectively distinguish sea fog from low stratus. The models, however, were less effective at detecting sub-cloud fog, with only 11.6% POD for fog below low stratus, and 57.4% POD for fog below mid–high-level clouds. In light of this, future research should focus on improving sub-cloud fog detection by considering cloud layers.

Retrieving leaf area index from FY-3D MERSI-II data using a sensor-adaptive algorithm

ABSTRACT Leaf area index (LAI), officially listed as one of essential climate variables, quantifies the structure and amount of vegetation and characterizes the interaction between vegetation and climate. The advanced MEdium Resolution Spectral Imager (MERSI-II) onboard FengYun-3D (FY-3D) can provide twice-daily global observations of earth at a spatial resolution of 250 m. Therefore, it has great potential for promoting the improvement of global LAI products and boosting the development of earth system modelling. However, the existing methods are mostly sensor-specific, they can not be directly applied to MERSI-II data for LAI generation. In this paper, we proposed a sensor-adaptive approach for LAI estimation based on MERSI-II observations. This method is composed of an LAI retrieval look-up table based on radiative transfer theory, which was calibrated by a global optimizing algorithm to adapt MERSI-II characteristics, and a backup algorithm training neural networks with MERSI-II Normalized Difference Vegetation Index (NDVI) and Moderate Resolution Imaging Spectroradiometer (MODIS) LAI. We evaluated the MERSI-II LAI retrievals by intercomparison to MODIS products over mainland China and direct validation using ground-based upscaled LAI reference maps. The assessments demonstrate that (1) the derived MERSI-II LAI products agree well with the MODIS benchmarks in both spatial and temporal respects; (2) compared to the MODIS LAI, the MERSI-II LAI retrievals with less gaps show great potential in improving spatiotemporal continuity; (3) validation versus in-situ measurements reveals acceptable accuracy of the MERSI-II LAI products with a R2 of 0.85 and RMSE of 0.82; (4) the proposed parametric optimization strategy could successfully transplant MODIS LAI algorithm to MERSI-II data.

Temperature and Relative Humidity Profile Retrieval from Fengyun-3D/VASS in the Arctic Region Using Neural Networks

In this study, a new technique is proposed to retrieve temperature and relative humidity profiles under clear sky conditions in the Arctic region based on the artificial neural network (ANN) algorithm using Fengyun-3D (FY-3D) vertical atmospheric sounder suit (VASS: HIRAS, MWTS-II, and MWHS-II) observations. This technology combines infrared (IR) and microwave (MW) observations to improve retrieval accuracy in the middle and low troposphere by reducing the sensitivity of the neural networks (NNs) to cloud coverage. The approach was compared against other methods available in the literature on retrieving profiles only from FY-3D/HIRAS data. Furthermore, its retrieval performance was tested by comparing the NNs’ prediction accuracy versus the corresponding FY-3D/VASS and Aqua/AIRS L2 products. The results showed that: (1) NNs retrieval accuracy is higher during the warm season and over the ocean; (2) the retrieval accuracy of NNs has been significantly improved compared with satellite L2 products; (3) referring to radiosonde observations, the retrieval accuracy of NNs below 600 hPa is effectively improved by adding the information of the MW channel, especially on land where cloud clearing is more difficult. The root mean square error (RMSE) of temperature and relative humidity in the cold season were reduced by 0.3 K and 2%, respectively. The advanced NNs proposed herein offer a more stable retrieval performance compared with NNs built only by FY-3D/HIRAS data. The study results indicated the potential value in time and space domain of the NN algorithm in retrieving temperature and relative humidity profiles of the Arctic region from FY-3D/VASS observations under clear-sky conditions. All in all, this work enhances our knowledge towards improving operational use of FY-3D satellite data in the Arctic region.

Evaluation of oceanic precipitable water vapor products from Microwave Radiation Imager (MWRI) onboard the Fengyun-3D satellite

Precipitable water vapor (PWV) products from the Microwave Radiation Imager (MWRI) onboard the Fengyun-3D (FY-3D), a second-generation polar orbit meteorological satellite of China, have a high spatiotemporal resolution and play an increasingly key role in the atmosphere and climate analysis. We focus on the FY-3D/MWRI PWV products in the global ocean area and comprehensively evaluate their quality by using the coastal radiosonde data from the Integrated Global Radiosonde Archive (IGRA) along with the PWV reanalysis products from the fifth generation of European Centre for Medium-Range Weather Forecasts Atmospheric Reanalysis (ERA5). With respect to the radiosonde data, the correlation coefficient is 0.97, the root mean squared error (RMSE) is 2.71 mm, and the mean bias and mean relative bias is 0.36 mm and 2.57% respectively. The discrepancy between FY-3D/MWRI and radiosonde PWV is smallest at high latitudes and is biggest in the ocean around Africa. The seasonal variation characteristics perform best in spring and worst in summer. We also find that FY-3D/MWRI PWV products are in consistent with the ERA5 PWV, with a RMSE of 2.39 mm and 2.85 mm for daily and monthly data, respectively. Among globally 551,213 daily oceanic PWV grid data, there are 11.3% of matched points, which are mainly located between 15°S and 30°N latitudes, has a RMSE larger than 3 mm, while 68.5% of points have a RMSE of less than 2 mm, meeting the requirements of climate research.

A high-precision aerosol retrieval algorithm for FY-3D MERSI-II images

The Medium Resolution Spectral Imager-II (MERSI-II) onboard the recently launched Chinese Fengyun-3D (FY-3D) satellite has great capability in detecting global aerosols as it includes aerosol bands similar to Moderate Resolution Imaging Spectroradiometer (MODIS). However, to date, aerosol retrieval based on MERSI-II is still limited to dark target regions and there is no official aerosol products for the MERSI-II. This study focuses on developing a high-precision algorithm to retrieve aerosol optical depth (AOD) suitable for entire land areas (except snow/ice and inland waters) based on MERSI-II measurements. Considering both the accuracy and retrieval efficiency, a new cost function is constructed based on (1) the fact that the AOD (550 nm) retrieved independently from different bands should be the same with the correct aerosol model, and (2) the assumption that the aerosol types are the same in the 5 × 5 km spatial range. The retrieval method based on the new cost function is nearly 50 times faster than most current methods using iterative calculations. To extend the application adaption of the FY-3D MERSI-II AOD retrieval and avoid the errors caused by the surface Lambertian hypothesis, a bidirectional reflectance distribution function (BRDF) database is built using MODIS products. Eight candidate aerosol models in different natural zones of China are constructed based on AERONET aerosol products from 2010 − 2021. The new method is applied to MERSI-II images over China and validated against ground-based measurements at 14 stations from 2020 to 2021. MODIS aerosol products from three operational algorithms are also used for comparison purposes. The results show that MERSI-II AOD retrievals agree well with the ground-based measurements with correlation coefficient (R), root mean square error (RMSE), and relative mean bias (RMB) of 0.913, 0.123, and 0.955, respectively. In addition, 72.19 % of AOD matchups fall within the expected error (EE) envelopes. The MERSI-II retrievals show higher accuracy than that of MODIS dark target (DT) and deep blue (DB) products and comparable accuracy of the MODIS Multi-Angle Implementation of Atmospheric Correction (MAIAC) product. MERSI-II AOD also shows higher stability in terms of spatial and temporal and better performance under heavy aerosol loading conditions than MODIS products. A good AOD agreement with R from 0.777 to 0.863 between MERSI-II and MODIS products is found over the land of China. The new method showing high retrieval efficiency and accuracy has great potential to be operationally applied on AOD retrieval for MERSI-II.

A Study on the Retrieval of Ozone Profiles Using FY-3D/HIRAS Infrared Hyperspectral Data

Atmospheric ozone is a pollutant gas that has an important influence on the process of atmospheric radiation transmission and climate change. The Fengyun-3D (FY-3D) satellite Hyperspectral Infrared Atmospheric Sounder (HIRAS) has better spectral performance than other remote sensing payloads. Its observation radiation data contains abundant atmospheric vertical information, which can be used for ozone retrieval, but there are no ozone profile business products being generated at present. Therefore, for the mainland of Hong Kong, based on HIRAS infrared hyperspectral observation data, we used the traditional one-dimensional variational (1D-VAR) physical retrieval algorithm, combined with the radiative transfer model for TOVS (RTTOV), and selected the spectrum channel according to the optimal sensitive profile algorithm. The artificial neural network (ANN) algorithm was used to optimize the prior profiles, and the atmospheric ozone profile retrieval system was established. Finally, a set of ozone profile retrieval schemes suitable for FY-3D/HIRAS were summarized. We used ERA5 reanalysis data and World Ozone and Ultraviolet Radiation Data Centre (WOUDC) data to determine true values. The retrieval results were compared with Global Forecast System (GFS) forecast data, Ozone Mapping and Profile Suite (OMPS) ozone products, and Atmospheric Infrared Sounder (AIRS) ozone products. The results show that our ozone profile retrieval scheme makes up for the shortcomings of the conventional physical methods in some atmospheric pressure levels. The overall root-mean-square error (RMSE) of the ozone from the ground to the top of the stratosphere is within 30% on average, which was better than that for the GFS forecast data; the retrieval accuracy RMSE (%) was less than 20% in the pressure layer with the highest ozone concentration (15–25 hPa), which is better than that of OMPS ozone products and AIRS ozone products. The retrieval results prove that FY3D/HIRAS observation data allow ozone profile retrieval. This paper provides a reference for generating independent HIRAS ozone profile product data sets in business, and provides support for the subsequent application of Fengyun-3 series meteorological satellites in atmospheric parameter remote sensing.

A Simplified Coastline Inflection Method for Correcting Geolocation Errors in FengYun-3D Microwave Radiation Imager Images

Passive microwave (PMW) sensors are popularly applied to Earth observations. However, the satellite PMW radiometer data sometimes have non-negligible errors in geolocation. Coastline inflection methods (CIMs) are widely used to improve geolocation errors of PMW images. However, they commonly require accuracy satellite flight parameters, which are difficult to obtain by users. In this study, a simplified coastline inflection method (SCIM) is proposed to correct the geolocation errors without demanding for the satellite flight parameters. SCIM is applied to improve geolocation errors of FengYun-3D (FY-3D) Microwave Radiation Imager (MWRI) brightness temperature images from 2018 and 2019. It reduces the geolocation errors of MWRI images to 0.15 pixels in the along-track and cross-track direction. This means reductions of 75% and 86% of the geolocation errors, respectively. The mean brightness temperature differences between the ascending and descending MWRI images are reduced by 34%, demonstrating the improved geolocation accuracy of SCIM. The corrected images are also used to estimate Arctic sea ice concentration (SIC). By comparing with SICs retrieved from the un-corrected images, the root mean square error (RMSE) and mean absolute error (MAE) of the SICs from the corrected images are reduced from 13.7% to 10.2% and 8.9% to 6.9%, respectively. The mean correlation coefficient (R) increases from 0.91 to 0.95. All these results indicate that SCIM can reduce geolocation errors of satellite-based PMW images significantly. As SCIM is very simple and easy to be applied, it could be a useful method for users of PMW images.

Impacts of 3DEnVar-Based FY-3D MWHS-2 Radiance Assimilation on Numerical Simulations of Landfalling Typhoon Ampil (2018)

The module for assimilating radiance data of the Microwave Humidity Sounder-2 (MWHS-2) onboard the Feng Yun 3D (FY-3D) satellite is built in the Weather Research and Forecasting (WRF) model data assimilation (WRFDA) system. The CONV, 3DVar, and EnVar experiments are conducted to investigate the impact of assimilating the new humidity sounder based on Typhoon Ampil (2018). Both the 3DVar and EnVar experiments assimilate FY-3D MWHS-2 radiance data on top of the conventional data, while the CONV experiment only applies conventional data. In the EnVar experiment, notable geopotential height increment is observed around the typhoon, leading the typhoon to move northeast. In addition, the moisture field is improved to some extent. Finally, from the analysis of the dynamic field of the typhoon, it can be found that the EnVar experiment can adjust the dynamic structure of the typhoon. Furthermore, the assimilation of FY-3D MWHS-2 radiance data reduces the forecast error of the typhoon track and intensity. Additionally, the precipitation skill is improved in terms of rainfall pattern and the verification score. This improvement in the precipitation may be closely related to the features of the circulation structure concerning the evolution of the typhoon. The improved prediction of the position and intensity of rainbands in the FY-3D MWHS-2 radiance data assimilation experiment corresponds to a better prediction of typhoon structure.

Retrieving the Motion of Beaufort Sea Ice Using Brightness Temperature Data from FY-3D Microwave Radiometer Imager.

Sea ice is an important marine phenomenon in the Arctic region, and it is of great importance to study the motion of Arctic sea ice in the present day when its melting is accelerated by global warming. This study proposes a method to retrieve the motion of sea ice based on the maximum cross-correlation (MCC) and the successive correction method (SCM). The proposed method can apply different scales of search ranges to template matching according to the location of sea ice in the Arctic area. In addition, the data assimilation method can assign different weights to different data. We used 36.5 GHz and 89 GHz brightness temperature (Tb) data from the microwave radiometer imager (MWRI) aboard the Fengyun-3D (FY-3D) satellite, for the first time in the literature, to retrieve the sea ice motion in the Beaufort Sea from January to April 2019. The retrieved sea ice motion results were in good agreement with those obtained from the motion of the buoys. Compared with the data from the buoys, the root mean-squared error (RMSE) of the sea ice motion retrieved from FY-3D/MWRI Tb data was 1.1418 cm/s in the zonal direction and 1.0481 cm/s in the meridional direction, and the mean absolute error (MAE) between them was 0.7166 cm/s in the zonal direction and 0.6777 cm/s in the meridional direction. The RMSE between the sea ice motion obtained from the National Snow and Ice Data Center (NSIDC) and the motion of the buoys was 0.9515 cm/s in the zonal direction and 0.67003 cm/s in the meridional direction, and the MAE between them was 0.6576 cm/s in the zonal direction and 0.4922 cm/s in the meridional direction. The RMSE of daily average velocity from the FY-3D/MWRI results and NSIDC data product was 2.2726 cm/s in zonal and 1.9270 cm/s in meridional, and the MAE was 1.5103 cm/s in zonal and 1.1071 cm/s in zonal. The density of the merged data was higher than that obtained from a single polarization or frequency in this paper. The results indicate that FY-3D/MWRI Tb data can retrieve the sea ice motion successfully.