HESS Opinions "Should we apply bias correction to global and regional climate model data?"

Despite considerable progress in recent years, output of both global and regional circulation models is still afflicted with biases to a degree that precludes its direct use, especially in climate change impact studies. This is well known, and to overcome this problem, bias correction (BC; i.e. the correction of model output towards observations in a post-processing step) has now become a standard procedure in climate change impact studies. In this paper we argue that BC is currently often used in an invalid way: it is added to the GCM/RCM model chain without sufficient proof that the consistency of the latter (i.e. the agreement between model dynamics/model output and our judgement) as well as the generality of its applicability increases. BC methods often impair the advantages of circulation models by altering spatiotemporal field consistency, relations among variables and by violating conservation principles. Currently used BC methods largely neglect feedback mechanisms, and it is unclear whether they are time-invariant under climate change conditions. Applying BC increases agreement of climate model output with observations in hindcasts and hence narrows the uncertainty range of simulations and predictions without, however, providing a satisfactory physical justification. This is in most cases not transparent to the end user.We argue that this hides rather than reduces uncertainty, which may lead to avoidable forejudging of end users and decision makers. We present here a brief overview of state-of-the-art bias correction methods, discuss the related assumptions and implications, draw conclusions on the validity of bias correction and propose ways to cope with biased output of circulation models in the short term and how to reduce the bias in the long term. The most promising strategy for improved future global and regional circulation model simulations is the increase in model resolution to the convection-permitting scale in combination with ensemble predictions based on sophisticated approaches for ensemble perturbation. With this article, we advocate communicating the entire uncertainty range associated with climate change predictions openly and hope to stimulate a lively discussion on bias correction among the atmospheric and hydrological community and end users of climate change impact studies

Simulation of semi-arid biomass plantations and irrigation using the WRF-NOAH model – a comparison with observations from Israel

A 10 × 10 km irrigated biomass plantation was simulated in an arid region of Israel to simulate diurnal energy balances during the summer of 2012 (JJA). The goal is to examine daytime horizontal flux gradients between plantation and desert. Simulations were carried out within the coupled WRF-NOAH atmosphere/land surface model. MODIS land surface data was adjusted by prescribing tailored land surface and soil/plant parameters, and by adding a controllable sub-surface irrigation scheme to NOAH. Two model cases studies were compared – Impact and Control. Impact simulates the irrigated plantation. Control simulates the existing land surface, where the predominant land surface is bare desert soil. Central to the study is parameter validation against land surface observations from a desert site and from a 400 ha Simmondsia chinensis (jojoba) plantation. Control was validated with desert observations, and Impact with Jojoba observations. Model evapotranspiration was validated with two Penman–Monteith estimates based on the observations. Control simulates daytime desert conditions with a maximum deviation for surface 2 m air temperatures (T2) of 0.2 °C, vapour pressure deficit (VPD) of 0.25 hPa, wind speed (U) of 0.5 m s−1, surface radiation (Rn) of 25 W m−2, soil heat flux (G) of 30 W m−2 and 5 cm soil temperatures (ST5) of 1.5 °C. Impact simulates irrigated vegetation conditions with a maximum deviation for T2 of 1–1.5 °C, VPD of 0.5 hPa, U of 0.5 m s−1, Rn of 50 W m−5, G of 40 W m−2 and ST5 of 2 °C. Latent heat curves in Impact correspond closely with Penman–Monteith estimates, and magnitudes of 160 W m−2 over the plantation are usual. Sensible heat fluxes, are around 450 W m−2 and are at least 100–110 W m−2 higher than the surrounding desert. This surplus is driven by reduced albedo and high surface resistance, and demonstrates that high evaporation rates may not occur over Jojoba if irrigation is optimized. Furthermore, increased daytime T2 over plantations highlight the need for hourly as well as daily mean statistics. Daily mean statistics alone may imply an overall cooling effect due to surplus nocturnal cooling, when in fact a daytime warming effect is observed

Continuous high-resolution midlatitude-belt simulations for July–August 2013 with WRF

