Assessment and enhancement of MERRA land surface hydrology estimates

The Modern-Era Retrospective Analysis for Research and Applications (MERRA) is a state-of-the-art reanalysis that provides, in addition to atmospheric fields, global estimates of soil moisture, latent heat flux, snow, and runoff for 1979 present. This study introduces a supplemental and improved set of land surface hydrological fields ("MERRA-Land") generated by rerunning a revised version of the land component of the MERRA system. Specifically, the MERRA-Land estimates benefit from corrections to the precipitation forcing with the Global Precipitation Climatology Project pentad product (version 2.1) and from revised parameter values in the rainfall interception model, changes that effectively correct for known limitations in the MERRA surface meteorological forcings. The skill (defined as the correlation coefficient of the anomaly time series) in land surface hydrological fields from MERRA and MERRA-Land is assessed here against observations and compared to the skill of the state-of-the-art ECMWF Re-Analysis-Interim (ERA-I). MERRA-Land and ERA-I root zone soil moisture skills (against in situ observations at 85 U.S. stations) are comparable and significantly greater than that of MERRA. Throughout the Northern Hemisphere, MERRA and MERRA-Land agree reasonably well with in situ snow depth measurements (from 583 stations) and with snow water equivalent from an independent analysis. Runoff skill (against naturalized stream flow observations from 18 U.S. basins) of MERRA and MERRA-Land is typically higher than that of ERA-I. With a few exceptions, the MERRA-Land data appear more accurate than the original MERRA estimates and are thus recommended for those interested in using MERRA output for land surface hydrological studies

Final Report of the DAUFIN project

