Satellite-derived Surface Soil Moisture Research Articles

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.

Soil Moisture Assimilation Improves Terrestrial Biosphere Model GPP Responses to Sub-Annual Drought at Continental Scale

Due to the substantial gross exchange fluxes with the atmosphere, the terrestrial carbon cycle plays a significant role in the global carbon budget. Drought commonly affects terrestrial carbon absorption negatively. Terrestrial biosphere models exhibit significant uncertainties in capturing the carbon flux response to drought, which have an impact on estimates of the global carbon budget. Through plant physiological processes, soil moisture tightly regulates the carbon cycle in the environment. Therefore, accurate observations of soil moisture may enhance the modeling of carbon fluxes in a model–data fusion framework. We employ the Carbon Cycle Data Assimilation System (CCDAS) to assimilate 36-year satellite-derived surface soil moisture observations in combination with flask samples of atmospheric CO2 concentrations. We find that, compared to the default model, the performance of optimized net ecosystem productivity (NEP) and gross primary productivity (GPP) has increased with the RMSEs reduced by 1.62 gC/m2/month and 10.84 gC/m2/month, which indicates the added value of the ESA-CCI soil moisture observations as a constraint on the terrestrial carbon cycle. Additionally, the combination of soil moisture and CO2 concentration in this study improves the representation of inter-annual variability of terrestrial carbon fluxes as well as the atmospheric CO2 growth rate. We thereby investigate the ability of the optimized GPP in responding to drought by comparing continentally aggregated GPP with the drought index. The assimilation of surface soil moisture has been shown to efficiently capture the influences of the sub-annual (≤9 months drought durations) and large-scale (e.g., regional to continental scales) droughts on GPP. This study highlights the significant potential of satellite soil moisture for constraining inter-annual models of the terrestrial biosphere’s carbon cycle and for illustrating how GPP responds to drought at a continental scale.

Benchmarking downscaled satellite-based soil moisture products using sparse, point-scale ground observations

While strides have been made in their accuracy and availability, the overall utility of satellite-derived surface soil moisture (SM) datasets derived from passive microwave radiometry is still reduced by their relatively coarse spatial resolution (typically >30 km). In response to this shortcoming, many independent satellite-based SM downscaling approaches have been introduced recently. However, owing to limitations in the spatial sampling characteristics of existing SM ground-monitoring networks, it has proven difficult to obtain reliable reference SM observations at the target downscaling resolution for these approaches (typically 1 to 10 km). As a result, the objective evaluation of SM downscaling approaches is often challenging and/or limited to very localized conditions. Here, we introduce and evaluate a point-scale downscaling (PSD) benchmarking strategy whereby spatially sparse, long-term, point-scale SM observations available from existing ground-based SM networks are utilized for the objective benchmarking of downscaled satellite-based SM products. First, we define criteria that must be met for a given SM downscaling strategy to add either temporal accuracy or spatial skill relative to its coarse-resolution SM baseline. Next, we illustrate, both analytically and numerically, that such criteria can be accurately evaluated using sparse, point-scale SM observations available from existing ground-based SM networks. Finally, we apply our new PSD benchmarking approach to evaluate existing fine-scale SM products. Results demonstrate that the PSD approach, in concert with existing ground-based network data, can be leveraged to robustly evaluate SM downscaling approaches.

Evaluation of Satellite-Derived Surface Soil Moisture Products over Agricultural Regions of Canada

