Land Surface Data Assimilation Research Articles

Evaluation of Daily and Hourly Performance of Multi-Source Satellite Precipitation Products in China’s Nine Water Resource Regions

Satellite precipitation products (SPPs) are of great significance for water resource management and utilization in China; however, they suffer from considerable uncertainty. While numerous researchers have evaluated the accuracy of various SPPs, further investigation is needed to assess their performance across China’s nine major water resource regions. This study used the latest precipitation dataset of the China Meteorological Administration’s Land Surface Data Assimilation System (CLDAS-V2.0) as the benchmark and evaluated the performance of six SPPs—GSMaP, PERSIANN, CMORPH, CHIRPS, GPM IMERG, and TRMM—using six indices: correlation coefficient (CC), root mean square error (RMSE), mean absolute error (MAE), probability of detection (POD), false alarm rate (FAR), and critical success index (CSI), at both daily and hourly scales across China’s nine water resource regions. The conclusions of this study are as follows: (1) The performance of the six SPPs was generally weaker in the west than in the east, with the Continental Basin (CB) exhibiting the poorest performance, followed by the Southwest Basin (SB). (2) At the hourly scale, the performance of the six SPPs was weaker compared to the daily scale, particularly in the high-altitude CB and the high-latitude Songhua and Liaohe River Basin (SLRB), where observing light precipitation and snowfall presents significant challenges. (3) GSMaP, CMORPH, and GPM IMERG demonstrated superior overall performance compared to CHIRPS, PERISANN, and TRMM. (4) CMORPH was found to be better suited for application in drought-prone areas, showcasing optimal performance in the CB and SB. GSMaP excelled in humid regions, displaying the best overall performance in the remaining seven basins. GPM IMERG serves as a complementary precipitation data source for the first two.

Comparison of Met Office regional model soil moisture with COSMOS‐UK field‐scale in situ observations

AbstractThe UK Met Office state‐of‐the‐art, deterministic, convection‐permitting, coupled land‐atmosphere, regional weather forecasting system, known as the UKV or UK Variable resolution model (Tang et al. Meteorological Applications, 2013; 20:417–426), has been operational since 2015. Science updates are regularly made to the UKV land surface data assimilation scheme when those updates improve predictions of screen temperature and humidity, since these quantities have a direct impact on atmospheric states and weather forecasts. Less attention has been paid to whether UKV soil moisture analyses are close to independent, in‐situ soil moisture observations, partly because it is difficult to make meaningful comparisons between 1.5 km2 gridded model outputs and traditional point sensor measurements. Soil moisture is recognized to be important when hydrological forecasts for runoff and rivers are required. This is because soil moisture controls the extent to which rainfall can infiltrate the soil, and the amount of surface runoff affects the timing of peak river flows (Ward & Robinson, Principles of Hydrology. McGraw‐Hill Publishing Company; 2000; Singh et al. Water Resources Research, 2021, 57, e2020WR028827). Gómez et al. (Remote Sensing, 2020; 12:3691) report benefits to river flow forecasts when using soil moisture data assimilation in the UKV system instead of a daily downscaled product from the Met Office global model. The Met Office measures soil temperature and soil moisture at Cardington (Osborne & Weedon, Journal of Hydrometeorology, 2021, 22:279–295); there is no other UK Met Office site at which soil moisture is measured. In this study, we use field‐scale (~200 m radius) soil moisture measurements from the UK Centre for Ecology and Hydrology's (UKCEH's) COSMOS‐UK network to provide independent verification and analysis of UKV soil moisture during summer 2018, an unusually dry period in the United Kingdom. We find that the match to COSMOS‐UK soil moisture observations is generally good, and that changes made to the land data assimilation approach during a recent operational upgrade had a generally beneficial impact on UKV soil moisture analyses under very dry conditions.

Using a Reanalysis-Driven Land Surface Model for Initialization of a Numerical Weather Prediction System

