Soil Moisture And Ocean Salinity Research Articles

Soil Moisture and Sea Surface Salinity Derived from Satellite-Borne Sensors

AbstractThe monitoring of soil moisture and sea surface salinity over the Earth has been profoundly enhanced during the last thirteen years due to a new generation of satellite sensors. L-band radiometry is currently the only technology providing direct measurements of soil moisture, insensitive to surface roughness and distribution of elements in the soil, and the only technology the only technology for measuring that allows us to measure sea surface salinity from space. The Soil Moisture and Ocean Salinity (SMOS) and Soil Moisture Active Passive (SMAP) satellite missions resolve global and local variability with a spatial resolution of approximately 43 km, a swath width close to 1000 km, and a sampling time, for each mission, of at least twice every 3 days. These resolutions and samplings can be increased by either merging data from the two sensors, and with complementary information gathered from other passive or active sensors, or with in situ information at higher spatial resolution. Numerous scientific studies based on the use of this new type of measurement have led to a better understanding and constraint of the processes governing the variability of the water cycle, ocean circulation and the Earth's climate. The continuity of measurements, and the increased spatial and radiometric resolution is critical for fulfilling scientific needs. Future L-band radiometry missions currently being planned in Europe (the Copernicus Imaging Microwave Radiometer), and in China (the Ocean Salinity mission) should provide better constraints on auxiliary parameters by combining multiple frequencies, but they will not have improved spatial resolution beyond SMOS and SMAP. The temporal continuity with SMOS and SMAP will likely not be ensured. In parallel, new concepts are being developed to increase spatial resolution of both land and ocean parameters.

Estimation of High-Resolution Soil Moisture in Canadian Croplands Using Deep Neural Network with Sentinel-1 and Sentinel-2 Images

Soil moisture (SM) is a crucial hydrologic factor that affects the global cycle of energy, carbon, and water, as well as plant growth and crop yield; therefore, an accurate estimate of SM is important for both the global environment and agriculture. Satellite-based SM data have been provided by the National Aeronautics and Space Administration (NASA)’s Soil Moisture Active Passive (SMAP) and the European Space Agency (ESA)’s Soil Moisture and Ocean Salinity (SMOS) satellite missions, but these data are based on passive microwave sensors, which have limited spatial resolution. Thus, detailed observations and analyses of the local distribution of SM are limited. The recent emergence of deep learning techniques, such as rectified linear unit (ReLU) and dropout, has produced effective solutions to complex problems. Deep neural networks (DNNs) have been used to accurately estimate hydrologic factors, such as SM and evapotranspiration, but studies of SM estimates derived from the joint use of DNN and high-resolution satellite data, such as Sentinel-1 and Sentinel-2, are lacking. In this study, we aim to estimate high-resolution SM at 30 m resolution, which is important for local-scale SM monitoring in croplands. We used a variety of input data, such as radar factors, optical factors, and vegetation indices, which can be extracted from Sentinel-1 and -2, terrain information (e.g., elevation), and crop information (e.g., cover type and month), and developed an integrated SM model across various crop surfaces by using these input data and DNN (which can learn the complexity and nonlinearity of the various data). The study was performed in the agricultural areas of Manitoba and Saskatchewan, Canada, and the in situ SM data for these areas were obtained from the Agriculture and Agri-Food Canada (AAFC) Real-time In Situ Soil Monitoring for Agriculture (RISMA) network. We conducted various experiments with several hyperparameters that affected the performance of the DNN-based model and ultimately obtained a high-performing SM model. The optimal SM model had a root-mean-square error (RMSE) of 0.0416 m3/m3 and a correlation coefficient (CC) of 0.9226. This model’s estimates showed better agreement with in situ SM than the SMAP 9 km SM. The accuracy of the model was high when the daily precipitation was zero or very low and also during the vegetation growth stage. However, its accuracy decreased when precipitation or the vitality of the vegetation were high. This suggests that precipitation affects surface erosion and water layer formation, and vegetation adds complexity to the SM estimate. Nevertheless, the distribution of SM estimated by our model generally reflected the local soil characteristics. This work will aid in drought and flood prevention and mitigation, and serve as a tool for assessing the potential growth of crops according to SM conditions.

Integration of microwave satellite soil moisture products in the contextual surface temperature-vegetation index models for spatially continuous evapotranspiration estimation

