Towards Consistent Soil Moisture Records from China’s FengYun-3 Microwave Observations

: Soil moisture plays an essential role in the land-atmosphere interface. It has become necessary to develop quality large-scale soil moisture data from satellite observations for relevant applications in climate, hydrology, agriculture, etc. Speciﬁcally, microwave-based observations provide more consistent land surface records because they are unhindered by cloud conditions. The recent microwave radiometers onboard FY-3B, FY-3C and FY-3D satellites launched by China’s Meteorological Administration (CMA) extend the number of available microwave observations, covering late 2011 up until the present. These microwave observations have the potential to provide consistent global soil moisture records to date, ﬁlling the data gaps where soil moisture estimates are missing in the existing records. Along these lines, we studied the FY-3C to understand its added value due to its unique time of observation in a day (ascending: 22:15, descending: 10:15) absent from the existing satellite soil moisture records. Here, we used the triple collocation technique to optimize a benchmark retrieval model of land surface temperature (LST) tailored to the observation time of FY3C, by evaluating various soil moisture scenarios obtained with different bias-imposed LSTs from 2014 to 2016. The globally optimized LST was used as an input for the land parameter retrieval model (LPRM) algorithm to obtain optimized global soil moisture estimates. The obtained FY-3C soil moisture observations were evaluated with global in situ and reanalysis datasets relative to FY3B soil moisture products to understand their differences and consistencies. We found that the RMSEs of their anomalies were mostly concentrated between 0.05 and 0.15 m 3 m − 3 , and correlation coefﬁcients were between 0.4 and 0.7. The results showed that the FY-3C ascending data could better capture soil moisture dynamics than the FY-3B estimates. Both products were found to consistently complement the skill of each other over space and time globally. Finally, a linear combination approach that maximizes temporal correlations merged the ascending and descending soil moisture observations separately. The results indicated that superior soil moisture estimates are obtained from the combined product, which provides more reliable global soil moisture records both day and night. Therefore, this study aims to show that there is merit to the combined usage of the two FY-3 products, which will be extended to the FY-3D, to ﬁll the gap in existing long-term global satellite soil moisture records. show that the largest value (4.13 K) falls in the horizontal polarization of the X-band and Ka-band, demonstrating that the two satellites can be relied on for consistent long-term daily time series of SM and LST dynamics.


Introduction
The introduction of soil moisture (SM), also known as soil water content, is the per unit volume expressed as the dimensionless ratio of soil and water, which forms only a fraction of the world's freshwater resources [1,2]. SM plays an indispensable role in the global water cycle, including hydrological processes such as precipitation, runoff, infiltration, evapotranspiration [3][4][5], and energy and biogeochemical cycles [6,7]. SM is an carrying an active L-band radar sensor and a passive L-band microwave radiometer sensor for observing SM at higher resolutions than previous microwave satellite observations [4].
In the last decade, passive microwave observations from China's FY satellites have shown increasing importance in the long-term SM records for global and regional applications [40]. The first sensor of this series was the Microwave Radiation Imager (MWRI) onboard the FY-3B satellite, which had a local equatorial crossing time of 01:40 for its descending local time and 13:40 for its ascending local time. The FY-3B MWRI brightness observations have been used to estimate SM, LST, and VOD since its operation in July 2011, and has reliably demonstrated its skill to capture SM dynamics across the global and regional scales [38,[41][42][43][44]. Additionally, it has also been successfully used to develop improvement schemes for both SM [34] and LST retrievals [10]. Following FY-3B, the FY-3C satellite was launched in September 2013 with a MWRI onboard. However, unlike FY-3B, the FY-3C crosses the equator at 10:15 and 22:15 for descending and ascending local time, respectively. This provides unique microwave observations which are not available from other sensors. Thus far, very few studies have focused on the scientific and operational benefits of passive observations around this time, such as with FY-3C.
Due to the strong diurnal cycle of LST variability, it is essential to obtain its observations at different overpass times and properly understand its diurnal dynamics [45,46]. LST is fundamentally nonlinearly linked to soil moisture in the relative algorithms [47,48], and satellite observations at such an overpass time would significantly contribute to understanding land parameter changes. The most recent MWRI observations come from the FY-3D satellite, which became operational in 2020 to extend the legacy of its predecessor [49,50] of FY-3B. In addition to having a similar overpass time, observations from these two sensors also have an overlapping period coverage, which ensures continuity. In developing long-term soil moisture records, satellite observations from these three sensors become essential to the existing ones. This study investigates the potential of global FY-3C SM to understand its added value to long-term SM data. In this study, we used the triple collocation analysis (TCA) method to optimize LST retrievals from the FY-3C TB, which were then used as inputs in the LPRM to obtain the global SM for both ascending and descending paths. Next, the FY-3C SM retrievals were compared with FY-3B, which already has an established benchmark LST retrieval model [47]. Finally, their strengths were merged into one framework to assess their long-term global SM development contributions.
The structure of the paper is as follows: Section 2 presents the datasets used, and introduces the FY-3C SM retrieval scheme and the data merging method. In Section 3, the SM of FY-3C and FY-3B are analyzed and evaluated relative to ground observations and a reanalysis product; their relative contributions in the merged daily SM records are assessed within a merging scheme. Section 4 presents the discussions and conclusions.  [45], and the X-band (10.65 GHz) vertical and horizontal polarization TB is used to retrieve SM. The similarities that exist in the FY-3B and FY-3C, in terms of instrument specification and errors and their spatial representations, allow us to isolate and focus more on the impact of their observation times instead. In addition, the differences between these two serve as the backdrop for assessing the added value of FY-3C here, since they are both similarly impacted by systematic differences with the ground and model datasets.

