Temperature Effects in AMSR2 Soil Moisture Products and Development of a Removal Method Using Data at Ascending and Descending Overpasses

: Soil moisture is among the most essential variables in hydrology and earth science. Many satellite missions, such as AMSR-E/2, have been launched to observe it in broader spatial coverage to overcome the shortage of in situ observations. However, the satellite soil moisture products have been reported to comprise errors caused by the so-called “temperature effects” widely observed in dielectrically measured in situ volumetric soil water content (SWC). In this work, we confirmed the existence of these errors in AMSR2 soil moisture products. A new algorithm was developed to remove these errors using satellite data at ascending and descending overpasses. The application of this algorithm to both satellite and in situ data of SWC and soil temperature at the Mongolia site shows that the difference between SWC values at ascending and descending overpasses caused by temperature effects is effectively removed. We assess the impact of this removal method on satellite data by comparing it with in situ data, utilizing metrics such as the correlation coefficient and other widely adopted evaluation methods. It is shown that the difference between the original and corrected in situ SWC is much smaller than that between AMSR2 and in situ SWC, either corrected or not. The results indicate that the metric values between the corrected AMSR2 and in situ SWC, after removing apparent differences caused by temperature effects, slightly improved compared to those between the original AMSR2 and in situ SWC. Though these findings imply that the removed errors may not be the most dominant, considering the current significant difference between AMSR2 and in situ SWC, the removal makes the ascending and descending data have close characteristics. It may allow using data at both ascending and descending overpasses and double the temporal resolution of AMSR2 SWC data.