Abstract Realistic initialization of the land surface is important to produce accurate NWP forecasts. Therefore, making use of available observations is essential when estimating the surface state. In this work, sequential land surface data assimilation of soil variables is replaced with an offline cycling method. To obtain the best possible initial state for the lower boundary of the NWP system, the land surface model is rerun between forecasts with an analyzed atmospheric forcing. We found a relative reduction of 2-m temperature root-mean-square errors and mean errors of 6% and 12%, respectively, and 4.5% and 11% for 2-m specific humidity. During a convective event, the system was able to produce useful (fractions skill score greater than the uniform forecast) forecasts [above 30 mm (12 h)−1] down to a 100-km length scale where the reference failed to do so below 200 km. The different precipitation forcing caused differences in soil moisture fields that persisted for several weeks and consequently impacted the surface fluxes of heat and moisture and the forecasts of screen level parameters. The experiments also indicate diurnal- and weather-dependent variations of the forecast errors that give valuable insight on the role of initial land surface conditions and the land–atmosphere interactions in southern Scandinavia.

Quick estimation of parameters for the land surface data assimilation system and its influence based on the extended Kalman filter and automatic differentiation

Quick estimation of parameters for the land surface data assimilation system and its influence based on the extended Kalman filter and automatic differentiation

Evaluation of Long Time-Series Soil Moisture Products Using Extended Triple Collocation and In Situ Measurements in China

Currently, satellite-based soil moisture (SM) products and land surface model assimilation techniques are widely utilized. However, the presence of systematic errors in the observation process, algorithmic discrepancies between products, and variations in spatial and temporal scales result in diverse accuracy characteristics and applicability. This study evaluates three prominent SM products in China, namely, the Essential Climate Variable Soil Moisture (ECV), the European Centre for Medium-Range Weather Forecasts’ Fifth-Generation Land Surface Reanalysis Data (ERA5-Land), and the Global Land Surface Data Assimilation System (GLDAS). The evaluation was conducted using extended triple collocation (ETC) analysis and in situ validation methods at a monthly scale from 2000 to 2020. The ETC analysis results show that among the three products, GLDAS exhibits the highest correlation coefficient (CC) and the lowest standard deviation of error (ESD), indicating its superior performance in China. ECV and ERA5-Land follow, with slightly lower performance. In the in situ validation results, ERA5-Land displays the highest correlation, capturing the temporal trend of the ground SM well. Comparatively, in terms of overall accuracy, ECV performs the best, with a slightly smaller mean error (ME) and root mean square error (RMSE) than GLDAS, and ERA5-Land has the lowest accuracy. The discrepancy between the in situ validation results and ETC analysis can be attributed to the limited number of sites and their representativeness errors. Notably, ERA5-Land exhibits a highly consistent trend of interannual fluctuations between ESD and precipitation. Furthermore, a strong association is observed between the ME and RMSE of ECV and GLDAS and precipitation. These findings serve as valuable references for future SM studies in China.

Prediction of apple first flowering date using daily land surface temperature spatio-temporal reconstruction and machine learning

The first flowering date (FFD) is a critical phenological parameter closely related to the apple yield, so the accurate prediction of the FFD is important for precise orchard production management. Existing methods to predict the FFD are mostly based on air temperature (Ta) measured at meteorological stations, but to great differences in meteorological variations and the ecological conditions, these methods cannot present the differences of FFD under complex meteorological conditions and provide spatially continuous FFD information at the level of a region. Therefore, we propose a method to predict spatially continuous apple FFD from remote sensing land surface temperature (LST) based on flowering prediction model. Firstly, the missing LST data were reconstructed by using spatio-temporal reconstruction (STR) approach developed. Next, new air temperature (NAT) data were generated by using the daily Ta estimation (DTE) model and the reconstructed LST. Finally, apple FFD was predicted by the NAT data and the apple flowering prediction model established based on random forest (RF) algorithm and the phenology sequential model, and the prediction accuracy was verified by comparison with the independently measured apple FFD. The LST reconstructed by using the STR approach has mean absolute error (MAE) ranging from 0.51 to 0.68 °C, and root mean square error (RMSE) ranging from 1.07 to 1.21 °C. The MAE between the NAT data and the High-Resolution Land Surface Data Assimilation System (HR-CLDAS) meteorological data ranges from 2.15 to 3.23 °C, and the RMSE ranges from 2.81 to 4.27 °C. In addition, the determination coefficient (R2) and RMSE between the predicted and measured FFD is 0.72 and 2.96 days, respectively. These results demonstrate that the developed method maximizes the potential of MODIS LST in predicting spatially continuous apple FFD, which is valuable for flower and fruit thinning, to defend against frost disasters, and in general for refined orchard production management.

