Soil Moisture And Ocean Salinity Research Articles

Simulated Sea Surface Salinity Data from a 1/48° Ocean Model

To study the validation process for sea surface salinity (SSS) we have generated one year (November 2011- October 2012) of simulated satellite and in situ “ground truth” data. This was done using the ECCO (Estimating the Circulation and Climate of the Oceans) 1/48° simulation, the highest resolution global ocean model currently available. The ground tracks of three satellites, Aquarius, SMAP (Soil Moisture Active Passive) and SMOS (Soil Moisture and Ocean Salinity) were extracted and used to sample the model with a gaussian weighting similar to that of the spaceborne sensor ground footprint. This produced simulated level 2 (L2) data. Simulated level 3 (L3) data were then produced by averaging L2 data onto a regular grid. The model was sampled to produce simulated Argo and tropical mooring SSS datasets. The Argo data were combined into a simulated gridded monthly 1° Argo product. The simulated data produced from this effort have been used to study sampling errors, matchups, subfootprint variability and the validation process for SSS at L2 and L3.

Influence of Supraglacial Lakes on Accuracy of Inversion of Greenland Ice Sheet Surface Melt Data in Different Passive Microwave Bands

The occurrence of Supraglacial Lakes (SGLs) may influence the signals acquired with microwave radiometers, which may result in a degree of uncertainty when employing microwave radiometer data for the detection of surface melt. Accurate monitoring of surface melting requires a reasonable assessment of this uncertainty. However, there is a scarcity of research in this field. Therefore, in this study, we computed surface melt in the vicinity of Automatic Weather Stations (AWSs) by employing Defense Meteorological Satellite Program (DMSP) Ka-band data and Soil Moisture and Ocean Salinity (SMOS) satellite L-band data and extracted SGL pixels by utilizing Sentinel-2 data. A comparison between surface melt results derived from AWS air temperature estimates and those obtained with remote sensing inversion in the two different bands was conducted for sites below the mean snowline elevation during the summers of 2016 to 2020. Compared with sites with no SGLs, the commission error (CO) of DMSP morning and evening data at sites where these water bodies were present increased by 36% and 30%, respectively, and the number of days with CO increased by 12 and 3 days, respectively. The omission error (OM) of SMOS morning and evening data increased by 33% and 32%, respectively, and the number of days with OM increased by 17 and 21 days, respectively. Identifying the source of error is a prerequisite for the improvement of surface melt algorithms, for which this study provides a basis.

Contribution of satellite sea surface salinity to the estimation of liquid freshwater content in the Beaufort Sea

Abstract. The hydrography of the Arctic Ocean has experienced profound changes over the last 2 decades. The sea ice extent has declined by more than 10 % per decade, and its liquid freshwater content has increased mainly due to glaciers and sea ice melting. Further, new satellite retrievals of sea surface salinity (SSS) in the Arctic might contribute to better characterizing the freshwater changes in cold regions. Ocean salinity and freshwater content are intimately related such that an increase (decrease) in one entails a decrease (increase) in the other. In this work, we evaluate the freshwater content in the Beaufort Gyre using surface salinity measurements from the satellite radiometric mission Soil Moisture and Ocean Salinity (SMOS) and TOPAZ4b reanalysis salinity at depth, estimating the freshwater content from 2011 to 2019 and validating the results with in situ measurements. The results highlight the underestimation of the freshwater content using reanalysis data in the Beaufort Sea and a clear improvement in the freshwater content estimation when adding satellite sea surface salinity measurements in the mixed layer. The improvements are significant, with up to a 70 % reduction in bias in areas near the ice melting. Our research demonstrates how remotely sensed salinity can assist us in better monitoring the changes in the Arctic freshwater content and understanding key processes related to salinity variations that cause density differences with potential to influence the global circulation system that regulates Earth's climate.

Intercomparisons and Evaluations of Satellite-Derived Arctic Sea Ice Thickness Products

