# Assimilation of Satellite-Derived Soil Moisture and Brightness Temperature in Land Surface Models: A Review

Reza Khandan, Wigneron J.-P., Stefania Bonafoni, Arastoo Pour Biazar, Mehdi Gholamnia

<figure><img src="image_1.png" /></figure>

## To cite this version:

Reza Khandan, Wigneron J.-P., Stefania Bonafoni, Arastoo Pour Biazar, Mehdi Gholamnia. Assimilation of Satellite-Derived Soil Moisture and Brightness Temperature in Land Surface Models: A Review. Remote Sensing, 2022, 14 (3), pp.770. 10.3390/rs14030770. hal-03609946

HAL Id: hal-03609946

https://hal.inrae.fr/hal-03609946

Submitted on 16 Mar 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_3.png" /></figure>

Review

# Assimilation of Satellite-Derived Soil Moisture and Brightness Temperature in Land Surface Models: A Review

Reza Khandan <fnref n="1" />*, Jean-Pierre Wigneron <fnref n="2" />, Stefania Bonafoni <fnref n="3" /> ID, Arastoo Pour Biazar <fnref n="4" /> ID
and Mehdi Gholamnia <fnref n="5" /> ID

<figure><img src="image_4.png" /></figure>

check for
updates

**Citation:** Khandan, R.;
Wigneron, J.-P.; Bonafoni, S.;
Biazar, A.P.; Gholamnia, M.
Assimilation of Satellite-Derived Soil
Moisture and Brightness
Temperature in Land Surface Models:
A Review. *Remote Sens.* **2022**, *14*, 770.
https://doi.org/10.3390/rs14030770

Academic Editor: Prashant K.
Srivastava

Received: 26 December 2021

Accepted: 3 February 2022

Published: 7 February 2022

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_5.png" /></figure>

**Copyright:** © 2022 by the authors.
Licensee MDPI, Basel, Switzerland.
This article is an open access article
distributed under the terms and
conditions of the Creative Commons
Attribution (CC BY) license (https://
creativecommons.org/licenses/by/
4.0/).

**Abstract:** The correction of Soil Moisture (SM) estimates in Land Surface Models (LSMs) is considered essential for improving the performance of numerical weather forecasting and hydrologic models used in weather and climate studies. Along with surface screen-level variables, the satellite data, including Brightness Temperature (BT) from passive microwave sensors, and retrieved SM from active, passive, or combined active–passive sensor products have been used as two critical inputs in improvements of the LSM. The present study reviewed the current status in correcting LSM SM estimates, evaluating the results with in situ measurements. Based on findings from previous studies, a detailed analysis of related issues in the assimilation of SM in LSM, including bias correction of satellite data, applied LSMs and in situ observations, input data from various satellite sensors, sources of errors, calibration (both LSM and radiative transfer model), are discussed. Moreover, assimilation approaches are compared, and considerations for assimilation implementation are presented. A quantitative representation of results from the literature review, including ranges and variability of improvements in LSMs due to assimilation, are analyzed for both surface and root zone SM. A direction for future studies is then presented.

**Keywords:** Soil Moisture (SM); assimilation; Land Surface Model (LSM); Radiative Transfer Model (RTM); surface; root zone

# 1. Introduction

Soil Moisture (SM) is a crucial variable in the partitioning of water (into infiltration and runoff) and energy (into sensible and latent heat flux) and has been considered in many atmospheric and hydrological studies [1,2]. SM influences the carbon cycling by photosynthesis and transpiration processes [3] and affects short and medium-range forecasts in meteorology [4]. The variability in SM is due to the spatial distribution of many factors such as rainfall, the spatial variation of the soil texture, vegetation, and topography. Volumetric SM measurements at different depths are limited to sparse ground measurements, which cannot provide required spatial and temporal resolution data for numerical weather initialization. To solve this problem, Land Surface Models (LSMs) are used to model the near-surface (~0–5 cm) and root zone (~5 to 100) SM (hereafter referred to as SSM and RSM) behavior through physical and hydrological laws. However, the estimation of SM by LSMs is limited by forcing data, built-in physical rules, and parameterizations.

In order to reduce the mentioned uncertainties, the assimilation of surface observations in LSMs was proposed. Generally, the assimilation procedure in LSMs consists of three main

<footnote n="1">¹ Department of GIS and Remote Sensing, Faculty of Geography, University of Tehran, Teheran 14155-6619, Iran</footnote>

<footnote n="2">² INRAE, nUMR1391 ISPA, F-33140 Villenave d'Ornon, France; jean-pierre.wigneron@inra.fr</footnote>

<footnote n="3">³ Department of Engineering, University of Perugia, Via G. Duranti 93, 06125 Perugia, Italy;
stefania.bonafoni@unipg.it</footnote>

<footnote n="4">⁴ Atmospheric Science Department, University of Alabama in Huntsville, Huntsville, AL 35805, USA;
biazar@nsstc.uah.edu</footnote>

<footnote n="5">5 Department of Civil Engineering, Islamic Azad University, Sanandaj Branch, Sanandaj P.O. Box 618, Iran;
mehdi_gholamnia@ut.ac.ir</footnote>

<footnote>* Correspondence: rs.reza_khandan@ut.ac.ir</footnote>
components: (i) observation, (ii) modeling (LSM) for background and forecast generation of SM, and (iii) correction of background SM in the LSM based on observations. As stated in [5], the initial assimilation experiments based on the optimal interpolation method used screen-level temperature and relative humidity as assimilation inputs. However, this approach had some built-in limitations, such as applied fix coefficients, the accumulation of errors induced by the model, forcing data in the root zone, and sparse data sampling [5-7].

The emergence of satellite data both in the infrared and microwave spectrum has made it possible to measure the Brightness Temperature (BT) and to estimate volumetric SMs at the near surface (SSM). Such observations can be applied as an external source for the correction of biases in SM initialization in LSMs. For this purpose, microwave imaging sensors (both passive and active) have been more desirable due to their low attenuation by clouds and atmospheric aerosols [8]. The passive sensors include: Soil Moisture Active and Passive (SMAP) radiometer (L-band 1.41 GHz) [9]; Soil Moisture and Ocean Salinity (SMOS) (L-band 1.4 GHz) [10]; Advanced Microwave Scanning Radiometer-Earth Observing System (AMSR-E) (C, X, K, Ka, and W bands at six frequencies: 6.93, 10.65, 18.7, 23.8, 36.5, and 89.0 GHz) [11]; Advanced Microwave Scanning Radiometer 2 (AMSR2) (C, X, K, Ka, and W bands at seven frequencies: 6.93, 7.3, 10.65, 18.7, 23.8, 36.5, and 89 GHz) [12-14]; Scanning Multichannel Microwave Radiometer (SMMR) (C, X, K, and Ka bands at five frequencies: 6.6, 10.7, 18.0, 21.0, and 37.0 GHz) [15]; and Tropical Rainfall Measuring Mission's (TRMM) Microwave Imager (TMI) (X, K, Ka, and W bands at five frequencies 10.7, 19.4, 21.3, 37.0, and 85.5 GHz) [16]. The active sensors include: Advanced Scatterometer (ASCAT) (C-band 5.255 GHz) [13]; Sentinel-1 (C-band 5.405 GHz) [14]; and ERS Scatterometer (C-band 5.32 GHz) [17]. The merged products of these sensors, Climate Change Initiative (CCI) [18] and Soil Moisture Operational Products System (SMOPS) [19], can also be applied for such assimilation studies.

A successful assimilation procedure depends on the assimilation approach, LSM physical laws and parameterizations, the Radiative Transfer Model (RTM) parameterization (as an observation operator for BT), and the observations and errors in both the model and observations [20,21]. The theoretical issues about the assimilation of both satellite BT and SM were discussed in [20,22,23]. As miscellaneous papers used different satellite data, LSMs, data preparation, forcing data, assimilation approaches (with different settings), observation operators (RTM), and different locations, it is helpful to normalize and summarize their results for future studies.

In order to search and select the related papers, the following procedure was conducted. (A) We searched the keywords "Assimilation", "Land surface model", "Soil moisture", "Brightness temperature", "satellite", and "Remote sensing" in the Google Scholar search engine to the end of 2019. (B) The found papers were sorted based on "relevance" and were from "any type" provided by Google Scholar, excluding conference papers. (C) After that, we read the abstract of retrieved papers one by one. Then, we selected the papers that used satellite-derived SM or BT for assimilation in LSMs (even if it was not the main procedure of the research). (D) We then performed the same procedure for different orders of keywords. (E) Finally, we downloaded selected papers and classified them based on year of publication, and we chose to consider the papers from 1994 to 2019 (this being the final paper). (F) For boosting the search of related papers, when reading each selected paper: (1) we noted the used references in the paper and downloaded them (in the range 1994-2019) and (2) we searched those papers that cited that paper (through Google Scholar). Then, we updated our archive. (G) After reading all papers and extracting their results, we started to filter them; in this stage, we removed those papers that (i) used synthetic brightness temperatures for assimilation, (ii) did not provide the same evaluation metrics (or did not provide quantitative results), (iii) did not use in situ measurements for the validation of an assimilation procedure. Finally, we started to analyze the results. Examples of ground station networks found in the reviewed papers (from 2005 to 2019) are shown in Figure 1. The horizontal axis indicates the number of papers that used a particular network, and the vertical axis shows the various networks applied in assimilation studies.
The ISMN network's sub-networks are shown in blue, and other individual networks are depicted with different colors. The acronyms used for the networks are described in the Abbreviations.

<figure><img src="image_6.png" /><figcaption>**Figure 1.** The application frequency of SM networks used in assimilation studies from the selected papers in the time interval 2005–2019. The acronyms are explained in the Abbreviations. ISMN sub-networks are shown in blue and other networks have different colors; \* ISMN in the vertical axis includes papers that used the ISMN network but did not specify the subcategories of ISMN.</figcaption></figure>

The review paper is organized as follows. In Section 2, the current LSM SM improvement procedures are presented, and the uncertainties in LSMs, satellite BT, and SM for assimilation are reported. The correction of bias between the LSM and satellite is also reviewed. This section also refers to the applied assimilation methods regarding their characteristics and limitations. Section 3 presents the outcomes of studies based on LSM classification, and Section 4 identifies some of the issues with current approaches assimilating SM and BT in LSMs and offers recommendations for improvement.

## 2. Procedures for the Improvement of SM in LSMs

Based on the reviewed papers, when using both pieces of satellite data inputs (BT and retrieved SM), the correction of LSM SM estimates can be divided into the following categories: (A) some papers followed the correction of LSM climatology by calibration (optimization of the fixed parameters) of the LSM, (B) others conducted the correction of forcing data (required data for initialization of LSMs) by assimilation. In Figure 2, the procedures for both calibration (approaches 1 and 2 in Figure 2) and assimilation (approaches 3 and 4 in Figure 2) in the LSM are presented. In approach 1, the LSM SM estimates are applied to simulate the BTs (through RTM), and then an optimization method is used to decrease the difference between the simulated and observed BTs, and hence to optimize the critical parameters in the LSM (e.g., [24]). In approach 2, the optimization of LSM parameters is conducted by comparing the estimated LSM SM to satellite retrieved SM products (e.g., [25]). In approach 3, RTM is calibrated by adjusting LSM SM to reduce the difference between the observed and simulated BTs (e.g., [26–28]). Systems such as the Coupled Land and Vegetation Data Assimilation System (CLVDAS) were designed (as a dual-pass assimilation system), which simultaneously calibrates the parameters and assimilates the BT in the LSM [29]. In approach 4, satellite SM is used through assimilation to correct the LSM SM estimates. In both approaches (calibration and assimilation), RTM calibration can be conducted.
<figure><img src="image_7.png" /><figcaption>**Figure 2.** The current procedures for calibration of LSM and RTM and assimilation approaches.</figcaption></figure>

## 2.1. Land Surface Models

Over the last decade, 2 m relative humidity, air temperature (from ground stations), and BT and SM products (from satellites) have been assimilated in different LSMs to improve SM estimates. Figure 3 shows different LSMs (total number: 73) applied in SM assimilations across the years in reviewed papers. According to Figure 3, the number of studies has increased since 2007, and most of these were conducted between 2013 and 2017. Most of the assimilation studies were implemented in Noah, CLSM, and the Community Land Model (CLM).

<figure><img src="image_8.png" /><figcaption>**Figure 3.** SM assimilation studies classified by LSMs between 2005 and 2019 (legend acronym explanation in Table 1).</figcaption></figure>

The applied LSMs along with their default soil layers are shown in Table 1.
<table><caption>**Table 1.** The layering assumptions in LSMs; **: depends on land cover, tall vegetation fraction 0.15 and 1, bare soil 0.15. *: Discretized to 5.</caption><thead><tr><th>LSM</th><th>Number and Depth of Soil Layers (m)</th><th>Reference</th></tr></thead><tbody><tr><td>Noah</td><td>4</td><td>0.1, 0.3, 0.6, 1.0</td><td>[30,31]</td></tr><tr><td>Australian Water Resources Assessment system (AWRA-L)</td><td>3</td><td>0.1, 0.2, and 8</td><td>[32]</td></tr><tr><td>Catchment Land Surface Model (CLSM)</td><td>2</td><td>0.02, 1.0</td><td>[33]</td></tr><tr><td>Community Land Model (CLM)</td><td>10<br>0.204, 0.336, 0.553,<br>0.913, 1.506</td><td>0.018, 0.028, 0.045, 0.075, 0.124</td><td>[34]</td></tr><tr><td>Joint UK Land Environment Simulator (JULES)</td><td>4</td><td>0.1, 0.25, 0.65, 2.0</td><td>[35]</td></tr><tr><td>Interactions between Surface, Biosphere, and Atmosphere (ISBA.Diffusion)</td><td>14</td><td>0.01, 0.04, 0.1, 0.2, 0.4, 0.6, 0.8,<br>1.0, 1.5, 2.0, 3.0, 5.0, 8.0, 12.0</td><td>[36]</td></tr><tr><td>Simple Biosphere Model 2 (SiB2)</td><td>3</td><td>0.05, 0.20 and 2</td><td>[37]</td></tr><tr><td>Variable Infiltration Capacity model (VIC)</td><td>3</td><td>3; 0.1, 1.6, and 1.9</td><td>[38]</td></tr><tr><td>Soil, Vegetation, and Snow (SVS)</td><td>7</td><td>0.1, 0.2, 1, 2, and 3</td><td>[39]</td></tr><tr><td>World Wide Water (W3)</td><td>3</td><td>0.05, 1, and 10</td><td>[40]</td></tr><tr><td>Global Land Evaporation Amsterdam Model (GLEAM)</td><td>**</td><td>0.15, 1, and 2.5, and short</td><td>[41]</td></tr><tr><td>Distributed Hydrology Soil Vegetation Model (DHSVM)</td><td>4</td><td>0.01, 0.05, 0.4, and 0.8</td><td>[42]</td></tr><tr><td>HYDRUS</td><td>100</td><td>100 cm discretized to 1 cm layers</td><td>[43]</td></tr><tr><td>Modélisation Environmentale Communautaire (MEC)<br>Surface and Hydrology (MESH)</td><td>3</td><td>0.10, 0.35, and 0.41</td><td>[44]</td></tr><tr><td>SW-ET</td><td>*</td><td>0.05, 0.10, 0.20, and 0.50 and 1</td><td>[45]</td></tr><tr><td>Tiled ECMWF Scheme for Surface Exchanges over Land (TESSEL) with improved land surface hydrology (H-TESEL)</td><td>4</td><td>0.07, 0.21, 0.72, and 1.89</td><td>[46]</td></tr><tr><td>Met Office Surface Exchange Scheme version 2 (MOSES 2)</td><td>4</td><td>0.1, 0.25, 0.65, and 2</td><td>[47]</td></tr></tbody></table>

## 2.2. Uncertainties in LSMs, RTM, and Satellite BT and SM Products

### 2.2.1. LSMs

The sources of errors in LSMs include a large number of their parameters influenced by the vegetation and soil state estimation [48], model physics and structure, forcing data, and boundary conditions [49]. The role of precipitation as forcing data for assimilation was proven to be significant (discussed further). Figure 4 shows the frequency of applied forcing data for the assimilation of SM in LSMs in the reviewed papers, showing that Modern-Era Retrospective analysis for Research and Applications (MERRA) [50] was used the most. Blyverket et al. [51] compared the performance of MERRA-2 and NLDAS (with precipitation correction) over the contiguous United States (US). They suggested that the assimilation procedure has higher improvement skills with MERRA-2; the correction of precipitation in NLDAS may improve the results of SM estimation, and hence reduce the assimilation skill [52].
<figure><img src="image_9.png" /><figcaption>**Figure 4.** The application frequency (count) of different forcing data for assimilating SM in LSMs in the reviewed papers.</figcaption></figure>