The contextual surface temperature-vegetation index (TVX) models have been widely used for the retrieval of soil moisture (SM) and evapotranspiration (ET). One of the key premises for their application is to determine quantitatively the theoretical boundaries of this contextual TVX space. Usually these theoretical boundaries are determined based on land surface energy balance principle. Although sound in physical mechanism, the complex parameterization involved has hampered its remote sensing (RS) applications. Our recent studies show that ground-based SM observations can be used for a quick and continuous retrieval of these theoretical boundaries. The main disadvantage lies in its reliance on in-situ SM observations. Fortunately the development of microwave RS technology has made it possible to monitor SM at both regional and global scale. Under this background, a practical framework based on the contextual TVX models was proposed in this paper for continuous ET estimation through the combination of microwave-based SM products and MODIS (Moderate Resolution Imaging Spectroradiometer) products. Daily actual ET over the Southern Great Plains of the United States of America was calculated as the product of net radiation and evaporative fraction (EF). The focus of this framework was to investigate how optical and microwave RS can be combined to achieve a quick and continuous retrieval of EF independently of ground observations. Specifically, EF was retrieved from the contextual TVX models, the theoretical boundaries of which were calibrated at the annual scale based on an optimization scheme. Assuming the seasonal variation of the theoretical dry edge follows a cosine function of solar zenith, the essence of this optimization scheme is to search for the optimal amplitude parameter that maximizes the correlation between microwave-based SM observations and soil moisture index retrieved from the contextual TVX models. Two microwave-based SM products including SMAP (Soil Moisture Active Passive) and SMOS (Soil Moisture and Ocean Salinity) were adopted to demonstrate this optimization scheme. The correlation coefficient, root mean square error and bias of EF estimates calibrated by SMAP were 0.700, 0.146 and 0.047, and those achieved by SMOS were 0.663, 0.174 and 0.168, respectively. The optimization scheme has not only made it possible to conduct a continuous monitoring of ET based entirely on RS observations but also achieved the downscaling of microwave satellite SM products. Thus the combination of optical and microwave RS holds great potential for coupled and continuous estimation of SM and ET over large scale.

Soil moisture at 30 m from multiple satellite datasets fused by random forest

As a key variable used to characterize the climate process between the land surface and atmosphere, the surface soil moisture (SM) plays an irreplaceable role in the fields of hydrology, meteorology and agriculture. However, the relatively coarse spatial resolution of SM products currently limit the application in water resource management at the field scale. In this study, we proposed a high-spatial-resolution SM retrieval framework based on the random forest algorithm (RF-SM) to integrate in situ SM dataset from in the International Soil Moisture Network (ISMN), Landsat 8 optical and thermal observations, soil properties from SoilGrids V2.0, meteorological variables from the fifth generation of the European ReAnalysis (ERA5) dataset and four coarse-scale SM products including the Soil Moisture Active/Passive (SMAP), the Soil Moisture and Ocean Salinity (SMOS), European Space Agency Climate Change Initiative (ESA CCI) and U.S. National Climate Assessment Land Data Assimilation System (NCA-LDAS). Compared to three other machine learning (ML) algorithms [extremely randomized trees (Extra-Trees), gradient boosting regression tree (GBRT) and extreme gradient boosting (XGBoost)], the random forest (RF) algorithm exhibited the best performance against a subset of 100 validation sites with a Kling–Gupta efficiency (KGE) of 0.69 and root-mean-square-error (RMSE) of 0.063 m3/m3. In terms of different land cover types and typical sites, RF-SM also showed a better accuracy than any of the individual SM product. Finally, the retrieval framework was applied to map the 30-m resolution SM spatial distributions in five substudy areas in the U.S. The results suggest that it is feasible to retrieve accurate SM at a 30-m spatial resolution from multiple satellite datasets based on the RF algorithm, which has important practical significance for agricultural drought monitoring at the field scale.

A Bayesian machine learning method to explain the error characteristics of global-scale soil moisture products

Estimating accurate surface soil moisture (SM) dynamics from space, and knowing the error characteristics of these estimates, is of great importance for the application of satellite-based SM data throughout many Earth Science/Environmental Engineering disciplines. Here, we introduce the Bayesian inference approach to analyze the error characteristics of widely used passive and active microwave satellite-derived SM data sets, at different overpass times, acquired from the Soil Moisture Active Passive (SMAP), Soil Moisture and Ocean Salinity (SMOS), and Advanced Scatterometer (ASCAT) missions. In particular, we apply Bayesian hierarchical modeling (BHM) and triple collocation analysis (TCA) to investigate the relative importance of different environmental factors and human activities on the accuracy of satellite-based data.To start, we compare the BHM-based sensitivity analysis method to the classic multiple regression models using a frequentist approach, which includes complete pooling and no-pooling models that have been widely used for sensitivity analysis in the field of remote sensing and demonstrate the BHM's adaptability and great potential for providing insight into sensitivity analysis that can be used by various remote sensing research communities.Next, we conduct an uncertainty analysis on BHM's model parameters using a full range of uncertainties to assess the association of various environmental factors with the accuracy of satellite-derived SM data. We focus on investigating human-induced error sources such as disturbed surface soil layers caused by irrigation activities on microwave satellite systems, naturally introduced error sources such as vegetation and soil organic matter, and errors related to the disregard of SM retrieval algorithmic assumptions - such as the thermal equilibrium passive microwave systems. Based on the BHM-based sensitivity analysis, we find that assessments of SM data quality with a single variable should be avoided, since numerous other factors simultaneously influence their quality. As such, this provides a useful framework for applying Bayesian theory to the investigation of the error characteristics of satellite-based SM data and other time-varying geophysical variables.

