Search Results
You are looking at 1 - 10 of 28 items for
- Author or Editor: Florian Pappenberger x
- Refine by Access: All Content x
Abstract
In this paper the suitability of ECMWF forecasts for hydrological applications is evaluated. This study focuses on three spatial scales: the upper Danube (which is upstream of Bratislava, Slovakia), the entire Danube catchment, and the whole of Europe. Two variables, 2-m temperature and total precipitation, are analyzed. The analysis shows that precipitation forecasts follow largely in pattern the observations. The timing of the peaks between forecasted and observed precipitation and temperature is good although precipitation amounts are often underestimated. The catchment scale influences the skill scores significantly. Small catchments exhibit a larger variance as well as larger extremes. A water balance analysis suggest a 10% underestimation by the ensemble mean and an overestimation by the high-resolution forecast over the past few years. Precipitation and temperature predictions are skillful up to days 5–7. Forecasts accumulated over a longer time frame are largely more skillful than forecasts accumulated over short time periods.
Abstract
In this paper the suitability of ECMWF forecasts for hydrological applications is evaluated. This study focuses on three spatial scales: the upper Danube (which is upstream of Bratislava, Slovakia), the entire Danube catchment, and the whole of Europe. Two variables, 2-m temperature and total precipitation, are analyzed. The analysis shows that precipitation forecasts follow largely in pattern the observations. The timing of the peaks between forecasted and observed precipitation and temperature is good although precipitation amounts are often underestimated. The catchment scale influences the skill scores significantly. Small catchments exhibit a larger variance as well as larger extremes. A water balance analysis suggest a 10% underestimation by the ensemble mean and an overestimation by the high-resolution forecast over the past few years. Precipitation and temperature predictions are skillful up to days 5–7. Forecasts accumulated over a longer time frame are largely more skillful than forecasts accumulated over short time periods.
Abstract
The land surface forms an important component of Earth system models and interacts nonlinearly with other parts such as ocean and atmosphere. To capture the complex and heterogeneous hydrology of the land surface, land surface models include a large number of parameters impacting the coupling to other components of the Earth system model.
Focusing on ECMWF’s land surface model Hydrology Tiled ECMWF Scheme of Surface Exchanges over Land (HTESSEL), the authors present in this study a comprehensive parameter sensitivity evaluation using multiple observational datasets in Europe. The authors select six poorly constrained effective parameters (surface runoff effective depth, skin conductivity, minimum stomatal resistance, maximum interception, soil moisture stress function shape, and total soil depth) and explore their sensitivity to model outputs such as soil moisture, evapotranspiration, and runoff using uncoupled simulations and coupled seasonal forecasts. Additionally, the authors investigate the possibility to construct ensembles from the multiple land surface parameters.
In the uncoupled runs the authors find that minimum stomatal resistance and total soil depth have the most influence on model performance. Forecast skill scores are moreover sensitive to the same parameters as HTESSEL performance in the uncoupled analysis. The authors demonstrate the robustness of these findings by comparing multiple best-performing parameter sets and multiple randomly chosen parameter sets. The authors find better temperature and precipitation forecast skill with the best-performing parameter perturbations demonstrating representativeness of model performance across uncoupled (and hence less computationally demanding) and coupled settings.
Finally, the authors construct ensemble forecasts from ensemble members derived with different best-performing parameterizations of HTESSEL. This incorporation of parameter uncertainty in the ensemble generation yields an increase in forecast skill, even beyond the skill of the default system.
Abstract
The land surface forms an important component of Earth system models and interacts nonlinearly with other parts such as ocean and atmosphere. To capture the complex and heterogeneous hydrology of the land surface, land surface models include a large number of parameters impacting the coupling to other components of the Earth system model.
Focusing on ECMWF’s land surface model Hydrology Tiled ECMWF Scheme of Surface Exchanges over Land (HTESSEL), the authors present in this study a comprehensive parameter sensitivity evaluation using multiple observational datasets in Europe. The authors select six poorly constrained effective parameters (surface runoff effective depth, skin conductivity, minimum stomatal resistance, maximum interception, soil moisture stress function shape, and total soil depth) and explore their sensitivity to model outputs such as soil moisture, evapotranspiration, and runoff using uncoupled simulations and coupled seasonal forecasts. Additionally, the authors investigate the possibility to construct ensembles from the multiple land surface parameters.
In the uncoupled runs the authors find that minimum stomatal resistance and total soil depth have the most influence on model performance. Forecast skill scores are moreover sensitive to the same parameters as HTESSEL performance in the uncoupled analysis. The authors demonstrate the robustness of these findings by comparing multiple best-performing parameter sets and multiple randomly chosen parameter sets. The authors find better temperature and precipitation forecast skill with the best-performing parameter perturbations demonstrating representativeness of model performance across uncoupled (and hence less computationally demanding) and coupled settings.
Finally, the authors construct ensemble forecasts from ensemble members derived with different best-performing parameterizations of HTESSEL. This incorporation of parameter uncertainty in the ensemble generation yields an increase in forecast skill, even beyond the skill of the default system.
Abstract
This paper presents an approach to postprocess ensemble forecasts for the discrete and bounded weather variable of total cloud cover. Two methods for discrete statistical postprocessing of ensemble predictions are tested: the first approach is based on multinomial logistic regression and the second involves a proportional odds logistic regression model. Applying them to total cloud cover raw ensemble forecasts from the European Centre for Medium-Range Weather Forecasts improves forecast skill significantly. Based on stationwise postprocessing of raw ensemble total cloud cover forecasts for a global set of 3330 stations over the period from 2007 to early 2014, the more parsimonious proportional odds logistic regression model proved to slightly outperform the multinomial logistic regression model.
Abstract
This paper presents an approach to postprocess ensemble forecasts for the discrete and bounded weather variable of total cloud cover. Two methods for discrete statistical postprocessing of ensemble predictions are tested: the first approach is based on multinomial logistic regression and the second involves a proportional odds logistic regression model. Applying them to total cloud cover raw ensemble forecasts from the European Centre for Medium-Range Weather Forecasts improves forecast skill significantly. Based on stationwise postprocessing of raw ensemble total cloud cover forecasts for a global set of 3330 stations over the period from 2007 to early 2014, the more parsimonious proportional odds logistic regression model proved to slightly outperform the multinomial logistic regression model.
Abstract
Heat waves represent a threat to human health and excess mortality is one of the associated negative effects. A health-based definition for heat waves is therefore relevant, especially for early warning purposes, and it is here investigated via the universal thermal climate index (UTCI). The UTCI is a bioclimate index elaborated via an advanced model of human thermoregulation that estimates the thermal stress induced by air temperature, wind speed, moisture, and radiation on the human physiology. Using France as a test bed, the UTCI was computed from meteorological reanalysis data to assess the thermal stress conditions associated with heat-attributable excess mortality in five cities. UTCI values at different climatological percentiles were defined and evaluated in their ability to identify periods of excess mortality (PEMs) over 24 years. Using verification metrics such as the probability of detection (POD), the false alarm ratio (FAR), and the frequency bias (FB), daily minimum and maximum heat stress levels equal to or above corresponding UTCI 95th percentiles (15° ± 2°C and 34.5° ± 1.5°C, respectively) for 3 consecutive days are demonstrated to correlate to PEMs with the highest sensitivity and specificity (0.69 ≤ POD ≤ 1, 0.19 ≤ FAR ≤ 0.46, 1 ≤ FB ≤ 1.48) than minimum, maximum, and mean heat stress level singularly and other bioclimatological percentiles. This finding confirms the detrimental effect of prolonged, unusually high heat stress at day- and nighttime and suggests the UTCI 95th percentile as a health-meaningful threshold for a potential heat-health watch warning system.
Abstract
Heat waves represent a threat to human health and excess mortality is one of the associated negative effects. A health-based definition for heat waves is therefore relevant, especially for early warning purposes, and it is here investigated via the universal thermal climate index (UTCI). The UTCI is a bioclimate index elaborated via an advanced model of human thermoregulation that estimates the thermal stress induced by air temperature, wind speed, moisture, and radiation on the human physiology. Using France as a test bed, the UTCI was computed from meteorological reanalysis data to assess the thermal stress conditions associated with heat-attributable excess mortality in five cities. UTCI values at different climatological percentiles were defined and evaluated in their ability to identify periods of excess mortality (PEMs) over 24 years. Using verification metrics such as the probability of detection (POD), the false alarm ratio (FAR), and the frequency bias (FB), daily minimum and maximum heat stress levels equal to or above corresponding UTCI 95th percentiles (15° ± 2°C and 34.5° ± 1.5°C, respectively) for 3 consecutive days are demonstrated to correlate to PEMs with the highest sensitivity and specificity (0.69 ≤ POD ≤ 1, 0.19 ≤ FAR ≤ 0.46, 1 ≤ FB ≤ 1.48) than minimum, maximum, and mean heat stress level singularly and other bioclimatological percentiles. This finding confirms the detrimental effect of prolonged, unusually high heat stress at day- and nighttime and suggests the UTCI 95th percentile as a health-meaningful threshold for a potential heat-health watch warning system.
Abstract
Climate trends have been observed over the recent decades in many parts of the world, but current global climate models (GCMs) for seasonal climate forecasting often fail to capture these trends. As a result, model forecasts may be biased above or below the trendline. In our previous research, we developed a trend-aware forecast postprocessing method to overcome this problem. The method was demonstrated to be effective for embedding observed trends into seasonal temperature forecasts. In this study, we further develop the method for postprocessing GCM seasonal precipitation forecasts. We introduce new formulation and evaluation features to cater for special characteristics of precipitation amounts, such as having a zero lower bound and highly positive skewness. We apply the improved method to calibrate ECMWF SEAS5 forecasts of seasonal precipitation for Australia. Our evaluation shows that the calibrated forecasts reproduce observed trends over the hindcast period of 36 years. In some regions where observed trends are statistically significant, forecast skill is greatly improved by embedding trends into the forecasts. In most regions, the calibrated forecasts outperform the raw forecasts in terms of bias, skill, and reliability. Wider applications of the new trend-aware postprocessing method are expected to boost user confidence in seasonal precipitation forecasts.
Abstract
Climate trends have been observed over the recent decades in many parts of the world, but current global climate models (GCMs) for seasonal climate forecasting often fail to capture these trends. As a result, model forecasts may be biased above or below the trendline. In our previous research, we developed a trend-aware forecast postprocessing method to overcome this problem. The method was demonstrated to be effective for embedding observed trends into seasonal temperature forecasts. In this study, we further develop the method for postprocessing GCM seasonal precipitation forecasts. We introduce new formulation and evaluation features to cater for special characteristics of precipitation amounts, such as having a zero lower bound and highly positive skewness. We apply the improved method to calibrate ECMWF SEAS5 forecasts of seasonal precipitation for Australia. Our evaluation shows that the calibrated forecasts reproduce observed trends over the hindcast period of 36 years. In some regions where observed trends are statistically significant, forecast skill is greatly improved by embedding trends into the forecasts. In most regions, the calibrated forecasts outperform the raw forecasts in terms of bias, skill, and reliability. Wider applications of the new trend-aware postprocessing method are expected to boost user confidence in seasonal precipitation forecasts.
Abstract
In the absence of a dynamical fire model that could link the emissions to the weather dynamics and the availability of fuel, atmospheric composition models, such as the European Copernicus Atmosphere Monitoring Services (CAMS), often assume persistence, meaning that constituents produced by the biomass burning process during the first day are assumed constant for the whole length of the forecast integration (5 days for CAMS). While this assumption is simple and practical, it can produce unrealistic predictions of aerosol concentration due to an excessive contribution from biomass burning. This paper introduces a time-dependent factor
Abstract
In the absence of a dynamical fire model that could link the emissions to the weather dynamics and the availability of fuel, atmospheric composition models, such as the European Copernicus Atmosphere Monitoring Services (CAMS), often assume persistence, meaning that constituents produced by the biomass burning process during the first day are assumed constant for the whole length of the forecast integration (5 days for CAMS). While this assumption is simple and practical, it can produce unrealistic predictions of aerosol concentration due to an excessive contribution from biomass burning. This paper introduces a time-dependent factor
Abstract
Early awareness of extreme precipitation can provide the time necessary to make adequate event preparations. At the European Centre for Medium-Range Weather Forecasts (ECMWF), one tool that condenses the forecast information from the Integrated Forecasting System ensemble (ENS) is the extreme forecast index (EFI), an index that highlights regions that are forecast to have potentially anomalous weather conditions compared to the local climate. This paper builds on previous findings by undertaking a global verification throughout the medium-range forecast horizon (out to 15 days) on the ability of the EFI for water vapor transport [integrated vapor transport (IVT)] and precipitation to capture extreme observed precipitation. Using the ECMWF ENS for winters 2015/16 and 2016/17 and daily surface precipitation observations, the relative operating characteristic is used to show that the IVT EFI is more skillful than the precipitation EFI in forecast week 2 over Europe and western North America. It is the large-scale nature of the IVT, its higher predictability, and its relationship with extreme precipitation that result in its potential usefulness in these regions, which, in turn, could provide earlier awareness of extreme precipitation. Conversely, at shorter lead times the precipitation EFI is more useful, although the IVT EFI can provide synoptic-scale understanding. For the whole globe, the extratropical Northern Hemisphere, the tropics, and North America, the precipitation EFI is more useful throughout the medium range, suggesting that precipitation processes not captured in the IVT are important (e.g., tropical convection). Following these results, the operational implementation of the IVT EFI is currently being planned.
Abstract
Early awareness of extreme precipitation can provide the time necessary to make adequate event preparations. At the European Centre for Medium-Range Weather Forecasts (ECMWF), one tool that condenses the forecast information from the Integrated Forecasting System ensemble (ENS) is the extreme forecast index (EFI), an index that highlights regions that are forecast to have potentially anomalous weather conditions compared to the local climate. This paper builds on previous findings by undertaking a global verification throughout the medium-range forecast horizon (out to 15 days) on the ability of the EFI for water vapor transport [integrated vapor transport (IVT)] and precipitation to capture extreme observed precipitation. Using the ECMWF ENS for winters 2015/16 and 2016/17 and daily surface precipitation observations, the relative operating characteristic is used to show that the IVT EFI is more skillful than the precipitation EFI in forecast week 2 over Europe and western North America. It is the large-scale nature of the IVT, its higher predictability, and its relationship with extreme precipitation that result in its potential usefulness in these regions, which, in turn, could provide earlier awareness of extreme precipitation. Conversely, at shorter lead times the precipitation EFI is more useful, although the IVT EFI can provide synoptic-scale understanding. For the whole globe, the extratropical Northern Hemisphere, the tropics, and North America, the precipitation EFI is more useful throughout the medium range, suggesting that precipitation processes not captured in the IVT are important (e.g., tropical convection). Following these results, the operational implementation of the IVT EFI is currently being planned.
Abstract
Flood simulation models and hazard maps are only as good as the underlying data against which they are calibrated and tested. However, extreme flood events are by definition rare, so the observational data of flood inundation extent are limited in both quality and quantity. The relative importance of these observational uncertainties has increased now that computing power and accurate lidar scans make it possible to run high-resolution 2D models to simulate floods in urban areas. However, the value of these simulations is limited by the uncertainty in the true extent of the flood. This paper addresses that challenge by analyzing a point dataset of maximum water extent from a flood event on the River Eden at Carlisle, United Kingdom, in January 2005. The observation dataset is based on a collection of wrack and water marks from two postevent surveys. A smoothing algorithm for identifying, quantifying, and reducing localized inconsistencies in the dataset is proposed and evaluated showing positive results. The proposed smoothing algorithm can be applied in order to improve flood inundation modeling assessment and the determination of risk zones on the floodplain.
Abstract
Flood simulation models and hazard maps are only as good as the underlying data against which they are calibrated and tested. However, extreme flood events are by definition rare, so the observational data of flood inundation extent are limited in both quality and quantity. The relative importance of these observational uncertainties has increased now that computing power and accurate lidar scans make it possible to run high-resolution 2D models to simulate floods in urban areas. However, the value of these simulations is limited by the uncertainty in the true extent of the flood. This paper addresses that challenge by analyzing a point dataset of maximum water extent from a flood event on the River Eden at Carlisle, United Kingdom, in January 2005. The observation dataset is based on a collection of wrack and water marks from two postevent surveys. A smoothing algorithm for identifying, quantifying, and reducing localized inconsistencies in the dataset is proposed and evaluated showing positive results. The proposed smoothing algorithm can be applied in order to improve flood inundation modeling assessment and the determination of risk zones on the floodplain.
Abstract
Land surface models (LSMs) have traditionally been designed to focus on providing lower-boundary conditions to the atmosphere with less focus on hydrological processes. State-of-the-art application of LSMs includes a land data assimilation system (LDAS), which incorporates available land surface observations to provide an improved realism of surface conditions. While improved representations of the surface variables (such as soil moisture and snow depth) make LDAS an essential component of any numerical weather prediction (NWP) system, the related increments remove or add water, potentially having a negative impact on the simulated hydrological cycle by opening the water budget. This paper focuses on evaluating how well global NWP configurations are able to support hydrological applications, in addition to the traditional weather forecasting. River discharge simulations from two climatological reanalyses are compared: one “online” set, which includes land–atmosphere coupling and LDAS with an open water budget, and an “offline” set with a closed water budget and no LDAS. It was found that while the online version of the model largely improves temperature and snow depth conditions, it causes poorer representation of peak river flow, particularly in snowmelt-dominated areas in the high latitudes. Without addressing such issues there will never be confidence in using LSMs for hydrological forecasting applications across the globe. This type of analysis should be used to diagnose where improvements need to be made; considering the whole Earth system in the data assimilation and coupling developments is critical for moving toward the goal of holistic Earth system approaches.
Abstract
Land surface models (LSMs) have traditionally been designed to focus on providing lower-boundary conditions to the atmosphere with less focus on hydrological processes. State-of-the-art application of LSMs includes a land data assimilation system (LDAS), which incorporates available land surface observations to provide an improved realism of surface conditions. While improved representations of the surface variables (such as soil moisture and snow depth) make LDAS an essential component of any numerical weather prediction (NWP) system, the related increments remove or add water, potentially having a negative impact on the simulated hydrological cycle by opening the water budget. This paper focuses on evaluating how well global NWP configurations are able to support hydrological applications, in addition to the traditional weather forecasting. River discharge simulations from two climatological reanalyses are compared: one “online” set, which includes land–atmosphere coupling and LDAS with an open water budget, and an “offline” set with a closed water budget and no LDAS. It was found that while the online version of the model largely improves temperature and snow depth conditions, it causes poorer representation of peak river flow, particularly in snowmelt-dominated areas in the high latitudes. Without addressing such issues there will never be confidence in using LSMs for hydrological forecasting applications across the globe. This type of analysis should be used to diagnose where improvements need to be made; considering the whole Earth system in the data assimilation and coupling developments is critical for moving toward the goal of holistic Earth system approaches.
Abstract
Seasonal streamflow prediction skill can derive from catchment initial hydrological conditions (IHCs) and from the future seasonal climate forecasts (SCFs) used to produce the hydrological forecasts. Although much effort has gone into producing state-of-the-art seasonal streamflow forecasts from improving IHCs and SCFs, these developments are expensive and time consuming and the forecasting skill is still limited in most parts of the world. Hence, sensitivity analyses are crucial to funnel the resources into useful modeling and forecasting developments. It is in this context that a sensitivity analysis technique, the variational ensemble streamflow prediction assessment (VESPA) approach, was recently introduced. VESPA can be used to quantify the expected improvements in seasonal streamflow forecast skill as a result of realistic improvements in its predictability sources (i.e., the IHCs and the SCFs)—termed “skill elasticity”—and to indicate where efforts should be targeted. The VESPA approach is, however, computationally expensive, relying on multiple hindcasts having varying levels of skill in IHCs and SCFs. This paper presents two approximations of the approach that are computationally inexpensive alternatives. These new methods were tested against the original VESPA results using 30 years of ensemble hindcasts for 18 catchments of the contiguous United States. The results suggest that one of the methods, end point blending, is an effective alternative for estimating the forecast skill elasticities yielded by the VESPA approach. The results also highlight the importance of the choice of verification score for a goal-oriented sensitivity analysis.
Abstract
Seasonal streamflow prediction skill can derive from catchment initial hydrological conditions (IHCs) and from the future seasonal climate forecasts (SCFs) used to produce the hydrological forecasts. Although much effort has gone into producing state-of-the-art seasonal streamflow forecasts from improving IHCs and SCFs, these developments are expensive and time consuming and the forecasting skill is still limited in most parts of the world. Hence, sensitivity analyses are crucial to funnel the resources into useful modeling and forecasting developments. It is in this context that a sensitivity analysis technique, the variational ensemble streamflow prediction assessment (VESPA) approach, was recently introduced. VESPA can be used to quantify the expected improvements in seasonal streamflow forecast skill as a result of realistic improvements in its predictability sources (i.e., the IHCs and the SCFs)—termed “skill elasticity”—and to indicate where efforts should be targeted. The VESPA approach is, however, computationally expensive, relying on multiple hindcasts having varying levels of skill in IHCs and SCFs. This paper presents two approximations of the approach that are computationally inexpensive alternatives. These new methods were tested against the original VESPA results using 30 years of ensemble hindcasts for 18 catchments of the contiguous United States. The results suggest that one of the methods, end point blending, is an effective alternative for estimating the forecast skill elasticities yielded by the VESPA approach. The results also highlight the importance of the choice of verification score for a goal-oriented sensitivity analysis.