### 2.2.2. RTM

RTMs include a large number of parameters that are highly variable and sometimes unavailable [53] as well as showing dependency to land cover [28]. See Section 2.3 for more details.

### 2.2.3. Satellite Data

The inherent uncertainties in both BT and SM products that may influence assimilation are:

1. SM retrieval approaches. Different strategies and simplifications applied for parameterization of RTMs and soil dielectric mixing models in SM retrieval. As stated in [54,55], most of the SM retrieval algorithms used T-ω zero-order RTM as the baseline. The major difference between these approaches originates from their parameterizations and estimation of vegetation optical depth. The modeling of roughness, soil, and canopy temperatures, vegetation structure and its optical depth, and atmospheric effects are the most significant factors in SM retrieval. Choosing between effective single or multiple scattering albedo, single or double polarization, static or dynamic roughness coefficients, and single- or multi-angle satellite observations, the applied frequency and ancillary data make the algorithms complex and define the quality of SM products. Moreover, as discussed in [54], the selection of mixing dielectric modeling is another critical issue that needs to be regarded when using SM products and evaluating assimilation results.

2. The depth of satellite-derived SM. The depth of measurement is not definite and is a function of moisture, temperature, and soil texture. The introduced depths vary spatially and temporally [23,56,57].

3. The time of satellite orbital pass. Some studies have investigated the accuracy of the SM products of a sensor at different acquisition times. For example, Martens et al. [58] and Yee et al. [59] showed that ascending SMOS data (about 6:00 a.m.) were more accurate than the descending pass (about 6:00 p.m.). However, Jing et al. [60] maintained that the descending pass of SMOS performed better than the ascending pass. In [17,61], it was stated that ASCAT SM data from the descending pass (9:30 a.m.) have more accuracy than the ascending one (9:30 p.m.). Jing et al. [60] noted that
the AMSR-E ascending (1:30 p.m.) product was better than the AMSR-E descending (1:30 a.m.) product. Yee et al. [59] stated that AMSR2 X-band products have better performance in evening overpasses (1:30 pm) in contrast with morning overpasses (1:30 a.m.). As stated in [62], it is expected that night time or early morning SM products be more accurate due to: (a) assumed temperature equilibrium between soil and vegetation canopy in the SM retrieval modeling [63] and (b) expected minimum of Faraday rotations during the night [64,65]. However, as indicated above in some studies, other revisit times showed more accuracy, which can be attributed to other factors [62]. These are: (a) Radio Frequency Interference (RFI) contamination in some time periods and (b) diurnal variability of surface conditions such as water in vegetation. Therefore, it is expected that the accuracy of SM products in different acquisition times be checked in the study area before conducting an assimilation study.

4. Land cover. The vegetation coverage and type affect the SM retrieval, especially at higher frequencies. In [17,61], it was reported that X and C bands are opaque in dense forests and shrubs with green vegetation fractions of more than 0.5. Other land covers, such as frozen ground, snow cover, water body areas (e.g., flooded areas, rivers, wetlands, lakes, precipitation at the time of satellite passes), steep topography (e.g., rocks), urban area, and heavily forested areas [66] would decrease the quality flag of data recording and, as a result, reduce the available data for assimilation.

5. Temperature limitation. Due to the freezing temperature, currently, no reliable satellite-derived SM product is available for low-temperature latitudes and snow covers [67].

6. Temporal resolution. The higher the temporal resolution, the more accurate SM temporal dynamics (e.g., caused by irrigation) can be retrieved. Yin et al. [68] proposed that SMOPS products (as combined products) have a higher temporal resolution, which can capture SM dynamics with more detail.

7. Spatial resolution. The low spatial resolution of passive microwave sensors, the mismatch between the LSM and satellite spatial resolutions, and the unmolded sub-pixel heterogeneity (mixed land covers mentioned in item 3, land cover, within a pixel) could increase the errors. Downscaling the SM products can be a solution, but the accuracy of downscaling depends on applied ancillary data, downscaling the method and accuracy of in situ measurements, which needs to be analyzed [69]. Therefore, in some cases, it may bring forth lower results. For example, Lievens et al. [70] compared the performance of downscaled SM (by MODIS thermal data) with coarse resolution and showed that downscaling could not provide better results.

8. Polarization. Vertical and horizontal polarization might not retrieve the same BTs; vertical polarization BTs are less sensitive to vegetation heterogeneity and roughness than horizontal polarization while they are more sensitive to SM [4,70,71]. Tian et al. [26] showed that the horizontal polarization BTs from AMSR-E have lower temporal variability than the vertical, which can make them less suitable for assimilation experiments. The role of incidence angle could make the problem more complex; Lievens et al. [70] showed that horizontal polarization with incidence angle less than 42.5 degrees could provide better results than in combination with vertical polarization and multi-angle observations (with SMOS sensor). They stated that it could be attributed to the correlation between dual-polarization and multi-angle observations. The higher sensitivity of horizontal polarization to SM and the lower sensitivity of vertical polarization to vegetation and roughness are issues that require more research to find out the weight of each one of these characteristics.

9. Dynamic drying rate. The drying rate of estimated SMs in LSMs could be different from satellite-derived SMs. Shellito et al. [72] compared the drying rates of SMAP and those simulated by Noah and showed that the drying SSM from SMAP is faster than Noah simulations. They also stated that when SM content (Noah and SMAP SM contents are linearly related) is high, potential evaporation is high, vegetation cover is low, and SMAP drying is the fastest. They reported that the effect of vegetation on the Noah model is simplified and is less than SMs retrieved by SMAP.
To find out the biases in satellite products, several studies compared them with ground networks. Although most of them stated that SMAP, SMOS-IC (V105 and V2), and CCI (based on LSM corrections) have slightly better performance [73–75], it should be considered that SMAP product biases differ based on soil and vegetation conditions. Figure 5 shows the applied satellite data in assimilation studies. As Figure 5 shows, several microwave data and an algorithm (ALEXI, Atmosphere-Land Exchange Inverse, based on the thermal infrared band and used for comparison with microwave data) were applied for the assimilation studies, and the AMSR-E was used the most. The increasing trend in the application of such data requires that the best products are identified based on the spatio-temporal conditions of the study area and the characteristics of applied evaluation methods [76].

<figure><img src="image_10.png" /><figcaption>**Figure 5.** The application frequency (count) of satellite/sensor data for SM assimilation studies in the reviewed papers.</figcaption></figure>

## 2.3. Calibration of LSM and RTM Parameters

As LSM and RTM parameter values are assumed to be fixed, their resultant error in the assimilation process is inevitable. It is essential that the biases between RTM and LSM (through simulated BT) and satellite data decrease [77,78]. Therefore, these parameters in both models should be calibrated by comparison with satellite-observed BT. In parameter calibration (optimization), several issues need to be considered, including the applied methods for optimization, the selection of parameters to be optimized, uncertainties in calibration, time window setting, and the number of used microwave frequencies in the calibration process.

The applied methods for the optimization of parameters in LSM and RTM are:
(1) Shuffled Complex Evolution (SCE-UA); this minimizes a cost function composed of observed and simulated BT (by LSM or RTM), wavelength, and polarization of observation and then optimizes the parameters in their value domain using complex evolution algorithm [26,79–81]. (2) Ensemble Proper Orthogonal Decomposition-based Parameter (EnPOD_P); this constructs an ensemble of state variables (by the Monte Carlo method), including the perturbed states (simulated BTs) and perturbed parameters at each time of satellite BT, and finally optimizes a cost function and obtains the parameters [49,53,82].
(3) Particle Swarm Optimization; this uses multi-angular SMOS BTs and tries to minimize a
cost function composed of the penalty terms for bias between simulated BTs and SMOS
BTs, the variability (temporal standard deviation) and the parameters by Particle Swarm
Optimization algorithm [83]. (4) Particle Filtering; this simultaneously estimates SM and
some LSM parameters using state augmentation method [84]. (5) Markov Chain Monte
Carlo (MCMC); this uses MCMC algorithm to estimate RTM parameters by minimizing
the differences between simulated BTs and multi-angle SMOS BTs [21]. Tian et al. [82]
showed that the EnPOD_P considerably outperforms the SCE-UA due to the simultaneous
optimization of the model states and parameters. Furthermore, De Lannoy et al. [21]
indicated that the parameter values estimated through Particle Swarm Optimization were
in close agreement with those obtained from the MCMC simulation. Further research is
essential to compare the performance of these optimization methods for both LSMs and
RTMs in different conditions.

The selection of parameters for calibration is another crucial issue that is used in
both RTM and LSM. Generally, the sensitive parameters in RTM to BT (microwave emis-
sivity) are mainly from vegetation and surface roughness parameters, soil structure,
and texture [77,85,86]. Based on the conducted sensitivity analysis in observation operators
(Land Emissivity Model [87], Community Microwave Emission Model [88], and T-ω) [89],
these parameters include the following. (1) For vegetation: leaf thickness, affecting vegeta-
tion optical thickness [53,90], vegetation structure parameter, affecting vegetation optical
thickness [53,91], the radius of dense medium scatters, fraction volume of dense medium
scatters [90], and b parameter, affecting vegetation optical depth in the T-ω model [91].
(2) For soil: surface roughness parameter (HR), affecting soil emissivity [92,93], charac-
terized as the standard deviation of surface height [4,53,90,94], correlation length of the
surface [4,29,78,95], soil bulk volume density [49,53] sand percentage, clay percentage,
and soil porosity. Montzka et al. [96] expressed that as the sensitivity of the T-ω model to
roughness parameters is less than vegetation optical depth, roughness parameters showed
higher uncertainties. Sawada et al. [95] found that texture parameters (clay and sand) are
less sensitive in contrast with roughness parameters.

In LSM, the most sensitive parameters to SM also need to be considered. Sawada et al. [29] and Sawada et al. [95] conducted a sensitivity analysis and selected the most sensitive parameters (both hydrological and ecological). These parameters include saturated hydraulic conductivity, saturation water content (porosity), parameters in van Genuchten's water retention curve, the maximum rubisco capacity of top leaf, the normal turnover rate of leaves, the empirical parameter of root biomass fraction function, point of water-related stress loss, and point of temperature-related stress loss. After choosing the sensitive parameters, the works of [29,78,97] optimized the selected parameters by comparing simulated BTs and satellite-observed BTs. Koster et al. [27] applied a different method, where a time-invariant parameter was determined for calibration of LSM SM estimates, and through comparison with SMAP SM observations, the parameter value was determined.

Many factors should be considered during calibration. In RTM calibration, Zhao et al. [77] maintained that, due to the limited number of iterations in the optimization process, it is possible that calibrated parameters are not efficiently optimized. This can produce degraded results in contrast to non-calibrated parameters. Moreover, another issue is the consideration of land cover in calibration. For example, De Lannoy et al. [28] showed that the calibration of the RTM revealed higher errors in BT predictions over cropland. In LSM calibration, Yang et al. [97] recognized that the estimated parameter values might not be the truth values of soil parameters because, firstly, dense vegetation could reduce the model calibration effect; in areas with dense vegetation, microwave signals provide more information about vegetation rather than soil state. This may reduce the accuracy of the soil model calibration. Secondly, estimated parameter values may be biased due to inappropriate parameterizations or simplified physical assumptions in an LSM, such as the partitioning between runoff generation and SM content and the relationship between surface temperature and SM. Furthermore, current RTMs and LSMs are usually calibrated by passive microwave satellites with pixel sizes of more than 10 km. The issues of spatial
matching between models, satellite data, and heterogeneity within each pixel could cause
errors in the calibration process.

The time window for both RTM and LSM calibration is another issue that needs to be considered. For RTM, no study analyzed the effect of time window on the calibration process in the analyzed period of this study, but for LSM, [97] showed that a time window of 1 year for the model calibration might be reasonable for grasslands. More studies are needed to analyze the effect of a time window.

In most studies for RTM and LSM calibration, one satellite frequency was applied, but Sawada et al. [29] applied a cost function that used three frequencies to provide more stability as it could provide more robust results.

2.4. Bias Correction in SM Assimilation

The systematic error or bias removal between the LSM and satellite-derived SM should be considered for the assimilation procedure. The applied strategies to deal with these biases are shown in Figure 6.

<figure><img src="image_11.png" /><figcaption>**Figure 6.** The applied strategies to deal with bias between LSM estimates and satellite-derived SM.</figcaption></figure>

Several online and a priori methods have been introduced for bias correction between the LSM, in situ measurements, and satellite data [98–106]. In online (dynamic) bias correction, the biases are corrected in the assimilation process. This could be used when transient biases appear and they could not be identified in the climatology of both model and observation. In this case, the bias is attributed to either model [103,104] or observations [107].

The applied a priori bias correction methods (bias correction prior to assimilation) for satellite SM are Triple Collocation (TC), linear scaling (linear regression), and Cumulative Distribution Function (CDF). In TC, it is possible to use three independent datasets to remove the bias. Yilmaz et al. [101] evaluated the performance of least squares regression, variance matching, and TC. They showed that TC had better results compared to others.

Most of the studies reviewed here used the CDF technique for bias removal in SM
assimilation. In CDF matching, the CDFs of both observation and model at each grid cell are
extracted, and the observation values are bias-corrected to the model climatology through
the mean, variance, and possibly higher-order moments; Lievens et al. [108] emphasized the
mean as the most accurate rescaling method. A disadvantage of this method is that it needs
a long data record to provide a suitable sample size. Reichle et al. [109] proposed the replace-
ment of spatial sampling with temporal sampling to reduce the required time scale to one
year. Regarding the required long time series for CDF matching, Mahfouf [110] proposed
considering all points in the study area (global rescaling). Such a method was compared
with localized recalling (point-to-point matching between model and satellite observations)
in [111,112]. Kolassa et al. [111] showed that the global approach had higher skills in im-
proving SM estimation by having a higher impact of observation on assimilation procedure.
They stated that this approach is subjected to these conditions: (1) a good observation error
characterization; (2) rigorous observation quality control; (3) potential component-wise
assimilation to better isolate the reliable satellite information. Despite these results, in the point-wise approach, there is more sensitivity to locations where observation is uncertain or the range of change is high (mountainous areas). Schneider et al. [112] also showed the influential role of domain spatial diversity in models where global rescaling could remove important observation data. Blankenship et al. [113] and Blankenship et al. [114] used soil type-based CDF matching and showed that this method had better performance than localized matching. It needs to be considered that the spatial resolution (or subcategories) of soil types is an important variable. Blankenship et al. [114] also showed that vegetation could not have significant results in their applied CDF matching in the warm season of 2011 in central and southeastern United States.

According to some recent studies, from a temporal point of view, CDF matching can be conducted in two ways: lumped and seasonal [98,115]. In the lumped approach (one year CDF curve for all points), it is assumed that the difference between the model and observations is stationary in time. In the seasonal (or monthly) approach, the biases are corrected through temporal stratification. In the case of unmolded processes such as irrigation, the seasonal method would be more effective but it requires an appropriate sample size [116–119]. Yin et al. [120] compared the performance of lumped CDF matching, monthly CDF matching, and linear matching of monthly mean and standard deviation of SMAP and model SM simulations. They showed that monthly CDF matching had the best results in the prediction of precipitation.

As it is possible that in some situations (such as localized bias correction) the inaccurate model climatology is forced toward satellite one, some studies proposed a correction of the model toward observations. This would be applicable when the model biases are spatially known.

Some studies used bias-blind assimilation. In this case, it is assumed that observations and model estimates are not biased. Lin et al. [121] believe that all available data sets for assimilation should be kept intact to use the useful information within the data. Rescaling the observation to the model climatology could remove valuable information [113]. This is more critical when the sources of biases are unknown in both observation and the model.

## 2.5. Applied Approaches for Assimilation of SM and BT in LSMs

### 2.5.1. Recruitment of Assimilation Methods

Since 2005, several assimilation techniques have been applied for the improvement of SM in LSMs. Table 2 lists and describes the techniques used in these studies, while Figure 7 shows the specific methods used in the studies reviewed here. This figure also shows the increase in the number of studies using assimilation techniques (in particular, from 2014 to 2016), with most of the studies using EnKF (Ensemble Kalman Filter), as based on reviewed papers.