Soil moisture is a critical indicator for climate change and agricultural drought, but its measurement is challenging due to large variability with land cover, soil type, time, space and depth. Satellite estimates of soil moisture are highly desirable and have become more widely available over the past decade. This study investigates and compares the performance of four surface soil moisture satellite datasets over Canada, namely, Soil Moisture and Ocean Salinity Level 3 (SMOS L3), versions 3.3 and 4.2 of European Space Agency Climate Change Initiative (ESA CCI) soil moisture product and a recent product called SMOS-INRA-CESBIO (SMOS-IC) that contains corrections designed to reduce several known sources of uncertainty in SMOS L3. These datasets were evaluated against in situ networks located in mostly agricultural regions of Canada for the period 2012 to 2014. Two statistical comparison methods were used, namely, metrics for mean soil moisture and median of metrics. The results suggest that, while both methods show similar comparisons for regional networks, over large networks, the median of metrics method is more representative of the overall correlation and variability and is therefore a more appropriate method for evaluating the performance of satellite products. Overall, the SMOS products have higher daily temporal correlations, but larger biases, against in situ soil moisture than the ESA CCI products, with SMOS-IC having higher correlations and smaller variability than SMOS L3. The SMOS products capture daily wetting and drying events better than the ESA CCI products, with the SMOS products capturing at least 75% of observed drying as compared to 55% for the ESA CCI products. Overall, for periods during which there are sufficient observations, both SMOS products are more suitable for agricultural applications over Canada than the ESA CCI products, even though SMOS-IC is able to capture soil moisture variability more accurately than SMOS L3.

Towards a Long-Term Reanalysis of Land Surface Variables over Western Africa: LDAS-Monde Applied over Burkina Faso from 2001 to 2018

This study focuses on the ability of the global Land Data Assimilation System, LDAS-Monde, to improve the representation of land surface variables (LSVs) over Burkina-Faso through the joint assimilation of satellite derived surface soil moisture (SSM) and leaf area index (LAI) from January 2001 to June 2018. The LDAS-Monde offline system is forced by the latest European Centre for Medium-Range Weather Forecasts (ECMWF) atmospheric reanalysis ERA5 as well as ERA-Interim former reanalysis, leading to reanalyses of LSVs at 0.25° × 0.25° and 0.50° × 0.50° spatial resolution, respectively. Within LDAS-Monde, SSM and LAI observations from the Copernicus Global Land Service (CGLS) are assimilated with a simplified extended Kalman filter (SEKF) using the CO2-responsive version of the ISBA (Interactions between Soil, Biosphere, and Atmosphere) land surface model (LSM). First, it is shown that ERA5 better represents precipitation and incoming solar radiation than ERA-Interim former reanalysis from ECMWF based on in situ data. Results of four experiments are then compared: Open-loop simulation (i.e., no assimilation) and analysis (i.e., joint assimilation of SSM and LAI) forced by either ERA5 or ERA-Interim. After jointly assimilating SSM and LAI, it is noticed that the assimilation is able to impact soil moisture in the first top soil layers (the first 20 cm), and also in deeper soil layers (from 20 cm to 60 cm and below), as reflected by the structure of the SEKF Jacobians. The added value of using ERA5 reanalysis over ERA-Interim when used in LDAS-Monde is highlighted. The assimilation is able to improve the simulation of both SSM and LAI: The analyses add skill to both configurations, indicating the healthy behavior of LDAS-Monde. For LAI in particular, the southern region of the domain (dominated by a Sudan-Guinean climate) highlights a strong impact of the assimilation compared to the other two sub-regions of Burkina-Faso (dominated by Sahelian and Sudan-Sahelian climates). In the southern part of the domain, differences between the model and the observations are the largest, prior to any assimilation. These differences are linked to the model failing to represent the behavior of some specific vegetation species, which are known to put on leaves before the first rains of the season. The LDAS-Monde analysis is very efficient at compensating for this model weakness. Evapotranspiration estimates from the Global Land Evaporation Amsterdam Model (GLEAM) project as well as upscaled carbon uptake from the FLUXCOM project and sun-induced fluorescence from the Global Ozone Monitoring Experiment-2 (GOME-2) are used in the evaluation process, again demonstrating improvements in the representation of evapotranspiration and gross primary production after assimilation.

An Evaluation of the EnKF vs. EnOI and the Assimilation of SMAP, SMOS and ESA CCI Soil Moisture Data over the Contiguous US

