Framed within the Copernicus Climate Change Service (C3S) of the European Commission, the European Centre for Medium-Range Weather Forecasts (ECMWF) is producing an enhanced global dataset for the land component of the fifth generation of European ReAnalysis (ERA5), hereafter referred to as ERA5-Land. Once completed, the period covered will span from 1950 to the present, with continuous updates to support land monitoring applications. ERA5-Land describes the evolution of the water and energy cycles over land in a consistent manner over the production period, which, among others, could be used to analyse trends and anomalies. This is achieved through global high-resolution numerical integrations of the ECMWF land surface model driven by the downscaled meteorological forcing from the ERA5 climate reanalysis, including an elevation correction for the thermodynamic near-surface state. ERA5-Land shares with ERA5 most of the parameterizations that guarantees the use of the state-of-the-art land surface modelling applied to numerical weather prediction (NWP) models. A main advantage of ERA5-Land compared to ERA5 and the older ERA-Interim is the horizontal resolution, which is enhanced globally to 9 km compared to 31 km (ERA5) or 80 km (ERA-Interim), whereas the temporal resolution is hourly as in ERA5. Evaluation against independent in situ observations and global model or satellite-based reference datasets shows the added value of ERA5-Land in the description of the hydrological cycle, in particular with enhanced soil moisture and lake description, and an overall better agreement of river discharge estimations with available observations. However, ERA5-Land snow depth fields present a mixed performance when compared to those of ERA5, depending on geographical location and altitude. The description of the energy cycle shows comparable results with ERA5. Nevertheless, ERA5-Land reduces the global averaged root mean square error of the skin temperature, taking as reference MODIS data, mainly due to the contribution of coastal points where spatial resolution is important. Since January 2020, the ERA5-Land period available has extended from January 1981 to the near present, with a 2- to 3-month delay with respect to real time. The segment prior to 1981 is in production, aiming for a release of the whole dataset in summer/autumn 2021. The high spatial and temporal resolution of ERA5-Land, its extended period, and the consistency of the fields produced makes it a valuable dataset to support hydrological studies, to initialize NWP and climate models, and to support diverse applications dealing with water resource, land, and environmental management. The full ERA5-Land hourly (Muñoz-Sabater, 2019a) and monthly (Muñoz-Sabater, 2019b) averaged datasets presented in this paper are available through the C3S Climate Data Store at https://doi.org/10.24381/cds.e2161bac and https://doi.org/10.24381/cds.68d2bb30, respectively.

The Global Land Evaporation Amsterdam Model (GLEAM) is a set of algorithms dedicated to the estimation of terrestrial evaporation and root-zone soil moisture from satellite data. Ever since its development in 2011, the model has been regularly revised, aiming at the optimal incorporation of new satellite-observed geophysical variables, and improving the representation of physical processes. In this study, the next version of this model (v3) is presented. Key changes relative to the previous version include (1) a revised formulation of the evaporative stress, (2) an optimized drainage algorithm, and (3) a new soil moisture data assimilation system. GLEAM v3 is used to produce three new data sets of terrestrial evaporation and root-zone soil moisture, including a 36-year data set spanning 1980–2015, referred to as v3a (based on satellite-observed soil moisture, vegetation optical depth and snow-water equivalent, reanalysis air temperature and radiation, and a multi-source precipitation product), and two satellite-based data sets. The latter share most of their forcing, except for the vegetation optical depth and soil moisture, which are based on observations from different passive and active C- and L-band microwave sensors (European Space Agency Climate Change Initiative, ESA CCI) for the v3b data set (spanning 2003–2015) and observations from the Soil Moisture and Ocean Salinity (SMOS) satellite in the v3c data set (spanning 2011–2015). Here, these three data sets are described in detail, compared against analogous data sets generated using the previous version of GLEAM (v2), and validated against measurements from 91 eddy-covariance towers and 2325 soil moisture sensors across a broad range of ecosystems. Results indicate that the quality of the v3 soil moisture is consistently better than the one from v2: average correlations against in situ surface soil moisture measurements increase from 0.61 to 0.64 in the case of the v3a data set and the representation of soil moisture in the second layer improves as well, with correlations increasing from 0.47 to 0.53. Similar improvements are observed for the v3b and c data sets. Despite regional differences, the quality of the evaporation fluxes remains overall similar to the one obtained using the previous version of GLEAM, with average correlations against eddy-covariance measurements ranging between 0.78 and 0.81 for the different data sets. These global data sets of terrestrial evaporation and root-zone soil moisture are now openly available at www.GLEAM.eu and may be used for large-scale hydrological applications, climate studies, or research on land–atmosphere feedbacks.

This paper outlines a new strategy to derive evaporation from satellite observations. The approach uses a variety of satellite-sensor products to estimate daily evaporation at a global scale and 0.25 degree spatial resolution. Central to this methodology is the use of the Priestley and Taylor (PT) evaporation model. The minimalistic PT equation combines a small number of inputs, the majority of which can be detected from space. This reduces the number of variables that need to be modelled. Key distinguishing features of the approach are the use of microwave-derived soil moisture, land surface temperature and vegetation density, as well as the detailed estimation of rainfall interception loss. The modelled evaporation is validated against one year of eddy covariance measurements from 43 stations. The estimated annual totals correlate well with the stations' annual cumulative evaporation (R=0.80, N=43) and present a low average bias (−5%). The validation of the daily time series at each individual station shows good model performance in all vegetation types and climate conditions with an average correlation coefficient of R=0.83, still lower than the R=0.90 found in the validation of the monthly time series. The first global map of annual evaporation developed through this methodology is also presented.

We present Multi-Source Weighted-Ensemble Precipitation, version 2 (MSWEP V2), a gridded precipitation P dataset spanning 1979–2017. MSWEP V2 is unique in several aspects: i) full global coverage (all land and oceans); ii) high spatial (0.1°) and temporal (3 hourly) resolution; iii) optimal merging of P estimates based on gauges [WorldClim, Global Historical Climatology Network-Daily (GHCN-D), Global Summary of the Day (GSOD), Global Precipitation Climatology Centre (GPCC), and others], satellites [Climate Prediction Center morphing technique (CMORPH), Gridded Satellite (GridSat), Global Satellite Mapping of Precipitation (GSMaP), and Tropical Rainfall Measuring Mission (TRMM) Multisatellite Precipitation Analysis (TMPA) 3B42RT)], and reanalyses [European Centre for Medium-Range Weather Forecasts (ECMWF) interim reanalysis (ERA-Interim) and Japanese 55-year Reanalysis (JRA-55)]; iv) distributional bias corrections, mainly to improve the P frequency; v) correction of systematic terrestrial P biases using river discharge Q observations from 13,762 stations across the globe; vi) incorporation of daily observations from 76,747 gauges worldwide; and vii) correction for regional differences in gauge reporting times. MSWEP V2 compares substantially better with Stage IV gauge–radar P data than other state-of-the-art P datasets for the United States, demonstrating the effectiveness of the MSWEP V2 methodology. Global comparisons suggest that MSWEP V2 exhibits more realistic spatial patterns in mean, magnitude, and frequency. Long-term mean P estimates for the global, land, and ocean domains based on MSWEP V2 are 955, 781, and 1,025 mm yr⁻¹, respectively. Other P datasets consistently underestimate P amounts in mountainous regions. Using MSWEP V2, P was estimated to occur 15.5%, 12.3%, and 16.9% of the time on average for the global, land, and ocean domains, respectively. MSWEP V2 provides unique opportunities to explore spatiotemporal variations in P, improve our understanding of hydrological processes and their parameterization, and enhance hydrological model performance.