Estimating Maize Yield in the Black Soil Region of Northeast China Using Land Surface Data Assimilation: Integrating a Crop Model and Remote Sensing.

Accurate yield estimation at the regional scale has always been a persistent challenge in the agricultural sector. With the vigorous emergence of remote sensing land surface observations in recent decades, data assimilation methodology has become an effective means to promote the accuracy and efficiency of yield estimation by integrating regional data and point-scale crop models. This paper focuses on the black soil area of Northeast China, a national strategic grain production base, applying the AquaCrop crop growth model to simulate the fractional vegetation cover (FVC) and maize yield from 2000 to 2020 and then forming a reliable FVC optimization dataset based on an ensemble Kalman filter (EnKF) assimilation algorithm with remote sensing products. Using the random forest model, the regression relationship between FVC and yield was established from the long-term time series data, which is crucial to achieve better yield estimation through the optimized FVC. The major findings include the following: (1) The R2 of the assimilated FVC and maize yield can reach 0.557. (2) When compared with the local statistical yield, our method reduced the mean absolute error (MAE) from 1.164 ton/ha (based on GLASS FVC products) to 1.004 ton/ha (based on the calibrated AquaCrop model) and then to 0.888 ton/ha (the result after assimilation). The above results show that we have proposed a yield estimation method to provide accurate yield estimations by combining data assimilation and machine learning. This study provided deep insights into understanding the variations in FVC and revealed the spatially explicit yield prediction ability from the time series land surface parameters, which has significant potential for optimizing water and soil resource management.

How does precipitation data influence the land surface data assimilation for drought monitoring?

Droughts are among the costliest natural hazards that occur annually worldwide. Their socioeconomic impacts are significant and widespread, affecting the sustainable development of human societies. This study investigates the influence of different forcing precipitation data in driving Land Surface Models (LSMs) and characterizing drought conditions. Here, we utilize our recently developed LSM data assimilation system for probabilistically monitoring drought over the Contiguous United States (CONUS). The Noah-MP LSM model is forced with two widely used precipitation data including IMERG (Integrated Multi-satellitE Retrievals for GPM) and NLDAS (North American Land Data Assimilation System). Soil moisture and evapotranspiration are known to have a strong relationship in the land-atmospheric interaction processes. Unlike other studies that attempted the individual assimilation of these variables, here we propose a multivariate data assimilation framework. Therefore, in both modeling scenarios, the data assimilation approach is used to integrate remotely sensed MODIS (Moderate Resolution Imaging Spectroradiometer) evapotranspiration and SMAP (Soil Moisture Active Passive) soil moisture observations into the Noah-MP LSM. The results of this study indicate that the source of precipitation data has a significant impact on the performance of LSM data assimilation system for drought monitoring. The findings revealed that NLDAS and IMERG precipitation can result in a significant difference in identifying drought severity depending on the region and time of the year. Furthermore, our analysis indicates that regardless of the precipitation forcing data product used in the land surface data assimilation system, our modeling framework can effectively detect the drought impacts on crop yield. Additionally, we calculated the drought probability based on the ensemble of soil moisture percentiles and found that there exist temporal and spatial discrepancies in drought probability maps generated from the NLDAS and IMERG precipitation forcings.

Improving the prediction of monsoon depressions by assimilating ASCAT soil moisture in NCUM-R modeling system