ERA5 LST and Soil Moisture
LST and SM data from ERA5 reanalysis are used. ERA5 is the fifth-generation atmospheric reanalysis product of the European Centre for Medium-Range Weather Forecast (ECMWF) [53]. Compared with its predecessor, ERA-Interim (ERAI), the improvements within the ERA5 framework include several newly reprocessed observations and high spatial and temporal resolutions [54]. ERA5 assimilation uses a collective 4DVar data assimilation scheme with a model spatial resolution of 31 km × 31 km, and a 6-hourly temporal resolution [55]. The skin temperature is the theoretical temperature required to satisfy the surface energy balance. It represents the temperature of the uppermost surface layer. In this study, we used the skin temperature to represent the LST. Here, ERA5 skin temperature products are obtained at 10:00 and 22:00 from 2014 to 2016, which comes at a gridded spatial resolution of 0.25 • × 0.25 • . Evaluation studies have demonstrated the reliability of the ERA5 skin temperature [56].
ERA5-Land is a reanalysis dataset providing consistent land variables over several decades at an enhanced resolution compared to ERA5. It was produced by replaying the land component of ERA5 s global atmospheric reanalysis framework [57]. The ERA5-Land SM products come with four soil layer depths (0-7 cm, 7-28 cm, 28-100 cm, and 100-289 cm) at hourly intervals. This study uses the ERA5-Land volumetric soil water layer 1 (0-7 cm), resampled at 0.25 • × 0.25 • resolution. Several studies have shown that the ERA5-Land SM can capture soil moisture dynamics from ground observations [58].

Noah Soil Moisture
The Global land data assimilation system (GLDAS) data are an assimilation product based on satellite, land surface model, and in situ data jointly released by NASA's Goddard Space Flight Center (GSFC) and NOAA's National Center for environmental prediction (NCEP) [59]. These high-quality global land surface datasets are widely used in weather and climate forecasting, water cycle research, and water resources application. Compared with other models, GLDAS Noah2.1 has the advantages of a stable driving field, advanced model, and long time series. In this study, the soil moisture data from 2014 to 2016 were selected and processed by unit conversion to match the satellite datasets. Both global and regional evaluation studies have also shown that Noah2.1 significantly improves over its older versions [43,58,60].