Currently, Arctic sea ice thickness (SIT) data with extensive spatiotemporal coverage primarily comes from satellite observations, including CryoSat-2, Soil Moisture and Ocean Salinity (SMOS), and the Ice, Cloud, and Land Elevation Satellite-2 (ICESat-2). The studies of the intercomparison and evaluation of multi-source satellite products in recent years are limited. In this study, three latest version products of ICESat-2, CryoSat-2, and CS2SMOS (a merged product of CryoSat-2 and SMOS) were examined from October to April, between 2018 and 2022. Three types of observation including airborne data from the Operation IceBridge (OIB) and IceBird, and in situ data from Beaufort Gyre Exploration Project (BGEP) are selected as the reference in the evaluation. The intercomparison results show that the mean SIT is generally largest in ICESat-2, second largest in CryoSat-2, and smallest in CS2SMOS. The SIT in CryoSat-2 is closer to the SIT in ICESat-2. The thickness displayed by the three satellite products starts to increase at different freezing months, varying between October and November. The three satellite products demonstrated the strongest agreements in SIT in the Beaufort Sea and Central Arctic regions, and exhibited the most distinct differences in the Barents Sea. In the evaluation with OIB data, three satellite-derived SIT were generally underestimated and CS2SMOS demonstrates the closest match. The evaluation using IceBird data indicates an underestimation for all satellites, with CryoSat-2 showing the best agreement. In the assessment with BGEP data, ICESat-2 displayed a more pronounced degree of overestimation or underestimation compared to the other two satellites, and CS2SMOS exhibited the optimal agreement. Based on the comprehensive consideration, CS2SMOS demonstrated the best performance with the airborne and in situ observational data, followed by CryoSat-2 and ICESat-2. The intercomparison and evaluation results of satellite products can contribute to a further understanding of the accuracies and uncertainties of the latest version SIT retrieval and the appropriate selection and utilization of satellite products.

Validation of downscaled 1‐km SMOS and SMAP soil moisture data in 2010–2021

AbstractSoil moisture (SM) is an important component for many applications in agriculture, hydrology, meteorology, and ecology. In past decades, passive/active microwave sensors onboard Earth observation satellites are utilized to obtain SM estimates from radiometer or radar observations. In this study, the Soil Moisture and Ocean Salinity (SMOS) Level 3 daily SM retrievals at 25‐km spatial resolution between 2010 and 2021 were downscaled through an apparent thermal inertia principle‐based algorithm. The 1‐km downscaled SMOS SM retrievals were validated by in situ measurements from 635 sites of 19 SM networks in the world, which were acquired from the International Soil Moisture Network and Texas Soil Observation Network. Additionally, the validation results of the SMOS SM products were compared with those of the Soil Moisture Active Passive (SMAP) global Level 2 enhanced SM products at 1‐km downscaled and original 9‐km resolution in 2015–2021. It shows that the downscaled SMOS SM data have an overall improved accuracy and outperform the coarse‐resolution 25‐km data, with a lower unbiased Root Mean Squared Error of 0.114 m3/m3 on average.

Land cover and management effects on ecosystem resistance to drought stress

Abstract. Drought events are projected to become more severe and frequent across many regions in the future, but their impacts will likely differ among ecosystems depending on their ability to maintain functioning during droughts, i.e., ecosystem resistance. Plant species have diverse strategies to cope with drought. As a result, divergent responses of different vegetation types for similar levels of drought severity have been observed. It remains unclear whether such divergence can be explained by different drought duration; co-occurring compounding effects, e.g., heat stress or memory effects; management practices; etc. Here, we provide a global synthesis of vegetation resistance to drought and heat using different proxies for vegetation condition, namely the vegetation optical depth (L-VOD) data from the ESA Soil Moisture and Ocean Salinity (SMOS) passive L-band microwave mission and enhanced vegetation index (EVI) and kernel normalized difference vegetation index (kNDVI) from NASA MODIS. Due to its longer wavelength, L-VOD has the advantage over more commonly used vegetation indices (such as kNDVI, EVI) in that it provides different information on vegetation structure and biomass and suffers from less saturation over dense forests. We apply a linear model accounting for drought and temperature effects to characterize ecosystem resistance by their sensitivity to drought duration and temperature anomalies. We analyze how ecosystem resistance varies with land cover across the globe and investigate the potential effects of forest management and crop irrigation. We compare estimates of ecosystem resistance to drought and heat as retrieved from L-VOD, kNDVI, and EVI products. We find that regions with higher forest fraction show stronger ecosystem resistance to extreme droughts than cropland for all three vegetation proxies. L-VOD indicates that primary forests tend to be more resistant to drought events than secondary forests when controlling for the differences in background climate, but this cannot be detected in EVI and kNDVI. The difference is possibly related to EVI and kNDVI saturation in dense forests. In tropical primary evergreen broadleaf forests, old-growth trees tend to be more resistant to drought than young trees from L-VOD and kNDVI. Irrigation increases the drought resistance of cropland substantially. Forest harvest decreases the drought resistance of forests. Our results suggest that ecosystem resistance can be better monitored using L-VOD in dense forests and highlight the role of forest cover, forest management, and irrigation in determining ecosystem resistance to droughts.

Automated retrieval of sea surface salinity from satellite data and investigation of the spatio-temporal variation of ocean salinity across the West and Central African coasts