The present study investigates the impact of assimilation of satellite derived soil wetness data from Advanced Scatterometer (ASCAT) on simulation of inland moving monsoon depressions (MDs) over Indian region using high resolution NCUM-R forecasting system. The simplified Extended Kalman Filter (sEKF) based Land Surface Data Assimilation (LSDA) system is used to create land surface soil moisture initial conditions for NCUM-R by assimilating ASCAT soil moisture (SM) along with screen level observations. Two numerical experiments, namely CTL (incorporating only screen level observations in LSDA) and ASCAT (assimilating both ASCAT SM and screen level observations in LSDA) are carried out to update the model soil moisture initial conditions. Identical atmospheric initial condition prepared considering the 4DVAR data assimilation technique is used in both experiments.The results of this study show that the assimilation of remotely sensed SM has a beneficial impact on MD simulation, its movement, structure and associated precipitation. The analysis increment of fluxes from ASCAT experiment clearly shows properly coupled positive feedback between land and atmosphere. SM in the vicinity of MD is considerably improved in the ASCAT analysis as compared to the CTL analysis and is well correlated with the observed rainfall over the regions. The root mean square error of top layer SM is decreased by 5 to 15% in ASCAT analyses compared to CTL analyses. The land-atmosphere coupling feedback is investigated by using Dirmeyer's indices, it is noticed that the precipitation modulated by the latent heat flux is higher in ASCAT simulation than in the CTL simulation. This improvement is due to the proper representation of surface conditions by assimilating the satellite SM observations into LSDA. The contribution of individual terms of vorticity equation and total vorticity tendency is considerably well represented in ASCAT experiment as compared to the CTL experiment during the various phases of MD. Overall, this study indicates that the use of ASCAT SM in the preparation of soil moisture initial conditions helps to improve the forecast of MDs by the regional NWP system.

Temporal and Spatial Variation Characteristics of Soil Moisture in Spring in the Arid Regions of Northwest China in the Past 60s years

Soil moisture (SM) is one of the important parameters in the process of land-atmosphere interactions. The spatial-temporal distribution of SM plays a significant role in weather and climate research. In this study, based on the monthly SM datasets from GLDAS (Global Land Surface Data Assimilation System), the temporal and spatial changes of shallow SM are discussed, and the applicability of five domestic models from Coupled Model Intercomparison Project 6 (CMIP6) is also evaluated in the Northwest China. The results show that: 1) The shallow SM (0 - 10 cm) in spring in Northwest China was generally low during 1948-2015, and the low value areas were mainly located in the Tarim Basin in Xinjiang and the Gobi and desert areas in western Inner Mongolia. In most parts of Northwest China, SM had an increasing trend in spring, this implies it became wetter in recent 60 years; 2) There are larger difference between the five models for simulating SM. Except for BCC-ESM1, all the four models (including BCC-CSM2-MR, FGOALS-f3-L, FGOALS-g3, TaiESM1) can basically simulate the spatial distribution and trend of SM in spring, all the spatial correlation coefficients between the four models and GLDAS data pass the 99% significant level; 3) After multi-ensemble mean, the simulation performance can be obviously improved, the spatial correlation can reach about 0.55 in spring, the spatial trend is much closer to the GLDAS.

Reconstruction of Hourly All-Weather Land Surface Temperature by Integrating Reanalysis Data and Thermal Infrared Data From Geostationary Satellites (RTG)

Thermal infrared (TIR) land surface temperature (LST) products derived from geostationary satellites have a high temporal resolution in a diurnal cycle, but they have many missing values under cloudy-sky conditions. Therefore, it is pressing to obtain all-weather LST (AW LST) with a high temporal resolution by filling the gap of TIR LST. In this study, a method integrating reanalysis data and TIR data from geostationary satellites (RTG) was proposed for reconstructing hourly AW LST. Then, taking the Tibetan Plateau, which is a focus of climate change as a case, RTG was applied to the Chinese Fengyun-4A (FY-4A) TIR LST and China Land Surface Data Assimilation System (CLDAS) data. Validation based on the <italic xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">in-situ</i> LST shows that the accuracy of the AW LST is better than the FY-4A LST and CLDAS LST under clear-sky, cloudy-sky, and all-weather conditions. The mean RMSEs are 3.02 K for clear-sky conditions, 3.94 K for cloudy-sky conditions, and 3.57 K for all-weather conditions. Uncertainty and coarse resolution of the original FY-4A and CLDAS data affect the accuracy of the obtained AW LST. The results of the LST time series comparison also show that the reconstructed AW LST is consistent with <italic xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">in-situ</i> LST. The reconstructed AW LST also has good image quality and provides reliable spatial patterns. RTG is practical in obtaining high temporal resolution AW LST from the Chinese FY-4A to satisfy related applications. It can also be extended to other geostationary satellites and reanalysis datasets.

Toward the Estimation of All-Weather Daytime Downward Longwave Radiation over the Tibetan Plateau