Introduction
Though soil moisture is a very tiny fraction of the total water budget (accounting for about 0.005% of whole water storage and 0.15% of the freshwater in the world), it is one of the most crucial variables for studies in the water and energy exchanges between the Earth's surface and atmosphere [1][2][3][4].Due to the high interaction with the atmosphere and the ground surface, soil moisture plays a significant role in developing weather patterns, including heat fluxes [5] and precipitation [6], and also in predicting extreme climate events such as drought [7,8].Soil moisture is like a bridge between the air and the solid earth, in which water can return to the atmosphere via evaporation, infiltrate to the lower soil layers, percolate, and become groundwater or join to baseflow [9,10].Furthermore, the energy exchange across the Earth's surface is concurrent with these processes to create balance in the water cycle.Some studies reported that soil moisture is vital in climate change research [11,12].Many studies have been and are currently being conducted to observe soil moisture [13][14][15][16].
As a numeric metric to express the soil moisture, volumetric soil water content (SWC) is the most widely used.Though the in situ SWC observation is known as the most precise soil moisture data, it suffers a noticeable limitation due to the spatial coverage.Although the number of soil moisture networks is increasing [17,18], more is needed to satisfy the demand.Meanwhile, extreme weather and climate events have become more frequent during the past decade [19,20].There was an urgent need to establish techniques that satisfy soil moisture's temporal and spatial measurements.Hence, the microwave remote sensing technique for soil moisture has been proposed to meet the need [21].With the launching of two dedicated satellite missions, SMOS [22] and SMAP [23], along with the products from other projects such as AMSR-E/2 [24][25][26], numerous opportunities have been opened for global hydrometeorological studies.
Nonetheless, satellite-based soil moisture faces another obstacle related to accuracy.Although many efforts have been made to improve the quality of these estimations, the accuracy in some areas still needs improvement [27].It is widely recognized that satellite soil water content data at ascending and descending overpasses have systematic differences in the values [28][29][30][31][32].According to Hoang and Lu [33], this discrepancy is mainly caused by the difference in soil temperatures at ascending and descending overpasses, the so-called 'temperature effects' (hereafter referred to as TEs).Though these effects were observed in dielectrically measured in situ soil moisture several decades ago [34][35][36] and various methods have been proposed to make corrections [35,37], few studies on satellite soil moisture have been carried out.
It has been reported that the TEs in satellite soil moisture are similar to their impact on in situ data and have the same physical base [28,[37][38][39].Hoang and Lu [33] checked the existence of TEs in SMAP soil moisture product and showed that the temperature dependency of the dielectric permittivity [38,40] can have significant effects on surface brightness temperature from which AMSR2 SWC is retrieved.
During the last few decades, several studies have been initiated to propose the temperature effect removal algorithms.The current algorithms can be sorted mainly into three groups, namely, data-driven approach [37,41], empirical [42,43], and mixture model-based methods [35].However, none of these methods is dedicated to temperature correction in satellite soil moisture data.Though Hoang and Lu [33] also tried to apply Kapilaratne and Lu's method [37,41] to remove the TEs in the SMAP SWC product, these methods require in situ data as a reference or high temporal resolution of SWC and soil temperature to calculate daily amplitudes.Modifications are necessary to apply to satellite data, which are only observed once or twice a day.
Based on Kapilaratne and Lu's concept [37,41], a new method suitable for sunsynchronized satellite observation is developed in this study.This method is designed to remove the TEs for AMSR2 soil moisture using only the data at ascending and descending overpasses (hereafter referred to as ascending and descending data).This method is then evaluated by comparison of the differences between the soil moisture at ascending and descending overpasses before and after this removal.Also, the impact of this removal on AMSR2 soil moisture product is investigated.

The Review of Studies to Remove Temperature Effects
Several studies tried to determine the total TEs on soil moisture, yet the influences of diurnal changes in soil temperature on dielectric measurement of SWC still need to be eliminated.Through a detailed analysis of the data, Lu et al. [37] clarified that the diurnal amplitude of the in situ SWC is proportional to the product of the diurnal amplitude of soil temperature and the daily mean SWC, as shown in Equation (1).
where A θ (m 3 /m 3 ) and A T ( o C) represent the daily amplitude of in situ SWC and soil temperature, respectively; θ m (m 3 /m 3 ) stands for the mean in situ SWC of a day and α (m 3 /m 3 / o C) is the temperature correction coefficient for in situ SWC.In the study by Kapilaratne and Lu [41], it was found that the α value varies significantly across specific regions.This variability can be attributed to the fact that soil moisture demonstrates diverse characteristics influenced by factors such as climate, soil properties, and land coverage, resulting in different dependencies on temperature.Lu et al. [37] also introduced an expansion form for Equation (1) as a general form by assuming the apparent diurnal changes of SWC caused by the changes in soil temperature from the reference temperature.
where T ( o C) and T re f ( o C) are the actual soil temperature and reference temperature at which the sensor calibration curve was created; θ (m 3 /m 3 ) and θ re f (m 3 /m 3 ) stand for measured and actual in situ SWC at T and T re f , respectively; and θ c (m 3 /m 3 ) represents the in situ SWC at temperature T c ( o C), which is between T and T re f .Usually, T re f can be set to 20 • C. To make SWC correction, namely, to calculate the actual in situ SWC, Equation ( 3) is derived by introducing approximation θ c = θ re f .
Lu et al. [37] also confirmed that 20 • C soil temperature change can lead to a 16 percent relative error in the actual in situ SWC at the Mongolia site.In order to avoid the effects of rainfall, Lu et al. [37] excluded the data of the days with significant daily rainfall amount and their subsequent days.This exclusion limits its application to soil moisture sites without rainfall observation.In 2017, Kapilaratne and Lu [41] automated their method to eliminate the requirement of in situ rainfall information.
However, estimating the temperature correction coefficient requires the daily amplitude of SWC and soil temperature, which are often difficult to obtain from satellite data.Even sun-synchronized satellites such as SMAP and AMSR2 only have two observations a day at most.Furthermore, some satellites cannot provide direct observation of the soil temperature.Hoang and Lu [33] confirmed the existence of TEs in SMAP soil moisture product and showed the possibility of making SWC correction using the in situ soil temperature and the temperature correction coefficient α estimated from in situ data.However, this solution, again, faces the limitation regarding the spatial coverage.Considering the differences in observation technology and the microwave frequencies, a new method using satellite data is expected.

The New Method to Remove TEs from AMRS2 Soil Moisture Product
Theoretically, the amplitudes in Equation (1) can be substituted by differences between data at AMSR2 ascending and descending overpasses.However, the difference in AMSR2 SWC may include changes caused by soil temperature and other physical phenomena such as evapotranspiration and soil water redistribution [33,41].Separating these two will be crucial for developing a new TE removal algorithm.
This study investigates the relationship among three consecutive AMSR2 observations: one at descending overpass and its previous and following ascending overpasses.This sequence forms what we refer to as the previous Ascending-Descending-following Ascending relationship (ADA triple hereafter).The ADA triples of soil temperature and SWC are formed by the ratios of them at two adjacent ascending overpasses to their values at descending overpass and are calculated as follows where R expresses ratios, θ (m 3 /m 3 ) and T ( o C) denote SWC and temperature, and subscript X =A p , A f or D represents the previous ascending overpass, the following ascending overpass and the descending overpass, respectively.The mean ratios of all ADA triples calculated from Equation ( 4) are displayed in Figure 1.They form two triangles (ADA triangle hereafter), one for temperature and another for SWC.In Figure 1, R θ and R T stand for the ratios of SWC and temperature of each overpass relative to their values at descending overpass, respectively.A m signifies the average value of two adjacent ascending overpasses (the previous Ascending, Ap, and the following Ascending, Af).Following this, R θ,Am and R T,Am denote the ratios between the A m values and descending data of SWC and temperature, respectively.It is observable that the temperature ADA triangle is almost bilaterally symmetrical.The temperatures at previous and following ascending overpasses are almost identical.The SWC ADA triangle shows a decreasing trend from the previous ascending and the following ascending overpasses.Considering the TEs in SWCs at two ascending overpasses should be almost the same, the difference between SWCs at the previous ascending overpass and the following ascending overpass can be a decrease other than TEs.The average of the SWCs, θ Am , at two ascending overpasses and the SWC at descending overpass may have very close actual SWCs.Their difference is mainly caused by the difference in temperatures at ascending and descending overpasses.Based on the above analysis, we can obtain ∆θ = αθ D,re f ∆T (5) where θ D,re f (m 3 /m 3 ) expresses the actual SWC at the descending overpass; θ (m 3 /m 3 ) and T ( o C) stand for the soil water content and temperature; the subscript A and D mean ascending and descending; small subscripts p, f and m represent the words "previous", "following" and "mean", respectively; for instance, θ Am (m 3 /m 3 ) denotes the mean SWC of two adjacent ascending overpasses.
By slightly changing Kapilaratne and Lu's algorithm [37,41], the temperature correction coefficient α can be estimated using ascending and descending data only.Because θ D,re f is unknown and needs to be calculated from α, θ D and T re f , we estimate the α as follows.
Step 0: set initial θ D,re f = (θ Am + θ D )/2; Step 1: apply the recursive regression method [41] to estimate α; By repeating steps 1 and 2, the α value will converge very quickly.In many cases, the second estimation of α is accurate enough.In step 1, linear regression and outlier removal are carried out recursively.The outlier removal removes the data points outside the confidence interval which depends on pre-assigned confidence level, γ.In this study, it is set to 0.01.Kapilaratne and Lu [41] showed this removal can effectively remove large changes caused by other factors, e.g., rainfall.
Just like the process of removing temperature effects from in situ data, the temperature correction coefficient (α value) varies significantly across distinct regions.This variance arises from the inclusion of the vegetation factor as a crucial component within the Radiative Transfer Model utilized for satellite soil moisture retrieval, primarily depicted by land cover [44].Consequently, the α value necessitates adjustment tailored to the characteristics of individual regions.Although the precise temporal and spatial distribution of α remains unclear, generating a global map of α can be automated by employing the ADA triangle method with AMSR2 soil moisture dataset.

The Tools Using for Evaluating New TE Removal Method on AMSR2 Soil Water Content
We evaluate the impact of temperature effect removal from three aspects.First and foremost, the visual confirmation will be conducted to examine the difference between ascending and descending AMSR2 soil moisture data of before and after the removal.Secondly, the quantitative evaluation of the impact of temperature effect removal on in situ and AMSR2 SWC datasets will be carried out.In this approach, both original and corrected AMSR2 SWC datasets are compared with original and corrected in situ SWC datasets and quantitatively evaluated by using the following metrics: where Y 1 and Y 2 are two variables to be compared; µ 1 , µ 1 , σ 1 and σ 2 are their mean values and standard deviations; t and N t are sequential numbers of time and the total number of the data; ME, MAE, R, RMSE and ubRMSE are mean error, mean absolute error, correlation coefficient, root mean square error and unbiased root mean square error, respectively.They are widely used to evaluate satellite soil moisture products [28,32,33].These metrics are affected by all errors.The large errors are dominant, especially in RMSE and ubRMSE.They are used to evaluate the satellite soil moisture product's total performance compared to the in situ ones.
In order to avoid the effects of significant errors caused by factors other than temperature, the following equation is used to pick up the median value of absolute errors (Med AE) for all data, including the outliers.Last but not least, we will assess the influence of the TE removal basing on the difference between ascending and descending SWC data quantitatively.In this case, the absolute difference (AD) between θ Am and θ D is analyzed as the following equation: To exclude the impacts of the outlier, Equation (13) was applied to pick up the median value for all AD values (MedAD) of Equation (12).

Data Description
In this study, we chose the Mongolia site as a representative location to assess the effectiveness of the ADA triangle method.As a validation site maintained by The Japan Aerospace Exploration Agency (JAXA) science team, its data are widely used for calibrating and validating AMSR2 soil moisture product [30,31].The confirmation and removal of TEs in AMSR2 soil moisture product focused on a 50 × 50 km specific target area (45.75°N, 106.9°E and 46.2°N, 106.25°E), corresponding to an AMSR2 soil moisture single pixel (refer to Figure 2) that is known as the Core Validation Site (CVS) in this study.This target area was within a more extensive study area on the Mongolian Plateau, spanning a region of 1.1°× 1.1° [30].The study area is dominated by the cold, arid climatic (BSk) area (as sorted in the Köppen-Geiger climate classification [45]).According to Yamanaka [46], this area exhibits a fully desert climate characterized by a high Budyko's radiative dryness index (greater than three) and low annual precipitation.Though short grass and small shrubs grow from late spring to mid-autumn, the predominant land cover is bare soil.Land use mainly consists of rangeland, interspersed with sparsely distributed ephemeral ponds and small towns.Soil types are predominantly sandy loam or sandy silt loam with abundant gravel.It is noted that soil-hydraulic properties play a significant role in regulating the soil moisture patterns in the Mongolia CVS area [30].The topography of this site is primarily flat, with an elevation from 1300 m to 1500 m [30].

In Situ Data
Based on data availability and quality of all months at all stations, we chose nine stations from four AWSs and twelve ASSHs for analysis from 5 September 2016 to 31 October 2019.These selected stations have successfully observed the essential elements of meteorology and soil moisture for the Coordinated Enhanced Observing Period (CEOP) and several projects, e.g., the ADEOS-II Mongolian Plateau Experiment (AMPEX) for ground truth and the Mongol AMSR/AMSR-E/ALOS Validation Experiment (MAVEX).During the installation process, particular attention was paid to ensuring that each station represents the surrounding geomorphological and vegetation conditions [30,31].The Global Change Observation Mission (GCOM) website of JAXA provides the in situ database for this study (https://suzaku.eorc.jaxa.jp/GCOM/insitu/index.html,accessed on 31 July 2023).For analysis, we accessed this website and downloaded the necessary data.The details of selected stations at the Mongolia site are displayed in Table 1 and Figure 2. As shown in Table 1, this site comprises two automatic weather stations (AWSs) and seven automatic stations for soil hydrology (ASSHs).The AWSs monitor important meteorological parameters and soil moisture readings at 30 min intervals at depths of 3 cm, 10 cm, and 40 cm for all stations.The ASSHs, on the other hand, measure soil moisture and soil temperature every two hours at depths of 3 and 10 cm.Before installation, all sensors used at AWSs are calibrated and verified against a reference marker and the Japanese Meteorological Agency standard in a laboratory setting [30].Soil moisture sensors are the TDR sensors inserted horizontally.Further details of AWS and ASSH networks can be found in [30].Considering the penetration depth of the microwave, this study used soil water content (θ InSitu ) and temperature (T S ) observed at 3-centimeter depth.We resampled the data of AWSs to generate a bi-hourly dataset.Then, the areal average values over the CVS are calculated from this dataset by using the Thiessen polygon method as follows: where N s is the number of stations; V is the variable to be averaged; V is the average value; and w is the Thiessen weight, and The CVS and Thiessen polygons are shown in Figure 2 and the Thiessen weights are displayed in Table 1.This study will compare the average SWC with the satellite data over the CVS.Over the study period, the longterm averages of V and S V of the SWC are 0.0624 m 3 /m 3 and 0.0160 m 3 /m 3 , respectively.The variability among nine stations and their deviation from the areal average values is illustrated in Figure 3 over a period of five days, showcasing both temporal and spatial fluctuations in the Mongolia CVS dataset.Figure 3a displays the soil water content values, while Figure 3b represents the soil temperature readings.On the left side of Figure 3 are the observed values, while on the right side are their corresponding areal averages and standard deviations.Also, a pseudo-satellite dataset was made by picking up data closest to the ascending and descending overpasses of AMSR2 (13:30 and 1:30 local time, respectively).This is a test bed for developing the algorithm to remove TEs.In many cases, especially the cases comparing with AMSR2 data, this pseudo-satellite dataset is simply referred to as in situ dataset.

Satellite Data
AMSR2 level 2 (L2) soil moisture content (SMC) product was used in this study.AMSR2, part of the Global Change Observation Mission 1-Water (GCOM-W1), was launched by the Japan Aerospace Exploration Agency (JAXA) in May 2012.Starting from August 2012, data from AMSR2 became accessible [24,32].They can be accessed through the website https://gportal.jaxa.jp/gpr/?lang=en(accessed on 30 July 2023).The JAXA Earth Observation Research Center (EORC) provides soil moisture product derived from AMSR2.This product is retrieved from brightness temperature data observed at ascending and descending orbits [29].The ascending and descending overpasses are about 13:30 and 1:30 local time, respectively.The spatial resolution is about 50 kilometers.The soil water content from the JAXA SMC product is referred to as AMSR2 SWC (θ AMSR2 ).
In order to confirm the possibility of removing TEs using satellite data only, the AMSR2 level 2 (L2) land surface temperature (LST) product is also used.This is called AMSR2 surface temperature (T SF ).

The Existence of Temperature Effects
Figure 4a shows bi-hourly in situ soil temperature and SWC, and AMSR2 SWC and AMSR2 surface temperature from 24 June 2017 to 28 June 2017, picked up from the whole study period.The in situ SWC, namely the TDR sensor reading and AMSR2 SWC, shows diurnal changes almost synchronized with the soil temperature.These diurnal changes differ from the ones caused by evaporation [47] because evaporation usually reduces the daytime SWC.Kapilaratne and Lu [48] also showed an opposite diurnal change in SWC to the potential evaporation.This phenomenon implies that TEs exist in both the in situ SWC and AMSR2 SWC and are more dominant than changes caused by evaporation.Figure 4b shows the ADA triangles made from the pseudo-satellite dataset made from the whole in situ dataset, and Figure 4c is the ADA triangles made from AMSR2 SWC and surface temperature.The two ADA triangles of temperature are almost bilaterally symmetrical.The ADA triangle of AMSR2 surface temperature has a more significant difference between ascending and descending overpasses than in situ soil temperature.This trend may reflect the difference between these two temperatures: the in situ temperature is observed at a depth of 3 cm and AMSR2 surface temperature is land surface temperature derived from AMSR2 brightness temperature [49].The two SWC ADA triangles show decreasing trends of SWCs from the previous ascending overpass to the following one.Considering the temperatures at the two adjacent ascending times are almost the same, this declination may not be caused by temperature.Recognizing the difference between the averaged ascending SWC and the descending SWC as TEs is reasonable.TEs in in situ SWC are more significant than those in AMSR2 SWC.Based on these observations, applying the removal algorithm to these datasets is highly possible.

Results of the Removal Algorithm
Figure 5 shows the results of the removal algorithm applied to the in situ dataset and AMSR2 dataset.Both are the results of the second round application of the recursive regression.The blue dots are original data; the red circles are data points outside the upper and lower bounds of the 99% confidence interval and are excluded as outliers in the recursive regression [41].Among many possible reasons for these outliers, rainfall may be the most obvious one.A meaningful rainfall before the descending overpass may increase θ D and make the data below the confidence interval's lower bound.On the other hand, the rainfall before the following ascending overpass will increase θ Am and make the data lie above the average of the upper bound.The temperature effect coefficients are 0.0060 and 0.0054, respectively.As mentioned in the previous subsection, the TEs in in situ data are more significant than those in AMSR2.These values mean that a 20 • C temperature change may cause 12.0 percent and 10.8 percent apparent relative error in in situ and AMSR2 SWC, respectively.Figure 6 displays part of the time series of the in situ SWC and AMSR2 SWC (θ InSitu and θ AMSR2 ) before and after correction (indicated by a superscript prime), and ADA triangles of in situ data and AMSR2 data without outliers.In Figure 6a, the SWCs at ascending overpass from 24 July 2023 to 28 July 2017 become closer to that at descending overpass after removing the TEs.The green lines in Figure 6b,c show the ADA triangles of the in situ SWC and the AMSR2 SWC, respectively, after applying the TE removal.It is observable that the TEs in the original ADA triangles are reduced, revealing that the distances between θ Am and θ D have become smaller.Moreover, the differences between the previous and following ascending overpasses are nearly unchanged.This result means the new removal algorithm can remove TEs and conserve the SWC declination caused by other factors.) datasets are derived.These datasets will undergo analysis using five metrics outlined in Section 2.3 to examine the relationship between four pairs as follows: 1.
Corrected AMSR2 SWC (θ ′ AMSR2 ) and corrected in situ SWC (θ ′ InSitu ).Table 2 displays the values of the quantitative metrics of averaged data of the ascending data and descending overpasses in the following four pairs of the original and corrected AMSR2 and in situ SWC products from September 2016 to October 2019.The first pair is used to check the impact of TE removal, namely data correction on in situ data and the ground truth in satellite SWC evaluation.The metric values showed that this pair has a much smaller difference than the other three pairs of AMSR2 and in situ SWC data.This result implies that the direct impacts of TE removal on soil moisture products will be limited.The other three pairs give metric values similar to other evaluation studies (for example, [28,32,33]).Their metric values show relatively significant differences between AMSR2 and in situ SWC data, indicating the need to improve the retrieval algorithm.
However, there are still differences worth noting.The metric values of the averaged ascending data in the second pair are much worse than those of the data at the descending overpass.This means that the descending retrievals would be more stable than the ascending because the effects of variations in both the spatial and profile variability of AMSR2 surface temperature are relatively small [28].In the other two pairs, the corrected in situ data are compared with AMSR2 data before and after correction.Given the relatively small difference between original and corrected in situ SWCs, these two pairs show values of all metrics similar to those of the original in situ and AMSR2 SWC.
The fourth pair, corrected AMSR2 and in situ SWC, is considered a potential pair for the evaluation of satellite soil moisture products because errors caused by TEs are apparent and do not physically exist.Though the difference is insignificant, it is compared to the pair of the original AMSR2 and the in situ SWC in detail.After correction, ME values are improved for both ascending and descending data.These values indicate the overestimation of AMSR2 SWC reduced by 5.4 percent and 9.2 percent for ascending and descending data, respectively.For the other four metrics, MAE, R, RMSE and ubRMSE, the values are improved by 4.1 percent, 20.7 percent, 3.7 percent and 3.5 percent for ascending data, and worsened by 4.6 percent, 2.4 percent, 3.8 percent and 4.3 percent for descending data.For the most commonly used metric, R, the improvement in the ascending data is more significant than the worsening in the descending data.
In addition to the averaged ascending and descending data, Table 3 also shows the MedAE values of the previous and following ascending data.Compared to MAEs, all MedAEs are much smaller because the effects of significant absolute errors are excluded.Like the other five metrics, MedAEs are improved in ascending data and worsened in descending data.In order to evaluate the impact of temperature effect removal on the difference between ascending and descending SWC, the absolute difference between θ Am and θ D is analyzed using Equation ( 12). Figure 7 shows the accumulated distribution functions of the absolute difference between θ Am and θ D of the in situ (Figure 7a) and AMSR2 SWC (Figure 7b).Clearly, the low part of the accumulated distribution functions moved left after the TE removal.For both the in situ and AMSR2 SWC, more than 80% of the absolute differences become smaller.Equation ( 13) is further utilized to determine the median value across all datasets.Employing this approach helps mitigate the influence of outliers mainly attributed to rainfall factors, enabling a more comprehensive assessment of the efficacy of TE removal on data affected by TEs.For the in situ SWC, the median value of the absolute difference between ascending and descending data was reduced from 0.0072 to 0.0031.For AMSR2 SWC, the same value was reduced from 0.0072 to 0.0036.Both reductions are almost half.

Discussion
This study confirmed that temperature effects exist in AMSR2 SWC product as well as in dielectrically measured in situ SWC [37] and SMAP SWC product [33] by visual inspection of time series of SWC and temperature and analysis of ADA triangles.The ascending and descending AMSR2 SWC data show diurnal fluctuation synchronized with AMSR2 surface temperature.The ascending SWC data in the daytime are more significant affected by the surface temperature than the descending nighttime ones.Similar to dielectrically measured in situ SWC and SMAP SWC products, this phenomenon can be considered as an apparent diurnal change caused by temperature effects because it cannot be reasonably explained by other physical factors, such as evapotranspiration [47,48].
The analysis of ADA triangles makes it possible to identify the contribution of temperature effects in the differences between ascending and descending SWC data.Then, a new temperature effect removal algorithm suitable for satellite SWC products, which only have two observations at ascending and descending overpasses, is designed following Kapilaratne and Lu's concepts [37,41].The application of this new algorithm proves itself capable of reducing the differences between ascending and descending SWC data in both the AMSR2 product and the in situ SWC dataset using observations at satellite overpasses.Furthermore, the ADA triangles after TE removal represent very close slopes between the previous and following ascending data, implying that the new algorithm mainly removes TEs while conserving the SWC declination caused by other factors.This algorithm can effectively remove TEs using data observed at the ascending and descending overpasses only.
The impacts of the newly developed TE removal algorithm on AMSR2 and in situ SWC are then evaluated using corrected and original data.Besides the widely used mean error, mean absolute error, correlation coefficient, root mean square error and unbiased root mean square error, the median absolute error is introduced to evaluate the differences between original and corrected SWC data without the effects of significant errors caused by factors other than temperature.It is shown that the differences between in situ and AMSR2 SWC, corrected or not, are much larger than those between original and corrected in situ SWC data.Though the relative importance of TE removal may increase while AMSR2 SWC improves and its difference with the in situ SWC narrows, the impact of TE removal is currently limited.
Though the difference is insignificant, a detailed comparison between the original SWC pair and the corrected SWC pair is made because the latter is considered a potential pair capable of excluding apparent errors caused by TEs in evaluating satellite soil moisture products.After correction, ME values showing the retrieval bias are improved for ascending and descending data.The other five metrics, MedAE, MAE, R, RMSE and ubRMSE, present a significant improvement in the ascending data and a slight worsening in the descending data.
Finally, the absolute difference between a descending SWC and the average of its two adjacent ascending SWCs is analyzed.It shows that the low part of the accumulated distribution functions of the absolute difference of the in situ and AMSR2 SWC moved leftward after TE removal.The median absolute differences of absolute differences in AMSR2 and in situ SWC are almost half.More than 80 percent of data points reduced the absolute difference between ascending and descending data.This reduction denotes that the ascending and descending data of AMSR2 and the in situ SWC become more homogenous after TE removal.This result may allow using both ascending and descending data and double the temporal resolution of AMSR2 SWC data.

Conclusions
The findings of this study mark a significant milestone in our understanding of the Advanced Microwave Scanning Radiometer 2 (AMSR2) soil moisture products.Firstly, our research confirms the presence of TEs in these datasets, aligning with previous studies that have identified TEs in SMAP soil water content data.Additionally, we have successfully devised a novel method for removing these effects from AMSR2 soil water content, marking the first instance of such a solution being developed.Our results indicate a modest improvement in the relationship between AMSR2 and in situ soil water content following the application of our correction method.While it is worth noting that the errors attributed to TEs may not be the predominant factor contributing to the observed differences between AMSR2 and in situ soil water content, our correction method nonetheless brings the ascending and descending data into closer alignment.Consequently, this enhancement allows for the combined utilization of both ascending and descending data, effectively doubling the temporal resolution of AMSR2 soil water content data for future research endeavors.

Figure 1 .
Figure1.Ratios of temperature (R T ) and SWC (R θ ) in three categories, Ap/D: ratio of value at previous ascending overpass to that at descending overpass, D/D: ratio of value at descending overpass to itself, Af/D: ratio of value at following ascending overpass to that at descending overpass.The red diamond is the average ratio of temperature of all ADA triples and the blue circle is the average ratio of SWC of ADA triples.The red square shows (R T,Am ), the average ratios of temperature in category Ap/D and Af/D; and the blue square shows (R θ,Am ), the average ratios of SWC in category Ap/D and Af/D.The yellow double-head arrow line indicates the differences between two ratios of SWC in categories Ap/D and Af/D which may be caused by factors other than temperature, and the blue double-head arrow line indicates the difference caused by temperature effects (TEs).

Figure 2 .
Figure 2. Study area in Mongolia.The green and red dots are automatic stations for soil hydrology (ASSHs) and automatic weather stations (AWSs).

Figure 3 .
Figure 3.The variety among nine stations in the Mongolia CVS, presenting their respective average values alongside the standard deviation over a five-day period.Panel (a) displays the soil water content values, while panel (b) showcases the soil temperature readings.

Figure 4 .
Figure 4.The bi-hourly in situ soil temperature, in situ SWC and AMSR2 SWC and AMSR2 LST at ascending and descending overpasses averaged over CVS in Mongolia (a).The blue dashed line and orange dotted line represent the in situ SWC (θ InSitu ) and soil temperature (T S ) at 3 cm depth, respectively.The red and green stars stand for AMSR2 LST at ascending (T SF,A ) and descending (T SF,D ) overpasses, and the red and green pluses are AMSR2 SWC at ascending (θ AMSR2,A ) and descending overpasses (θ AMSR2,D ), respectively.(b,c) The ADA triangles of temperature and SWC of in situ and satellite data.The red lines are the temperature and the blue lines are the SWC in both.

Figure 5 .
Figure 5. Estimation of temperature effect coefficient from in situ data (a) and AMSR2 data (b).The blue dots are the original data and the red circles are the outliers determined in the recursive regression.

3. 3 .
Evaluation of the Impact of Temperature Effect Removal 3.3.1.Visual Confirmation of the Impact of Temperature Effect Removal on the Difference between the Data from Ascending and Descending Overpasses

Figure 6 .
Figure 6.The bi-hourly in situ soil temperature, in situ SWC and the original and corrected AMSR2 SWC and AMSR2 LST at ascending and descending overpasses averaged over the CVS in Mongolia (a).In addition to Figure 4a, a red line representing the corrected in situ SWC (θ ′ InSitu ) at 3 cm depth and the red and green circles showing the corrected AMSR2 SWC at ascending (θ ′ AMSR2,A ) and descending (θ ′ AMSR2,D ) overpasses are added.The green ADA triangles of the corrected in situ and the AMSR2 SWC are also added in both (b,c).

3. 3 . 2 .
Quantitative Evaluation of the Impact of Temperature Effect Removal on In Situ and AMSR2 SWC Datasets By applying the newly developed temperature removal algorithm, corrected in situ and AMSR2 SWC (θ ′ InSitu and θ ′ AMSR2

Figure 7 .
Figure 7. Accumulated absolute difference between θ Am and θ D of in situ SWC (a) and AMSR2 SWC (b).The blue and red lines are the data before and after the correction.

Table 2 .
Values of the quantitative metrics, mean error, mean absolute error, correlation coefficient, root mean square error and unbiased root mean square error of the averaged ascending data and descending data in four pairs of the original and corrected AMSR2 and in situ SWC products from September 2016 to October 2019.

Table 3 .
Values of median absolute errors of the averaged ascending data and descending data in four pairs of the original and corrected AMSR2 and in situ SWC products from September 2016 to October 2019.,Af , A m and D are the same as those defined in Equation (5).3.3.3.Quantitative Evaluation of the Impact of Temperature Effect Removal on the Difference between Ascending and Descending SWC