Satellite and In Situ Soil Moisture
In 2010, the European Space Agency (ESA) listed soil moisture as an essential climate variable (ECV) under the Climate Change Initiative (CCI). Under this project, multiple satellite microwave SM observations were merged into one framework to develop a longterm global SM product to research the impact of climate change on the global water cycle (simply referred to as CCI SM). As a result, three merged CCI SM products have been developed: one from passive microwave observations, one from active microwave observations, and another where the active and passive are combined (the combined) [61]. This study uses the active CCI SM as an input in TCA [62], which comes at a spatial resolution, which is 0.25 • × 0.25 • within the selected period of 2014-2016.
Apart from the CCI SM, we have also used the FY-3B LPRM SM retrievals in our study. Parinussa et al. [51] used the LPRM to retrieve global SM data from the MWRI onboard and have since undergone systematic evaluations [38,40] and improvements [34]. Parinussa et al. [38] compared the quality of SM data retrieved from the fifth version of the LPRM algorithm with that of the sixth version in China. The results showed that the performance of the updated algorithm was improved in the whole study area. These studies have also demonstrated the skill of FY-3B to capture reliable SM variability both globally and regionally, even under complex conditions. This study uses global SM data retrieved from the sixth version of the LPRM algorithm [28,38].
Ideally, the go-to choice for all applicability of soil moisture would be ground observations. Ground observations are generally used to evaluate alternate gridded SM products at regional to global applications. However, a significant challenge to their use is their insufficient availability, and their point-scale representation does not allow for use beyond very localized (spatial-wise) applications. The International Soil Moisture Network (ISMN) is the largest surface in situ SM sharing database globally. It was established by the Vienna University of Science and Technology in 2009. The SM is collected from different regions and climate conditions across the globe and uses quality control procedures to filter and label the dataset [60,63]. The point-scale representation of the station datasets and the area-averaged gridded satellite observations introduces mismatched uncertainties in the comparisons [64]. As a result, we prioritize the comparisons from the correlation analysis and RMSEs based on the anomalies which are less impacted by systematic differences between the datasets. In this study, we used in situ data at 02:00 and 14:00, suitable for FY-3B satellite local time, and in situ data at 10:00 and 22:00, suitable for FY-3C satellite local time for 2014 to 2015. An overview of the networks used is presented in Table 1. Previous studies have noted that the quality of SM datasets from both models [40] and satellites [47] varies as a vegetation density function. Here, the Normalized Difference Vegetation Index (NDVI) obtained from the Moderate Resolution Imaging Spectro-radiometer (MODIS) is used to quantify the impact of vegetation density on the skill of SM retrievals. We used monthly NDVI images (MOD13C) aggregated from its native 0.05 • × 0.05 • grid into a global 0.25 • × 0.25 • grid to compare against the root mean square error from a verification technique and the correlations with the reference dataset.

Land Parameter Retrieval Model
The land parameter retrieval model (LPRM), based on a radiative transfer equation and the vegetation optical depth derived from the Microwave Polarization Difference Index (MPDI) [75], is one of the most widely used soil moisture retrieval algorithms for a range of microwave frequencies, which includes C, X, and Ka-bands [18]. It simultaneously retrieves SM and vegetation optical depth for each location. A unique advantage of the LPRM is the minimized use of ancillary external data such as model-based LST, SM ground observations, or vegetation information inputs commonly used in other retrieval methods, thus reducing the external errors [17,47]. In this study, a recently updated version of the LPRM, e.g., LPRMv06 [28], is used to retrieve SM from FY-3C TB.

Triple Collocation Analysis
An objective and quantitative assessment of the accuracy of various remote sensing products can make them a reliable source of information for important scientific research in earth system science and climate change and expand their application scope [76]. Therefore, it is important to evaluate the skill of retrieved satellite SM, ideally over entire areal extents. The TCA estimates three linearly related mutually uncorrelated sources' relative root mean square error patterns [77]. TCA was first used to evaluate wind and wave height in oceanography [78], and later used to estimate error variance of remote sensing SM data. Since then, several studies have applied TCA to estimate relative errors for various SM estimates [79][80][81][82].
A brief description of the TCA is as follows: three estimates, θ FY3C , θ CCI , and θ ERA5_land , which represent SM data from FY-3C, CCI active, and ERA5_land, respectively. According to Scipal et al. [83], the error variance formula of SM is as follows: , σ 2 CCI , and σ 2 FY3C represent the error variance of ERA5, CCI active, and FY-3C data, respectively, and <·> represents the time series mean, based on at least 100 overlapping observation points.