Evaluation and Application of SMRT Model for L-Band Brightness Temperature Simulation in Arctic Sea Ice

Using L-band microwave radiative transfer theory to retrieve ice and snow parameters is one of the focuses of Arctic research. At present, due to limitations of frequency and substrates, few operational microwave radiative transfer models can be used to simulate L-band brightness temperature (TB) in Arctic sea ice. The snow microwave radiative transfer (SMRT) model, developed with the support of the European Space Agency in 2018, has been used to simulate high-frequency TB in polar regions and has obtained good results, but no studies have shown whether it can be used appropriately in the L-band. Therefore, in this study, we systematically evaluate the ability of the SMRT model to simulate L-band TB in the Arctic sea ice and snow environment, and we show that the results are significantly optimized by improving the simulation method. In this paper, we first consider the thermal insulation effect of snow by adding the thermodynamic equation, then use a reasonable salinity profile formula for multi-layer model simulation to solve the problem of excessive L-band penetration in the SMRT single-layer model, and finally add ice lead correction to resolve the large influence it has on the results. The improved SMRT model is evaluated using Operation IceBridge (OIB) data from 2012 to 2015 and compared with the snow-corrected classical L-band radiative transfer model for Arctic sea ice proposed in 2010 (KA2010). The results show that the SMRT model has better simulation results, and the correlation coefficient (R) between SMRT-simulated TB and Soil Moisture and Ocean Salinity (SMOS) satellite TB is 0.65, and the RMSE is 3.11 K. Finally, the SMRT model with the improved simulation method is applied to the whole Arctic from November 2014 to April 2015, and the simulated R is 0.63, and the RMSE is 5.22 K. The results show that the SMRT multi-layer model is feasible for simulating L-band TB in the Arctic sea ice and snow environment, which provides a basis for the retrieval of Arctic parameters.

Assessment of 9 km SMAP soil moisture: Evidence of narrowing the gap between satellite retrievals and model-based reanalysis

A number of global surface soil moisture (SM) datasets have been retrieved from the L-band frequency Soil Moisture Active Passive (SMAP) and the Soil Moisture and Ocean Salinity (SMOS) missions to study the terrestrial water, energy, and carbon cycles. This paper presents the performance of the recently developed 9 km global SMAP product (hereafter SMAP-INRAE-BORDEAUX, SMAP-IB9). The product retrieves SM from the 9 km SMAP radiometric products using the forward model (L-MEB, L-band Microwave Emission of the Biosphere) of SMOS INRA-CESBIO (SMOS-IC) and SMOS L2 algorithms. We inter-compared SMAP-IB9 with two other products with a similar grid resolution (∼10 km): the SMAP Enhanced Level-3 SM dataset (SMAP-E) and the enhanced global dataset for the land component of the fifth generation of European reanalysis (ERA5-Land) with the main objective of assessing the discrepancy in accuracy between remotely sensed and model SM datasets. We found that ERA5-Land and SMAP-IB9 SM had the overall highest correlations (R = 0.62(±0.15) for ERA5-Land vs. 0.60(±0.17) for SMAP-IB9 and 0.50(±0.15) for SMAP-E) by comparing with the International Soil Moisture Network (ISMN) in-situ measurements from 22 networks. ERA5-Land showed better performances in the forest areas where SMAP-IB9 and SMAP-E still showed high potential in detecting the time variations of the observed SM, particularly in terms of median correlation values (0.62(±0.18) for SMAP-IB9 vs. 0.66(±0.16) for ERA5-and). The discrepancy in R between satellite and model SM products that were reported in some past studies has decreased to statistically insignificant levels over time. For instance, in the non-forest areas, we found that the latest versions of the SMAP SM products (SMAP-E and SMAP-IB9) had relatively comparable performances with ERA5-Land with regard to median ubRMSE (0.07(±0.02) m3/m3 for both SMAP-E and ERA5-Land) and R (0.59(±0.16) for SMAP-IB9 vs. 0.61(±0.15) for ERA5-Land), respectively.