A number of studies have shown that assimilation of satellite derived soil moisture using the ensemble Kalman Filter (EnKF) can improve soil moisture estimates, particularly for the surface zone. However, the EnKF is computationally expensive since an ensemble of model integrations have to be propagated forward in time. Here, assimilating satellite soil moisture data from the Soil Moisture Active Passive (SMAP) mission, we compare the EnKF with the computationally cheaper ensemble Optimal Interpolation (EnOI) method over the contiguous United States (CONUS). The background error–covariance in the EnOI is sampled in two ways: (i) by using the stochastic spread from an ensemble open-loop run, and (ii) sampling from the model spinup climatology. Our results indicate that the EnKF is only marginally superior to one version of the EnOI. Furthermore, the assimilation of SMAP data using the EnKF and EnOI is found to improve the surface zone correlation with in situ observations at a 95 % significance level. The EnKF assimilation of SMAP data is also found to improve root-zone correlation with independent in situ data at the same significance level; however this improvement is dependent on which in situ network we are validating against. We evaluate how the quality of the atmospheric forcing affects the analysis results by prescribing the land surface data assimilation system with either observation corrected or model derived precipitation. Surface zone correlation skill increases for the analysis using both the corrected and model derived precipitation, but only the latter shows an improvement at the 95 % significance level. The study also suggests that assimilation of satellite derived surface soil moisture using the EnOI can correct random errors in the atmospheric forcing and give an analysed surface soil moisture close to that of an open-loop run using observation derived precipitation. Importantly, this shows that estimates of soil moisture could be improved using a combination of assimilating SMAP using the computationally cheap EnOI while using model derived precipitation as forcing. Finally, we assimilate three different Level-2 satellite derived soil moisture products from the European Space Agency Climate Change Initiative (ESA CCI), SMAP and SMOS (Soil Moisture and Ocean Salinity) using the EnOI, and then compare the relative performance of the three resulting analyses against in situ soil moisture observations. In this comparison, we find that all three analyses offer improvements over an open-loop run when comparing to in situ observations. The assimilation of SMAP data is found to perform marginally better than the assimilation of SMOS data, while assimilation of the ESA CCI data shows the smallest improvement of the three analysis products.

Evaluation of a MetOp ASCAT-Derived Surface Soil Moisture Product in Tundra Environments.

Satellite‐derived surface soil moisture data are available for the Arctic, but detailed validation is still lacking. Previous studies have shown low correlations between in situ and modeled data. It is hypothesized that soil temperature variations after soil thaw impact MetOp ASCAT satellite‐derived surface soil moisture (SSM) measurements in wet tundra environments, as C band backscatter is sensitive to changes in dielectric properties. We compare in situ measurements of water content within the active layer at four sites across the Arctic in Alaska (Barrow, Sagwon, Toolik) and Siberia (Tiksi), taken in the spring after thawing and in autumn prior to freezing. In addition to the long‐term measurement fields, where sensors are installed deeper in the ground, we designed a monitoring setup for measuring moisture very close to the surface in the Lena River Delta, Siberia. The volumetric water content (VWC) and soil temperature sensors were placed in the moss organic layer in order to account for the limited penetration depth of the radar signal. ASCAT SSM variations are generally very small, in line with the low variability of in situ VWC. Short‐term changes after complete thawing of the upper organic layer, however, seem to be mostly influenced by soil temperature. Correlations between SSM and in situ VWC are generally very low, or even negative. Mean standard deviation matching results in a comparably high root‐mean‐square error (on average 11%) for predictions of VWC. Further investigations and measurement networks are needed to clarify factors causing temporal variation of C band backscatter in tundra regions.

LDAS-Monde Sequential Assimilation of Satellite Derived Observations Applied to the Contiguous US: An ERA-5 Driven Reanalysis of the Land Surface Variables