The Linear Data Combination Technique
The merging step in the study is based on the approach proposed by Kim et al. [84] to linearly combine two satellite SM retrievals by maximizing their correlations to a reference. Recently, it was used to develop a long-term SM product using six satellite SM products [40]. For the two sets of soil moisture data θ 1 and θ 2 in the same time window, the weighted factor w (0 to 1) is applied to linearly combine them into the fusion data θ c at each point: The optimal weighting factor w is calculated as follows: where R is the temporal correlation coefficient between θ 1 , θ 2 , and reference data (Noah2.1). σ 1 and σ 2 are the standard errors of SM data θ 1 and θ 2 . According to the definition of R and Equation (3), the temporal correlation coefficient (R) between θ c and reference data (θ ref ) can be expressed as a function of w. The following is the optimization of the function: where µ C and µ ref represent the average value of fusion data and reference data, respectively, and σ C and σ ref represent the standard error of θ C and θ ref , respectively. It is important to note that since the SM values obtained by different satellites have different ranges and systematic differences, a normalization method is necessary to unify the range of SM values.
Here, the normalization approach used is as follows: where θ p represents the original SM data and σ p and σ ref represent the standard errors of θ p and θ ref , respectively.
The merging approach is separately applied to the ascending and descending observations of the SM retrievals from both satellites to assess how they complementarily depict soil moisture daily variability in one framework. The Noah2.1 SM is chosen as the reference because it is free from interdependency with the FY-3B and FY-3C products aside from its reported high skill to capture SM dynamics [59].
To do this, the two satellite retrievals (parent products) are first normalized to the reference to reduce systematic differences between them. From here, the pixel-wise optimal weights are obtained for pixels where more than 100 days of observations are available. Regions that fail to meet this threshold of available days are masked out. Finally, the parent products are merged to obtain combined ascending/descending FY-3 SM products.

Optimizing Solution
LST plays a unique role in the radiative transfer equation and influences the final quality of soil moisture anomalies. In the LPRM and other widely used land models, the relationship between SM and LST is fundamentally nonlinear [47]. As a result, optimizing the LST inputs from the Ka-band observations is necessary to obtain optimal SM solutions. Passive microwave observations can be an alternative to measure LST using the Ka-band (36.50 GHz) vertical polarization, which balances a reduced sensitivity to soil surface characteristics with a relatively high atmospheric transmissivity. While there are benchmark LST retrieval equations for FY-3B Ka-bands at 1:30 a.m. and p.m., there are none for the overpass time of FY-3C observations. Therefore, it is necessary to obtain a benchmark LST retrieval equation for the FY-3C Ka-band, which will be used as an input routine to retrieve soil moisture. Holmes et al. [45] showed a linear relationship between Ka-band brightness temperatures and LST measurements. Therefore, LST retrievals from FY-3C are obtained by routinely linking its TB with LST estimates. Following Holmes et al. [45] and Parinussa et al. [47], the study systematically optimizes a retrieval approach.
Ideally, LST ground observations would be the first choice to calibrate such an equation; however, the lack of a consistent global spatial coverage of LST ground observations makes this difficult. Parinussa et al. [47] showed that reanalysis LST products are potentially useful due to their high quality and spatiotemporal consistency. Thus, we rely on ERA5 skin temperature to obtain LST from Ka-bands of both the descending and ascending paths of FY-3C. After an initial descending LST retrieval formula (un-calibrated) is obtained, additive biases from −10 K to −1 K (selection of the range for the additive biases are informed by an initial analysis) are added to it at an increment of 1 K. Based on preliminary analysis, additive biases from −5 K to 5 K are also added to the un-calibrated ascending LST at 1 K increments.
In the following step, we used CCI active and ERA5-Land to evaluate the quality of the retrieved SM data under different bias scenarios based on the RMSE (TCA) method. For each LST bias scenario, we used the LPRM to retrieve SM. We then optimized the slope and offset of the LST retrieval formula by selecting the best scenario (lowest RMSE) for each pixel. Finally, by using the optimized LST as input and evaluating with the in situ data from ISMN we obtained an optimized FY-3C SM data (calibrated). Figure 1 describes the multi-step optimization process with descending data as an example.