A climatology of thermodynamic vs. dynamic Arctic wintertime sea ice thickness effects during the CryoSat-2 era

Abstract. Thermodynamic and dynamic sea ice thickness processes are affected by differing mechanisms in a changing climate. Independent observational datasets of each are essential for model validation and accurate projections of future sea ice conditions. Here, we present a monthly, Arctic-basin-wide, and 25 km resolution Eulerian estimation of thermodynamic and dynamic effects on wintertime sea ice thickness from 2010–2021. Estimates of thermodynamic growth rate are determined by coupling passive microwave-retrieved snow–ice interface temperatures to a simple sea ice thermodynamic model, total growth is calculated from a weekly Alfred Wegener Institute (AWI) European Space Agency (ESA) CryoSat-2 and Soil Moisture and Ocean Salinity (SMOS) combination product (CS2SMOS), and dynamic effects are calculated as their difference. The dynamic effects are further separated into advection and residual effects using a sea ice motion dataset. Our results show new detail in these fields and, when summed to a basin-wide or regional scale, are in line with previous studies. Across the Arctic, dynamic effects are negative and about one-fourth the magnitude of thermodynamic growth. Thermodynamic growth varies from less than 0.1 m per month in the central Arctic to greater than 0.3 m per month in the seasonal ice zones. High positive dynamic effects of greater than 0.1 m per month, twice that of thermodynamic growth or more in some areas, are found north of the Canadian Arctic Archipelago, where the Transpolar Drift and Beaufort Gyre deposit ice. Strong negative dynamic effects of less than −0.2 m per month are found where the Transpolar Drift originates, nearly equal to and opposite the thermodynamic effects in these regions. Monthly results compare well with a recent study of the dynamic and thermodynamic effects on sea ice thickness along the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) drift track during the winter of 2019–2020. Couplets of deformation and advection effects with opposite signs are common across the Arctic, with positive advection effects and negative deformation effects found in the Beaufort Sea and negative advection effects and positive deformation effects found in most other regions. The seasonal cycle shows residual deformation effects and overall dynamic effects increasing as the winter season progresses.

Polar firn properties in Greenland and Antarctica and related effects on microwave brightness temperatures

Abstract. In studying the mass balance of polar ice sheets, fluctuations in firn density near the surface is a major uncertainty. In this paper, we explore these variations at locations on the Greenland Ice Sheet and at the Dome C location in Antarctica. Borehole in situ measurements, snow radar echoes, microwave brightness temperatures, and modeling results from the Community Firn Model (CFM) are used. It is shown that firn density profiles can be represented using three processes: “long-scale” and “short-scale” density variations and “refrozen layers”. Consistency with this description is observed in the dynamic range of airborne 0.5–2 GHz brightness temperatures and snow radar echo peaks in measurements performed in Greenland in 2017. Based on these insights, a new analytical partially coherent model is implemented to explain the microwave brightness temperatures using the three-scale description of the firn. Short- and long-scale firn processes are modeled as a 3D continuous random medium with finite vertical and horizontal correlation lengths as opposed to past 1D randomly layered medium descriptions. Refrozen layers are described as deterministic sheets with planar interfaces, with the number of refrozen-layer interfaces determined by radar observations. Firn density and correlation length parameters used in forward modeling to match measured 0.5–2 GHz brightness temperatures in Greenland show consistency with similar parameters in CFM predictions. Model predictions also are in good agreement with multi-angle 1.4 GHz vertically and horizontally polarized brightness temperature measured by the Soil Moisture and Ocean Salinity (SMOS) satellite at Dome C, Antarctica. This work shows that co-located active and passive microwave measurements can be used to infer polar firn properties that can be compared with predictions of the CFM. In particular, 0.5–2 GHz brightness temperature measurements are shown to be sensitive to long-scale firn density fluctuations with density standard deviations in the range of 0.01–0.06 g cm−3 and vertical correlation lengths of 6–20 cm.

A global analysis of water storage variations from remotely sensed soil moisture and daily satellite gravimetry