<table><caption>**Table 2.** Assimilation methods applied in the reviewed papers (from 2005 to 2019).</caption><thead><tr><td></td><td>Assimilation Methods</td><td>Description</td><td>References</td></tr></thead><tbody><tr><td>1</td><td>Kalman Filter (KF)</td><td>The state vector and covariance matrix of error are computed at each time step and are then forwarded in time by a linear model.</td><td>[122]</td></tr><tr><td>2</td><td>Extended Kalman Filter (EKF)</td><td>The Kalman Filter model and observation operators are replaced by non-linear equations to solve the inherent linear problem in the Kalman filter approach.</td><td>[123]</td></tr><tr><td>4</td><td>Ensemble Adjustment Kalman Filter (EAKF)</td><td>A deterministic configuration of EnKF where the observations are not perturbed; uses a linear operator for updating the prior ensemble.</td><td>[124]</td></tr><tr><td>5</td><td>Ensemble-Based Kalman Smoother (EBKS)</td><td>Time is included in EnKF, which causes the smoothing.</td><td>[125]</td></tr><tr><td>6</td><td>Ensemble based Four-dimensional Variational (En4DVAR)</td><td>A hybrid method where background error covariance is estimated from ensemble forecasts.</td><td>[126]</td></tr><tr><td>7</td><td>Ensemble Kalman Filter (EnKF)</td><td>The model states and error statistics are generated by perturbations in the model initial conditions (causing a set of ensembles) and are forwarded by a non-linear model without the need for linearization.</td><td>[127]</td></tr></tbody></table>
<table><caption>**Table 2. Cont.**</caption><thead><tr><td></td><td>Assimilation Methods</td><td>Description</td><td>References</td></tr></thead><tbody><tr><td>8</td><td>Ensemble Optimal Interpolation (EnOI)</td><td>A variant of EnKF that uses the prescribed background error-covariance.</td><td>[128]</td></tr><tr><td>9</td><td>Ensemble Square Root Filter (ESRF)</td><td>A deterministic configuration of EnKF where the observations are not perturbed; uses square-root of the forecast or analysis error covariance matrices.</td><td>[129]</td></tr><tr><td>10</td><td>Ensemble Transform Kalman Filter (ETKF)</td><td>A deterministic configuration of EnKF where the observations are not perturbed; uses an ensemble transformation matrix to compute analysis perturbation matrix.</td><td>[130]</td></tr><tr><td>11</td><td>Evolutionary (Evol)</td><td>Relies on evolutionary algorithms to provide evolved members for analysis step in assimilation.</td><td>[131]</td></tr><tr><td>12</td><td>Particle Filter (PF)</td><td>The probability density function of model SM estimates is represented by particles; a set of weighted trajectories (based on observations) are defined between background and analysis steps.</td><td>[132]</td></tr><tr><td>15</td><td>Newtonian Nudging (NNudg)</td><td>The model is dynamically relaxed toward the observations, and therefore the model is considered as a weak constraint and observations are considered as perfect.</td><td>[133,134]</td></tr><tr><td>16</td><td>Optimal Interpolation (OI)</td><td>Based on the known errors, observations are weighted and the analysis is estimated based on the gain matrix to obtain the analysis.</td><td>[135]</td></tr><tr><td>17</td><td>Variational (3D and 4D) (Var)</td><td>Model state is estimated by minimization of cost function, which explains the misfit between the model simulations and observations. If the cost function is conducted in a time interval, it is converted to 4D Var type.</td><td>[136]</td></tr></tbody></table>

<figure><img src="image_12.png" /><figcaption>**Figure 7.** The application frequency of each assimilation method for the improvement of SM in LSMs in each year from 2005 to 2019.</figcaption></figure>

### 2.5.2. Comparison of Assimilation Methods in Improvement of SM Estimates

Over the last decade, several studies have compared the performance of different assimilation methods on the improvement of SM in LSMs. De Rosnay et al. [5] replaced the Optimal Interpolation (OI) method with EKF using 2 m air temperature and relative humidity observations. They tested the possible replacement of ASCAT SM. They showed that EKF had superior performance over OI (as it accounts for the nonlinear relationship of meteorological forcing and SM). Dumedah et al. [137] used SMOS SM data and found that the surface estimation accuracy of SM for EnKF and Evolutionary Data Assimilation (EDA) assimilation methods increased, and EDA had better results than EnKF when
evaluating with SMOS level 2 SMs. EDA could converge the model parameters and improve the SM estimation. Blyverket et al. [51] compared the performance of EnKF, EnOI, and climatological background error (EnOI-Clim) assimilation methods using Climate Change Initiative (ESA CCI) by assimilating SM retrievals from SMAP and SMOS over the contiguous US. They showed that, although EnKF had slightly better performance, EnOI is recommended because it is computationally cheaper. The disadvantage of EnOI is that it does not use flow-dependent ensemble spread. Comparisons of the EKF and EnKF showed that EnKF had slightly better results in the southeastern United States using synthetic observations [138] and over an experimental site in southwestern France [139]. Bi et al. [140] compared the performance of EnKF, the standard PF, and the Improved PF (IPF) using brightness temperatures from the AMSR-E into the VIC LSM in the NaQu network region at the Tibetan Plateau. The IPF improved the estimation along with requiring fewer particles (explaining the probability density function of the model by a set of random samples) than EnKF and the standard PF.

Yan et al. [141] analyzed the estimation of SM and hydrologic properties by PF (used AMSR-E) in the Oklahoma Little Washita Watershed. The assimilation methods were PF with commonly used sampling importance resampling (PF-SIR) and Markov chain Monte Carlo sampling (PFMCMC). The applied assimilation methods were overestimated in the dry season and underestimated in the wet season. PF-MCMC showed less error than PF-SIR simulations.

### 2.5.3. Optimal Ensemble Size in EnKF and PF

To conduct EnKF and its variants, it is necessary to determine the optimal ensemble members for calculating the covariance matrices. The determination of ensemble size has been studied in [142–144]. The number of ensemble members and particles that were considered in the analyzed papers are shown in Figure 8 (based on available data). For EnKF (and its variants), most studies used 12 ensembles. In the En4DVAR approach, the numbers vary between 30 and 100. For the PF method, the minimum number was 30 and the maximum was 1000. However, the optimal ensemble size may vary depending on the study and, more importantly, the computation cost.

<figure><img src="image_13.png" /><figcaption>**Figure 8.** Number of ensembles chosen for EnKF, En4DVAR, and PF methods in the reviewed papers (2005–2019).</figcaption></figure>
### 3. The Results of Assimilating BT and SM in LSMs

In the following section, the results of assimilations (for both BT and SM) are categorized based on LSMs considering the papers reviewed in this work. Figures 9 and 10 compare the performance of non-assimilation cases with respect to assimilation cases based on bias (Figure 9A), correlation (Figure 9B), RMSE (Figure 9C), and unbiased RMSE (uRMSE) (Figure 9D) computed between the SM estimated in LSM and ground SM networks for SSM and RSM, respectively. Non-assimilation cases include those bias-corrected/non-bias corrected. In Figure 9A, which shows the bias for NOAH, Jules, and VIC, the SM bias range varies between −0.096 and 0.13 (m3 m-3) for non-assimilation studies and between −0.077 and 0.069 (m3 m-3) for assimilation studies. The NOAH has the greater variability in both cases, with assimilation cases having lower biases. In Figure 9B, the performances of NOAH, CLSM, CLM, Jules, ISBA, SiB2, and VIC are shown for a correlation index. The correlation varies between −0.243 and 0.95 in non-assimilations and between 0.343 and 0.93 for assimilations. The mean of assimilation (0.69) cases is higher than non-assimilation (0.64) cases.

<figure><img src="image_14.png" /><figcaption>**Figure 9.** Assimilation results for SSM: (A) bias (m3 m-3), (B) correlation, (C) RMSE (m3 m-3), and (D) uRMSE (m3 m-3); the metrics are computed based on estimated SM in LSMs and ground networks.</figcaption></figure>

In Figure 9C, the RMSE in NOAH, CLM, Jules, SiB2, and VIC varies between 0.016 and 0.13 (m3 m-3) in non-assimilation cases and between 0.01 and 0.15 (m3 m-3) for assimilation cases. The mean of RMSE (0.058 m3 m-3) is lower in assimilation cases with respect to non-assimilations (0.069 m3 m-3). In Figure 9D, the uRMSE in NOAH and CLSM varies between 0.018 and 0.079 (m3 m-3) in non-assimilations and between 0.012 and 0.075 (m3 m-3) in assimilations. The uRMSE mean of assimilations (0.041 m3 m-3) is lower than non-assimilations (0.046 m3 m-3).

Figure 10 is the same as Figure 9 but regarding RSM. In Figure 10A, the performances of NOAH and Jules in the bias index are shown. The range for non-assimilation cases varies between 0.025 and 0.1 (m3 m-3) and between −0.017 and 0.04 (m3 m-3) for assimilation cases. The mean for assimilation cases (0.005 m3 m-3) is lower compared with non-assimilation cases (0.054 m3 m-3). In Figure 10B, the correlations for NOAH, CLSM, Jules, ISBA, and SiB2 are shown. The ranges vary between 0.25 and 0.89 in non-assimilations and between 0.35 and 0.89 in assimilations. In Figure 10C, the RMSE is shown for NOAH, CLM, Jules, and SiB2, varying between 0.022 and 0.111 (m3 m-3) for non-assimilations and between 0.024 and 0.109 (m3 m-3) for assimilations. In Figure 10D, the uRMSE is shown for NOAH and CLSM. The ranges vary between 0.027 and 0.053 (m3 m-3) for non-assimilations
and between 0.025 and 0.048 (m3 m-3) for assimilations. The uRMSE mean of assimilations (0.036 m3 m-3) is lower than that of non-assimilations (0.040 m3 m-3). Such results (both in SSM and RSM) are shown to demonstrate the quantitative range of error indices in the reviewed papers. These ranges are a function of many factors, such as the accuracy of in situ measurements, study area, the applied sensors, assimilation methods, etc.

<figure><img src="image_15.png" /><figcaption>**Figure 10.** Assimilation results for RSM: (A) bias (m3 m-3), (B) correlation, (C) RMSE (m3 m-3), and (D) uRMSE (m3 m-3); the metrics are computed based on estimated SM in LSMs and ground networks.</figcaption></figure>

To show the effect of the assimilation procedures, the differences between assimilation and non-assimilation cases are reported (sorted based on correlation) for SSM (Figure 11). Some authors conducted multiple approaches, which are shown as slots with different colors in each bar; for example, [29] conducted two assimilation scenarios with different configurations. As Figure 11 shows, most of assimilation approaches improved the SM estimation. In some papers [29,78,97], some configurations resulted in negative correlation differences where only one station was used, which cannot be considered a robust result. Cases providing correlation improvements of more than 0.12 are presented in Table 3.

Most of these scenarios used AMSR-E and SMOS in their assimilation studies. Overall, when BT and SM are used as input, no significant priority can be found between them. The highest improvement was observed in [53] with CLM LSM and the En4DVAR approach in China. Most configurations used EnKF for conducting the assimilation.

Based on Figure 11, RMSE improvements for assimilation cases vary between -0.002 and -0.07 (m3 m-3) (excluding positive values of retrogression). Table 4 shows RMSE improvements in SSM with RMSEs less than -0.024 (m3 m-3). These studies applied AMSR-E and SMOS BT and SM products. Out of ten configurations, five scenarios used the EnKF approach for assimilation.

From Figure 11, bias varies between -0.003 and -0.09 (m3 m-3) (excluding positive values of retrogression). Table 5 demonstrates the bias improvements in SSM with biases less than -0.023 (m3 m-3). These studies used AMSR-E and SMOS and EnKF, which is the most frequent assimilation approach.

In Figure 12, the improvements in correlation, RMSE, and bias are shown for RSM. In most cases, the correlation was improved. The correlation improvements vary between 0.005 and 0.184 (excluding negative values). In [29], a significant negative correlation was found, which is the subtraction of assimilation from open loop (the open loop is initialized with 35 days of assimilation results). It shows that when the Sib2 LSM was initialized with
35 days of assimilation in this study, better results were obtained, and the open-loop run in this configuration had better performance than assimilation in the whole study period. Regarding RMSE, most of the studies reported an improvement between -0.007 and -0.04 (m3 m-3) (excluding positive values). All available studies show bias improvements in RSM by assimilation. In Tables 6–8, the configuration of high score results are shown.

<figure><img src="image_16.png" /><figcaption>**Figure 11.** The difference between metrics in assimilation and non-assimilation scenarios for SSM in terms of correlation, bias, and RMSE; the y axis shows the authors of papers, and the three columns show the metrics.</figcaption></figure>

<figure><img src="image_17.png" /><figcaption>**Figure 12.** The difference between metrics in assimilation and non-assimilation scenarios for RSM: bias (m3 m-3), correlation, RMSE (m3 m-3), and uRMSE (m3 m-3); the y axis shows the authors of papers, and the three columns show the metrics.</figcaption></figure>
<table><caption>**Table 3.** Correlation improvements of more than 0.12 for SSM for papers reviewed in this study.</caption><thead><tr><th>Authors</th><th>LSM</th><th>Sensors</th><th>In Situ Data</th><th>Study Area</th><th>RTM</th><th>Forcing</th><th>Bias Correction</th><th>Approach</th><th>Ensembles</th></tr></thead><tbody><tr><td>[53]</td><td>CLM</td><td>BT-AMSR-E (6.9 GHz) V</td><td>226 field sites</td><td>China</td><td>QH, LandEM, and CMEM</td><td>Fengyun-2C and NCEP-NCAR</td><td>-</td><td>En4DVAR</td><td>60</td></tr><tr><td>[97]</td><td>SiB2</td><td>BT-AMSR-E (6.9, 10.7, and 18.7 GHz) V</td><td>CTP-Naqu</td><td>Tibetan and Mongolian Plateau</td><td>Described in the paper</td><td>CMFD and GLDAS</td><td>-</td><td>LSM calibration</td><td>-</td></tr><tr><td>[121]</td><td>NOAH</td><td>SM-SMOS (S1)</td><td>SCAN</td><td>US</td><td>-</td><td>FNL</td><td>Blind-bias</td><td>1D-Var</td><td>-</td></tr><tr><td>[121]</td><td>NOAH</td><td>SM-SMOS (S2)</td><td>SCAN</td><td>US</td><td>-</td><td>FNL</td><td>Blind-bias</td><td>1D-Var</td><td>-</td></tr><tr><td>[84]</td><td>SiB2</td><td>BT-AMSR-E (18.7 GHz and 6.9 GHz) V</td><td>CEOP/Mongolia site</td><td>Mongolia</td><td>Q-h model</td><td>GLDAS</td><td>-</td><td>PF</td><td>500</td></tr><tr><td>[108]</td><td>VIC</td><td>SM-SMOS-Level 3 CATDS</td><td>OzNet</td><td>Australia</td><td>-</td><td>MERRA</td><td>CDF</td><td>EnKF</td><td>32</td></tr><tr><td>[70]</td><td>VIC</td><td>SM-SMOS-Level 3 CATDS (S1)</td><td>OzNet</td><td>Australia</td><td>CMEM</td><td>MERRA</td><td>CDF</td><td>EnKF</td><td>32</td></tr><tr><td>[70]</td><td>VIC</td><td>BT-SMOS-TB (42.5 degree) H (S2)</td><td>OzNet</td><td>Australia</td><td>CMEM</td><td>MERRA</td><td>CDF</td><td>EnKF</td><td>32</td></tr><tr><td>[70]</td><td>VIC</td><td>BT-SMOS-TB (42.5 degree) HV (S3)</td><td>OzNet</td><td>Australia</td><td>CMEM</td><td>MERRA</td><td>CDF</td><td>EnKF</td><td>32</td></tr><tr><td>[114]</td><td>NOAH</td><td>SM-SMOS-Level 2 (S1)</td><td>NASMD</td><td>US</td><td>-</td><td>NLDAS-2</td><td>-</td><td>EnKF</td><td>32</td></tr><tr><td>[114]</td><td>NOAH</td><td>SM-SMOS-Level 2 (S2)</td><td>NASMD</td><td>US</td><td>-</td><td>NLDAS-2</td><td>CDF (uniform)</td><td>EnKF</td><td>32</td></tr><tr><td>[114]</td><td>NOAH</td><td>SM-SMOS-Level 2 (S3)</td><td>NASMD</td><td>US</td><td>-</td><td>NLDAS-2</td><td>CDF (based on land cover)</td><td>EnKF</td><td>32</td></tr><tr><td>[145]</td><td>CLSM</td><td>Sentinel back scatter and SMAP TB</td><td>SMAPCVS, SCAN, USCRN SMOSMANIA, and OzNet</td><td>Eastern USA and western Europe</td><td>Tau-omega</td><td>GEOS-5</td><td>RTM calibration</td><td>3D-EnKF</td><td>24</td></tr></tbody></table>