Comparison and Analysis of the FY-3B/3C Brightness Temperature
In order to develop continuous and consistent SM datasets from the FY-3B and FY-3C satellite observations, it is for their TB observations to be compared to determine whether a mutual calibration is needed. Parinussa et al. [85] showed the necessity of such mutual calibration, which has led to the consistency between SM retrieved from WindSat and AMSR-E. Here, the X-band data, mainly needed for SM retrieval, and Ka-band data, for LST retrieval, are selected and analyzed using the range of 60 • S-60 • N ( Figure 2). Overall, there is a very good consistency between the two satellite observations. Nonetheless, it can be seen from the figure that the consistency of the descending TB data of the two sensors is slightly higher than that of the ascending data. For instance, in the case of Ka-band horizontal polarization, the coefficient of determination (R2) of TB descending path is 0.98, while the R2 of the ascending path is 0.96. Additionally, the R2 of X-band data is slightly higher than Ka-band data. Earlier studies have shown that the open water signals cause the scattered spread at the low ends of the Ka-band, with larger variation in observations close to midday, as shown by the FY3B ascending data in Figure 2 [10,86]. The RMSD calculation results (here, differences (RMSD) instead of errors (RMSE) are used because the comparison is just a relative difference between the two observations; RMSEs are used here in the case of TCA and comparisons with in situ, which are taken as ground truth) between the FY-3B and FY-3C TB data shown in Figure 2 show that the largest value (4.13 K) falls in the horizontal polarization of the X-band and Ka-band, demonstrating that the two satellites can be relied on for consistent long-term daily time series of SM and LST dynamics.  Figure 3 presents the optimization for the LST bias scenario with the lowest SM bias. We applied TCA to all the SM scenarios to evaluate their relative qualities and obtain an optimal solution. A pixel-wise correlation analysis was also performed with each bias scenario and ERA5-Land SM as an added verification step. Figure 3a,c represent the lowest RMSE of SM based on different LST bias scenarios. Figure 3b,d represent the highest correlation of SM based on different LST bias scenarios. The bias scenarios of the descending LST are from −10 K (blue) to 0 K (gray), and the bias scenarios of the ascending LST are from −5 K (blue) to 5 K (red). Thus, the results in Figure 3 show which LST bias scenario per pixel produces the highest quality of soil moisture for the pixel. Our global optimization results show low-bias scenarios around the equator, whereas higher latitudes have higher-bias scenarios. Figure 3a shows the spatial distribution map of the bias scenario with the least TCA RMSE for SM in the FY-3C descending observations. The results show that the low-bias scenarios (low or high biases here refer to their absolute magnitudes) are generally found at locations at or near the equator. In contrast, high-bias scenarios are found near higher latitudes. These findings are consistent with the previous studies where SM retrievals from LPRM were found to have large biases in high latitude regions due to the region's wet climate [87]. Hagan et al. [10] later capitalized on this to reduce LST biases by considering open water fraction in the grids. It is obvious from the results that no bias scenario presents an optimal solution for accurate SM estimates. The results also show that the LST scenario with the lowest bias does not necessarily produce the best quality of SM, which confirms the nonlinear relationship between SM and LST in the radiative transfer equation. Figure 3b shows the SM bias scenarios with the highest correlation coefficient with ERA5-Land SM. The results agree with the TCA RMSE in Figure 3a, which indicates that RMSE (TCA) is feasible as the basis for the optimization. A similar optimization routine was carried out with the ascending observations; however, optimal solutions were found within a bias scenario range from −5 K to 5 K. These significant differences between the bias ranges of the two acquisition times have been shown and elaborated by Lei et al. [88], who attributed the difference to how evapotranspiration is represented at the time of data acquisition. Figure 3c shows the optimal bias scenarios based on TCA, while Figure 3d shows the optimal solutions based on correlations with ERA5-Land. As noted with the descending results, the optimal solution does not lie with the 0 K scenarios or one particular scenario but within a spectrum of scenarios. Positive biases were found over the high latitudes and in the tropical regions. In contrast, the negative biases were found in the mid-latitudes, where the TCA and correlation-based optimal maps are not as high as the descending orbits results, suggesting that the accuracy for optimal solutions from descending orbits would be more significant than that from the ascending orbits.