Abstract. Water storage changes in the soil can be observed on a global scale with different types of satellite remote sensing. While active or passive microwave sensors are limited to the upper few centimeters of the soil, satellite gravimetry can detect changes in the terrestrial water storage (TWS) in an integrative way, but it cannot distinguish between storage variations in different compartments or soil depths. Jointly analyzing both data types promises novel insights into the dynamics of subsurface water storage and of related hydrological processes. In this study, we investigate the global relationship of (1) several satellite soil moisture products and (2) non-standard daily TWS data from the Gravity Recovery and Climate Experiment/Follow-On (GRACE/GRACE-FO) satellite gravimetry missions on different timescales. The six soil moisture products analyzed in this study differ in the post-processing and the considered soil depth. Level 3 surface soil moisture data sets of the Soil Moisture Active Passive (SMAP) and Soil Moisture and Ocean Salinity (SMOS) missions are compared to post-processed Level 4 data products (surface and root zone soil moisture) and the European Space Agency Climate Change Initiative (ESA CCI) multi-satellite product. On a common global 1∘ grid, we decompose all TWS and soil moisture data into seasonal to sub-monthly signal components and compare their spatial patterns and temporal variability. We find larger correlations between TWS and soil moisture for soil moisture products with deeper integration depths (root zone vs. surface layer) and for Level 4 data products. Even for high-pass filtered sub-monthly variations, significant correlations of up to 0.6 can be found in regions with a large, high-frequency storage variability. A time shift analysis of TWS versus soil moisture data reveals the differences in water storage dynamics with integration depth.

Ice Freeze‐Up and Break‐Up in Arctic Rivers Observed With Satellite L‐Band Passive Microwave Data From 2010 to 2020

AbstractThe timing of ice freeze‐up and break‐up in the Arctic may be responding to climate change. Passive microwave remote sensing is a powerful technique for monitoring this timing. We processed low‐frequency microwave time series from the European Space Agency Soil Moisture and Ocean Salinity (SMOS) mission for a set of 31 satellite gauging reaches (SGRs) above 65°N between 2010 and 2020 to determine timing of freeze‐up and break‐up and annual river ice durations. We found indication of progressive ice cover reduction over more than half of the monitored river reaches, with possibly the fastest rate occurring over northeast Russia. Some rivers in high‐latitude North America experienced a slight increase in ice cover. Across the data set, we observed an average 2.2 days shift toward later ice freeze‐up in autumn and an average 0.6 days shift toward earlier ice break‐up in spring, resulting in an average decrease of 3.4 days in ice duration between 2010 and 2020. River reaches with the longest duration of ice cover appeared to have experienced the fastest rate of decrease. A possible reduction of the time lag between air temperature rise or fall and corresponding river ice break‐up and freeze‐up was also observed. Yet results on variability are carefuly interpreted given the short length of the time series (2010–2020) and the low statistical confidence rates calculated for the decadal tendency. Still outcomes are consistent with increases in global and Arctic surface air temperature. Following these time series over the next decade using passive microwave satellite sensors can monitor ice cover duration in the Arctic and will further determine temporal and regional trends.

Assessing and comparing the subseasonal variations of summer soil moisture of satellite products over eastern China

AbstractUsing the 43‐station soil moisture (SM) at a depth of 0–10 cm over eastern China, the near‐surface (about 0–5 cm) SM products from satellite missions including the level‐2 product of Chinese FengYun 3C (FY3C), level‐2 neural network product of the European Soil Moisture and Ocean Salinity (SMOS) and level‐4 product of the US Soil Moisture Active Passive (SMAP) were evaluated and compared during summers of 2015–2018. Thus, features of those products can be identified for their further application of climate study over eastern China. Large diversity is found among the satellite and field SM in terms of spatial distribution. This disagreement may cause the poor performance of the relationship between spatial coefficients of variation (CV) and mean SM. Compared with the field SM, the FY3C SM has smaller bias than the other satellite products. FY3C also performs better than SMOS in terms of the root‐mean‐square error (RMSE) and unbiased RMSE (ubRMSE) but has the smallest correlation coefficient (R). The SMAP product is generally the best among the three products in terms of RMSE, ubRMSE and R. However, a good performance in those metrics does not guarantee the same results on various time scales. On subseasonal time scales, the R in FY3C is the smallest among the three products, and the SMAP product has the largest R, but its amplitudes of the subseasonal variations are much smaller than the field observations. This indicates that when the SMAP products are applied for the analysis on subseasonal SM–atmosphere interactions, the effects of SM may be underestimated. On 10–30 days and above 60 days, dry period tends to have large spatial CV but this phenomenon is weak in all the satellite products. On the other hand, dry area tends to have large temporal CV on the four time scales, and this relationship is the strongest in SMAP but the weakest in FY3C. Therefore, there are large uncertainties of variability among different satellite SM products on subseasonal time scales over eastern China. Besides seasonal and overall performance, more attention is called to the variations on different time scales.