The “-” in the table shows that this item was not discussed in the reviewed paper or the paper did not apply this item. “S” character in column of “sensors” refers to different configurations that have been applied in a paper.
<table><caption>**Table 4.** RMSE improvements in SSM with RMSEs of less than −0.024 (m3 m-3) for papers reviewed in this study.</caption><thead><tr><th>Authors</th><th>LSM</th><th>Sensors</th><th>In Situ Data</th><th>Study Area</th><th>RTM</th><th>Forcing</th><th>Bias Correction</th><th>Approach</th><th>Ensembles</th></tr></thead><tbody><tr><td>[146]</td><td>NOAH</td><td>SM-AMSR-E (S1)</td><td>One station</td><td>Tibet Plateau</td><td>-</td><td>-</td><td>-</td><td>Newtonian relaxation</td><td>-</td></tr><tr><td>[97]</td><td>SiB2</td><td>BT-AMSR-E 6.9, 10.7, and 18.7 GHz</td><td>CTP-Naqu</td><td>Tibetan and Mongolian Plateau</td><td>Described in the paper</td><td>CMFD and GLDAS</td><td>-</td><td>LSM calibration</td><td>-</td></tr><tr><td>[146]</td><td>NOAH</td><td>SM-AMSR-E (S2)</td><td>One station</td><td>Tibet Plateau</td><td>-</td><td>-</td><td>-</td><td></td><td></td></tr><tr><td>[147]</td><td>NOAH</td><td>SM-AMSR-E (S1)</td><td>ARS and one SCAN station and SGP</td><td>US</td><td>-</td><td>GDAS</td><td>Mass conservation constraint</td><td>1D-EnKF</td><td>-</td></tr><tr><td>[4]</td><td>SiB2</td><td>BT-AMSR-EE-V 6.925 and 10.65 GHz</td><td>One station</td><td>Tibetan Plateau</td><td>-</td><td>NCEP</td><td>-</td><td>EnKF</td><td>50</td></tr><tr><td>[114]</td><td>NOAH</td><td>SM-SMOS (S1)</td><td>NASMD</td><td>US</td><td>-</td><td>NLDAS-2</td><td>CDF-matching</td><td>EnKF</td><td>32</td></tr><tr><td>[114]</td><td>NOAH</td><td>SM-SMOS (S2)</td><td>NASMD</td><td>US</td><td>-</td><td>NLDAS-2</td><td>CDF-matching</td><td>EnKF</td><td>32</td></tr><tr><td>[114]</td><td>NOAH</td><td>SM-SMOS (S3)</td><td>NASMD</td><td>US</td><td>-</td><td>NLDAS-2</td><td>CDF-matching</td><td>EnKF</td><td>32</td></tr><tr><td>[121]</td><td>NOAH</td><td>SM-SMOS (S1)</td><td>SCAN</td><td>US</td><td>-</td><td>FNL</td><td>Blind-bias</td><td>1D-Var</td><td>-</td></tr><tr><td>[146]</td><td>NOAH</td><td>SM-AMSR-E (S3)</td><td>One station</td><td>Tibet Plateau</td><td>-</td><td>-</td><td>-</td><td>Newtonian relaxation</td><td>-</td></tr><tr><td>[147]</td><td>NOAH</td><td>SM-AMSR-E (S2)</td><td>ARS and one SCAN station and SGP</td><td>US</td><td>-</td><td>GDAS</td><td>Mass conservation constraint</td><td>1D-EnKF</td><td>-</td></tr><tr><td>[113]</td><td>NOAH</td><td>SM-SMAP</td><td>A gauge in Elora, ON, Canada</td><td>US</td><td>-</td><td></td><td>CDF-matching based on soil type</td><td>EnKF</td><td>12</td></tr><tr><td>[90]</td><td>CLM</td><td>BT-AMSR-E-V</td><td>2 stations</td><td>Two stations in China</td><td>LandEM</td><td>Ground observation-based</td><td>-</td><td>EnKF</td><td>-</td></tr><tr><td>[121]</td><td>NOAH</td><td>SM-SMOS (S2)</td><td>SCAN</td><td>US</td><td>-</td><td>FNL</td><td>Blind-bias</td><td>1D-Var</td><td>-</td></tr><tr><td>[84]</td><td>SiB2</td><td>BT-AMSR-E-V 18.7 GHz and 6.9 GHz</td><td>CEOP/Mongolia site</td><td>Mongolia</td><td>Q-h model</td><td>GLDAS</td><td></td><td>PF</td><td>500</td></tr></tbody></table>

The “-” in the table shows that this item was not discussed in the reviewed paper or the paper did not apply this item. “S” character in column of “sensors” refers to different configurations that have been applied in a paper.
<table><caption>**Table 5.** Bias improvements in SSM with biases of less than -0.023 (m3 m-3) for papers reviewed in this study.</caption><thead><tr><th>Authors</th><th>LSM</th><th>Sensors</th><th>In Situ Data</th><th>Study Area</th><th>RTM</th><th>Forcing</th><th>Bias Correction</th><th>Approach</th><th>Ensembles</th></tr></thead><tbody><tr><td>[148]</td><td>NOAH</td><td>SM-AMSR-E</td><td>USDA-ARS (S1)</td><td>US</td><td>-</td><td>NLDAS</td><td>After assimilation, the analysis is optionally bias-corrected</td><td>1D-EnKF</td><td>20</td></tr><tr><td>[147]</td><td>NOAH</td><td>SM-AMSR-E</td><td>ARS and one SCAN station and SGP (S1)</td><td>US</td><td>-</td><td>GDAS</td><td>Mass conservation constraint</td><td>1D-EnKF</td><td>-</td></tr><tr><td>[147]</td><td>NOAH</td><td>SM-AMSR-E</td><td>ARS and one SCAN station and SGP (S2)</td><td>US</td><td>-</td><td>GDAS</td><td>Mass conservation constraint</td><td>1D-EnKF</td><td>-</td></tr><tr><td>[121]</td><td>NOAH</td><td>SM-SMOS</td><td>SCAN (S1)</td><td>US</td><td>-</td><td>FNL</td><td>Blind-bias</td><td>1D-Var</td><td>-</td></tr><tr><td>[148]</td><td>NOAH</td><td>SM-AMSR-E</td><td>USDA-ARS (S2)</td><td>US</td><td>-</td><td>NLDAS</td><td>After assimilation, the analysis is optionally bias-corrected</td><td>3D-EnKF</td><td>20</td></tr><tr><td>[121]</td><td>NOAH</td><td>SM-SMOS</td><td>SCAN (S2)</td><td>US</td><td>-</td><td>FNL</td><td>Blind-bias</td><td>1D-Var</td><td>-</td></tr><tr><td>[114]</td><td>NOAH</td><td>SM-SMOS</td><td>NASMD (S1)</td><td>US</td><td>-</td><td>NLDAS-2</td><td>-</td><td>EnKF</td><td>32</td></tr><tr><td>[114]</td><td>NOAH</td><td>SM-SMOS</td><td>NASMD (S2)</td><td>US</td><td>-</td><td>NLDAS-2</td><td>CDF-matching</td><td>EnKF</td><td>32</td></tr><tr><td>[114]</td><td>NOAH</td><td>SM-SMOS</td><td>NASMD (S3)</td><td>US</td><td>-</td><td>NLDAS-2</td><td>CDF-matching</td><td>EnKF</td><td>32</td></tr></tbody></table>

The “-” in the table shows that this item was not discussed in the reviewed paper or the paper did not apply this item. “S” character in column of “sensors” refers to different configurations that have been applied in a paper.

<table><caption>**Table 6.** RSM correlation improvements greater than 0.09 for papers reviewed in this study.</caption><thead><tr><th>Authors</th><th>LSM</th><th>Sensors</th><th>In Situ Data</th><th>Study Area</th><th>RTM</th><th>Forcing</th><th>Bias Correction</th><th>Approach</th><th>Ensembles</th></tr></thead><tbody><tr><td>[121]</td><td>NOAH</td><td>SM-SMOS (S1)</td><td>SCAN</td><td>US</td><td>-</td><td>FNL</td><td>Blind-bias</td><td>1D-Var</td><td>-</td></tr><tr><td>[121]</td><td>NOAH</td><td>SM-SMOS (S1)</td><td>SCAN and CRN</td><td>US</td><td>-</td><td>FNL</td><td>Estimated the SM bias and removed it</td><td>1D-Var</td><td>-</td></tr><tr><td>[29]</td><td>SiB2</td><td>SM-AMSR2</td><td>One station</td><td>Australia</td><td>Described in the paper</td><td>Ground station and CLVDAS</td><td>-</td><td>Genetic Particle Filter</td><td>512</td></tr></tbody></table>

The “-” in the table shows that this item was not discussed in the reviewed paper or the paper did not apply this item. “S” character in column of “sensors” refers to different configurations that have been applied in a paper.
<table><caption>**Table 7.** RSM RMSE improvements higher than −0.04 (m3 m-3) for papers reviewed in this study.</caption><thead><tr><th>Authors</th><th>LSM</th><th>Sensors</th><th>In Situ Data</th><th>Study Area</th><th>RTM</th><th>Forcing</th><th>Bias Correction</th><th>Approach</th><th>Ensembles</th></tr></thead><tbody><tr><td>[147]</td><td>NOAH</td><td>SM-AMSR-E (S1)</td><td>ARS and one SCAN station and SGP</td><td>US</td><td>-</td><td>GDAS</td><td>Mass conservation constraint</td><td>1D-EnKF</td><td>-</td></tr><tr><td>[149]</td><td>Jules</td><td>SM-SMOS</td><td>OzNet</td><td>Australia</td><td>-</td><td>ACCESS-A</td><td>-</td><td>Evolutionary based on Non-Dominated Sorting Genetic</td><td>20</td></tr><tr><td>[147]</td><td>NOAH</td><td>SM-AMSR-E (S2)</td><td>ARS and one SCAN station and SGP</td><td>SM</td><td>US</td><td>GDAS</td><td>Mass conservation constraint</td><td>1D-EnKF</td><td>-</td></tr><tr><td>[29]</td><td>SiB2</td><td>SM-AMSR2</td><td>One station</td><td>Australia</td><td>Described in the paper</td><td>Ground station and CLVDAS</td><td>-</td><td>Genetic Particle Filter</td><td>512</td></tr><tr><td>[150]</td><td>SiB2</td><td>BT-TMI 10.7 GHz</td><td>MS3608</td><td>Tibetan Plateau</td><td>AIEM</td><td>Station measurement</td><td>-</td><td>EnKF</td><td>100</td></tr></tbody></table>

The “-” in the table shows that this item was not discussed in the reviewed paper or the paper did not apply this item. “S” character in column of “sensors” refers to different configurations that have been applied in a paper.

<table><caption>**Table 8.** RSM Bias improvements of more than −0.07 (m3 m-3) for papers reviewed in this study.</caption><thead><tr><th>Authors</th><th>LSM</th><th>Sensors</th><th>In Situ Data</th><th>Study Area</th><th>RTM</th><th>Forcing</th><th>Bias Correction</th><th>Approach</th><th>Ensembles</th></tr></thead><tbody><tr><td>[147]</td><td>NOAH</td><td>SM-AMSR-E (S1)</td><td>ARS and one SCAN station and SGP</td><td>US</td><td>-</td><td>GDAS</td><td>Mass conservation constraint</td><td>1D-EnKF</td><td>-</td></tr><tr><td>[147]</td><td>NOAH</td><td>SM-AMSR-E (S2)</td><td>ARS and one SCAN station and SGP</td><td>US</td><td>-</td><td>GDAS</td><td>Mass conservation constraint</td><td>1D-EnKF</td><td>-</td></tr><tr><td>[114]</td><td>NOAH</td><td>SM-SMOS (S1)</td><td>NASMD</td><td>US</td><td>-</td><td>NLDAS-2</td><td>-</td><td>EnKF</td><td>32</td></tr><tr><td>[114]</td><td>NOAH</td><td>SM-SMOS (S2)</td><td>NASMD</td><td>US</td><td>-</td><td>NLDAS-2</td><td>CDF-matching</td><td>EnKF</td><td>32</td></tr><tr><td>[121]</td><td>NOAH</td><td>SM-SMOS</td><td>SCAN</td><td>US</td><td>-</td><td>FNL</td><td>Blind-bias</td><td>1D-Var</td><td>-</td></tr><tr><td>[114]</td><td>NOAH</td><td>SM-SMOS (S3)</td><td>NASMD</td><td>US</td><td>-</td><td>NLDAS-2</td><td>CDF-matching</td><td>EnKF</td><td>32</td></tr><tr><td>[114]</td><td>NOAH</td><td>SM-SMOS (S4)</td><td>NASMD</td><td>US</td><td>-</td><td>NLDAS-2</td><td>CDF-matching</td><td>EnKF</td><td>32</td></tr></tbody></table>

The “-” in the table shows that this item was not discussed in the reviewed paper or the paper did not apply this item. “S” character in column of “sensors” refers to different configurations that have been applied in a paper.
Table 6 displays the configurations of scenarios for RSM with correlation improvements greater than 0.09. These studies applied SMOS and AMSR-E SM data and 1D-VAR and PF assimilation methods. Table 7 shows RMSE improvements higher than -0.04 (m3 m-3). AMSR-E, SMOS, and TMI data were applied and EnKF, GPF, and evolutionary methods were used for conducting the assimilations.

Table 8 shows the bias improvement in RSM for configurations with improvements of more than -0.07 (m3 m-3). They applied SMOS and AMSR-E SM and EnKF (most of them) and 1D-VAR as assimilation approaches.

Based on Figures 10 and 12 and experimental and physical conclusions in recent studies, the assimilation results show lower improvements in RSM. SM is highly dependent on the LSM structure (subsurface physics) and water balance modeling [151] along with the lag time (especially in wetter condition) that are required for the correction of RSM by assimilation [113].

# 4. Discussion

This study reviewed the studies that analyzed the improvement of LSM SM estimates by assimilation of satellite-derived SM and BTs. The details in assimilation and improvement of LSMs SM estimates through LSM and RTM calibration were presented and, finally, the results of reviewed studies were investigated. It should be noted that the assimilation results and the conclusions of the papers reviewed here could be spatially and temporally variable. It is necessary to analyze the outcome of high score assimilation configurations in analyzed studies in different spatial and temporal conditions. Based on the analyzed papers, along with the mentioned parameters (input data, bias correction, and assimilation details) in the last sections, the following considerations could be applied to gain better results in LSM SM correction in the future:

## a Joint Use of Calibration and Assimilation

The calibration of LSM and RTM is a crucial approach, along with assimilation, for reducing biases. The calibration directly affects model climatology and the assimilation corrects the forcing data (precipitation) for driving the model. During calibration, model parameters (both fixed and time-variant) could be problematic, and the application of both calibration and assimilation may reduce such uncertainties. Zhao et al. [77] proposed pre-calibration of RTM parameters for gaining better results in both SSM and RSM estimation through assimilation. Tian et al. [26] applied a dual-pass assimilation system to simultaneously update model states and RTM parameters. Koster et al. [27] emphasized the effects of model calibration and assimilation as two separate methods for correcting model parameterization and errors induced by forcing data. They proposed that the joint use of these approaches could have complementary effects and result in more accurate surface SM estimates. However, it should be noted that such a simultaneous estimation is challenging when rapid changes are happening on vegetation opacity and soil surface roughness; this might reduce the compensatory role of parameters on each other [96].

## b The Role of Precipitation, Vegetation, and Land Cover

Recent studies show that assimilation results are dependent on precipitation forcing, vegetation, and land cover. The role of precipitation, wetness, temperature, and vegetation were investigated on assimilation results. Lin et al. [152] and Liu et al. [52] suggested that the assimilation of both SM and precipitation could provide better results. They stated that as both methods have independent corrections, their combination can have the largest impact on SM skills. Moreover, assimilation could be more beneficial in places where precipitation corrections are not possible [28,153].