Downward longwave radiation (DLR) is a critical parameter for radiation balance, energy budget, and water cycle studies at regional and global scales. Accurate estimation of the all-weather DLR with a high temporal resolution is important for the estimation of the surface net radiation and evapotranspiration. However, most DLR products involve instantaneous DLR estimates based on polar orbiting satellite data under clear-sky conditions. To obtain an in-depth understanding of the performances of different models in the estimation of DLR over the Tibetan Plateau, which is a focus area of climate change study, this study tests eight methods for clear-sky conditions and six methods for cloudy conditions based on ground-measured data. It is found that the Dilley and O’Brien model and the Lhomme model are most suitable for clear-sky conditions and cloudy conditions, respectively. For the Dilley and O’Brien model, the average root mean square error (RMSE) of DLR under clear-sky conditions is approximately 22.5 W/m2 for nine ground sites; for the Lhomme model, the average RMSE is approximately 23.2 W/m2. Based on the estimated cloud fraction and meteorological data provided by the China Land Surface Data Assimilation System (CLDAS), hourly all-weather daytime DLR with a 0.0625° resolution over the Tibetan Plateau is estimated. Results demonstrate that the average RMSE of the estimated hourly all-weather DLR is approximately 26.4 W/m2. With the combined all-weather DLR model, the hourly all-weather daytime DLR dataset with a 0.0625° resolution from 2008 to 2016 over the Tibetan Plateau is generated. This dataset can contribute to studies associated with the radiation balance and energy budget, water cycle, and climate change over the Tibetan Plateau.

Detecting and Attributing Evapotranspiration Deviations Using Dynamical Downscaling and Convection-Permitting Modeling over the Tibetan Plateau

Terrestrial evapotranspiration (ET) over the Tibetan Plateau (TP) exerts considerable impacts on the local climate and the water cycle. However, the high-altitude, mountainous areas over the TP pose a challenge for field observations. To finely capture its ET characteristics, we employed dynamical downscaling modeling (DDM) with a 28 km resolution and convection-permitting modeling (CPM) with a 4 km resolution in a normal climatology year, 2014. The benchmark data were the surface energy balance–based global land ET dataset (EB). Other compared data included the Global Land-Surface Data Assimilation System (GLDAS) and two reanalysis datasets: ERA-Interim and ERA5. Results showed that EB exhibits a gradient from the southeastern to northwestern TP, which is in line with the precipitation pattern. GLDAS generally reproduces the annual mean magnitude and pattern but poorly represents the seasonal variations. DDM and CPM perform well in the monsoon season but underestimate ET in the non-monsoon season. The two reanalysis datasets greatly overestimate the ET in the monsoon season, but ERA-Interim performs well in the non-monsoon season. All five datasets underestimate the ET over tundra and snow/ice areas, both in the annual and seasonal means. ET deviations are dominated by precipitation deviations in the monsoon season and by surface net radiation deviations in the non-monsoon season.

Snow Depth Fusion Based on Machine Learning Methods for the Northern Hemisphere

In this study, a machine learning algorithm was introduced to fuse gridded snow depth datasets. The input variables of the machine learning method included geolocation (latitude and longitude), topographic data (elevation), gridded snow depth datasets and in situ observations. A total of 29,565 in situ observations were used to train and optimize the machine learning algorithm. A total of five gridded snow depth datasets—Advanced Microwave Scanning Radiometer for the Earth Observing System (AMSR-E) snow depth, Global Snow Monitoring for Climate Research (GlobSnow) snow depth, Long time series of daily snow depth over the Northern Hemisphere (NHSD) snow depth, ERA-Interim snow depth and Modern-Era Retrospective Analysis for Research and Applications, version 2 (MERRA-2) snow depth—were used as input variables. The first three snow depth datasets are retrieved from passive microwave brightness temperature or assimilation with in situ observations, while the last two are snow depth datasets obtained from meteorological reanalysis data with a land surface model and data assimilation system. Then, three machine learning methods, i.e., Artificial Neural Networks (ANN), Support Vector Regression (SVR), and Random Forest Regression (RFR), were used to produce a fused snow depth dataset from 2002 to 2004. The RFR model performed best and was thus used to produce a new snow depth product from the fusion of the five snow depth datasets and auxiliary data over the Northern Hemisphere from 2002 to 2011. The fused snow-depth product was verified at five well-known snow observation sites. The R2 of Sodankylä, Old Aspen, and Reynolds Mountains East were 0.88, 0.69, and 0.63, respectively. At the Swamp Angel Study Plot and Weissfluhjoch observation sites, which have an average snow depth exceeding 200 cm, the fused snow depth did not perform well. The spatial patterns of the average snow depth were analyzed seasonally, and the average snow depths of autumn, winter, and spring were 5.7, 25.8, and 21.5 cm, respectively. In the future, random forest regression will be used to produce a long time series of a fused snow depth dataset over the Northern Hemisphere or other specific regions.