Sea Surface Salinity (SSS) is an important variable that affects the physical properties of the oceanic environment. Satellite sensors have been developed to routinely acquire SSS concentrations. However, many in the broader scientific community do not engage with this data, because it requires relevant remote sensing skills and tools to extract and process into formats that can be readily utilized by non-remote sensing experts. To address this issue, we used Python, PostgreSQL database, and Power BI visualization tool, to develop an automated tool, that facilitates the retrieval, processing, analysis, and visualisation of data from the Optimally Interpolated Sea Surface Salinity (OISSS) dataset, comprising data from the Aquarius/SAC-D, Soil Moisture Active Passive (SMAP), and Soil Moisture and Ocean Salinity (SMOS) satellite missions. Using this tool, we extracted long term series of SSS from 2011 to March 2022. In addition, an interactive reporting dashboard was developed in this research to enable users to visualize and interrogate the underlying dataset. Analysis of this data revealed detailed spatiotemporal patterns of SSS variations in the study area across varying spatial and temporal scales. We were able to identify the cold and hot spots of SSS, as well as longitudinal (interand intra-annual) variations in SSS concentrations. We anticipate that the outcome of this research will boost further research in this area, contribute to future climatic change studies and stimulate the democratization of remote sensing data. Climatic changes have been found to have profound effects on ocean salinity. The data obtained in this research could also feed into machine learning models for predicting future impacts of climate change on salinity concentrations across the ocean. In addition, the ability to monitor the salinity on a regular basis can be deployed alongside other relevant variables to forecast future climatic changes. Due to the fact that salinity is an excellent indicator for water cycle changes, the ability to routinely monitor salinity at relatively higher spatial and temporal resolutions is of great significance to the local community, as it could be used to understand current and future water cycle regimes in the area. The water cycle has profound effect on various socio-economic activities including agriculture and food security, ecology, hydrology (water quality), and drought.

Flash Drought and Its Characteristics in Northeastern South America during 2004–2022 Using Satellite-Based Products

The term flash drought describes a special category of drought with rapid onset and strong intensity over the course of days or weeks. To characterize the impact of flash droughts on vegetation coverage, this study assessed the influence of soil water deficits on vegetation dynamics in the northeastern South America region by combining time series of vegetation index, rainfall, and soil moisture based on satellite products at a daily time scale. An 18-year analysis, from 2004 to 2022, of the Normalized Difference Vegetation Index (NDVI), Standard Precipitation Index (SPI), and surface soil moisture (SSM) was performed based on three different satellite remote sensing estimates: the spinning enhanced visible and infrared imager (SEVIRI) and the integrated multi-satellite retrievals algorithm (IMERG), and the soil moisture and ocean salinity (SMOS). The results revealed that flash drought events exerted dramatic impacts on terrestrial ecosystems in the study region during the first two decades of the 2000s, with changes in seasonal and regional vegetation dynamics. Further, the fixed-threshold values to characterize flash drought events were suggested as the timing when the water deficit was less than −1.0 units and vegetation index reached the value equal to +0.3 during five consecutive weeks or more, coupled with soil moisture rates below 40% percentile, leading to a strong region-wide drought throughout the entire region. Additionally, the results of linear least squares trend analyses revealed a negative trend in the pentad-SEVIRI radiance for the solar channel 1 within the semiarid ecosystems of the study region (i.e., the Caatinga biome) that was suggested as a reduction in clouds in the 18 years of the study. Developing combined threshold measures of flash drought based on satellite remote sensing may lead to an accurate assessment of flash drought mitigation.

Bivariate sea-ice assimilation for global-ocean analysis–reanalysis

Abstract. In the last decade, various satellite missions have been monitoring the status of the cryosphere and its evolution. Besides sea-ice concentration data, available since the 1980s, sea-ice thickness retrievals are now ready to be used in global operational prediction and global reanalysis systems. Nevertheless, while univariate algorithms are commonly used to constrain sea-ice area or volume, multivariate approaches have not yet been employed due to the highly non-Gaussian distribution of sea-ice variables together with the low accuracy of thickness observations. This study extends a 3DVar system, called OceanVar, which is routinely employed in the production of global/regional operational/reanalysis products, to process sea-ice variables. The tangent/adjoint versions of an anamorphosis operator are used to locally transform the sea-ice anomalies into Gaussian control variables and back, minimizing in the latter space. The benefit achieved by such a transformation is described. Several sensitivity experiments are carried out using a suite of diverse datasets. The sole assimilation of the CryoSat-2 provides a good spatial representation of thickness distribution but still overestimates the total volume that requires the inclusion of Soil Moisture and Ocean Salinity (SMOS) mission data to converge towards the observation estimates. The intermittent availability of thickness data can lead to potential jumps in the evolution of the volume and requires a dedicated tuning. The use of the merged L4 product CS2SMOS shows the best skill score when validated against independent measurements during the melting season when satellite data are not available. This new sea-ice module is meant to simplify the future coupling with ocean variables.