Increasing computational resources and the demands of impact modelers, stake holders, and society envision seasonal and climate simulations with the convection-permitting resolution. So far such a resolution is only achieved with a limited-area model whose results are impacted by zonal and meridional boundaries. Here, we present the setup of a latitude-belt domain that reduces disturbances originating from the western and eastern boundaries and therefore allows for studying the impact of model resolution and physical parameterization. The Weather Research and Forecasting (WRF) model coupled to the NOAH land–surface model was operated during July and August 2013 at two different horizontal resolutions, namely 0.03 (HIRES) and 0.12° (LOWRES). Both simulations were forced by the European Centre for Medium-Range Weather Forecasts (ECMWF) operational analysis data at the northern and southern domain boundaries, and the high-resolution Operational Sea Surface Temperature and Sea Ice Analysis (OSTIA) data at the sea surface.The simulations are compared to the operational ECMWF analysis for the representation of large-scale features. To analyze the simulated precipitation, the operational ECMWF forecast, the CPC MORPHing (CMORPH), and the ENSEMBLES gridded observation precipitation data set (E-OBS) were used as references.Analyzing pressure, geopotential height, wind, and temperature fields as well as precipitation revealed (1) a benefit from the higher resolution concerning the reduction of monthly biases, root mean square error, and an improved Pearson skill score, and (2) deficiencies in the physical parameterizations leading to notable biases in distinct regions like the polar Atlantic for the LOWRES simulation, the North Pacific, and Inner Mongolia for both resolutions.In summary, the application of a latitude belt on a convection-permitting resolution shows promising results that are beneficial for future seasonal forecasting

HESS Opinions “should we apply bias correction to global and regional climate model data?”

Despite considerable progress in recent years, output of both global and regional circulation models is still afflicted with biases to a degree that precludes its direct use, especially in climate change impact studies. This is well known, and to overcome this problem, bias correction (BC; i.e. the correction of model output towards observations in a post-processing step) has now become a standard procedure in climate change impact studies. In this paper we argue that BC is currently often used in an invalid way: it is added to the GCM/RCM model chain without sufficient proof that the consistency of the latter (i.e. the agreement between model dynamics/model output and our judgement) as well as the generality of its applicability increases. BC methods often impair the advantages of circulation models by altering spatiotemporal field consistency, relations among variables and by violating conservation principles. Currently used BC methods largely neglect feedback mechanisms, and it is unclear whether they are time-invariant under climate change conditions. Applying BC increases agreement of climate model output with observations in hindcasts and hence narrows the uncertainty range of simulations and predictions without, however, providing a satisfactory physical justification. This is in most cases not transparent to the end user.We argue that this hides rather than reduces uncertainty, which may lead to avoidable forejudging of end users and decision makers. We present here a brief overview of state-of-the-art bias correction methods, discuss the related assumptions and implications, draw conclusions on the validity of bias correction and propose ways to cope with biased output of circulation models in the short term and how to reduce the bias in the long term. The most promising strategy for improved future global and regional circulation model simulations is the increase in model resolution to the convection-permitting scale in combination with ensemble predictions based on sophisticated approaches for ensemble perturbation. With this article, we advocate communicating the entire uncertainty range associated with climate change predictions openly and hope to stimulate a lively discussion on bias correction among the atmospheric and hydrological community and end users of climate change impact studies

Interactions of Generated Weather Raster and Soil Profiles in Simulating Adaptive Crop Management and Consequent Yields for Five Major Crops throughout a Region in Southern Germany