Global Verification
As a final step, the LST of the optimal SM scenarios were used to obtain calibrated retrieval equations for both the ascending and descending TB of FY-3C by determining their optimal slope and offsets in the equations below. Equation (6) is the obtained calibrated LST retrieval equation for ascending TB, and Equation (7) These equations are, thus, the proposed benchmark equations for retrieving LST from FY-3C. To quantify the improvements from this calibration process, TCA is applied to the SM retrievals from both the calibrated and un-calibrated (Figure 4a,b). Additionally, the correlation between the anomalies of the two SM retrievals and the ERA5-Land SM anomalies are also estimated (Figure 4c,d). Beyond these NDVI ranges, satellite SM loses skill and sensitivity to rainfall with a low signal-to-noise ratio in bone-dry regions. Figure 4 shows that more significant improvements are realized from the descending products than ascending products through the calibration process. The evolution of the TCA RMSE results is consistent with the skill of passive microwave-based SM retrievals reported in previous studies. The SM observations lose skill with increasing vegetation density, as seen in the descending (Figure 4a) and ascending data (Figure 4b). This difference in the quality of the SM retrieval after calibration could be attributed to an assumption in LPRM, where the vegetation temperature is considered equivalent to the bare LST [45,88]. Here, the assumption works best when there is increased thermal equilibrium such that the temperature at the top of the canopy is equivalent to the temperature at the surface of the land. In the daytime, the transpiration process affects the canopy temperature, leading to the cooling of the canopy. However, the LST generally warms up [89]. Therefore, different LST inputs of different observation times will affect the retrieval quality. Further improvement of daytime observations will lead to significant progress, while nighttime observation quality may not be obvious. As a result, optimizing the un-calibrated retrieval realizes more obvious improvements in the descending observations (taken at 10:00) than ascending observations (taken at 22:00). The correlation analysis results for descending ( Figure 4c) and ascending data (Figure 4d) also demonstrate that the calibrated improvement is more significant in the descending observations. The lower correlations in NDVI < 0.4 regions result from the reduced sensitivities in the satellite observations due to insufficient precipitation in the regions [28,47]. Nonetheless, higher correlations are found over regions of low vegetation densities in the ascending product and reduced skill with increasing vegetation density (NDVI > 0.4) for both products, thus validating the assumption's impact, as discussed above.
Next, an independent evaluation of the calibrated datasets with the ISMN in situ SM observations is also carried out. Here, the skills of the calibrated ascending and descending FY-3C SM retrievals are assessed to understand their strengths and weaknesses. The relevant information of the selected SM networks is described in detail in Table 1. Figure 5 shows the RMSE and correlation analysis of the two satellite products with the in situ datasets at their corresponding times. The results are presented for each network, which comprises multiple stations. The RMSE values are mostly concentrated between 0.05-0.15 m 3 m −3 , and the correlation coefficient values are mostly within 0.4-0.7. Generally, the ascending mean errors (ranging between 0.03-0.05 m 3 m −3 ) are lower than the descending mean errors (ranging between 0.04-0.05 m 3 m −3 ) for all the networks except SMOSMANIA. Additionally, the mean correlations are also higher for the ascending than the descending datasets. As mentioned above, this may be because ascending observations (nighttime) are less affected by vegetation transpiration than descending orbital observations (daytime); therefore, the data from ascending observations (nighttime) are more accurate. Figure 5 shows that the AMMA-CATCH SM in Africa shows the lowest errors with high correlations for both the ascending and descending data. Because this SM network is in the arid region with a low vegetation coverage, it has less influence on the SM retrieval. Given that this region is significantly challenged by lack of observation data; it behooves both scientific research and application of these datasets to have such good qualities. On the other hand, the lowest performance of the datasets is found with the SMOSMANIA network located in the south of France, which has the highest errors and lowest correlations. This SM network is located near open water, affecting the retrieval quality. Additionally, high vegetation coverage in the region also contributes to the reduced qualities observed. Nonetheless, the results show that the calibrated FY-3C SM datasets are reliable on the whole, having relatively low errors (Figure 5a) and good temporal correlations (Figure 5b). Furthermore, the highly varying qualities between the networks also demonstrate a spatial dependency. This is explored in detail in the next section. ; (c) correlation for descending data; (d) correlation for ascending data across different NDVI scenarios. NDVI < 0.1 and NDVI > 0.85 are not included because there is no sensitivity in the satellite product due to low signal-to-noise ratio and soil moisture losing skill, respectively.