C‐LAEF: Convection‐permitting Limited‐Area Ensemble Forecasting system

AbstractC‐LAEF (Convection‐permitting Limited‐Area Ensemble Forecasting) has been developed at the Austrian national weather service ZAMG (Zentralanstalt für Meteorologie und Geodynamik) and has been running operationally at the European Centre for Medium Range Weather Forecasts (ECMWF) supercomputer since November 2019. It includes (a) an ensemble 3D variational blending technique to deal with atmospheric initial uncertainties, (b) an ensemble of land surface data assimilation to account for uncertainties in the initial land surface conditions, (c) a hybrid stochastic physics perturbation scheme to treat model uncertainties in the different physics parametrization schemes, and (d) a coupling with the global ensemble system IFS‐ENS (Integrated Forecasting System‐ENSemble) to consider uncertainties in the lateral boundary conditions. C‐LAEF has a horizontal resolution of 2.5 km and consists of 16 perturbed members plus one unperturbed control run. It runs four times a day and provides probabilistic forecasts up to 60 hr on a domain covering the whole Alpine area. This article describes the C‐LAEF system in detail and evaluates the relative contributions of the different perturbation techniques. The ensemble variational blending technique and the ensemble surface data assimilation provide additional spread in the first forecast hours, while the hybrid stochastically perturbed parametrization scheme improves the performance of C‐LAEF during the whole forecast range. The performance of C‐LAEF is evaluated extensively over one summer and one winter period and compared with its mesoscale counterpart ALADIN‐LAEF (Aire Limitee Adaptation dynamique Developpement InterNational – Limited‐Area Ensemble Forecasting) and the IFS‐ENS. State‐of‐the‐art probabilistic measures indicate that C‐LAEF is able to outperform ALADIN‐LAEF for all considered upper‐air variables. At the surface, C‐LAEF outperforms ALADIN‐LAEF and IFS‐ENS according to most conventional measures, particularly for wind and precipitation. C‐LAEF benefits from the higher resolution and the explicit treatment of deep convection and can provide more accurate probabilistic information for weather warnings.

How Does Precipitation Data Influence the Land Surface Data Assimilation Drought Monitoring?

How Does Precipitation Data Influence the Land Surface Data Assimilation Drought Monitoring?

The Met Office Operational Soil Moisture Analysis System

In this study, the current Met Office operational land surface data assimilation system used to produce soil moisture analyses is presented. The main aim of including Land Surface Data Assimilation (LSDA) in both the global and regional systems is to improve forecasts of surface air temperature and humidity. Results from trials assimilating pseudo-observations of 1.5 m air temperature and specific humidity and satellite-derived soil wetness (ASCAT) observations are analysed. The pre-processing of all the observations is described, including the definition and construction of the pseudo-observations. The benefits of using both observations together to produce improved forecasts of surface air temperature and humidity are outlined both in the winter and summer seasons. The benefits of using active LSDA are quantified by the root mean squared error, which is computed using both surface observations and European Centre for Medium-Range Weather Forecasts (ECMWF) analyses as truth. For the global model trials, results are presented separately for the Northern (NH) and Southern (SH) hemispheres. When compared against ground-truth, LSDA in winter NH appears neutral, but in the SH it is the assimilation of ASCAT that contributes to approximately a 2% improvement in temperatures at lead times beyond 48 h. In NH summer, the ASCAT soil wetness observations degrade the forecasts against observations by about 1%, but including the screen level pseudo-observations provides a compensating benefit. In contrast, in the SH, the positive effect comes from including the ASCAT soil wetness observations, and when both observations types are assimilated there is a compensating effect. Finally, we demonstrate substantial improvements to hydrological prediction when using land surface data assimilation in the regional model. Using the Nash-Sutcliffe Efficiency (NSE) metric as an aggregated measure of river flow simulation skill relative to observations, we find that NSE was improved at 106 of 143 UK river gauge locations considered after LSDA was introduced. The number of gauge comparisons where NSE exceeded 0.5 is also increased from 17 to 28 with LSDA.