Land data assimilation system (LDAS)-Monde, an offline land data assimilation system with global capacity, is applied over the CONtiguous US (CONUS) domain to enhance monitoring accuracy for water and energy states and fluxes. LDAS-Monde ingests satellite-derived surface soil moisture (SSM) and leaf area index (LAI) estimates to constrain the interactions between soil, biosphere, and atmosphere (ISBA) land surface model (LSM) coupled with the CNRM (Centre National de Recherches Météorologiques) version of the total runoff integrating pathways (CTRIP) continental hydrological system (ISBA-CTRIP). LDAS-Monde is forced by the ERA-5 atmospheric reanalysis from the European Center for Medium Range Weather Forecast (ECMWF) from 2010 to 2016 leading to a seven-year, quarter degree spatial resolution offline reanalysis of land surface variables (LSVs) over CONUS. The impact of assimilating LAI and SSM into LDAS-Monde is assessed over North America, by comparison to satellite-driven model estimates of land evapotranspiration from the Global Land Evaporation Amsterdam Model (GLEAM) project, and upscaled ground-based observations of gross primary productivity from the FLUXCOM project. Taking advantage of the relatively dense data networks over CONUS, we have also evaluated the impact of the assimilation against in situ measurements of soil moisture from the USCRN (US Climate Reference Network), together with river discharges from the United States Geological Survey (USGS) and the Global Runoff Data Centre (GRDC). Those data sets highlight the added value of assimilating satellite derived observations compared with an open-loop simulation (i.e., no assimilation). It is shown that LDAS-Monde has the ability not only to monitor land surface variables but also to forecast them, by providing improved initial conditions, which impacts persist through time. LDAS-Monde reanalysis also has the potential to be used to monitor extreme events like agricultural drought. Finally, limitations related to LDAS-Monde and current satellite-derived observations are exposed as well as several insights on how to use alternative datasets to analyze soil moisture and vegetation state.

Using Satellite-Derived Vegetation Products to Evaluate LDAS-Monde over the Euro-Mediterranean Area

Within a global Land Data Assimilation System (LDAS-Monde), satellite-derived Surface Soil Moisture (SSM) and Leaf Area Index (LAI) products are jointly assimilated with a focus on the Euro-Mediterranean region at 0.5∘ resolution between 2007 and 2015 to improve the monitoring quality of land surface variables. These products are assimilated in the CO2 responsive version of ISBA (Interactions between Soil, Biosphere and Atmosphere) land surface model, which is able to represent the vegetation processes including the functional relationship between stomatal aperture and photosynthesis, plant growth and mortality (ISBA-A-gs). This study shows the positive impact on SSM and LAI simulations through assimilating their satellite-derived counterparts into the model. Using independent flux estimates related to vegetation dynamics (evapotranspiration, Sun-Induced Fluorescence (SIF) and Gross Primary Productivity (GPP)), it is also shown that simulated water and CO2 fluxes are improved with the assimilation. These vegetation products tend to have higher root-mean-square deviations in summer when their values are also at their highest, representing 20–35% of their absolute values. Moreover, the connection between SIF and GPP is investigated, showing a linear relationship depending on the vegetation type with correlation coefficient values larger than 0.8, which is further improved by the assimilation.

Building a Flood-Warning Framework for Ungauged Locations Using Low Resolution, Open-Access Remotely Sensed Surface Soil Moisture, Precipitation, Soil, and Topographic Information

Soil moisture (SM) plays an important role in determining the antecedent condition of a watershed, while topographic attributes define how and where SM and rainfall interact to create floods. Based on this principle, we present a method to identify flood risk at a location in a watershed by using remotely sensed SM and open-access information on rainfall, soil properties, and topography. The method consists of three hydrologic modules that represent the generation, transfer, and accumulation of direct runoff. To simplify the modeling and provide timely warnings, the flood risk is ascertained based on frequency of exceedance, with warnings issued if above a specified threshold. The simplicity of the method is highlighted by the use of only three parameters for each watershed of interest, with effective regionalization allowing use in ungauged watersheds. For this proof-of-concept study, the proposed model was calibrated and tested for 65 hydrologic reference stations in the Murray–Darling Basin in Australia over a 35-year study period by using satellite-derived surface SM. The three model parameters were first estimated using the first ten-year data and then the model performance was evaluated through flood threshold exceedance analyses over the remaining 25-year study period. The results for estimated parameters and skill scores showed promise. The three model parameters can be regionalized as a function of watershed characteristics, and/or representative values estimated from neighboring watersheds, allowing use in ungauged basins everywhere.