As for vegetation cover, Draper et al. [154] showed that significant improvements were found over croplands and mixed cover for the root zone. Furthermore, Tian et al. [155] found better results over croplands, grasslands, and mixed cover for SSM. It was proposed that SM improvement in dry to moderate wetness conditions was higher in some studies [116,153]. Additionally, [156] reported that due to the presence of frozen soil (low
temperature), high vegetation cover and noise in dry regions resulted in a reduction in the accuracy of satellite retrieval, and larger improvements could be observed in areas with moderate vegetation cover.

### c The Joint Assimilation of Satellite Data

It is vital to use the best data; Gevaert et al. [153] stated that the joint assimilation of L-and C-band (or X-band) retrievals had higher results in contrast with individually assimilating C-band (or X-band) retrievals in AWRA-L LSM. Draper et al. [154] used ASCAT and AMSR-E SM products individually and in combination. They stated that when these two active and passive data are used together, better results are obtained. Yin et al. [68] used SMOPS SM products and showed that they had superior results compared to ASCATA, ASCATB, SMAP, SMOS, and AMSR2. However, Kolassa et al. [157] conducted assimilation experiments from active and passive SM products in CLSM LSM. They used the retrieved SM from both active (ASCAT) and passive (AMSR-E) platforms and found similarly improved results in the surface from both experiments in contrast with individual products. As shown in Tables 3–8 (based on available data), the highest improvements may be gained by the SMOS and AMSR-E sensors.

The application of other satellite products can also be a beneficial solution. Gavahi et al. [158] applied MODIS evaporation (MODIS16 ET) along with SMOPS SM product in an assimilation study and found improved results compared to a univariate assimilation procedure. According to the abovementioned studies, the joint use of satellite products may provide better results in the future.

### d Selection of Assimilation Method

The selection of the assimilation method is highly dependent on the nature of the problem and observation data [159]. According to Moradkhani et al. [159], the following factors should be considered. (1) The linearity and non-linearity of the LSM and updating procedure. As the SM probability distribution in LSMs is realized as a non-linear process, it is expected that a non-linear assimilation assumption could provide more plausible results. (2) Filtering and smoothing updating. It can be more beneficial in SM estimation to include backward observations simultaneously (4D-VAR). (3) The deterministic (KF, variation EKF) and ensemble-based assimilation procedures (EnKF and PF). The deterministic approaches could result in poor results in complex LSMs. (4) The computation cost. The required runtime in frequently used assimilation methods (Variational, EnKF, and PF) could be a challenge, especially for operational applications.

As shown in Section 3, most of the reviewed studies with the highest improvements used EnKF and variational approaches, respectively. It should be noted that in EnKF, the assumption of the non-Gaussian error structure in state variables and observations and the linear updating procedure can have negative impacts on assimilation performance [160–163]. Although variational methods are not bound to solve state-space modeling, they need the derivation of an adjoin model (linearized form of the model) to optimize the error minimization formulas. Moreover, the assumption of time-invariant model covariance could result in inherent errors in SM improvements [163]. The PF can be a solution when skewed model ensembles and non-Gaussian error structures are observed. In contrary to EnKF, where model state variables are updated at each time step, in PF, the probability distributions of state variables and observations evolve through time. Although PF was proposed to overcome the suboptimal problems of EnKF, it suffers from computational burden and an inherent problem called “particle degeneracy” when the weights of particles are not considerable, which causes inaccurate posteriori distribution estimates. To solve this problem, several methods were proposed (the details can be tracked in [162]).

As shown in Section 2.5.2, some studies compared the performance of assimilation methods regarding the improvement of SM in LSMs. The results from these studies may support the theoretical conclusions above that EKF (dealing with non-linearity) had better performance than OI [5] and that EnKF had better performance than EKF [138,139,164]
and EnOI (it is cheaper but does not use the flow-dependent ensembles) [51]. EDA and PF (proposed IPF) had better results than EnKF [137,165].

Although the internal setting in the reviewed papers in this study are different and the number of studies using PF is low, in Table 3 (correlation difference in SSM), PF had better results than EnKF, but in Table 4 (RMSE difference in SSM), the PF had lower results than EnKF. This can be an issue for future studies in which evaluation metrics for PF have better performance than EnKF.

## e Selection of Assimilation Period

Draper et al. [165] stated that periods longer than 31 days are essential for achieving better results from assimilation in the root zone; SM anomalies can take between 1 and 3 months [166,167]. Draper et al. [168] showed that the sub-seasonal scale had the largest impact on SM estimation. It is necessary that future studies investigate the role of the assimilation period while considering diverse spatial and temporal conditions.

## f Satellite BTs and SMs

Kolassa, et al. [111] (2017) showed that, on average, using both BT and SM retrieval in assimilation has similar results (excluding local effects). De Lannoy et al. [169] showed that BT assimilation has slightly better results than SM products. They stated that this may show that BT could contain more information than the filtered information within SM retrievals. However, the results from our analysis in Section 3 showed no significantly greater improvements in BT with respect to the SM products.

## g The Spatial Resolution of Satellite Data

According to Peng et al. [170], the application of high spatial resolution SM products (less than 1 km) is highly demanding in the hydrometeorological community. Promising projects such as the Sat'elite Argentino de Observacion' COn Microondas (SAOCOM) mission, the NASAISRO Synthetic Aperture Radar (NISAR), the Radar Observing System for Europe L (ROSE-L), and the Tandem-L satellites could provide such required resolutions, but at the present time, ESA Sentinel-1 European Radar Observatory can provide the highest resolution. Moreover, recently proposed downscaling algorithms could be applied to acquire better results in mid to low spatial scales using atmospheric and geophysical data from high spatial resolution remote-sensing data from other platforms [69,171].

**Author Contributions:** Conceptualization, R.K.; methodology, S.B.; software, M.G.; validation, J.-P.W., A.P.B. and S.B.; formal analysis, A.P.B.; investigation, J.-P.W.; resources, R.K.; data curation, S.B.; writing—original draft preparation, R.K.; writing—review and editing, R.K.; visualization, M.G.; supervision, A.P.B.; project administration, S.B. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research was not funded by any source.

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** Not applicable.

**Conflicts of Interest:** The authors declare no conflict of interest.

## Abbreviations

The full form of applied acronyms in this study are shown in the following:

AACES Australian Airborne Cal/Val Experiments for SMOS-AACES

AARD Alberta Agriculture and Rural Development

ACCESS-A Australian Community Climate Earth-System Simulato

AGCD Australian Gridded Climate Data

AGDMN Alberta Ground Drought Monitoring Network
<table>
  <tr>
    <td>ALADIN</td>
    <td>Me'te'o-France's Aire Limite'e Adaptation Dynamique de'veloppement InterNational</td>
  </tr>
  <tr>
    <td>ALEXI</td>
    <td>Atmosphere Land Exchange Inverse</td>
  </tr>
  <tr>
    <td>AMSR2</td>
    <td>Advanced Microwave Scanning Radiometer 2</td>
  </tr>
  <tr>
    <td>AMSR-E</td>
    <td>Advanced Microwave Scanning Radiometer-Earth Observing System</td>
  </tr>
  <tr>
    <td>ARS</td>
    <td>Agricultural Research Service</td>
  </tr>
  <tr>
    <td>ASCAT</td>
    <td>Advanced Scatterometer</td>
  </tr>
  <tr>
    <td>AWDN</td>
    <td>Automated Weather Data Network</td>
  </tr>
  <tr>
    <td>AWRA-L</td>
    <td>Australian Water Resources Assessment landscape hydrology model</td>
  </tr>
  <tr>
    <td>BAWAP</td>
    <td>Bureau of Meteorology-Australian Water Availability Project</td>
  </tr>
  <tr>
    <td>CAM 4.0</td>
    <td>Community Atmosphere Mode</td>
  </tr>
  <tr>
    <td>caPA</td>
    <td>Canadian Precipitation Analysis</td>
  </tr>
  <tr>
    <td>CCI</td>
    <td>Climate Change Initiative</td>
  </tr>
  <tr>
    <td>CERES</td>
    <td>Clouds and the Earth's Radiant Energy System</td>
  </tr>
  <tr>
    <td>CLDAS</td>
    <td>CMA Land Data Assimilation System</td>
  </tr>
  <tr>
    <td>CLM</td>
    <td>Community Land Model</td>
  </tr>
  <tr>
    <td>CLSM</td>
    <td>Catchment Land Surface Model</td>
  </tr>
  <tr>
    <td>CMFD</td>
    <td>China Meteorological Forcing Dataset</td>
  </tr>
  <tr>
    <td>COSMOZ</td>
    <td>Cosmic ray Soil Moisture Observing System</td>
  </tr>
  <tr>
    <td>DHSVM</td>
    <td>Distributed Hydrology Soil Vegetation Model (DHSVM)</td>
  </tr>
  <tr>
    <td>EAKF</td>
    <td>Ensemble Adjustment Kalman Filter</td>
  </tr>
  <tr>
    <td>EBKS</td>
    <td>Ensemble-Based Kalman smoother</td>
  </tr>
  <tr>
    <td>ECCC</td>
    <td>Environment and Climate Change Canada</td>
  </tr>
  <tr>
    <td>ECMWF</td>
    <td>European Centre for Medium-Range Weather Forecasts (ECMWF)</td>
  </tr>
  <tr>
    <td>EDA</td>
    <td>Evolutionary Data Assimilation</td>
  </tr>
  <tr>
    <td>EKF</td>
    <td>Extende kalman Filter</td>
  </tr>
  <tr>
    <td>En4DVAR</td>
    <td>Ensemble Based Four-Dimensional Variational</td>
  </tr>
  <tr>
    <td>EnKF</td>
    <td>Ensemble Kalman Filter</td>
  </tr>
  <tr>
    <td>EnOI</td>
    <td>Ensemble Optimal Interpolation</td>
  </tr>
  <tr>
    <td>ERA</td>
    <td>ECMWF Re-Analysis</td>
  </tr>
  <tr>
    <td>ESRF</td>
    <td>Ensemble Square Root Filter</td>
  </tr>
  <tr>
    <td>ETKF</td>
    <td>Ensemble Transform Kalman Filter</td>
  </tr>
  <tr>
    <td>EuropeNet</td>
    <td>European Network</td>
  </tr>
  <tr>
    <td>Evol</td>
    <td>Evolutionary</td>
  </tr>
  <tr>
    <td>FCRN</td>
    <td>Fluxnet-Canada Research Network (FCRN)</td>
  </tr>
  <tr>
    <td>FluxNet (OzFlux)</td>
    <td>Australian and New Zealand Flux Research and Monitoring Network</td>
  </tr>
  <tr>
    <td>FNL</td>
    <td>NCEP Final Analysis</td>
  </tr>
  <tr>
    <td>GDAS</td>
    <td>Global Data Assimilation System</td>
  </tr>
  <tr>
    <td>GEM</td>
    <td>Global Environmental Multi-scale Model</td>
  </tr>
  <tr>
    <td>GEOS-5</td>
    <td>Goddard Earth Observing System Model</td>
  </tr>
  <tr>
    <td>GLDAS</td>
    <td>Global Land Data Assimilation System</td>
  </tr>
  <tr>
    <td>GLEAM</td>
    <td>Global Land Evaporation Amsterdam Model</td>
  </tr>
  <tr>
    <td>GPCP</td>
    <td>Global Precipitation Climatology Project</td>
  </tr>
  <tr>
    <td>GPF</td>
    <td>Genetic Particle Filter</td>
  </tr>
  <tr>
    <td>GRACE</td>
    <td>Gravity Recovery and Climate Experiment</td>
  </tr>
  <tr>
    <td>GSI</td>
    <td>Gridpoint Statistical Interpolation System</td>
  </tr>
  <tr>
    <td>GSMDB</td>
    <td>Global Soil Moisture Data Bank</td>
  </tr>
  <tr>
    <td>HTESSEL</td>
    <td>Hydrology Tiled ECMWF scheme of Surface Exchanges over land</td>
  </tr>
  <tr>
    <td>ICN</td>
    <td>Illinois Climate Network</td>
  </tr>
  <tr>
    <td>IPF</td>
    <td>Improved PF</td>
  </tr>
  <tr>
    <td>ISBA</td>
    <td>Interactions between Soil, Biosphere, and Atmosphere</td>
  </tr>
  <tr>
    <td>ISCCP</td>
    <td>International Satellite Cloud Climatology Project</td>
  </tr>
  <tr>
    <td>ISMN</td>
    <td>International Soil Moisture Network</td>
  </tr>
  <tr>
    <td>JULES</td>
    <td>Joint UK Land Environment Simulator</td>
  </tr>
</table>
<table><tr><td>KF</td><td>Kalman Filter</td></tr><tr><td>MAWN</td><td>Michigan Automated Weather Network</td></tr><tr><td>MERRA</td><td>Modern-Era Retrospective analysis for Research and Applications</td></tr><tr><td>MESH</td><td>Modélisation Environmentale-Surface et Hydrologie</td></tr><tr><td>Mesonet</td><td>Mesoscale and network</td></tr><tr><td>MongoliaNet</td><td>Mongolia network</td></tr><tr><td>MSMMN</td><td>Murrumbidgee Soil Moisture Monitoring Network</td></tr><tr><td>MVIRI</td><td>Meteosat Visible Infra-Red Imager</td></tr><tr><td>NASMD</td><td>North American Soil Moisture Database</td></tr><tr><td>NCEP-NCAR</td><td>National Centers for Environmental Prediction-National Center for Atmospheric Research</td></tr><tr><td>NLDAS</td><td>North American Land Data Assimilation System</td></tr><tr><td>Nudg</td><td>Nudging</td></tr><tr><td>OI</td><td>Optimal Interpolation</td></tr><tr><td>OK Mesonet</td><td>Oklahoma Mesonet</td></tr><tr><td>Oznet</td><td>Australian monitoring network for soil moisture and micrometeorology</td></tr><tr><td>PF</td><td>Particle Filter</td></tr><tr><td>PFMCMC</td><td>Markov chain Monte Carlo sampling</td></tr><tr><td>PF-SIR</td><td>PF with commonly used sampling importance resampling</td></tr><tr><td>SAFRAN</td><td>Système d'Analyse Fournissant des Renseignements Atmosphériques à la Neige</td></tr><tr><td>SCAN</td><td>Soil Climate Analysis Network</td></tr><tr><td>SEKF</td><td>Simplified Extended Kalman Filter</td></tr><tr><td>SiB2</td><td>Simple Biosphere Model 2</td></tr><tr><td>SMAP</td><td>Soil Moisture Active and Passive</td></tr><tr><td>SMAPCVS</td><td>SMAP core validation sites</td></tr><tr><td>SMMR</td><td>Scanning Multichannel Microwave Radiometer</td></tr><tr><td>SMOPS</td><td>Soil Moisture Operational Products System</td></tr><tr><td>SMOSMANIA</td><td>Soil Moisture Observing System-Meteorological Automatic Network Integrated Application</td></tr><tr><td>SMTMN</td><td>Soil Moisture and Temperature Monitoring Network on the central Tibetan Plateau</td></tr><tr><td>SVS</td><td>Soil, Vegetation, and Snow</td></tr><tr><td>TESSEL</td><td>Tiled ECMWF Scheme for Surface Exchanges over Land</td></tr><tr><td>TMPA</td><td>TRMM Multisatellite Precipitation Analysis</td></tr><tr><td>TRMM</td><td>Tropical Rainfall Measuring Mission</td></tr><tr><td>USCRN</td><td>US Climate Reference Network</td></tr><tr><td>USDA-ARS</td><td>United State Agricultural Research Service</td></tr><tr><td>VIC</td><td>Variable Infiltration Capacity</td></tr><tr><td>W3</td><td>World Wide Water</td></tr><tr><td>WFDEI</td><td>WATCH Forcing Data methodology applied to ERA-Interim reanalysis data</td></tr></table>

## References

1. Loew, A. Impact of surface heterogeneity on surface soil moisture retrievals from passive microwave data at the regional scale: The Upper Danube case. *Remote Sens. Environ.* **2008**, *112*, 231–248. [CrossRef]

2. Seneviratne, S.I.; Corti, T.; Davin, E.L.; Hirschi, M.; Jaeger, E.B.; Lehner, I.; Orlowsky, B.; Teuling, A.J. Investigating soil moisture-climate interactions in a changing climate: A review. *Earth-Sci. Rev.* **2010**, *99*, 125–161. [CrossRef]