A Performance Analysis of Soil Dielectric Models over Organic Soils in Alaska for Passive Microwave Remote Sensing of Soil Moisture

Passive microwave remote sensing of soil moisture (SM) requires a physically based dielectric model that quantitatively converts the volumetric SM into the soil bulk dielectric constant. Mironov 2009 is the dielectric model used in the operational SM retrieval algorithms of the NASA Soil Moisture Active Passive (SMAP) and the ESA Soil Moisture and Ocean Salinity (SMOS) missions. However, Mironov 2009 suffers a challenge in deriving SM over organic soils, as it does not account for the impact of soil organic matter (SOM) on the soil bulk dielectric constant. To this end, we presented a comparative performance analysis of nine advanced soil dielectric models over organic soil in Alaska, four of which incorporate SOM. In the framework of the SMAP single-channel algorithm at vertical polarization (SCA-V), SM retrievals from different dielectric models were derived using an iterative optimization scheme. The skills of the different dielectric models over organic soils were reflected by the performance of their respective SM retrievals, which was measured by four conventional statistical metrics, calculated by comparing satellite-based SM time series with in-situ benchmarks. Overall, SM retrievals of organic-soil-based dielectric models tended to overestimate, while those from mineral-soil-based models displayed dry biases. All the models showed comparable values of unbiased root-mean-square error (ubRMSE) and Pearson Correlation (R), but Mironov 2019 exhibited a slight but consistent edge over the others. An integrated consideration of the model inputs, the physical basis, and the validated accuracy indicated that the separate use of Mironov 2009 and Mironov 2019 in the SMAP SCA-V for mineral soils (SOM <15%) and organic soils (SOM ≥15%) would be the preferred option.

Assessing the sensitivity of multi-frequency passive microwave vegetation optical depth to vegetation properties

Abstract. Vegetation attenuates the microwave emission from the land surface. The strength of this attenuation is quantified in models in terms of the parameter vegetation optical depth (VOD) and is influenced by the vegetation mass, structure, water content, and observation wavelength. Earth observation satellite sensors operating in the microwave frequencies are used for global VOD retrievals, enabling the monitoring of vegetation at large scales. VOD has been used to determine above-ground biomass, monitor phenology, or estimate vegetation water status. VOD can be also used for constraining land surface models or modelling wildfires at large scales. Several VOD products exist, differing by frequency/wavelength, sensor, and retrieval algorithm. Numerous studies present correlations or empirical functions between different VOD datasets and vegetation variables such as the normalized difference vegetation index, leaf area index, gross primary production, biomass, vegetation height, or vegetation water content. However, an assessment of the joint impact of land cover, vegetation biomass, leaf area, and moisture status on the VOD signal is challenging and has not yet been done. This study aims to interpret the VOD signal as a multi-variate function of several descriptive vegetation variables. The results will help to select VOD at the most suitable wavelength for specific applications and can guide the development of appropriate observation operators to integrate VOD with large-scale land surface models. Here we use VOD from the Land Parameter Retrieval Model (LPRM) in the Ku, X, and C bands from the harmonized Vegetation Optical Depth Climate Archive (VODCA) dataset and L-band VOD derived from Soil Moisture and Ocean Salinity (SMOS) and Soil Moisture Active Passive (SMAP) sensors. The leaf area index, live-fuel moisture content, above-ground biomass, and land cover are able to explain up to 93 % and 95 % of the variance (Nash–Sutcliffe model efficiency coefficient) in 8-daily and monthly VOD within a multi-variable random forest regression. Thereby, the regression reproduces spatial patterns of L-band VOD and spatial and temporal patterns of Ku-, X-, and C-band VOD. Analyses of accumulated local effects demonstrate that Ku-, X-, and C-band VOD are mostly sensitive to the leaf area index, and L-band VOD is most sensitive to above-ground biomass. However, for all VODs the global relationships with vegetation properties are non-monotonic and complex and differ with land cover type. This indicates that the use of simple global regressions to estimate single vegetation properties (e.g. above-ground biomass) from VOD is over-simplistic.

Assimilation of sea surface salinities from SMOS in an Arctic coupled ocean and sea ice reanalysis