Sequential assimilation of satellite-derived vegetation and soil moisture products using SURFEX_v8.0: LDAS-Monde assessment over the Euro-Mediterranean area

Abstract. In this study, a global land data assimilation system (LDAS-Monde) is applied over Europe and the Mediterranean basin to increase monitoring accuracy for land surface variables. LDAS-Monde is able to ingest information from satellite-derived surface soil moisture (SSM) and leaf area index (LAI) observations to constrain the interactions between soil–biosphere–atmosphere (ISBA, Interactions between Soil, Biosphere and Atmosphere) land surface model (LSM) coupled with the CNRM (Centre National de Recherches Météorologiques) version of the Total Runoff Integrating Pathways (ISBA-CTRIP) continental hydrological system. It makes use of the CO2-responsive version of ISBA which models leaf-scale physiological processes and plant growth. Transfer of water and heat in the soil rely on a multilayer diffusion scheme. SSM and LAI observations are assimilated using a simplified extended Kalman filter (SEKF), which uses finite differences from perturbed simulations to generate flow dependence between the observations and the model control variables. The latter include LAI and seven layers of soil (from 1 to 100 cm depth). A sensitivity test of the Jacobians over 2000–2012 exhibits effects related to both depth and season. It also suggests that observations of both LAI and SSM have an impact on the different control variables. From the assimilation of SSM, the LDAS is more effective in modifying soil moisture (SM) from the top layers of soil, as model sensitivity to SSM decreases with depth and has almost no impact from 60 cm downwards. From the assimilation of LAI, a strong impact on LAI itself is found. The LAI assimilation impact is more pronounced in SM layers that contain the highest fraction of roots (from 10 to 60 cm). The assimilation is more efficient in summer and autumn than in winter and spring. Results shows that the LDAS works well constraining the model to the observations and that stronger corrections are applied to LAI than to SM. A comprehensive evaluation of the assimilation impact is conducted using (i) agricultural statistics over France, (ii) river discharge observations, (iii) satellite-derived estimates of land evapotranspiration from the Global Land Evaporation Amsterdam Model (GLEAM) project and (iv) spatially gridded observation-based estimates of upscaled gross primary production and evapotranspiration from the FLUXNET network. Comparisons with those four datasets highlight neutral to highly positive improvement.

The effect of satellite-derived surface soil moisture and leaf area index land data assimilation on streamflow simulations over France

Abstract. This study evaluates the impact of assimilating surface soil moisture (SSM) and leaf area index (LAI) observations into a land surface model using the SAFRAN–ISBA–MODCOU (SIM) hydrological suite. SIM consists of three stages: (1) an atmospheric reanalysis (SAFRAN) over France, which forces (2) the three-layer ISBA land surface model, which then provides drainage and runoff inputs to (3) the MODCOU hydro-geological model. The drainage and runoff outputs from ISBA are validated by comparing the simulated river discharge from MODCOU with over 500 river-gauge observations over France and with a subset of stations with low-anthropogenic influence, over several years. This study makes use of the A-gs version of ISBA that allows for physiological processes. The atmospheric forcing for the ISBA-A-gs model underestimates direct shortwave and long-wave radiation by approximately 5 % averaged over France. The ISBA-A-gs model also substantially underestimates the grassland LAI compared with satellite retrievals during winter dormancy. These differences result in an underestimation (overestimation) of evapotranspiration (drainage and runoff). The excess runoff flowing into the rivers and aquifers contributes to an overestimation of the SIM river discharge. Two experiments attempted to resolve these problems: (i) a correction of the minimum LAI model parameter for grasslands and (ii) a bias-correction of the model radiative forcing. Two data assimilation experiments were also performed, which are designed to correct random errors in the initial conditions: (iii) the assimilation of LAI observations and (iv) the assimilation of SSM and LAI observations. The data assimilation for (iii) and (iv) was done with a simplified extended Kalman filter (SEKF), which uses finite differences in the observation operator Jacobians to relate the observations to the model variables. Experiments (i) and (ii) improved the median SIM Nash scores by about 9 % and 18 % respectively. Experiment (iii) reduced the LAI phase errors in ISBA-A-gs but had little impact on the discharge Nash efficiency of SIM. In contrast, experiment (iv) resulted in spurious increases in drainage and runoff, which degraded the median discharge Nash efficiency by about 7 %. The poor performance of the SEKF originates from the observation operator Jacobians. These Jacobians are dampened when the soil is saturated and when the vegetation is dormant, which leads to positive biases in drainage and/or runoff and to insufficient corrections during winter, respectively. Possible ways to improve the model are discussed, including a new multi-layer diffusion model and a more realistic response of photosynthesis to temperature in mountainous regions. The data assimilation should be advanced by accounting for model and forcing uncertainties.