Klimaanpassung und MitigationThe ability of bioeconomic simulation modelling to realistically predict agricultural adaptation is limited by the degree of detail in crucial model components. Model robustness must be tested before localized calibrations can be applied to regions of heterogenous environmental conditions. The agent-based model FARMACTOR was used to simulate the timing of field management actions (planting, harvest etc.) in response to environmental conditions, and consequent yields of winter wheat, barley and rapeseed, spring barley and silage maize as the predominant crops in a distinct region of Germany, by linking weather data and the crop growth simulation model EXPERT-N. The integrated models were calibrated to observed experimental data and official phenological observations and then run from 1990 to 2009, forced with climate data from ERA-interim Reanalyses data which was downscaled with the Weather and Research Forecast (WRF) model to a 12 X 12 km² grid. Variability in regional soils was replicated with 10 different soil profiles mapped at 1/25,000 scale. The nature of the forcing climate data dictates temporal aggregation for analysis, so that validity is examined by comparing mean simulated planting and harvest dates and yields to official records in the area. The mean predicted planting dates are very close to observations over the period, within a few days of observations, but show less variance. Harvest dates are accurately predicted as well, within one to two weeks, and the variances are closer to observations. Predicted winter wheat yields are well simulated in comparison to observed data, but maize yields are underestimated, while winter and spring barley and winter rapeseed yields are greater than observed district ("Landkreis") yields. The degree of variance in simulated yields is acceptable in wheat, winter barley and maize, but excessive in spring barley and winter rapeseed. Cross-sectional examination of yields shows that the different soil profiles are responsible for more yield variance than simulated weather cells in all crops. While the coupled models appear accurate in predicting crop management dates and physiological development, the inaccuracy in yields in all crops except winter wheat calls into question the reliability of the integrated models when applied, as is, outside of calibration conditions. That soil parameterization is responsible for more variance than generated weather is helpful in seeking to improve performance and encouraging in terms of the method of weather generation. Reliable extension of the coupled models to include all soils in an area together with artificial spatial climatic variability may require regionalized calibration to increase crop model stability

Afforestation impact on soil temperature in regional climate model simulations over Europe

In the context of the first phase of the Coordinated Regional Climate Downscaling Experiment in the European domain (EURO-CORDEX) flagship plot study on Land Use and Climate Across Scales (LUCAS), we investigate the biophysical impact of afforestation on the seasonal cycle of soil temperature over the European continent with an ensemble of 10 regional climate models. For this purpose, each ensemble member performed two idealized land cover experiments in which Europe is covered either by forests or grasslands. The multi-model mean exhibits a reduction of the annual amplitude of soil temperature (AAST) due to afforestation over all European regions, although this is not a robust feature among the models. In the Mediterranean, the spread of simulated AAST response to afforestation is between −4 and +2 ∘C at 1 m below the ground, while in Scandinavia the inter-model spread ranges from −7 to +1 ∘C. We show that the large range in the simulated AAST response is due to the representation of the summertime climate processes and is largely explained by inter-model differences in leaf area index (LAI), surface albedo, cloud fraction and soil moisture, when all combined into a multiple linear regression. The changes in these drivers essentially determine the ratio between the increased radiative energy at surface (due to lower albedo in forests) and the increased sum of turbulent heat fluxes (due to mixing-facilitating characteristics of forests), and consequently decide the changes in soil heating with afforestation in each model. Finally, we pair FLUXNET sites to compare the simulated results with observation-based evidence of the impact of forest on soil temperature. In line with models, observations indicate a summer ground cooling in forested areas compared to open lands. The vast majority of models agree with the sign of the observed reduction in AAST, although with a large variation in the magnitude of changes. Overall, we aspire to emphasize the biophysical effects of afforestation on soil temperature profile with this study, given that changes in the seasonal cycle of soil temperature potentially perturb crucial biochemical processes. Robust knowledge on biophysical impacts of afforestation on soil conditions and its feedbacks on local and regional climate is needed in support of effective land-based climate mitigation and adaption policies

Precipitation frequency in Med-CORDEX and EURO-CORDEX ensembles from 0.44° to convection-permitting resolution: impact of model resolution and convection representation

