ALBERT

Description: Hourly NOx concentrations were simulated for the city of Zurich, Switzerland, at 10 m resolution for the years 2013–2014. The simulations were generated with the nested mesoscale meteorology and microscale dispersion model system GRAMM/GRAL (versions v15.12/v14.8) by applying a catalogue-based approach. This approach was specifically designed to enable long-term city-wide building-resolving simulations with affordable computation costs. It relies on a discrete set of possible weather situations and corresponding steady-state flow and dispersion patterns that are pre-computed and then matched hourly with actual meteorological observations. The modelling system was comprehensively evaluated using eight sites continuously monitoring NOx concentrations and 65 passive samplers measuring NO2 concentrations on a 2-weekly basis all over the city. The system was demonstrated to fulfil the European Commission standards for air pollution modelling at nearly all sites. The average spatial distribution was very well represented, despite a general tendency to overestimating the observed concentrations, possibly due to a crude representation of traffic-induced turbulence. The temporal variability of concentrations explained by varying emissions and weather situations was accurately reproduced on different time scales. The seasonal cycle of concentrations, mostly driven by stronger vertical dispersion in summer than in winter, was very well captured in the two year simulation period. Short-term events, such as episodes of particularly high and low concentrations, were detected in most cases by the system, although some unrealistic pollution peaks were occasionally generated, pointing at some limitations of the steady-state approximation. The different patterns of the diurnal cycle of concentrations observed in the city were generally well captured as well. The evaluation confirmed the adequacy of the catalogue-based approach in the context of city scale air pollution modelling. The ability to reproduce not only the spatial gradients but also the hourly temporal variability over multiple years makes the model system particularly suitable for investigating individualized air pollution exposure in the city.