The potential of 2D Kalman filtering for soil moisture data assimilation

Abstract We examine the potential for parameterizing a two-dimensional (2D) land data assimilation system using spatial error auto-correlation statistics gleaned from a triple collocation analysis and the triplet of: (1) active microwave-, (2) passive microwave- and (3) land surface model-based surface soil moisture products. Results demonstrate that, while considerable spatial error auto-correlation exists in the errors for all three products, the inclusion of this information into a 2D assimilation system does not significantly improve the performance of the system relative to a one-dimensional (1D) baseline. This result is explained via an analytical evaluation of the impact of spatial error auto-correlation on the steady-state Kalman gain, which reveals that 2D filtering requires the existence of large auto-correlation differences (between the assimilation model and the assimilated observations) in order to enhance the analysis relative to a 1D filtering baseline. As a result, large error auto-correlations alone (in both the model or observations) are not sufficient to motivate the application of a 2D land assimilation system. These results have important consequences for the development of land data assimilation systems designed to ingest satellite derived surface soil moisture products for water resource and climate applications.

The Impact of Local Acquisition Time on the Accuracy of Microwave Surface Soil Moisture Retrievals over the Contiguous United States

Satellite-derived soil moisture products have become an important data source for the study of land surface processes and related applications. For satellites with sun-synchronous orbits, these products are typically derived separately for ascending and descending overpasses with different local acquisition times. Moreover, diurnal variations in land surface conditions, and the extent to which they are accurately characterized in retrieval algorithms, lead to distinct systematic and random error characteristics in ascending versus descending soil moisture products. Here, we apply two independent evaluation techniques (triple collocation and direct comparison against sparse ground-based observations) to quantify (correlation-based) accuracy differences in satellite-derived surface soil moisture acquired at different local acquisition times. The orbits from different satellites are separated into two overpass categories: AM (12:00 a.m. to 11:59 a.m. Local Solar Time) and PM (12:00 p.m. to 11:59 p.m. Local Solar Time). Results demonstrate how patterns in the accuracy of AM versus PM retrieval products obtained from a variety of active and passive microwave satellite sensors vary according to land cover and across satellite products with different local acquisition times.

Impact of Temporal Autocorrelation Mismatch on the Assimilation of Satellite-Derived Surface Soil Moisture Retrievals

Satellite-based surface soil moisture retrievals are commonly assimilated into ecohydrological models in order to obtain improved profile soil moisture estimates. However, differences in temporal autocorrelation structure between these retrievals and comparable model-based predictions can potentially undermine the efficiency of such assimilation. Here, we conduct a series of synthetic experiments to examine the magnitude of this problem and the potential for detecting the presence of retrieval/model autocorrelation differences using a simple diagnostic procedure. Our synthetic experiments are based on modifying the observation operator within a data assimilation system to artificially induce differences in temporal autocorrelation between assimilated surface soil moisture retrievals and comparable surface soil moisture estimates made by an off-line ecohydrological assimilation model. Results demonstrate that neglecting a mismatch in retrieval/assimilation model autocorrelation can reduce the benefit of surface soil moisture data assimilation. The impact is especially large for soil profiles with limited vertical coupling. However, the presence of this source of retrieval/model autocorrelation misfit is detectable using a simple diagnostic index derived from a time series of soil moisture retrievals and open loop model predictions. Using relatively short data sets (~2 years in length), the diagnostic is capable of identifying worst-case scenarios leading to the most significant degradation of assimilation results.