Abstract. In the Arctic, the sea surface salinity (SSS) plays a key role in processes related to water mixing and sea ice. However, the lack of salinity observations causes large uncertainties in Arctic Ocean forecasts and reanalysis. Recently the Soil Moisture and Ocean Salinity (SMOS) satellite mission was used by the Barcelona Expert Centre to develop an Arctic SSS product. In this study, we evaluate the impact of assimilating this data in a coupled ocean–ice data assimilation system. Using the deterministic ensemble Kalman filter from July to December 2016, two assimilation runs respectively assimilated two successive versions of the SMOS SSS product on top of a pre-existing reanalysis run. The runs were validated against independent in situ salinity profiles in the Arctic. The results show that the biases and the root-mean-squared differences (RMSD) of SSS are reduced by 10 % to 50 % depending on the area and highlight the importance of assimilating satellite salinity data. The time series of freshwater content (FWC) further shows that its seasonal cycle can be adjusted by assimilation of the SSS products, which is encouraging of the assimilation of SSS in a long-time reanalysis to better reproduce the Arctic water cycle.

Deep learning for daily spatiotemporally continuity of satellite surface soil Moisture over eastern China in summer

Except assimilated data of soil moisture (SM), there is always a circumstance that not every spatial point has valid value over a specific region in a day, and there are invalid values or gaps between days in a time series at a spatial point, especially the satellite and ground-based data. Thus, those data are not spatiotemporally seamless on daily time scale, and the seamless data of surface SM is fundamental for research on the land–atmosphere interactions on various time scales. Deep learning methods of multi-perceptron (MLP) and convolution (Conv) neural networks are applied for obtaining spatiotemporally seamless SM over eastern China during May-September of 2015–2020 to combine the advantages of ground-based and remote-sensed data. Firstly, taking atmospheric and land variables as inputs, MLP is trained against the satellite products, including the level-2 product of Chinese FengYun 3C (FY3C), level-3 product of the US Soil Moisture Active Passive (SMAP), and level-2 neural network product of the European Soil Moisture and Ocean Salinity (SMOS). After that, spatiotemporally seamless SM is predicted by MLP, which is stable during training and validation. Compared with the in-situ SM, MLP SM inherits the characteristics of the satellite products, though the quality is slightly improved. Secondly, Conv uses the MLP SM as inputs and is trained against in-situ SM. The unbiased errors are reduced notably that unbiased correlation increases and root mean square errors decreases. At last, from the pattern correlation between precipitation and SM anomalies, the error-reduced Conv SM can well reflect the relationships forced out by the major patterns of rainfall, but the satellite products cannot. In addition, during the evaluation, SMAP has the best quality over eastern China, and SMOS is generally better than FY3C. However, after processed by Conv, the differences of the quality among different satellites are negligible, except some differences in the mean patterns of Conv SM, which is mainly associated with the properties of the satellite products.

Towards soil moisture profile estimation in the root zone using L- and P-band radiometer observations: A coherent modelling approach

Precision irrigation management and crop water stress assessment rely on accurate estimation of root zone soil moisture. However, only the top 5 cm soil moisture can be estimated using the two current passive microwave satellite missions, Soil Moisture and Ocean Salinity (SMOS) and Soil Moisture Active Passive (SMAP), which operate at L-band (wavelength of ∼21 cm). Since the contributing depth of the soil to brightness temperature increases with observation wavelength, it is expected that a P-band (wavelength of ∼40 cm) radiometer could potentially provide soil moisture information from deeper layers of the soil profile. Moreover, by combining both L- and P- bands, it is hypothesized that the soil moisture profile can be estimated even beyond their individual observation depths. The aim of this study was to demonstrate the potential of combined L-band and P-band radiometer observations to estimate the soil moisture profile under flat bare soil using a stratified coherent forward model. Brightness temperature observations at L-band and P-band from a tower based experimental site across a dry (April 2019) and a wet (March 2020) period, covering different soil moisture profile shapes, were used in this study. Results from an initial synthetic study showed that the performance of a combined L-band and P-band approach was better than the performance of using either band individually, with an average depth over which reliable soil moisture profile information could be estimated (i.e. with a target root mean square error (RMSE) of less than 0.04 m3/m3) being 20 cm for linear and 15 cm for second-order polynomial functions. Other functions were also tested but found to have a poorer performance. Applying the method to the tower-based brightness temperature achieved an average estimation depth of 28 cm (20 cm) and 5 cm (5 cm) during the dry and wet periods respectively when using a second-order polynomial (linear) function. These findings highlight the opportunity of a satellite mission with L-band and P-band observations to accurately estimate the soil moisture profile to as deep as 30 cm globally.

Impact of the Antenna Spacing on the Brightness Temperature Maps Retrieved with a Synthetic Aperture Imaging Radiometer