Description: Most of the current large scale hydrological models do not contain a physically-based groundwater flow component. The main difficulties in large-scale groundwater modeling include the efficient representation of unsaturated zone flow, the characterization of dynamic groundwater-surface water interaction and the numerical stability while preserving complex physical processes and high resolution. To address these problems, we propose a highly-scalable coupled hydrologic and groundwater model (mHM#OGS) based on the integration of two open-source modeling codes: the mesoscale hydrologic Model (mHM) and the finite element simulator OpenGeoSys (OGS). mHM#OGS is coupled using a boundary condition-based coupling scheme that dynamically links the surface and subsurface parts. Nested time stepping allows smaller time steps for typically faster surface runoff routing in mHM and larger time steps for slower subsurface flow in OGS. mHM#OGS features the coupling interface which can transfer the groundwater recharge and river baseflow rate between mHM and OpenGeoSys. Verification of the coupled model was conducted using the time-series of observed streamflow and groundwater levels. Moreover, we force the transient model using groundwater recharge in two scenarios: (1) spatially variable recharge based on the mHM simulations, and (2) spatially homogeneous groundwater recharge. The modeling result in first scenario has a slightly higher correlation with groundwater head time-series, which further validates the plausibility of spatial groundwater recharge distribution calculated by mHM in the mesocale. The statistical analysis of model predictions shows a promising prediction ability of the model. The offline coupling method implemented here can reproduce reasonable groundwater head time series while keep a desired level of detail in the subsurface model structure with little surplus in computational cost. Our exemplary calculations show that the coupled model mHM#OGS can be a valuable tool to assess the effects of variability in land surface heterogeneity, meteorological, topographical forces and geological zonation on the groundwater flow dynamics.

Description: Hourly NOx concentrations were simulated for the city of Zürich, Switzerland, at 10 m resolution for the years 2013–2014. The simulations were generated with the nested mesoscale meteorology and micro-scale dispersion model system GRAMM–GRAL (versions v15.12 and v14.8) by applying a catalogue-based approach. This approach was specifically designed to enable long-term city-wide building-resolving simulations with affordable computation costs. It relies on a discrete set of possible weather situations and corresponding steady-state flow and dispersion patterns that are pre-computed and then matched hourly with actual meteorological observations. The modelling system was comprehensively evaluated using eight sites continuously monitoring NOx concentrations and 65 passive samplers measuring NO2 concentrations on a 2-weekly basis all over the city. The system was demonstrated to fulfil the European Commission standards for air pollution modelling at nearly all sites. The average spatial distribution was very well represented, despite a general tendency to overestimate the observed concentrations, possibly due to a crude representation of traffic-induced turbulence and to underestimated dispersion in the vicinity of buildings. The temporal variability of concentrations explained by varying emissions and weather situations was accurately reproduced on different timescales. The seasonal cycle of concentrations, mostly driven by stronger vertical dispersion in summer than in winter, was very well captured in the 2-year simulation period. Short-term events, such as episodes of particularly high and low concentrations, were detected in most cases by the system, although some unrealistic pollution peaks were occasionally generated, pointing at some limitations of the steady-state approximation. The different patterns of the diurnal cycle of concentrations observed in the city were generally well captured as well. The evaluation confirmed the adequacy of the catalogue-based approach in the context of city-scale air pollution modelling. The ability to reproduce not only the spatial gradients but also the hourly temporal variability over multiple years makes the model system particularly suitable for investigating individualized air pollution exposure in the city.

Description: Soil properties show high heterogeneity at different spatial scales and their correct characterization remains a crucial challenge over large areas. The aim of the study is to quantify the impact of different types of uncertainties that arise from the unresolved soil spatial variability on simulated hydrological states and fluxes. Three perturbation methods are presented for the characterization of uncertainties in soil properties. The methods are applied on the soil map of the upper Neckar catchment (Germany), as an example. The uncertainties are propagated through the distributed mesoscale hydrological model (mHM) to assess the impact on the simulated states and fluxes. The model outputs are analysed by aggregating the results at different spatial and temporal scales. These results show that the impact of the different uncertainties introduced in the original soil map is equivalent when the simulated model outputs are analysed at the model grid resolution (i.e. 500 m). However, several differences are identified by aggregating states and fluxes at different spatial scales (by subcatchments of different sizes or coarsening the grid resolution). Streamflow is only sensitive to the perturbation of long spatial structures while distributed states and fluxes (e.g. soil moisture and groundwater recharge) are only sensitive to the local noise introduced to the original soil properties. A clear identification of the temporal and spatial scale for which finer-resolution soil information is (or is not) relevant is unlikely to be universal. However, the comparison of the impacts on the different hydrological components can be used to prioritize the model improvements in specific applications, either by collecting new measurements or by calibration and data assimilation approaches. In conclusion, the study underlines the importance of a correct characterization of uncertainty in soil properties. With that, soil maps with additional information regarding the unresolved soil spatial variability would provide strong support to hydrological modelling applications.

Description: Travel-time distributions are a comprehensive tool for the characterization of hydrological system dynamics. Unlike the streamflow hydrograph, they describe the movement and storage of water within and throughout the hydrological system. Until recently, studies using such travel-time distributions have generally either been applied to lumped models or to real-world catchments using available time series, e.g., stable isotopes. Whereas the former are limited in their realism and lack information on the spatial arrangements of the relevant quantities, the latter are limited in their use of available data sets. In our study, we employ the spatially distributed mesoscale Hydrological Model (mHM) and apply it to a catchment in central Germany. Being able to draw on multiple large data sets for calibration and verification, we generate a large array of spatially distributed states and fluxes. These hydrological outputs are then used to compute the travel-time distributions for every grid cell in the modeling domain. A statistical analysis indicates the general soundness of the upscaling scheme employed in mHM and reveals precipitation, saturated soil moisture and potential evapotranspiration as important predictors for explaining the spatial heterogeneity of mean travel times. In addition, we demonstrate and discuss the high information content of mean travel times for characterization of internal hydrological processes.

Description: Long-term, high-resolution data about hydrologic fluxes and states are needed for many hydrological applications. Because continuous large-scale observations of such variables are not feasible, hydrologic or land surface models are applied to derive them. This study aims to analyze and provide a consistent high-resolution dataset of land surface variables over Germany, accounting for uncertainties caused by equifinal model parameters. The mesoscale Hydrological Model (mHM) is employed to derive an ensemble (100 members) of evapotranspiration, groundwater recharge, soil moisture, and runoff generated at high spatial and temporal resolutions (4 km and daily, respectively) for the period 1951–2010. The model is cross-evaluated against the observed daily streamflow in 222 basins, which are not used for model calibration. The mean (standard deviation) of the ensemble median Nash–Sutcliffe efficiency estimated for these basins is 0.68 (0.09) for daily streamflow simulations. The modeled evapotranspiration and soil moisture reasonably represent the observations from eddy covariance stations. Our analysis indicates the lowest parametric uncertainty for evapotranspiration, and the largest is observed for groundwater recharge. The uncertainty of the hydrologic variables varies over the course of a year, with the exception of evapotranspiration, which remains almost constant. This study emphasizes the role of accounting for the parametric uncertainty in model-derived hydrological datasets.

Description: There is growing evidence that climate change will alter water availability in Europe. Here, we investigate how hydrological low flows are affected under different levels of future global warming (i.e., 1.5, 2 and 3 K). The analysis is based on a multi-model ensemble of 45 hydrological simulations based on three RCPs (rcp2p6, rcp6p0, rcp8p5), five CMIP5 GCMs (GFDL-ESM2M, HadGEM2-ES, IPSL-CM5A-LR, MIROC-ESM-CHEM, NorESM1-M) and three state-of-the-art hydrological models (HMs: mHM, Noah-MP, and PCR-GLOBWB). High resolution model results are available at the unprecedented spatial resolution of 5 km across the pan-European domain at daily temporal resolution. Low river flow is described as the percentile of daily streamflow that is exceeded 90 % of the time. It is determined separately for each GCM/HM combinations and the warming scenarios. The results show that the change signal amplifies with increasing warming levels. Low flows decrease in the Mediterranean while they increase in the Alpine and Northern regions. In the Mediterranean, the level of warming amplifies the signal from −12 % under 1.5 K to −35 % under 3 K global warming largely due to the projected decreases in annual precipitation. In contrast, the signal is amplified from +22 % (1.5 K) to +45 % (3 K) in the Alpine region because of the reduced snow melt contribution. The changes in low flows are significant for regions with relatively large change signals and under higher levels of warming. Nevertheless, it is not possible to distinguish climate induced differences in low flows between 1.5 and 2 K warming because of the large variability inherent in the multi-model ensemble. The contribution by the GCMs to the uncertainty in the model results is generally higher than the one by the HMs. However, the uncertainty due to HMs cannot be neglected. In the Alpine and Northern region as well as the Mediterranean, the uncertainty contribution by the HMs is partly higher than those by the GCMs due to different representations of processes such as snow, soil moisture and evapotranspiration.

Description: Land surface and hydrologic models (LSMs/HMs) are used at diverse spatial resolutions ranging from catchment-scale (1–10 km) to global-scale (over 50 km) applications. Applying the same model structure at different spatial scales requires that the model estimates similar fluxes independent of the chosen resolution, i.e., fulfills a flux-matching condition across scales. An analysis of state-of-the-art LSMs and HMs reveals that most do not have consistent hydrologic parameter fields. Multiple experiments with the mHM, Noah-MP, PCR-GLOBWB, and WaterGAP models demonstrate the pitfalls of deficient parameterization practices currently used in most operational models, which are insufficient to satisfy the flux-matching condition. These examples demonstrate that J. Dooge's 1982 statement on the unsolved problem of parameterization in these models remains true. Based on a review of existing parameter regionalization techniques, we postulate that the multiscale parameter regionalization (MPR) technique offers a practical and robust method that provides consistent (seamless) parameter and flux fields across scales. Herein, we develop a general model protocol to describe how MPR can be applied to a particular model and present an example application using the PCR-GLOBWB model. Finally, we discuss potential advantages and limitations of MPR in obtaining the seamless prediction of hydrological fluxes and states across spatial scales.

Description: Land surface and hydrologic models (LSM/HM) are used at diverse spatial resolutions ranging from 1–10 km in catchment-scale applications to over 50 km in global-scale applications. Application of the same model structure at different spatial scales requires that the model estimates similar fluxes independent of the model resolution and fulfills a flux-matching condition across scales. An analysis of state-of-the-art LSMs and HMs reveals that most do not have consistent and realistic parameter fields for land surface geophysical properties. Multiple experiments with the mHM, Noah-MP, PCR-GLOBWB and WaterGAP models are conducted to demonstrate the pitfalls of poor parameterization practices currently used in most operational models, which are insufficient to satisfy the flux-matching condition. These examples demonstrate that J. Dooge's 1982 statement on the unsolved problem of parameterization in these models remains true. We provide a short review of existing parameter regionalization techniques and discuss a method for obtaining seamless hydrological predictions of water fluxes and states across multiple spatial resolutions. The multiscale parameter regionalization (MPR) technique is a practical and robust method that provides consistent (seamless) parameter and flux fields across scales. A general model protocol is presented to describe how MPR can be applied to a specific model, with an example of this application using the PCR-GLOBWB model. Applying MPR to PCR-GLOBWB substantially improves the flux-matching condition. Estimation of evapotranspiration without MPR at 5 arcmin and 30 arcmin spatial resolutions for the Rhine river basin results in a difference of approximately 29 %. Applying MPR reduce this difference to 9 %. For total soil water, the differences without and with MPR are 25 % and 7 %, respectively.