Recent studies using convection-permitting (CP) climate simulations have demonstrated a step-change in the representation of heavy rainfall and rainfall characteristics (frequency-intensity) compared to coarser resolution Global and Regional climate models. The goal of this study is to better understand what explains the weaker frequency of precipitation in the CP ensemble by assessing the triggering process of precipitation in the different ensembles of regional climate simulations available over Europe. We focus on the statistical relationship between tropospheric temperature, humidity and precipitation to understand how the frequency of precipitation over Europe and the Mediterranean is impacted by model resolution and the representation of convection (parameterized vs. explicit). We employ a multi-model data-set with three different resolutions (0.44°, 0.11° and 0.0275°) produced in the context of the MED-CORDEX, EURO-CORDEX and the CORDEX Flagship Pilot Study "Convective Phenomena over Europe and the Mediterranean" (FPSCONV). The multi-variate approach is applied to all model ensembles, and to several surface stations where the integrated water vapor (IWV) is derived from Global Positioning System (GPS) measurements. The results show that all model ensembles capture the temperature dependence of the critical value of IWV (IWVcv), above which an increase in precipitation frequency occurs, but the differences between the models in terms of the value of IWVcv, and the probability of its being exceeded, can be large at higher temperatures. The lower frequency of precipitation in convection-permitting simulations is not only explained by higher temperatures but also by a higher IWVcv necessary to trigger precipitation at similar temperatures, and a lower probability to exceed this critical value. The spread between models in simulating IWVcv and the probability of exceeding IWVcv is reduced over land in the ensemble of models with explicit convection, especially at high temperatures, when the convective fraction of total precipitation becomes more important and the influence of the representation of entrainment in models thus becomes more important. Over lowlands, both model resolution and convection representation affect precipitation triggering while over mountainous areas, resolution has the highest impact due to orography-induced triggering processes. Over the sea, since lifting is produced by large-scale convergence, the probability to exceed IWVcv does not depend on temperature, and the model resolution does not have a clear impact on the results

The first multi-model ensemble of regional climate simulations at kilometer-scale resolution. Part I: Evaluation of precipitation

Here we present the first multi-model ensemble of regional climate simulations at kilometer-scale horizontal grid spacing over a decade long period. A total of 23 simulations run with a horizontal grid spacing of ∼ 3 km, driven by ERA-Interim reanalysis, and performed by 22 European research groups are analysed. Six different regional climate models (RCMs) are represented in the ensemble. The simulations are compared against available high-resolution precipitation observations and coarse resolution (∼ 12 km) RCMs with parameterized convection. The model simulations and observations are compared with respect to mean precipitation, precipitation intensity and frequency, and heavy precipitation on daily and hourly timescales in different seasons. The results show that kilometer-scale models produce a more realistic representation of precipitation than the coarse resolution RCMs. The most significant improvements are found for heavy precipitation and precipitation frequency on both daily and hourly time scales in the summer season. In general, kilometer-scale models tend to produce more intense precipitation and reduced wet-hour frequency compared to coarse resolution models. On average, the multi-model mean shows a reduction of bias from ∼ −40 at 12 km to ∼ −3 at 3 km for heavy hourly precipitation in summer. Furthermore, the uncertainty ranges i.e. the variability between the models for wet hour frequency is reduced by half with the use of kilometer-scale models. Although differences between the model simulations at the kilometer-scale and observations still exist, it is evident that these simulations are superior to the coarse-resolution RCM simulations in the representing precipitation in the present-day climate, and thus offer a promising way forward for investigations of climate and climate change at local to regional scales. © 2021, The Author(s)

The first multi-model ensemble of regional climate simulations at kilometer-scale resolution, part I: evaluation of precipitation

Here we present the first multi-model ensemble of regional climate simulations at kilometer-scale horizontal grid spacing over a decade long period. A total of 23 simulations run with a horizontal grid spacing of ∼3 km, driven by ERA-Interim reanalysis, and performed by 22 European research groups are analysed. Six different regional climate models (RCMs) are represented in the ensemble. The simulations are compared against available high-resolution precipitation observations and coarse resolution (∼ 12 km) RCMs with parameterized convection. The model simulations and observations are compared with respect to mean precipitation, precipitation intensity and frequency, and heavy precipitation on daily and hourly timescales in different seasons. The results show that kilometer-scale models produce a more realistic representation of precipitation than the coarse resolution RCMs. The most significant improvements are found for heavy precipitation and precipitation frequency on both daily and hourly time scales in the summer season. In general, kilometer-scale models tend to produce more intense precipitation and reduced wet-hour frequency compared to coarse resolution models. On average, the multi-model mean shows a reduction of bias from ∼ −40% at 12 km to ∼ −3% at 3 km for heavy hourly precipitation in summer. Furthermore, the uncertainty ranges i.e. the variability between the models for wet hour frequency is reduced by half with the use of kilometer-scale models. Although differences between the model simulations at the kilometer-scale and observations still exist, it is evident that these simulations are superior to the coarse-resolution RCM simulations in the representing precipitation in the present-day climate, and thus offer a promising way forward for investigations of climate and climate change at local to regional scales