Study on the Variation Trend of Potential Evapotranspiration in the Three-River Headwaters Region in China Over the Past 20 years

To study the variation trend of potential evapotranspiration (PET) in the Three-River Headwaters (TRH) region of the Qinghai-Tibet Plateau in China, we use 2-m temperature and surface pressure observation data from 14 meteorological weather stations in the TRH region, and the surface PET is calculated by the Penman-Monteith formula. The global land surface data assimilation system (GLDAS) from 2000 to 2018 were used to compare and verify the accuracy and applicability of the calculated PET in the THR region. The results show that in the past 20 years, the PET of 14 weather stations in the TRH region has shown an increasing trend, with annual averaged growth rate of 1.4±1.2 mm·a-1, and the spatial distribution of the annual variation rate of PET has obvious difference. PET is higher in the eastern area of TRH region, and lower in the western area. The drought in this area increased from southeast to northwest, which was consistent with the spatial distribution of precipitation. The aridity index K has fluctuated and increased begore 2015, but there was a sudden change in 2018, and the aridity index K in the TRH began to decline after 2018, and the climate changed from dry to wet.

Improving Spatial Patterns Prior to Land Surface Data Assimilation via Model Calibration Using SMAP Surface Soil Moisture Data

AbstractPrior to their application in a land surface data assimilation (DA) system, remotely sensed (RS) soil moisture (SM) products are typically rescaled to remove systematic differences with respect to comparable SM estimates obtained from a land surface model (LSM). This preprocessing of RS SM data—commonly referred to as bias correction—implicitly discards spatial information present in the RS SM retrievals. Here, based on dense SM network measurements obtained within the Huai River Basin of China, we demonstrate that L‐band Soil Moisture Active Passive (SMAP) L3 surface SM retrievals provide a better representation of spatial variability in time‐averaged SM fields than a LSM—suggesting that, as typically applied, bias correction is neglecting important spatial information present in RS SM products. To maximally use RS information, we propose an alternative bias correction approach for SM DA that provides the benefits of traditional rescaling while also improving modeled SM spatial patterns. Specifically, the regional mean of SMAP SM is first corrected using modeled SM of the Variable Infiltration Capacity (VIC) model. Then, two key VIC soil parameters (EXPT and BULKD) are calibrated to match the relative subregional spatial variability captured by the SMAP SM product. Results show that our model calibration method successfully improves VIC SM spatial patterns and pixel‐wise time series. Furthermore, these SM improvements translate into enhanced VIC streamflow estimates. Overall, results suggest that the proposed bias correction framework can improve current land surface DA systems by maximally utilizing spatial information contained in RS SM products.

A 3\u2009km spatially and temporally consistent European daily soil moisture reanalysis from 2000 to 2015

High-resolution soil moisture (SM) information is essential to many regional applications in hydrological and climate sciences. Many global estimates of surface SM are provided by satellite sensors, but at coarse spatial resolutions (lower than 25 km), which are not suitable for regional hydrologic and agriculture applications. Here we present a 16 years (2000–2015) high-resolution spatially and temporally consistent surface soil moisture reanalysis (ESSMRA) dataset (3 km, daily) over Europe from a land surface data assimilation system. Coarse-resolution satellite derived soil moisture data were assimilated into the community land model (CLM3.5) using an ensemble Kalman filter scheme, producing a 3 km daily soil moisture reanalysis dataset. Validation against 112 in-situ soil moisture observations over Europe shows that ESSMRA captures the daily, inter-annual, intra-seasonal patterns well with RMSE varying from 0.04 to 0.06 m3m−3 and correlation values above 0.5 over 70% of stations. The dataset presented here provides long-term daily surface soil moisture at a high spatiotemporal resolution and will be beneficial for many hydrological applications over regional and continental scales.