3. van der Molen, M.K.; Dolman, A.J.; Ciais, P.; Eglin, T.; Gobron, N.; Law, B.E.; Meir, P.; Peters, W.; Phillips, O.L.; Reichstein, M.; et al. Drought and ecosystem carbon cycling. *Agric. For. Meteorol.* **2011**, *151*, 765–773. [CrossRef]

4. Rasmy, M.; Koike, T.; Boussetta, S.; Lu, H.; Li, X. Development of a satellite land data assimilation system coupled with a mesoscale model in the Tibetan Plateau. *IEEE Trans. Geosci. Remote Sens.* **2011**, *49*, 2847–2862. [CrossRef]

5. De Rosnay, P.; Drusch, M.; Vasiljevic, D.; Balsamo, G.; Albergel, C.; Isaksen, L. A simplified Extended Kalman Filter for the global operational soil moisture analysis at ECMWF. *Q. J. R. Meteorol. Soc.* **2013**, *139*, 1199–1213. [CrossRef]

6. Mahfouf, J.F.; Bergaoui, K.; Draper, C.; Bouyssel, F.; Taillefer, F.; Taseva, L. A comparison of two off-line soil analysis schemes for assimilation of screen level observations. *J. Geophys. Res. Atmos.* **2009**, *114*, D08105. [CrossRef]

7. de Rosnay, P.; Balsamo, G.; Albergel, C.; Muñoz-Sabater, J.; Isaksen, L. Initialisation of land surface variables for numerical weather prediction. *Surv. Geophys.* **2014**, *35*, 607–621. [CrossRef]

8. Petropoulos, G.P.; Ireland, G.; Barrett, B. Surface soil moisture retrievals from remote sensing: Current status, products & future trends. *Phys. Chem. Earth Parts A/B/C* **2015**, *83*, 36–56.
9. Entekhabi, D.; Njoku, E.G.; O'Neill, P.E.; Kellogg, K.H.; Crow, W.T.; Edelstein, W.N.; Entin, J.K.; Goodman, S.D.; Jackson, T.J.; Johnson, J.; et al. The soil moisture active passive (SMAP) mission. *Proc. IEEE* **2010**, *98*, 704–716. [CrossRef]

10. Kerr, Y.H.; Waldteufel, P.; Richaume, P.; Wigneron, J.P.; Ferrazzoli, P.; Mahmoodi, A.; Al Bitar, A.; Cabot, F.; Gruhier, C.; Juglea, S.E.; et al. The SMOS soil moisture retrieval algorithm. *IEEE Trans. Geosci. Remote Sens.* **2012**, *50*, 1384–1403. [CrossRef]

11. Njoku, E.G.; Jackson, T.J.; Lakshmi, V.; Chan, T.K.; Nghiem, S.V. Soil moisture retrieval from AMSR-E. *IEEE Trans. Geosci. Remote Sens.* **2003**, *41*, 215–229. [CrossRef]

12. Imaoka, K.; Kachi, M.; Kasahara, M.; Ito, N.; Nakagawa, K.; Oki, T. Instrument performance and calibration of AMSR-E and AMSR2. *Int. Arch. Photogramm. Remote Sens. Spat. Inf. Sci.* **2010**, *38*, 13–18.

13. Wagner, W.; Hahn, S.; Kidd, R.; Melzer, T.; Bartalis, Z.; Hasenauer, S.; Figa-Saldana, J.; De Rosnay, P.; Jann, A.; Schneider, S.; et al. The ASCAT soil moisture product: A review of its specifications, validation results, and emerging applications. *Meteorol. Z.* **2013**, *22*, 5–33. [CrossRef]

14. Paloscia, S.; Pettinato, S.; Santi, E.; Notarnicola, C.; Pasolli, L.; Reppucci, A. Soil moisture mapping using Sentinel-1 images: Algorithm and preliminary validation. *Remote Sens. Environ.* **2013**, *134*, 234–248. [CrossRef]

15. Njoku, E.G.; Stacey, J.; Barath, F.T. The Seasat scanning multichannel microwave radiometer (SMMR): Instrument description and performance. *IEEE J. Ocean. Eng.* **1980**, *5*, 100–115. [CrossRef]

16. Gao, H.; Wood, E.F.; Jackson, T.; Drusch, M.; Bindlish, R. Using TRMM/TMI to retrieve surface soil moisture over the southern United States from 1998 to 2002. *J. Hydrometeorol.* **2006**, *7*, 23–38. [CrossRef]

17. Wagner, W.; Lemoine, G.; Rott, H. A method for estimating soil moisture from ERS scatterometer and soil data. *Remote Sens. Environ.* **1999**, *70*, 191–207. [CrossRef]

18. Dorigo, W.; Wagner, W.; Albergel, C.; Albrecht, F.; Balsamo, G.; Brocca, L.; Chung, D.; Ertl, M.; Forkel, M.; Gruber, A.; et al. ESA CCI Soil Moisture for improved Earth system understanding: State-of-the art and future directions. *Remote Sens. Environ.* **2017**, *203*, 185–215. [CrossRef]

19. Liu, J.; Zhan, X.; Hain, C.; Yin, J.; Fang, L.; Li, Z.; Zhao, L. NOAA soil moisture operational product system (SMOPS) and its validations. In *Proceedings of the 2016 IEEE International Geoscience and Remote Sensing Symposium (IGARSS)*, Beijing, China, 10–15 July 2016; pp. 3477–3480.

20. De Lannoy, G.J.M.; de Rosnay, P.; Reichle, R.H. Soil moisture data assimilation. In *Handbook of Hydrometeorological Ensemble Forecasting*; Springer: Berlin/Heidelberg, Germany, 2016; pp. 1–43.

21. De Lannoy, G.J.; Reichle, R.H.; Vrugt, J.A. Uncertainty quantification of GEOS-5 L-band radiative transfer model parameters using Bayesian inference and SMOS observations. *Remote Sens. Environ.* **2014**, *148*, 146–157. [CrossRef]

22. Maggioni, V.; Houser, P.R. Soil moisture data assimilation. In *Data Assimilation for Atmospheric, Oceanic and Hydrologic Applications*; Springer: Berlin/Heidelberg, Germany, 2017; Volume 3, pp. 195–217.

23. Montzka, C. Soil Moisture Remote Sensing and Data Assimilation. In *Remote Sensing of Energy Fluxes and Soil Moisture Content*; CRC Press: Boca Raton, FL, USA, 2013; pp. 415–434.

24. Lu, H.; Koike, T.; Yang, K.; Hu, Z.; Xu, X.; Rasmy, M.; Kuria, D.; Tamagawa, K. Improving land surface soil moisture and energy flux simulations over the Tibetan plateau by the assimilation of the microwave remote sensing data and the GCM output into a land surface model. *Int. J. Appl. Earth Obs. Geoinf.* **2012**, *17*, 43–54. [CrossRef]

25. Zhao, L.; Yang, K.; Qin, J.; Chen, Y.; Tang, W.; Lu, H.; Yang, Z.-L. The scale-dependence of SMOS soil moisture accuracy and its improvement through land data assimilation in the central Tibetan Plateau. *Remote Sens. Environ.* **2014**, *152*, 345–355. [CrossRef]

26. Tian, X.; Xie, Z.; Dai, A.; Shi, C.; Jia, B.; Chen, F.; Yang, K. A dual-pass variational data assimilation framework for estimating soil moisture profiles from AMSR-E microwave brightness temperature. *J. Geophys. Res. Atmos.* **2009**, *114*, D16102. [CrossRef]

27. Koster, R.D.; Liu, Q.; Mahanama, S.P.; Reichle, R.H. Improved hydrological simulation using SMAP data: Relative impacts of model calibration and data assimilation. *J. Hydrometeorol.* **2018**, *19*, 727–741. [CrossRef]

28. De Lannoy, G.J.; Reichle, R.H. Global assimilation of multiangle and multipolarization SMOS brightness temperature observations into the GEOS-5 catchment land surface model for soil moisture estimation. *J. Hydrometeorol.* **2016**, *17*, 669–691. [CrossRef]

29. Sawada, Y.; Koike, T.; Walker, J.P. A land data assimilation system for simultaneous simulation of soil moisture and vegetation dynamics. *J. Geophys. Res. Atmos.* **2015**, *120*, 5910–5930. [CrossRef]

30. Mitchell, K. *The Community Noah Land-Surface Model (LSM)*; NOAA/NCEP/Environmental Modeling Center, 2005. Available online: https://ftp.emc.ncep.noaa.gov/mmb/gcp/ldas/noahlsm/ver_2.7.1/ (accessed on 2 February 2022).

31. Niu, G.Y.; Yang, Z.L.; Mitchell, K.E.; Chen, F.; Ek, M.B.; Barlage, M.; Kumar, A.; Manning, K.; Niyogi, D.; Rosero, E. The community Noah land surface model with multiparameterization options (Noah-MP): 1. Model description and evaluation with local-scale measurements. *J. Geophys. Res. Atmos.* **2011**, *116*, D12109. [CrossRef]

32. Van Dijk, A.; Warren, G. *The Australian Water Resources Assessment System*; Technical Report 3, Version 0.5; Australian Government Bureau of Meteorology: Melbourne, Australia, 2010.

33. Koster, R.D.; Suarez, M.J.; Ducharne, A.; Stieglitz, M.; Kumar, P. A catchment-based approach to modeling land surface processes in a general circulation model: 1. Model structure. *J. Geophys. Res. Atmos.* **2000**, *105*, 24809–24822. [CrossRef]

34. Lawrence, D.M.; Fisher, R.A.; Koven, C.D.; Oleson, K.W.; Swenson, S.C.; Bonan, G.; Collier, N.; Ghimire, B.; van Kampenhout, L.; Kennedy, D.; et al. The Community Land Model version 5: Description of new features, benchmarking, and impact of forcing uncertainty. *J. Adv. Model. Earth Syst.* **2019**, *11*, 4245–4287. [CrossRef]
35. Best, M.; Pryor, M.; Clark, D.; Rooney, G.; Essery, R.; Ménard, C.; Edwards, J.; Hendry, M.; Porson, A.; Gedney, N. The Joint UK Land Environment Simulator (JULES), model description—Part 1: Energy and water fluxes. *Geosci. Model Dev.* **2011**, *4*, 677–699.
[CrossRef]

36. Decharme, B.; Boone, A.; Delire, C.; Noilhan, J. Local evaluation of the Interaction between Soil Biosphere Atmosphere soil multilayer diffusion scheme using four pedotransfer functions. *J. Geophys. Res. Atmos.* **2011**, *116*, D20126. [CrossRef]

37. Sellers, P.; Randall, D.; Collatz, G.; Berry, J.; Field, C.; Dazlich, D.; Zhang, C.; Collelo, G.; Bounoua, L. A revised land surface parameterization (SiB2) for atmospheric GCMs. Part I: Model formulation. *J. Clim.* **1996**, *9*, 676–705. [CrossRef]

38. Liang, X.; Lettenmaier, D.P.; Wood, E.F.; Burges, S.J. A simple hydrologically based model of land surface water and energy fluxes for general circulation models. *J. Geophys. Res. Atmos.* **1994**, *99*, 14415–14428. [CrossRef]

39. Alavi, N.; Bélair, S.; Fortin, V.; Zhang, S.; Husain, S.Z.; Carrera, M.L.; Abrahamowicz, M. Warm season evaluation of soil moisture prediction in the Soil, Vegetation, and Snow (SVS) scheme. *J. Hydrometeorol.* **2016**, *17*, 2315–2332. [CrossRef]

40. Van Dijk, A.I.; Peña-Arancibia, J.L.; Wood, E.F.; Sheffield, J.; Beck, H.E. Global analysis of seasonal streamflow predictability using an ensemble prediction system and observations from 6192 small catchments worldwide. *Water Resour. Res.* **2013**, *49*, 2729–2746.
[CrossRef]

41. Miralles, D.; De Jeu, R.; Gash, J.; Holmes, T.; Dolman, A. Magnitude and variability of land evaporation and its components at the global scale. *Hydrol. Earth Syst. Sci.* **2011**, *15*, 967–981. [CrossRef]

42. Wigmosta, M.S.; Nijssen, B.; Storck, P.; Lettenmaier, D. The distributed hydrology soil vegetation model. In *Mathematical Models of Small Watershed Hydrology and Applications*; Water Resources Publications: Littleton, CO, USA, 2002; pp. 7–42.

43. Simunek, J.; Van Genuchten, M.T.; Sejna, M. The HYDRUS-1D software package for simulating the one-dimensional movement of water, heat, and multiple solutes in variably-saturated media. *Univ. Calif. Riverside Res. Rep.* **2005**, *3*, 1–240.

44. Pietroniro, A.; Fortin, V.; Kouwen, N.; Neal, C.; Turcotte, R.; Davison, B.; Verseghy, D.; Soulis, E.; Caldwell, R.; Evora, N.; et al. Development of the MESH modelling system for hydrological ensemble forecasting of the Laurentian Great Lakes at the regional scale. *Hydrol. Earth Syst. Sci.* **2007**, *11*, 1279–1294. [CrossRef]

45. Shuttleworth, W.J.; Wallace, J. Evaporation from sparse crops-an energy combination theory. *Q. J. R. Meteorol. Soc.* **1985**, *111*, 839–855.
[CrossRef]

46. Balsamo, G.; Beljaars, A.; Scipal, K.; Viterbo, P.; van den Hurk, B.; Hirschi, M.; Betts, A.K. A revised hydrology for the ECMWF model: Verification from field site to terrestrial water storage and impact in the Integrated Forecast System. *J. Hydrometeorol.* **2009**, *10*, 623–643. [CrossRef]

47. Essery, R.; Best, M.; Cox, P. MOSES 2.2 Technical Documentation (Hadley Centre Technical Note 30); Hadley Centre, Met Office: Exeter, UK, 2001.

48. Chen, F.; Dudhia, J. Coupling an advanced land surface-hydrology model with the Penn State-NCAR MM5 modeling system. Part I: Model implementation and sensitivity. *Mon. Weather Rev.* **2001**, *129*, 569–585. [CrossRef]

49. Tian, X.; Xie, Z.; Dai, A.; Jia, B.; Shi, C. A microwave land data assimilation system: Scheme and preliminary evaluation over China. *J. Geophys. Res. Atmos.* **2010**, *115*, D21113. [CrossRef]

50. Ramadevi, G.; Parvathi, S.; Kumaresan, A.; Vijayakumar, K. Modern-era retrospective analysis for research and applications. *Adv. Nat. Appl. Sci.* **2017**, *11*, 109–113.

51. Blyverket, J.; Hamer, P.D.; Bertino, L.; Albergel, C.; Fairbairn, D.; Lahoz, W.A. An Evaluation of the EnKF vs. EnOI and the Assimilation of SMAP, SMOS and ESA CCI Soil Moisture Data over the Contiguous US. *Remote Sens.* **2019**, *11*, 478. [CrossRef]

52. Liu, Q.; Reichle, R.H.; Bindlish, R.; Cosh, M.H.; Crow, W.T.; de Jeu, R.; De Lannoy, G.J.; Huffman, G.J.; Jackson, T.J. The contributions of precipitation and soil moisture observations to the skill of soil moisture estimates in a land data assimilation system. *J. Hydrometeorol.* **2011**, *12*, 750–765. [CrossRef]

53. Jia, B.; Tian, X.; Xie, Z.; Liu, J.; Shi, C. Assimilation of microwave brightness temperature in a land data assimilation system with multi-observation operators. *J. Geophys. Res. Atmos.* **2013**, *118*, 3972–3985. [CrossRef]

54. Wigneron, J.-P.; Jackson, T.; O’neill, P.; De Lannoy, G.; De Rosnay, P.; Walker, J.; Ferrazzoli, P.; Mironov, V.; Bircher, S.; Grant, J.; et al. Modelling the passive microwave signature from land surfaces: A review of recent results and application to the L-band SMOS & SMAP soil moisture retrieval algorithms. *Remote Sens. Environ.* **2017**, *192*, 238–262.

55. Karthikeyan, L.; Pan, M.; Wanders, N.; Kumar, D.N.; Wood, E.F. Four decades of microwave satellite soil moisture observations: Part 1. A review of retrieval algorithms. *Adv. Water Resour.* **2017**, *109*, 106–120. [CrossRef]