DAUFIN = Data Assimulation within Unifying Framework for Improved river basiN modeling (EC 5th framework Project

Aqueduct Metadata Document, Aqueduct Global Maps 2.0

This document describes the specific characteristics of the indicator data and calculations for the Aqueduct Water Risk Atlas Global Maps. Complete guidelines and processes for data collection, calculations, and mapping techniques are described fully in the Aqueduct Water Risk Framework. The Aqueduct Water Risk Atlas makes use of a Water Risk Framework, that includes 12 global indicators grouped into three categories of risk and one overall score

Multi-source data assimilation for physically based hydrological modeling of an experimental hillslope

Data assimilation has recently been the focus of much attention for integrated surface–subsurface hydrological models, whereby joint assimilation of water table, soil moisture, and river discharge measurements with the ensemble Kalman filter (EnKF) has been extensively applied. Although the EnKF has been specifically developed to deal with nonlinear models, integrated hydrological models based on the Richards equation still represent a challenge, due to strong nonlinearities that may significantly affect the filter performance. Thus, more studies are needed to investigate the capabilities of the EnKF to correct the system state and identify parameters in cases where the unsaturated zone dynamics are dominant, as well as to quantify possible tradeoffs associated with assimilation of multi-source data. Here, the CATHY (CATchment HYdrology) model is applied to reproduce the hydrological dynamics observed in an experimental two-layered hillslope, equipped with tensiometers, water content reflectometer probes, and tipping bucket flow gages to monitor the hillslope response to a series of artificial rainfall events. Pressure head, soil moisture, and subsurface outflow are assimilated with the EnKF in a number of scenarios and the challenges and issues arising from the assimilation of multi-source data in this real-world test case are discussed. Our results demonstrate that the EnKF is able to effectively correct states and parameters even in a real application characterized by strong nonlinearities. However, multi-source data assimilation may lead to significant tradeoffs: the assimilation of additional variables can lead to degradation of model predictions for other variables that are otherwise well reproduced. Furthermore, we show that integrated observations such as outflow discharge cannot compensate for the lack of well-distributed data in heterogeneous hillslopes.</p

Assessing the utility of geospatial technologies to investigate environmental change within lake systems

Over 50% of the world's population live within 3. km of rivers and lakes highlighting the on-going importance of freshwater resources to human health and societal well-being. Whilst covering c. 3.5% of the Earth's non-glaciated land mass, trends in the environmental quality of the world's standing waters (natural lakes and reservoirs) are poorly understood, at least in comparison with rivers, and so evaluation of their current condition and sensitivity to change are global priorities. Here it is argued that a geospatial approach harnessing existing global datasets, along with new generation remote sensing products, offers the basis to characterise trajectories of change in lake properties e.g., water quality, physical structure, hydrological regime and ecological behaviour. This approach furthermore provides the evidence base to understand the relative importance of climatic forcing and/or changing catchment processes, e.g. land cover and soil moisture data, which coupled with climate data provide the basis to model regional water balance and runoff estimates over time. Using examples derived primarily from the Danube Basin but also other parts of the World, we demonstrate the power of the approach and its utility to assess the sensitivity of lake systems to environmental change, and hence better manage these key resources in the future

Correcting the radar rainfall forcing of a hydrological model with data assimilation: application to flood forecasting in the Lez Catchment in Southern France

The present study explores the application of a data assimilation (DA) procedure to correct the radar rain- fall inputs of an event-based, distributed, parsimonious hy- drological model. An extended Kalman filter algorithm was built on top of a rainfall-runoff model in order to assimilate discharge observations at the catchment outlet. This work fo- cuses primarily on the uncertainty in the rainfall data and considers this as the principal source of error in the sim- ulated discharges, neglecting simplifications in the hydro- logical model structure and poor knowledge of catchment physics. The study site is the 114 km2 Lez catchment near Montpellier, France. This catchment is subject to heavy oro- graphic rainfall and characterised by a karstic geology, lead- ing to flash flooding events. The hydrological model uses a derived version of the SCS method, combined with a Lag and Route transfer function. Because the radar rainfall in- put to the model depends on geographical features and cloud structures, it is particularly uncertain and results in signifi- cant errors in the simulated discharges. This study seeks to demonstrate that a simple DA algorithm is capable of ren- dering radar rainfall suitable for hydrological forecasting. To test this hypothesis, the DA analysis was applied to estimate a constant hyetograph correction to each of 19 flood events. The analysis was carried in two different modes: by assimi- lating observations at all available time steps, referred to here as reanalysis mode, and by using only observations up to 3 h before the flood peak to mimic an operational environment, referred to as pseudo-forecast mode. In reanalysis mode, the resulting correction of the radar rainfall data was then com- pared to the mean field bias (MFB), a corrective coefficient determined using rain gauge measurements. It was shown that the radar rainfall corrected using DA leads to improved discharge simulations and Nash-Sutcliffe efficiency criteria compared to the MFB correction. In pseudo-forecast mode, the reduction of the uncertainty in the rainfall data leads to a reduction of the error in the simulated discharge, but un- certainty from the model parameterisation diminishes data assimilation efficiency. While the DA algorithm used is this study is effective in correcting uncertain radar rainfall, model uncertainty remains an important challenge for flood fore- casting within the Lez catchment

MSWEP : 3-hourly 0.25° global gridded precipitation (1979-2015) by merging gauge, satellite, and reanalysis data

Current global precipitation (P) datasets do not take full advantage of the complementary nature of satellite and reanalysis data. Here, we present Multi-Source Weighted-Ensemble Precipitation (MSWEP) version 1.1, a global P dataset for the period 1979-2015 with a 3hourly temporal and 0.25 degrees ffi spatial resolution, specifically designed for hydrological modeling. The design philosophy of MSWEP was to optimally merge the highest quality P data sources available as a function of timescale and location. The long-term mean of MSWEP was based on the CHPclim dataset but replaced with more accurate regional datasets where available. A correction for gauge under-catch and orographic effects was introduced by inferring catchment-average P from streamflow (Q) observations at 13 762 stations across the globe. The temporal variability of MSWEP was determined by weighted averaging of P anomalies from seven datasets; two based solely on interpolation of gauge observations (CPC Unified and GPCC), three on satellite remote sensing (CMORPH, GSMaP-MVK, and TMPA 3B42RT), and two on atmospheric model reanalysis (ERA-Interim and JRA-55). For each grid cell, the weight assigned to the gauge-based estimates was calculated from the gauge network density, while the weights assigned to the satellite-and reanalysis-based estimates were calculated from their comparative performance at the surrounding gauges. The quality of MSWEP was compared against four state-of-the-art gauge-adjusted P datasets (WFDEI-CRU, GPCP-1DD, TMPA 3B42, and CPC Unified) using independent P data from 125 FLUXNET tower stations around the globe. MSWEP obtained the highest daily correlation coefficient (R) among the five P datasets for 60.0% of the stations and a median R of 0.67 vs. 0.44-0.59 for the other datasets. We further evaluated the performance of MSWEP using hydrological modeling for 9011 catchments (< 50 000 km(2)) across the globe. Specifically, we calibrated the simple conceptual hydrological model HBV (Hydrologiska Byrans Vattenbalansavdelning) against daily Q observations with P from each of the different datasets. For the 1058 sparsely gauged catchments, representative of 83.9% of the global land surface (excluding Antarctica), MSWEP obtained a median calibration NSE of 0.52 vs. 0.29-0.39 for the other P datasets. MSWEP is available via http://www.gloh2o.org

Models of everywhere revisited: a technological perspective

The concept ‘models of everywhere’ was first introduced in the mid 2000s as a means of reasoning about the environmental science of a place, changing the nature of the underlying modelling process, from one in which general model structures are used to one in which modelling becomes a learning process about specific places, in particular capturing the idiosyncrasies of that place. At one level, this is a straightforward concept, but at another it is a rich multi-dimensional conceptual framework involving the following key dimensions: models of everywhere, models of everything and models at all times, being constantly re-evaluated against the most current evidence. This is a compelling approach with the potential to deal with epistemic uncertainties and nonlinearities. However, the approach has, as yet, not been fully utilised or explored. This paper examines the concept of models of everywhere in the light of recent advances in technology. The paper argues that, when first proposed, technology was a limiting factor but now, with advances in areas such as Internet of Things, cloud computing and data analytics, many of the barriers have been alleviated. Consequently, it is timely to look again at the concept of models of everywhere in practical conditions as part of a trans-disciplinary effort to tackle the remaining research questions. The paper concludes by identifying the key elements of a research agenda that should underpin such experimentation and deployment