Inter-Comparison of FY-3B and FY-3C Products
As mentioned above, the accuracy of SM retrievals has been found to vary as a function of vegetation distribution. Here, the SM anomalies of FY-3C are inter-compared with that of FY-3B to understand their consistency over different NDVI distributions based on the in situ observations, as shown in Figure 6. The results show that vegetation density impacts the errors of the ascending and descending SM estimates of both datasets. Figure 4a,b and Figure 6a,b also show that low vegetation areas have the lowest errors, while dense vegetation areas have the largest errors. Furthermore, most errors are found above the 1:1 line, indicating slightly larger in FY-3C than in FY-3B. In Figure 6c,d, on the other hand, the correlations are almost evenly generally scattered around the 1:1 line, showing that their temporal consistencies are very similar. In addition, higher correlations are found in the ascending datasets of FY-3C, since a higher density of the results are found above the 1:1 line in Figure 6d. Higher correlations are found in the dense vegetation areas within the descending FY-3C datasets, and higher correlations are also found in the low vegetation areas within the ascending FY-3C datasets. Generally, we expect that because of the observation time of FY-3C, we would find lower qualities in the dense vegetation areas than FY-3B due to higher thermal equilibrium in FY-3B at night. However, FY-3C appears to have higher correlations (Figure 6c), although larger RMSEs are also found there ( Figure 6a). As noted in the optimal maps from Figure 3a,b, the optimization would also improve the correlations of FY-3C, especially in the dense vegetation regions where thermal equilibrium is generally reduced due to transpiration. Higher correlations in the low vegetation areas of the ascending observations of FY-3C in Figure 6d are consistent with previous studies that also found higher correlations in nighttime observations [47]. Overall, Figure 6 distinctly demonstrates a very high consistency between the two satellite retrievals for both orbits, which is essential for developing long-term records from both observations. In aiming towards including the SM retrievals from FY-3B and FY-3C in long-term records, it is important to investigate how their strengths can be leveraged to describe daily SM dynamics. Additionally, it is helpful to identify regions where either retrieval can be relied on more to explain daily soil moisture variability. To do this, TCA is applied to the SM anomalies from both satellites observations and the regional differences are shown globally for both the ascending and descending products (Figure 7a,b). Here, red-colored regions imply lower errors in FY-3B, and blue-colored regions imply lower errors in FY-3C. As a routine in LPRM, very dense vegetation regions are masked out. Overall, the results show that the skill of FY-3B in capturing daily SM dynamics is better in descending products, while FY-3C is better in the ascending product. A strong contrast can be seen from Figure 7a,b, especially from the mid-latitudes in the northern hemisphere (NH) to the southern hemisphere (SH) that in regions where FY-3B has higher accuracies in descending products, the FY-3C picks that up in the ascending products. It is worth noting that the FY-3B is shown to dominate the high latitudes, eastern North America, and Australia in both the ascending and descending products from the TCA results. A similar accuracy assessment is done based on correlation analysis, as shown in Figure 7c,d. The results show a strong agreement with the TCA results. Higher accuracies are reported in FY-3C over high latitudes over the NH in the descending products, with higher accuracies realized in NH mid-latitudes all through the SH (Figure 7c). The opposite distribution of the accuracies between the two satellite observations is seen in the ascending products (Figure 7c). These results highlight how the two observation sources complement each other in capturing daily SM dynamics.

Evaluation of Combination Results
To demonstrate the added value of the combined product of merging the FY-3B and FY-3C ascending/descending observations, it is evaluated with in situ SM and compared with the parent products (FY-3B and FY-3C). In this study, the in situ data from 2014 to 2015 are selected for verification. The latitude and longitude positions of the chosen grid in the Figure 8 are 118 • W, 42 • N, and the corresponding SM site name is SNOTEL_DISATERPEAK. This site is located in the Midwest region of the United States. The timeseries of this site is more complete, which is easy to compare with the retrieval and combined data. Figure 8a,b shows the time series of a chosen pixel for both the descending and ascending products, respectively. The results demonstrate that the merged product has accurately captured the variability of in situ soil moisture. A clear increased performance is seen in the combined product compared to the parent products. Apart from capturing the seasonal variability, significant precipitation events seen in the in situ data are also well captured in the combined product, which verifies the validity of the parent products as well.
As a final verification step, all matching pixels in the three products: two parents and the combined, to the in situ observations, are compared based as a function of vegetation based on NDVI. This final evaluation aims to understand the relative strengths of the parent products that are leveraged in the combined product. As was done in the above analysis, NDVI > 0.1 to NDVI < 0.75 are used, since satellite SM loses skill beyond these boundaries [47]. In Figure 9, the blue and red lines represent the FY-3B and FY-3C SM products, respectively, while the black line presents the combined product. As was reported by Hagan et al. [40], improvements in the mean differences, here quantified with RMSEs with the in situ observations, are attributed to the normalization step. Figure 9a,b both show that the RMSEs are reduced in the combined products over all vegetation density scenarios. This partly shows the biases within the parent products attributed to a limitation of LPRM. Future studies could incorporate improvement schemes, such as improving the LST inputs as proposed by Hagan et al. [10] or applying a dynamic vegetation correction as suggested by Parinussa et al. [34], to reduce biases and errors within the retrievals from LPRM. Figure 9c,d show the added value of the merging process to maximize the correlations in the combined product. Significant improvements are realized in the combined product over most vegetation density scenarios, especially between NDVI > 0.2 and NDVI < 0.7. This attests to the importance of a co-usage of these products for both scientific research and applications. Both FY-3B and FY-3C clearly have individual merits that, when combined, have very promising uses. More importantly, the combined usage is essential to developing long-term soil moisture records. The results in Figure 9 also show that while FY-3B tends to have lower RMSE, FY-3C tends to have higher correlations in both the ascending and descending products. This is very important to note, especially when it comes to studying physical processes within the climate. Some processes would benefit from accurate means, and others would benefit from accurate temporal variability. Furthermore, this sheds some light on a possible benefit of SM observations at these overpass times of the two satellite observations.