56. Lv, S.; Zeng, Y.; Wen, J.; Zhao, H.; Su, Z. Estimation of penetration depth from soil effective temperature in microwave radiometry. *Remote Sens.* **2018**, *10*, 519. [CrossRef]

57. Usowicz, B.; Lukowski, M.; Lipiec, J. The SMOS-Derived Soil Water Extent and equivalent layer thickness facilitate determination of soil water resources. *Sci. Rep.* **2020**, *10*, 1–12. [CrossRef]

58. Martens, B.; Miralles, D.; Lievens, H.; Fernández-Prieto, D.; Verhoest, N.E. Improving terrestrial evaporation estimates over continental Australia through assimilation of SMOS soil moisture. *Int. J. Appl. Earth Obs. Geoinf.* **2016**, *48*, 146–162. [CrossRef]

59. Yee, M.S.; Walker, J.P.; Rüdiger, C.; Parinussa, R.M.; Koike, T.; Kerr, Y.H. A comparison of SMOS and AMSR2 soil moisture using representative sites of the OzNet monitoring network. *Remote Sens. Environ.* **2017**, *195*, 297–312. [CrossRef]

60. Jing, W.; Song, J.; Zhao, X. Evaluation of Multiple Satellite-Based Soil Moisture Products over Continental US Based on In Situ Measurements. *Water Resour. Manag.* **2018**, *32*, 3233–3246. [CrossRef]
61. Albergel, C.; Rüdiger, C.; Carrer, D.; Calvet, J.-C.; Fritz, N.; Naeimi, V.; Bartalis, Z.; Hasenauer, S. An evaluation of ASCAT surface soil moisture products with in-situ observations in Southwestern France. *Hydrol. Earth Syst. Sci.* **2009**, *13*, 115–124. [CrossRef]

62. Lei, F.; Crow, W.T.; Shen, H.; Parinussa, R.M.; Holmes, T.R. The impact of local acquisition time on the accuracy of microwave surface soil moisture retrievals over the contiguous United States. *Remote Sens.* **2015**, *7*, 13448–13465. [CrossRef]

63. O'Neill, P.; Bindlish, R.; Steven, C.; Julian, C.; Eni, N.; Tom, J. *Algorithm Theoretical Basis Document Level 2 & 3 Soil Moisture (Passive) Data Products*; Jet Propulsion Laboratory: Pasadena, CA, USA, 2020.

64. Yueh, S.H. Estimates of Faraday rotation with passive microwave polarimetry for microwave remote sensing of Earth surfaces. *IEEE Trans. Geosci. Remote Sens.* **2000**, *38*, 2434–2438. [CrossRef]

65. Beck, H.E.; Pan, M.; Miralles, D.G.; Reichle, R.H.; Dorigo, W.A.; Hahn, S.; Sheffield, J.; Karthikeyan, L.; Balsamo, G.; Parinussa, R.M.; et al. Evaluation of 18 satellite-and model-based soil moisture products using in situ measurements from 826 sensors. *Hydrol. Earth Syst. Sci.* **2021**, *25*, 17–40. [CrossRef]

66. Bosch, P.S.; Goodrich, D.; Prueger, J.; Palecki, M.; Small, E.; Zreda, M.; Calvet, J.; Crow, W. SMAP L2 Radiometer Half-Orbit 36 km EASE-Grid Soil Moisture, Version 6. Available online: https://nsidc.org/sites/nsidc.org/files/SPL2SMP-V006-UserGuide_0.pdf (accessed on 2 February 2022).

67. Zwieback, S.; Westermann, S.; Langer, M.; Boike, J.; Marsh, P.; Berg, A. Improving permafrost modeling by assimilating remotely sensed soil moisture. *Water Resour. Res.* **2019**, *55*, 1814–1832. [CrossRef]

68. Yin, J.; Zhan, X.; Liu, J.; Schull, M. An Intercomparison of Noah Model Skills With Benefits of Assimilating SMOS Blended and Individual Soil Moisture Retrievals. *Water Resour. Res.* **2019**, *55*, 2572–2592. [CrossRef]

69. Sabaghy, S.; Walker, J.P.; Renzullo, L.J.; Akbar, R.; Chan, S.; Chaubell, J.; Das, N.; Dunbar, R.S.; Entekhabi, D.; Gevaert, A.; et al. Comprehensive analysis of alternative downscaled soil moisture products. *Remote Sens. Environ.* **2020**, *239*, 111586. [CrossRef]

70. Lievens, H.; De Lannoy, G.; Al Bitar, A.; Drusch, M.; Dumedah, G.; Franssen, H.-J.H.; Kerr, Y.; Tomer, S.K.; Martens, B.; Merlin, O.; et al. Assimilation of SMOS soil moisture and brightness temperature products into a land surface model. *Remote Sens. Environ.* **2016**, *180*, 292–304. [CrossRef]

71. Fujii, H. Development of a Microwave Radiative Transfer Model for Vegetated Land Surface Based on Comprehensive In Situ Observations. Ph.D. Thesis, The University of Tokyo, Tokyo, Japan, 2005.

72. Shellito, P.J.; Small, E.E.; Livneh, B. Controls on surface soil drying rates observed by SMAP and simulated by the Noah land surface model. *Hydrol. Earth Syst. Sci.* **2018**, *22*, 1649–1663. [CrossRef]

73. Ma, H.; Zeng, J.; Chen, N.; Zhang, X.; Cosh, M.H.; Wang, W. Satellite surface soil moisture from SMAP, SMOS, AMSR2 and ESA CCI: A comprehensive assessment using global ground-based observations. *Remote Sens. Environ.* **2019**, *231*, 111215. [CrossRef]

74. Al-Yaari, A.; Wigneron, J.-P.; Dorigo, W.; Colliander, A.; Pellarin, T.; Hahn, S.; Mialon, A.; Richaume, P.; Fernandez-Moran, R.; Fan, L.; et al. Assessment and inter-comparison of recently developed/reprocessed microwave satellite soil moisture products using ISMN ground-based measurements. *Remote Sens. Environ.* **2019**, *224*, 289–303. [CrossRef]

75. Li, X.; Wigneron, J.-P.; Frappart, F.; Fan, L.; Wang, M.; Liu, X.; Al-Yaari, A.; Moisy, C. Development and validation of the SMOS-IC version 2 (V2) soil moisture product. In Proceedings of the IGARSS 2020–2020 IEEE International Geoscience and Remote Sensing Symposium, Waikoloa, HI, USA, 26 September–2 October 2020; pp. 4434–4437.

76. Gruber, A.; De Lannoy, G.; Albergel, C.; Al-Yaari, A.; Brocca, L.; Calvet, J.-C.; Colliander, A.; Cosh, M.; Crow, W.; Dorigo, W. Validation practices for satellite soil moisture retrievals: What are (the) errors? *Remote Sens. Environ.* **2020**, *244*, 111806. [CrossRef]

77. Zhao, L.; Yang, Z.-L.; Hoar, T.J. Global soil moisture estimation by assimilating AMSR-E brightness temperatures in a coupled CLM4-RTM-DART system. *J. Hydrometeorol.* **2016**, *17*, 2431–2454. [CrossRef]

78. Lu, H.; Yang, K.; Koike, T.; Zhao, L.; Qin, J. An improvement of the radiative transfer model component of a land data assimilation system and its validation on different land characteristics. *Remote Sens.* **2015**, *7*, 6358–6379. [CrossRef]

79. Yang, K.; Watanabe, T.; Koike, T.; Li, X.; Fujii, H.; Tamagawa, K.; Ma, Y.; Ishikawa, H. Auto-calibration system developed to assimilate AMSR-E data into a land surface model for estimating soil moisture and the surface energy budget. *J. Meteorol. Soc. Jpn. Ser. II* **2007**, *85*, 229–242. [CrossRef]

80. Yang, K.; Koike, T.; Kaihotsu, I.; Qin, J. Validation of a dual-pass microwave land data assimilation system for estimating surface soil moisture in semiarid regions. *J. Hydrometeorol.* **2009**, *10*, 780–793. [CrossRef]

81. Duan, Q.; Gupta, V.K.; Sorooshian, S. Shuffled complex evolution approach for effective and efficient global minimization. *J. Optim. Theory Appl.* **1993**, *76*, 501–521. [CrossRef]

82. Tian, X.; Xie, Z.; Dai, A. An ensemble conditional nonlinear optimal perturbation approach: Formulation and applications to parameter calibration. *Water Resour. Res.* **2010**, *46*, W09540. [CrossRef]

83. De Lannoy, G.J.; Reichle, R.H.; Pauwels, V.R. Global calibration of the GEOS-5 L-band microwave radiative transfer model over nonfrozen land using SMOS observations. *J. Hydrometeorol.* **2013**, *14*, 765–785. [CrossRef]

84. Qin, J.; Liang, S.; Yang, K.; Kaihotsu, I.; Liu, R.; Koike, T. Simultaneous estimation of both soil moisture and model parameters using particle filtering method through the assimilation of microwave signal. *J. Geophys. Res. Atmos.* **2009**, *114*, D15103. [CrossRef]

85. Tao, L.; Wang, G.; Chen, W.; Chen, X.; Li, J.; Cai, Q. Soil moisture retrieval from SAR and optical data using a combined model. *IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens.* **2019**, *12*, 637–647. [CrossRef]

86. Parrens, M.; Wigneron, J.-P.; Richaume, P.; Mialon, A.; Al Bitar, A.; Fernandez-Moran, R.; Al-Yaari, A.; Kerr, Y.H. Global-scale surface roughness effects at L-band as estimated from SMOS observations. *Remote Sens. Environ.* **2016**, *181*, 122–136. [CrossRef]

87. Weng, F.; Yan, B.; Grody, N.C. A microwave land emissivity model. *J. Geophys. Res. Atmos.* **2001**, *106*, 20115–20123. [CrossRef]
88. Drusch, M.; Holmes, T.; de Rosnay, P.; Balsamo, G. Comparing ERA-40-based L-band brightness temperatures with Skylab observations: A calibration/validation study using the Community Microwave Emission Model. *J. Hydrometeorol.* **2009**, *10*, 213–226. [CrossRef]

89. Mo, T.; Choudhury, B.; Schmugge, T.; Wang, J.; Jackson, T. A model for microwave emission from vegetation-covered fields. *J. Geophys. Res. Ocean.* **1982**, *87*, 11229–11237. [CrossRef]

90. Zhang, S.; Shi, J.; Dou, Y. A soil moisture assimilation scheme based on the microwave land emissivity model and the community land model. *Int. J. Remote Sens.* **2012**, *33*, 2770–2797. [CrossRef]

91. Wigneron, J.-P.; Pardé, M.; Waldteufel, P.; Chanzy, A.; Kerr, Y.; Schmidl, S.; Skou, N. Characterizing the dependence of vegetation model parameters on crop structure, incidence angle, and polarization at L-band. *IEEE Trans. Geosci. Remote Sens.* **2004**, *42*, 416–425. [CrossRef]

92. Wigneron, J.-P.; Laguerre, L.; Kerr, Y.H. A simple parameterization of the L-band microwave emission from rough agricultural soils. *IEEE Trans. Geosci. Remote Sens.* **2001**, *39*, 1697–1707. [CrossRef]

93. Wigneron, J.-P.; Kerr, Y.; Waldteufel, P.; Saleh, K.; Escorihuela, M.-J.; Richaume, P.; Ferrazzoli, P.; De Rosnay, P.; Gurney, R.; Calvet, J.-C.; et al. L-band Microwave Emission of the Biosphere (L-MEB) Model: Description and calibration against experimental data sets over crop fields. *Remote Sens. Environ.* **2007**, *107*, 639–655. [CrossRef]

94. Jia, B.; Xie, Z.; Tian, X.; Shi, C. A soil moisture assimilation scheme based on the ensemble Kalman filter using microwave brightness temperature. *Sci. China Ser. D Earth Sci.* **2009**, *52*, 1835. [CrossRef]

95. Sawada, Y.; Koike, T. Simultaneous estimation of both hydrological and ecological parameters in an ecohydrological model by assimilating microwave signal. *J. Geophys. Res. Atmos.* **2014**, *119*, 8839–8857. [CrossRef]

96. Montzka, C.; Grant, J.P.; Moradkhani, H.; Franssen, H.-J.H.; Weihermüller, L.; Drusch, M.; Vereecken, H. Estimation of radiative transfer parameters from L-band passive microwave brightness temperatures using advanced data assimilation. *Vadose Zone J.* **2013**, *12*, vzj2012.0040. [CrossRef]

97. Yang, K.; Zhu, L.; Chen, Y.; Zhao, L.; Qin, J.; Lu, H.; Tang, W.; Han, M.; Ding, B.; Fang, N. Land surface model calibration through microwave data assimilation for improving soil moisture simulations. *J. Hydrol.* **2016**, *533*, 266–276. [CrossRef]

98. Kumar, S.V.; Reichle, R.H.; Harrison, K.W.; Peters-Lidard, C.D.; Yatheendradas, S.; Santanello, J.A. A comparison of methods for a priori bias correction in soil moisture data assimilation. *Water Resour. Res.* **2012**, *48*, W03515. [CrossRef]

99. Afshar, M.; Yilmaz, M. The added utility of nonlinear methods compared to linear methods in rescaling soil moisture products. *Remote Sens. Environ.* **2017**, *196*, 224–237. [CrossRef]

100. Brocca, L.; Melone, F.; Moramarco, T.; Wagner, W.; Albergel, C. Scaling and filtering approaches for the use of satellite soil moisture observations. In *Remote Sensing of Energy Fluxes and Soil Moisture Content*; Routledge: Abingron-on-Thames, UK, 2013; pp. 411–426.

101. Yilmaz, M.T.; Crow, W.T. The optimality of potential rescaling approaches in land data assimilation. *J. Hydrometeorol.* **2013**, *14*, 650–660. [CrossRef]

102. Fertig, E.; Baek, S.-J.; Hunt, B.; Ott, E.; Szunyogh, I.; Aravéquia, J.; Kalnay, E.; Li, H.; Liu, J. Observation bias correction with an ensemble Kalman filter. *Tellus A Dyn. Meteorol. Oceanogr.* **2009**, *61*, 210–226. [CrossRef]

103. De Lannoy, G.J.; Houser, P.R.; Pauwels, V.R.; Verhoest, N.E. State and bias estimation for soil moisture profiles by an ensemble Kalman filter: Effect of assimilation depth and frequency. *Water Resour. Res.* **2007**, *43*, W06401. [CrossRef]

104. De Lannoy, G.J.; Reichle, R.H.; Houser, P.R.; Pauwels, V.R.; Verhoest, N.E. Correcting for forecast bias in soil moisture assimilation with the ensemble Kalman filter. *Water Resour. Res.* **2007**, *43*, W09410. [CrossRef]

105. Bosilovich, M.G.; Radakovich, J.D.; da Silva, A.; Todling, R.; Verter, F. Skin temperature analysis and bias correction in a coupled land-atmosphere data assimilation system. *J. Meteorol. Soc. Jpn. Ser. II* **2007**, *85*, 205–228. [CrossRef]

106. Reichle, R.H.; Kumar, S.V.; Mahanama, S.P.; Koster, R.D.; Liu, Q. Assimilation of satellite-derived skin temperature observations into land surface models. *J. Hydrometeorol.* **2010**, *11*, 1103–1122. [CrossRef]

107. Ridler, M.E.; Madsen, H.; Stisen, S.; Bircher, S.; Fensholt, R. Assimilation of SMOS-derived soil moisture in a fully integrated hydrological and soil-vegetation-atmosphere transfer model in Western Denmark. *Water Resour. Res.* **2014**, *50*, 8962–8981. [CrossRef]

108. Lievens, H.; Tomer, S.K.; Al Bitar, A.; De Lannoy, G.J.; Drusch, M.; Dumedah, G.; Franssen, H.-J.H.; Kerr, Y.H.; Martens, B.; Pan, M.; et al. SMOS soil moisture assimilation for improved hydrologic simulation in the Murray Darling Basin, Australia. *Remote Sens. Environ.* **2015**, *168*, 146–162. [CrossRef]

109. Reichle, R.H.; Koster, R.D. Bias reduction in short records of satellite soil moisture. *Geophys. Res. Lett.* **2004**, *31*, L19501. [CrossRef]