Suitability of modelled and remotely sensed essential climate variables for monitoring Euro-Mediterranean droughts

Abstract. Two new remotely sensed leaf area index (LAI) and surface soil moisture (SSM) satellite-derived products are compared with two sets of simulations of the ORganizing Carbon and Hydrology In Dynamic EcosystEms (ORCHIDEE) and Interactions between Soil, Biosphere and Atmosphere, CO2-reactive (ISBA-A-gs) land surface models. We analyse the interannual variability over the period 1991–2008. The leaf onset and the length of the vegetation growing period (LGP) are derived from both the satellite-derived LAI and modelled LAI. The LGP values produced by the photosynthesis-driven phenology model of ISBA-A-gs are closer to the satellite-derived LAI and LGP than those produced by ORCHIDEE. In the latter, the phenology is based on a growing degree day model for leaf onset, and on both climatic conditions and leaf life span for senescence. Further, the interannual variability of LAI is better captured by ISBA-A-gs than by ORCHIDEE. In order to investigate how recent droughts affected vegetation over the Euro-Mediterranean area, a case study addressing the summer 2003 drought is presented. It shows a relatively good agreement of the modelled LAI anomalies with the observations, but the two models underestimate plant regrowth in the autumn. A better representation of the root-zone soil moisture profile could improve the simulations of both models. The satellite-derived SSM is compared with SSM simulations of ISBA-A-gs only, as ORCHIDEE has no explicit representation of SSM. Overall, the ISBA-A-gs simulations of SSM agree well with the satellite-derived SSM and are used to detect regions where the satellite-derived product could be improved. Finally, a correspondence is found between the interannual variability of detrended SSM and LAI. The predictability of LAI is less pronounced using remote sensing observations than using simulated variables. However, consistent results are found in July for the croplands of the Ukraine and southern Russia.

Measurement and Simulation of Topographic Effects on Passive Microwave Remote Sensing Over Mountain Areas: A Case Study From the Tibetan Plateau

Knowledge about the surface soil water content is essential because it controls the surface water dynamics and land-atmosphere interaction. In high mountain areas in particular, soil surface water content controls infiltration and flood events. Although satellite-derived surface soil moisture data from passive microwave sensors are readily available for most regions globally, mountainous areas are often excluded from these data (or at least flagged as biased) due to the strong topographic influence on the retrieved signal. Even though a substantial volume of literature is available dealing with topographic effects on spaceborne brightness temperature, no systematic analysis has been reported. Therefore, we present a comprehensive analysis of topographic effects on brightness temperature at C-band using a two-step approach. First, a well-controlled field experiment is carried out using a mobile truck-mounted C-band radiometer to analyze the impact of geometric and adjacent effects on the radiometer signal. Additionally, a comprehensive radiative transfer model is developed accounting for both effects and tested on the ground-based data. Second, recorded Advanced Microwave Scanning Radiometer-Earth Observing System (AMSR-E) data over the Tibetan Plateau were used to analyze the error due to the impact of topography using the developed model. The results of the field experiment clearly show that the geometric effect of a single hill has a much larger impact on brightness temperature compared to the adjacent effect of multiple hills, whereby, due to the geometric effect, the bias is up to +20 K for horizontal and -13 K for vertical polarization. For the adjacent effect, the bias is less than 3 K for both polarizations. Additionally, the developed radio transfer model was able to reproduce both effects with high accuracy. For the AMSR-E data, the model shows that the brightness temperature recorded is biased in the same way as the ground-based measurements and that uncertainties induced by the wide existence of atypical mountain regions in the Tibetan Plateau will have a great impact on the retrieving error (maximum 30%). The largest impact on the retrieval error, on the other hand, is calculated for the soil moisture with a maximum relative error of 44%. The negligible impact can be attributed to false parameterization of the soil texture, soil surface temperature, and sky temperature. Finally, the overall absolute error in the estimated water content is quantified on average with 4%, whereby single pixels indicate a maximum absolute error of up to 16%. In conclusion, we show that recorded spaceborne brightness temperatures are highly biased by topographic effects in mountainous regions using a comprehensive radiative transfer model. Additionally, we suggest using this model to invert the effective surface emissivity of mountain areas for standard processing of higher level data products such as surface soil water content.