True global error maps for SMAP, SMOS, and ASCAT soil moisture data based on machine learning and triple collocation analysis

Quantifying the accuracy of the satellite-based soil moisture (SM) data is important for a number of key applications, such as: combining satellite-based SM products for long-term SM analyses, assimilating SM data into land surface models, and providing quality flags to mask bad quality SM data.A range of statistical methods have been proposed to estimate error statistics for large-scale SM datasets including the: instrumental variable (IV) method, triple collocation analysis (TCA), and quadruple collocation analysis (QCDA). While requiring only two input products, the IV method also imposes an additional assumption that one input product possesses serially uncorrelated errors - thus limiting its scope compared to TC. Likewise, QCDA requires four independent SM data products that are difficult to obtain and may not always be available for analysis. Nonetheless, TCA-based methods still cannot provide truly global error maps for satellite SM products due to the limited number of independent SM products and difficulties with baseline TCA assumptions. Moreover, temporal sampling requirements for TCA are often impractical because of low SM retrieval skill in forested and arid areas – as well as in regions prone to radio frequency interference.Here, we seek to fill significant spatial gaps in TCA results using machine learning (ML) and therefore provide spatially complete error maps for the satellite-based SM data products derived from the Soil Moisture Active Passive (SMAP), Soil Moisture and Ocean Salinity (SMOS), and Advanced Scatterometer (ASCAT) systems. Furthermore, we use SHapley Additive exPlanations (SHAP) values, a model-agnostic technique for interpreting ML models, to examine the impact of various environmental conditions on the quality of satellite-based SM retrievals.Globally, and across all three products, 72.0% of missing error information in a TCA-based analysis, due to either the lack of valid data or the inability of TCA to provide reliable results, can be reconstructed from the ensemble prediction mean of the ML models. Overall, we found that 22.7% (a.m.) and 34.2% (p.m.) of the Earth'sSM dynamics (between 60°S to 60°N) have not been investigated properly across all three satellite missions.

Estimation of Surface and Subsurface Soil Moisture Using Microwave Remote Sensing: A Typical Analysis

Accurate measurement and monitoring of surface and subsurface soil moisture is essential for understanding hydrological processes, crop growth modeling, crop water requirement, and climate studies. Accurate measurement of the soil moisture content (SMC) in the root zone is essential for precise irrigation authority and plant water stress evaluation. However, the existing passive microwave satellite missions, Soil Moisture and Ocean Salinity (SMOS) and Soil Moisture Active Passive (SMAP), that operate at L-band, can only estimate the top 5 cm of soil moisture. Microwave remote sensing has proven to be a valuable tool for non-invasive soil moisture estimation. This research aims to investigate and develop a methodology for estimating surface and subsurface soil moisture using microwave data from Sentinel-1. The study was conducted to establish the relationship between surface & the backscatter coefficient derived using the Sentinel-1 SAR microwave remote sensing satellite imagery, and relationship between surface and subsurface soil moisture at different depths, in the Godhra region. Two seasons namely summer (Zaid) and monsoon (Kharif) were taken into consideration to build up the relationship between surface soil moisture and co-polarization backscatter coefficient ( For the summer (Zaid) and monsoon (Kharif) seasons, the co-polarization backscatter coefficient ( and surface soil moisture (0-5, cm) were found to have a correlation in terms of R2 as 0.91 and 0.90, respectively. The study explores the relationship between microwave signals and surface soil moisture content (0-5, cm) and then the relationship between surface soil moisture and soil moisture at various depths were also modeled thereby contributing to improved soil moisture estimation techniques and applications. The value of the coefficient of determination (R2) of surface soil moisture (0-5, cm) to subsurface soil moisture at 6-20 cm, 21-40 cm, and 41-60 cm depths were found to be 0.60, 0.51, and 0.46, respectively, in the summer (Zaid) season. The value of the coefficient of determination (R2) of surface soil moisture (0-5, cm) to subsurface soil moisture at 6-20 cm, 21-40 cm, 41-60 cm, 61-80 cm, and 81-100 cm depths were found to be 0.83, 0.61, 0.51, 0.26, and 0.13, respectively. According to the study, it is observed that the relationship between co-polarization backscatter coefficient ( and soil moisture weakens as the depth of soil moisture increases. Overall, the regression models developed between the co-polarization backscatter coefficient ( and surface soil moisture showed very good results, whereas the regression models developed between the surface soil moisture and soil moisture at various depths showed reasonably acceptable results up to the depth of 60 cm. The findings in the present study suggest that Sentinel-1A C-band SAR data can be used to estimate surface soil moisture. It is also shown in this study that the surface soil moisture can be correlated with the subsurface soil moisture up to the depth of 60 cm, satisfactorily using regression equations.

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.