Conclusions
SM can regulate the exchange of water vapor and energy between the land and the atmosphere by changing surface fluxes. Therefore, SM data can help predict extreme events, such as droughts, floods, and heatwaves [3,13]. Microwave remote sensing technology is one of the most effective ways to monitor SM dynamics on a regional and global scale [90]. In this paper, we used China Meteorological Administration's FY-3C TB to solve for an optimal solution of SM by retrieving the LST input used in LPRM. The obtained soil moisture product was routinely evaluated using both global in situ and high-quality reanalysis SM products to understand its relative qualities. We also examined the complementarity of FY-3C and FY-3B satellite SM observations for obtaining consistent long-term SM records. Our conclusions are as follows.
The evaluation of the TB from the MWRI onboard FY-3B and FY-3C under the different polarization modes of X-and Ka-band suggests that the TB values of the two satellite platforms are consistent with correlation coefficients above 0.95. The consistency of TB data in the descending orbits is slightly higher than that in the ascending orbit. The X-band data showed a higher consistency than the Ka-band data in our comparison. Our overall results demonstrate the potential usefulness of the combined applications of FY-3B and FY-3C satellite observations. Unlike FY-3B, FY-3C comes at a different observation time of the day, implying the need to tune the SM retrieval algorithm, the LPRM, to fit the observation time of FY-3C. We adopted the LST retrieval approach proposed by Holmes et al. [45] to obtain the uncalibrated LST retrieval equations. However, because of the nonlinear linkage between LST and SM, optimizing the LST inputs was essential to obtain an optimal solution for SM. Following Parinussa et al. [47], the algorithm was optimized using the TCA method for both ascending and descending datasets of FY-3C. Comparison of the correlation coefficients and the RMSE (TCA) of un-calibrated and calibrated results showed obvious improvements in the descending (daytime) data, and a slight improvement in the ascending (nighttime) data, possibly due to increased thermal equilibrium of the nighttime SM and canopy temperature. Further evaluations of the improved FY-3C SM anomalies using the ISMN in situ data, we found the RMSE values were concentrated between 0.05-0.15 m 3 m −3 and the correlation coefficient between 0.4-0.7, although better skill was found in the ascending SM product. Additionally, the evaluations also showed that, similar to the FY-3B, the FY-3C SM improves skill with decreasing vegetation density.
In a final step, we aimed to unravel the complementarity of FY-3B and FY-3C SM by comparing their individual and combined potentials [84] relative to in situ observations. When we combined the daytime and nighttime SM products separately, the results showed superior skill to capture SM dynamics than the individual products, since they leveraged the best of each product. Each product had its unique advantage over space-time and time, demonstrating the FY-3C series can add value to existing SM records similar to the CCI.
This study has found that the FY-3C is useful in filling gaps in the global observational records of SM, LST, and VOD, for which there are little to no observations available. Furthermore, the FY-3C provides observations at a unique time that was initially unavailable, which makes it useful to understand land surface processes. Future studies are needed to extend this framework to the recently launched FY-3D satellite for global soil moisture observation to serve as a successor of the FY-3B and FY-3C satellites.
Supplementary Materials: The following supporting information can be downloaded at: https: //www.mdpi.com/article/10.3390/rs14051225/s1. Figure S1. A flowchart that sequentially describes the multi-step optimization procedure for FY-3C ascending LST. This procedure starts with the uncalibrated formula, followed by the additive bias scenarios used within the LPRM to obtain the calibrated LST equation.

Acknowledgments:
We are grateful to the teams at the National Aeronautics Space Administration, the European Space Agency, the European Center for Medium-Range Weather Forecast, TU Vienna, and the National Satellite Meteorological Centre of China for making their datasets available. We are also grateful to Robert Parinussa for his helpful suggestions throughout this study, as well as the anonymous reviewers for their helpful comments.

Conflicts of Interest:
The authors declare that they have no known commercial or financial competing interest that represents a conflict of interest in this paper.