Development of a Satellite Land Data Assimilation System Coupled With a Mesoscale Model in the Tibetan Plateau

Soil moisture is the central focus of land surface and atmospheric modeling because it controls surface water and energy fluxes and consequently affects land-atmosphere interactions. Although global or regional satellite-derived surface soil moisture data sets are readily available, knowledge about assimilating them into numerical weather prediction (NWP) models is limited. The methods of assimilating soil moisture products in NWP models have several limitations, and they cannot be applied in near-real-time applications. As a result, this paper focuses on the development of a system [a Land Data Assimilation System coupled with a mesoscale Atmospheric model (LDAS-A)] that couples satellite land data assimilation with a mesoscale model to physically introduce land surface heterogeneities into the mesoscale model. The LDAS-A consists of a sequential LDAS that directly assimilates the lower frequency passive microwave brightness temperatures, and therefore, its use is feasible for near-real-time NWP applications. The LDAS-A was validated for the Tibetan Plateau using surface, radiosonde, and satellite observations. The simulation results show that the LDAS-A effectively improved the land surface variables (i.e., surface soil moisture and skin temperature) compared with the no-assimilation case and that it has the potential to correct uncertainties resulting from model initialization, model-specific parameters, and model forcing on a wider scale. The improved land surface conditions in the LDAS-A improve the land-atmosphere feedback mechanism, and the assimilated results provide better prediction of atmospheric profiles (i.e., potential temperature and specific humidity) than the no-assimilation case when compared with radiosonde soundings. Improvements in solar radiation, in addition to soil moisture, are necessary to introduce realistic land-atmosphere interactions into a mesoscale model.

Impact of Accuracy, Spatial Availability, and Revisit Time of Satellite-Derived Surface Soil Moisture in a Multiscale Ensemble Data Assimilation System

This study evaluates the sensitivity of a multiscale ensemble assimilation system to different configurations of satellite soil moisture observations, namely the retrieval accuracy, spatial availability, and revisit time. We perform horizontally coupled assimilation experiments where pixels are updated not only by observations at the same location but also all in the study domain. Carrying out sensitivity studies within a multiscale assimilation system is a significant advancement over previous studies that used a 1-D assimilation framework where all horizontal grids are uncoupled. Twin experiments are performed with synthetic soil moisture retrievals. The hydrologic modeling system is forced with satellite estimated rainfall, and the assimilation performance is evaluated against model simulations using <i xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">in-situ</i> measured rainfall. The study shows that the assimilation performance is most sensitive to the spatial availability of soil moisture observations, then to revisit time and least sensitive to retrieval accuracy. The horizontally coupled assimilation system performs reasonably well even with large observation errors, and it is less sensitive to retrieval accuracy than the uncoupled system, as reported by previous studies. This suggests that more information may be extracted from satellite soil moisture observations using multiscale assimilation systems resulting in a potentially higher value of such satellite products.

Application of Satellite-Derived Surface Soil Moisture Data to Simulating Seasonal Precipitation by a Simple Soil Moisture Transfer Method

Abstract A simple algorithm is presented for transferring root-zone soil moisture from surface soil moisture data on a global scale. Analysis of offline soil moisture data shows that the climatological relationship between surface and root-zone soil moisture becomes linear when appropriate time lags are applied. The climatological relationship of root-zone soil moisture among different land surface models (LSMs) is also linear; therefore, the root-zone and surface soil moisture obtained from one LSM can be applied to another. The algorithm is then applied to the surface soil moisture observations made by the precipitation radar on board the Tropical Rainfall Measuring Mission precipitation radar (TRMM/PR), and the transferred root-zone soil moisture is input to a general circulation model (GCM) summer—June, July, August—precipitation simulation as the boundary condition. The approach is computationally efficient, and the simulation using the root-zone soil moisture by the transfer method is much better th...