For almost five years, phase 0 and phase A studies have been conducted by the French space agency for a second generation of the Soil Moisture and Ocean Salinity (SMOS) satellite devoted to a High Resolution (HR) follow-on mission. Within the frame of the preliminary results obtained with a candidate array for this SMOS-HR project, this contribution focuses on what could happen to any synthetic aperture imaging radiometer when the shortest spacing between the antennas of the interferometric array becomes smaller than a geometrical limit below which the synthesized field of view seems to be wider than the field of view seen by each elementary antenna. It is shown that in such a situation, the inversion of the complex visibilities becomes unstable in presence of noise and this instability is characterized by the undesirable presence of a phantom in the retrieved brightness temperature maps. The origin of this phantom is explained and a solution to cure the interferometric array from that issue is proposed and assessed with the aid of numerical simulations.

Passive Only Microwave Soil Moisture Retrieval in Indian Cropping Conditions: Model Parameterization and Validation

The present study carried out to parameterize the single channel soil moisture active passive (SMAP) passive soil moisture (SM) retrieval algorithm, over Indian conditions. The moderate resolution imaging spectroradiometer (MODIS) data products and soil texture data were used for an improved parameterization of the algorithm. The bias correction was applied to the MODIS leaf area index (LAI) for accurate computation of vegetation optical depth. The necessary vegetation and roughness parameter were calibrated through minimization of the error between model retrieved and ground measured SM. The value of root mean square error (RMSE) for retrieved SM was found as <inline-formula xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink"> <tex-math notation="LaTeX">$0.059\,\,m^{3}m^{-3}$ </tex-math></inline-formula> with bias and correlation coefficients of <inline-formula xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink"> <tex-math notation="LaTeX">$0.036\,\,m^{3}m^{-3}$ </tex-math></inline-formula> and 0.724 for ascending overpass, respectively, while a lower value was recorded (RMSE = <inline-formula xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink"> <tex-math notation="LaTeX">$0.059\,\,m^{3}m^{-3}$ </tex-math></inline-formula> , bias = <inline-formula xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink"> <tex-math notation="LaTeX">$0.024\,\,m^{3}m^{-3}$ </tex-math></inline-formula> , and correlation coefficients = 0.752) for descending overpass. The same method is also implemented on two other test sites in different regions of India to check the model robustness, which indicates that the current parameterization provides a better estimate of SM over croplands in India. The overall performance of new parameterized model is found as (RMSE = 0.052 and bias = 0.034) for ascending and descending (RMSE = 0.048 and bias = 0.026) satellite overpasses for all the three test sites. Additionally, the intercomparing of various operational SM products SMAP SM (L2_SM_P), Soil Moisture and Ocean Salinity (SMOS) SM (SMOS_L3_SM), and SMOS-IC data products was carried out with the SAC-ISRO PAN India SM network, which showed a significant RMSE, dry and wet biases over all three test sites as compared to the developed improved parameterized algorithm.

A Soil Moisture Retrieval Method for Reducing Topographic Effect: A Case Study on the Qinghai–Tibetan Plateau With SMOS Data

The topography can be very important for passive microwave remote sensing of soil moisture due to its complex influence on the emitted brightness temperature observed by a satellite microwave radiometer. In this study, a methodology of using the first brightness Stokes parameter (i.e., the sum of vertical and horizontal polarization brightness temperature) observed by the Soil Moisture and Ocean Salinity (SMOS) was proposed to improve the soil moisture retrieval under complex topographic conditions. The applicability of the proposed method is validated using in-situ soil moisture measurements collected at four networks (Pali, Naqu, Maqu and Wudaoliang) on the Qinghai-Tibetan Plateau. The results over Pali, which is a typical mountainous area, showed that soil moisture retrievals using the first brightness Stokes parameter are in better agreement with the in-situ measurements (the correlation coefficient R >0.75 and unbiased root mean square error < 0.04 m <sup xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">3</sup> /m <sup xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">3</sup> ) compared with that using the single-polarization brightness temperature. At the other three networks with relatively flatter terrains, soil moisture retrievals using the first brightness Stokes parameter are found to be comparable to the single-polarization retrievals. On the contrary, the maximum bias of the retrieved soil moisture caused by topographic effects exceeds 0.1 m <sup xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">3</sup> /m <sup xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">3</sup> when using vertical or horizontal polarization alone, which is far beyond the expected accuracy (0.04 m <sup xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">3</sup> /m <sup xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">3</sup> ) of SMOS satellite. In the regions on the Qinghai-Tibetan Plateau where the vegetation effect can be ignored, soil moisture retrieved using horizontal polarization brightness temperature is generally underestimated, overestimated when using vertical polarization brightness temperature. It is reasonable due to the polarization rotation effect (depolarization) caused by the topographic effects. It is concluded that the proposed method for soil moisture retrieval using the first brightness Stokes parameter has a great potential in reducing the influence of topographic effects.