110. Mahfouf, J.F. Assimilation of satellite-derived soil moisture from ASCAT in a limited-area NWP model. *Q. J. R. Meteorol. Soc.* *A J. Atmos. Sci. Appl. Meteorol. Phys. Oceanogr.* **2010**, *136*, 784–798. [CrossRef]

111. Kolassa, J.; Reichle, R.; Liu, Q.; Cosh, M.; Bosch, D.; Caldwell, T.; Colliander, A.; Holifield Collins, C.; Jackson, T.; Livingston, S.; et al. Data assimilation to extract soil moisture information from SMAP observations. *Remote Sens.* **2017**, *9*, 1179. [CrossRef]

112. Schneider, S.; Wang, Y.; Wagner, W.; Mahfouf, J.-F. Impact of ASCAT soil moisture assimilation on regional precipitation forecasts: A case study for Austria. *Mon. Weather Rev.* **2014**, *142*, 1525–1541. [CrossRef]

113. Blankenship, C.B.; Case, J.L.; Crosson, W.L.; Zavodsky, B.T. Correction of forcing-related spatial artifacts in a land surface model by satellite soil moisture data assimilation. *IEEE Geosci. Remote Sens. Lett.* **2018**, *15*, 498–502. [CrossRef]
114. Blankenship, C.B.; Case, J.L.; Zavodsky, B.T.; Crosson, W.L. Assimilation of SMOS retrievals in the Land Information System. *IEEE Trans. Geosci. Remote Sens.* **2016**, *54*, 6320–6332. [CrossRef]

115. Kumar, S.; Peters-Lidard, C.; Santanello, J.; Reichle, R.; Draper, C.; Koster, R.; Nearing, G.; Jasinski, M. Evaluating the utility of satellite soil moisture retrievals over irrigated areas and the ability of land data assimilation methods to correct for unmodeled processes. *Hydrol. Earth Syst. Sci.* **2015**, *19*, 4463–4478. [CrossRef]

116. Santanello, J.A., Jr.; Kumar, S.V.; Peters-Lidard, C.D.; Lawston, P.M. Impact of soil moisture assimilation on land surface model spinup and coupled land-atmosphere prediction. *J. Hydrometeorol.* **2016**, *17*, 517–540. [CrossRef]

117. Barbu, A.; Calvet, J.-C.; Mahfouf, J.-F.; Lafont, S. Integrating ASCAT surface soil moisture and GEOV1 leaf area index into the SURFEX modelling platform: A land data assimilation application over France. *Hydrol. Earth Syst. Sci.* **2014**, *18*, 173–192. [CrossRef]

118. Draper, C.; Mahfouf, J.F.; Walker, J.P. An EKF assimilation of AMSR-E soil moisture into the ISBA land surface scheme. *J. Geophys. Res. Atmos.* **2009**, *114*, D20104. [CrossRef]

119. Carrera, M.L.; Bilodeau, B.; Bélair, S.; Abrahamowicz, M.; Russell, A.; Wang, X. Assimilation of passive L-band microwave brightness temperatures in the Canadian Land Data Assimilation System: Impacts on short-range warm season Numerical Weather Prediction. *J. Hydrometeorol.* **2019**, *20*, 1053–1079. [CrossRef]

120. Yin, J.; Zhan, X. Impact of bias-correction methods on effectiveness of assimilating SMAP soil moisture data into NCEP global forecast system using the ensemble Kalman filter. *IEEE Geosci. Remote Sens. Lett.* **2018**, *15*, 659–663. [CrossRef]

121. Lin, L.F.; Ebtehaj, A.M.; Wang, J.; Bras, R.L. Soil moisture background error covariance and data assimilation in a coupled land-atmosphere model. *Water Resour. Res.* **2017**, *53*, 1309–1335. [CrossRef]

122. Kalman, R.E. A new approach to linear filtering and prediction problems. *Trans. ASME J. Basic Eng.* D **1960**, *82*, 35–45. [CrossRef]

123. Evensen, G. Using the extended Kalman filter with a multilayer quasi-geostrophic ocean model. *J. Geophys. Res. Ocean.* **1992**, *97*, 17905–17924. [CrossRef]

124. Anderson, J.L. An ensemble adjustment Kalman filter for data assimilation. *Mon. Weather Rev.* **2001**, *129*, 2884–2903. [CrossRef]

125. Evensen, G.; Van Leeuwen, P.J. An ensemble Kalman smoother for nonlinear dynamics. *Mon. Weather Rev.* **2000**, *128*, 1852–1867. [CrossRef]

126. Liu, C.; Xiao, Q.; Wang, B. An ensemble-based four-dimensional variational data assimilation scheme. Part I: Technical formulation and preliminary test. *Mon. Weather Rev.* **2008**, *136*, 3363–3373. [CrossRef]

127. Evensen, G. The ensemble Kalman filter: Theoretical formulation and practical implementation. *Ocean Dyn.* **2003**, *53*, 343–367. [CrossRef]

128. Oke, P.R.; Allen, J.S.; Miller, R.N.; Egbert, G.D.; Kosro, P.M. Assimilation of surface velocity data into a primitive equation coastal ocean model. *J. Geophys. Res. Ocean.* **2002**, *107*, 5–1–5–25. [CrossRef]

129. Sakov, P.; Oke, P.R. Implications of the form of the ensemble transformation in the ensemble square root filters. *Mon. Weather Rev.* **2008**, *136*, 1042–1053. [CrossRef]

130. Bishop, C.H.; Etherton, B.J.; Majumdar, S.J. Adaptive sampling with the ensemble transform Kalman filter. Part I: Theoretical aspects. *Mon. Weather Rev.* **2001**, *129*, 420–436. [CrossRef]

131. Dumedah, G. Formulation of the evolutionary-based data assimilation, and its implementation in hydrological forecasting. *Water Resour. Manag.* **2012**, *26*, 3853–3870. [CrossRef]

132. Arulampalam, M.S.; Maskell, S.; Gordon, N.; Clapp, T. A tutorial on particle filters for online nonlinear/non-Gaussian Bayesian tracking. *IEEE Trans. Signal Process.* **2002**, *50*, 174–188. [CrossRef]

133. Hoke, J.E.; Anthes, R.A. The initialization of numerical models by a dynamic-initialization technique. *Mon. Weather Rev.* **1976**, *104*, 1551–1556. [CrossRef]

134. Houser, P.R.; Shuttleworth, W.J.; Famiglietti, J.S.; Gupta, H.V.; Syed, K.H.; Goodrich, D.C. Integration of soil moisture remote sensing and hydrologic modeling using data assimilation. *Water Resour. Res.* **1998**, *34*, 3405–3420. [CrossRef]

135. Giard, D.; Bazile, E. Implementation of a new assimilation scheme for soil and surface variables in a global NWP model. *Mon. Weather Rev.* **2000**, *128*, 997–1015. [CrossRef]

136. Sasaki, Y. An objective analysis based on the variational method. *J. Meteorol. Soc. Jpn. Ser. II* **1958**, *36*, 77–88. [CrossRef]

137. Dumedah, G.; Walker, J.P. Evaluation of model parameter convergence when using data assimilation for soil moisture estimation. *J. Hydrometeorol.* **2014**, *15*, 359–375. [CrossRef]

138. Reichle, R.H.; Walker, J.P.; Koster, R.D.; Houser, P.R. Extended versus ensemble Kalman filtering for land data assimilation. *J. Hydrometeorol.* **2002**, *3*, 728–740. [CrossRef]

139. Sabater, J.M.; Jarlan, L.; Calvet, J.-C.; Bouyssel, F.; De Rosnay, P. From near-surface to root-zone soil moisture using different assimilation techniques. *J. Hydrometeorol.* **2007**, *8*, 194–206. [CrossRef]

140. Bi, H.; Ma, J.; Wang, F. An improved particle filter algorithm based on ensemble Kalman filter and Markov chain Monte Carlo method. *IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens.* **2014**, *8*, 447–459. [CrossRef]

141. Yan, H.; DeChant, C.M.; Moradkhani, H. Improving soil moisture profile prediction with the particle filter-Markov chain Monte Carlo method. *IEEE Trans. Geosci. Remote Sens.* **2015**, *53*, 6134–6147. [CrossRef]

142. Yin, J.; Zhan, X.; Zheng, Y.; Hain, C.R.; Liu, J.; Fang, L. Optimal ensemble size of ensemble Kalman filter in sequential soil moisture data assimilation. *Geophys. Res. Lett.* **2015**, *42*, 6710–6715. [CrossRef]
143. Reichle, R.H.; McLaughlin, D.B.; Entekhabi, D. Hydrologic data assimilation with the ensemble Kalman filter. *Mon. Weather Rev.* **2002**, *130*, 103–114. [CrossRef]

144. Kumar, S.V.; Reichle, R.H.; Peters-Lidard, C.D.; Koster, R.D.; Zhan, X.; Crow, W.T.; Eylander, J.B.; Houser, P.R. A land surface data assimilation framework using the land information system: Description and applications. *Adv. Water Resour.* **2008**, *31*, 1419–1432. [CrossRef]

145. Lievens, H.; Reichle, R.H.; Liu, Q.; De Lannoy, G.; Dunbar, R.S.; Kim, S.; Das, N.N.; Cosh, M.; Walker, J.P.; Wagner, W. Joint Sentinel-1 and SMAP data assimilation to improve soil moisture estimates. *Geophys. Res. Lett.* **2017**, *44*, 6145–6153. [CrossRef]

146. Shi, X.; Wen, J.; Wang, L.; Zhang, T.; Tian, H.; Wang, X.; Liu, R.; Zhang, J. Regional soil moisture retrievals and simulations from assimilation of satellite microwave brightness temperature observations. *Environ. Earth Sci.* **2010**, *61*, 1289–1299. [CrossRef]

147. Li, B.; Toll, D.; Zhan, X.; Cosgrove, B. Improving estimated soil moisture fields through assimilation of AMSR-E soil moisture retrievals with an ensemble Kalman filter and a mass conservation constraint. *Hydrol. Earth Syst. Sci.* **2012**, *16*, 105–119. [CrossRef]

148. Sahoo, A.K.; De Lannoy, G.J.; Reichle, R.H.; Houser, P.R. Assimilation and downscaling of satellite observed soil moisture over the Little River Experimental Watershed in Georgia, USA. *Adv. Water Resour.* **2013**, *52*, 19–33. [CrossRef]

149. Dumedah, G.; Walker, J.P.; Merlin, O. Root-zone soil moisture estimation from assimilation of downscaled Soil Moisture and Ocean Salinity data. *Adv. Water Resour.* **2015**, *84*, 14–22. [CrossRef]

150. Huang, C.; Li, X.; Lu, L.; Gu, J. Experiments of one-dimensional soil moisture assimilation system based on ensemble Kalman filter. *Remote Sens. Environ.* **2008**, *112*, 888–900. [CrossRef]

151. Kumar, S.V.; Reichle, R.H.; Koster, R.D.; Crow, W.T.; Peters-Lidard, C.D. Role of subsurface physics in the assimilation of surface soil moisture observations. *J. Hydrometeorol.* **2009**, *10*, 1534–1547. [CrossRef]

152. Lin, L.-F.; Ebtehaj, A.M.; Flores, A.N.; Bastola, S.; Bras, R.L. Combined assimilation of satellite precipitation and soil moisture: A case study using trmm and smos data. *Mon. Weather Rev.* **2017**, *145*, 4997–5014. [CrossRef]

153. Gevaert, A.I.; Renzullo, L.J.; Van Dijk, A.I.J.M.; van der Woerd, H.J.; Weerts, A.H.; De Jeu, R.A. Joint assimilation of soil moisture retrieved from multiple passive microwave frequencies increases robustness of soil moisture state estimation. *Hydrol. Earth Syst. Sci.* **2018**, *22*, 4605–4619. [CrossRef]

154. Draper, C.; Reichle, R.; De Lannoy, G.; Liu, Q. Assimilation of passive and active microwave soil moisture retrievals. *Geophys. Res. Lett.* **2012**, *39*, L04401. [CrossRef]

155. Tian, S.; Renzullo, L.J.; van Dijk, A.I.; Tregoning, P.; Walker, J.P. Global joint assimilation of GRACE and SMOS for improved estimation of root-zone soil moisture and vegetation response. *Hydrol. Earth Syst. Sci.* **2019**, *23*, 1067–1081. [CrossRef]

156. Yin, J.; Zhan, X.; Zheng, Y.; Liu, J.; Fang, L.; Hain, C.R. Enhancing model skill by assimilating SMOPS blended soil moisture product into Noah land surface model. *J. Hydrometeorol.* **2015**, *16*, 917–931. [CrossRef]

157. Kolassa, J.; Reichle, R.; Draper, C.S. Merging active and passive microwave observations in soil moisture data assimilation. *Remote Sens. Environ.* **2017**, *191*, 117–130. [CrossRef]

158. Gavahi, K.; Abbaszadeh, P.; Moradkhani, H.; Zhan, X.; Hain, C. Multivariate assimilation of remotely sensed soil moisture and evapotranspiration for drought monitoring. *J. Hydrometeorol.* **2020**, *21*, 2293–2308. [CrossRef]

159. Moradkhani, H.; Nearing, G.S.; Abbaszadeh, P.; Pathiraja, S. Fundamentals of data assimilation and theoretical advances. In *Handbook of Hydrometeorological Ensemble Forecasting*; Springer: Berlin/Heidelberg, Germany, 2019; pp. 675–699.

160. Crow, W.T.; Wood, E.F. The assimilation of remotely sensed soil brightness temperature imagery into a land surface model using ensemble Kalman filtering: A case study based on ESTAR measurements during SGP97. *Adv. Water Resour.* **2003**, *26*, 137–149. [CrossRef]

161. Dechant, C.; Moradkhani, H. Radiance data assimilation for operational snow and streamflow forecasting. *Adv. Water Resour.* **2011**, *34*, 351–364. [CrossRef]

162. Abbaszadeh, P.; Moradkhani, H.; Yan, H. Enhancing hydrologic data assimilation by evolutionary particle filter and Markov chain Monte Carlo. *Adv. Water Resour.* **2018**, *111*, 192–204. [CrossRef]

163. Moradkhani, H. Hydrologic remote sensing and land surface data assimilation. *Sensors* **2008**, *8*, 2986–3004. [CrossRef] [PubMed]

164. Fairbairn, D.; Barbu, A.; Mahfouf, J.-F.; Calvet, J.-C.; Gelati, E. Comparing the ensemble and extended Kalman filters for in situ soil moisture assimilation with contrasting conditions. *Hydrol. Earth Syst. Sci.* **2015**, *19*, 4811–4830. [CrossRef]

165. Draper, C.; Reichle, R.H. Assimilation of satellite soil moisture for improved atmospheric reanalyses. *Mon. Weather Rev.* **2019**, *147*, 2163–2188. [CrossRef]

166. Entin, J.K.; Robock, A.; Vinnikov, K.Y.; Hollinger, S.E.; Liu, S.; Namkhai, A. Temporal and spatial scales of observed soil moisture variations in the extratropics. *J. Geophys. Res. Atmos.* **2000**, *105*, 11865–11877. [CrossRef]

167. Vinnikov, K.Y.; Robock, A.; Speranskaya, N.A.; Schlosser, C.A. Scales of temporal and spatial variability of midlatitude soil moisture. *J. Geophys. Res. Atmos.* **1996**, *101*, 7163–7174. [CrossRef]

168. Draper, C.; Reichle, R. The impact of near-surface soil moisture assimilation at subseasonal, seasonal, and inter-annual timescales. *Hydrol. Earth Syst. Sci.* **2015**, *19*, 4831–4844. [CrossRef]

169. De Lannoy, G.; Reichle, R. Assimilation of SMOS brightness temperatures or soil moisture retrievals into a land surface model. *Hydrol. Earth Syst. Sci.* **2016**, *20*, 4895–4911. [CrossRef]

170. Peng, J.; Albergel, C.; Balenzano, A.; Brocca, L.; Cartus, O.; Cosh, M.H.; Crow, W.T.; Dabrowska-Zielinska, K.; Dadson, S.; Davidson, M.W.; et al. A roadmap for high-resolution satellite soil moisture applications–confronting product characteristics with user requirements. *Remote Sens. Environ.* **2020**, *252*, 112162. [CrossRef]
171. Abbaszadeh, P.; Moradkhani, H.; Zhan, X. Downscaling SMAP radiometer soil moisture over the CONUS using an ensemble learning method. *Water Resour. Res.* **2019**, *55*, 324–344. [CrossRef]