About the Scenario Service - Hydrology

These pages present the development of sub surface water, lakes and other water bodies as a consequence of the changing climate of the 21st century. The service builds on scenarios from multiple regional climate models forced by several different global climate models. The climate data are then used to force the hydrological model S-HYPE, in order to describe hydrological processes. The results are summarised in 262 geographical areas of the service.

In climate model calculations, information on future changes in the atmosphere are used. Th climate models describe the connections between processes in the atmosphere, on the surface and in the oceans. The results of climate model simulations are refined for each of the counties in Sweden, and cover the period 1951-2100. The climate model simulations are a part of the international research programme CORDEX (cordex.org) and Copernicus Climate Change Service (climate.copernicus.eu).

The climate models suffer from systematic errors, so-called bias, in the forcing variables temperature and precipitation. Before impact modelling with S-HYPE, these parameters are therefore bias adjusted to fit the calibration of the hydrological model.
The hydrological model describes water’s path through the ground, lakes and other water bodies, and how much water that is returned to the atmosphere through evaporation or return to the oceans. The S-HYPE models contains about 40 000 sub-catchments to simulated Swedish hydrology, but the results are summarised to 262 named geographical areas in the service.

Climate models

Climate models are used to calculate future climate. These are 3-dimensional representations of the atmosphere, the land surface, oceans, lakes and ice. In the model the atmosphere is divided into a 3-dimensional grid, across the surface and upwards. To get good results processes outside the entire atmosphere must be considered, that is around the globe and up in the air. Such models are called global climate models. The time evolution of different meteorological parameters are calculated in each cell of the grid.

Climate models create enormous amounts of information and requires therefore large amounts of computer power, which means that the 3-dimensional grid must be restrained. The grid in a global climate model is sparse, which gives less details on the regional scale. The grid is put on a smaller region in the regional model, for example Europe. In that way more detail can be achieved for a smaller region without using too much computer power.

The results from a global climate model governs what happens outside the model domain in the regional model. In that way the regional model also considers what happens outside the domain. In this service several regional climate models forced with several different global models are used. The regional models cover Europe and have a resolution (size of the grid boxes) of around 12.5x12.5 km.

Illustration som beskriver ett gridnät på en jordglob ner till över Europa till lokalt
Figure. A global model is required to simulate the global climate. By simulating only a part of the Earth in a regional model more detailed information can be achieved, although data from a global model is needed. Subsequent modelling may also be needed, for calculations of indicators, ocean currents or modelling of connectivity etc.

High resolution climate models require large resources, but more details can be achieved by so called downscaling. The higher resolution yields more data points which gives more details, climatological processes are also better described. This applies especially to short lived or local events. Precipitation in general and extreme precipitation in particular is an example of what can be improved by higher resolution. Precipitation may vary largely in space and time. That is better described by a model of high resolution (Rummukainen, 2010; Rummukainen, 2016).

Emissions Scenarios

Model simulations of the climate are based on emissions scenarios or radiative scenarios. Emissions scenarios are assumptions of future emissions of greenhouse gases. They are based on assumptions of future development of world economy, population growth, globalisation, transitions to environmentally friendly technology amongst others.

The amount of future emissions is affected by all these factors, which in turn affects the greenhouse effect. A measure of how the greenhouse effect is changing in the future is radiative forcing, measured in power per square metre (W/m2). The more emissions of greenhouse gases the more radiative forcing. Such scenarios are called RCP scenarios (Representative Concentration Pathways (Moss et al., 2010; van Vuuren et al., 2011)).
In this analysis three scenarios are used:

  • RCP2.6: Powerful climate policy makes emissions of greenhouse gases to culminate in year 2020, radiative forcing reaches 2.6 W/m² year 2100 (used in IPCC, AR5). This scenario is the one closest to the ambitions from the Paris Agreement.
  • RCP4.5: Strategies for reduced greenhouse gas emissions yield a stabilisation of the radiative forcing at 4.5 W/m² before year 2100 (used in IPCC, AR5).
  • RCP8.5: Increasing greenhouse gas emissions makes radiative forcing to reach 8.5 W/m² year 2100 (used in IPCC, AR5). 

Climate scenarios

A climate scenario is a combination of emissions or radiative scenario, global climate model, regional climate model and a time period of choice. See the table in the section on ensembles for more information on the climate models used.

The regional models start somewhere between 1951 and 1971 depending on the model. It is possible that the results deviate from observations already in the beginning. This is because the global models used not exactly represent the climate of today in the same way in all simulations. The period 1971-2000 is used as reference. Results for the future shows anomalies compared to the average of 1971-2000.

Since the results are provided as a grid, so called gridded data, it can be difficult to directly compare model results and observations. Observations give the conditions at a specific point, whereas the models give the average for the entire grid cell. The description of observed climate in this service is based on gridded data and can thus more easily be compared to climate model data. 

Reference data

The hydrological model is calibrated to produce accurate river discharge when forced with the PTHBV meteorological gridded data set (Johansson and Chen, 2003). PTHBV is constructed from several hundred station observations of precipitation and temperature across Sweden. The data is produced at a 4x4 km2 grid. Besides the use of PTHBV to calibrate the S-HYPE model, it is also used as reference for the bias adjustment of climate models.

Bias adjustment

Complex and intertwined processes in Earth’s climate system are not always straightforward to implement in climate models due to resolution, or simply by a lack of knowledge about some processes. The climate models will therefore inevitably introduce systematic deviations, or bias, of the modelled climate, in comparison to observations.

Bias will differ between variables, regions of the world, and models. When comparing a historical and a future projection, moderate bias is mostly not causing issues to the assessment of climate change. However, bias can prevent an indicator from being reasonably described when it is based on absolute thresholds, such as frost days, tropical nights, or precipitation above a defined intensity. Further, impact models are often sensitive to bias when they are calibrated with observational data. For example, the storage of water in frozen form as ice and snow in hydrological models requires a good description of temperature.

To prevent bias from generating misleading results in threshold based indicators, a so-called bias adjustment is first applied to the data. The data in the service has been bias adjusted using the MIdAS method (MultI-scale bias AdjuStment; Berg et al., SMHI Klimatologirapport nr 63). MIdAS is developed at SMHI and has been shown to perform well, also compared to international state-of-the-art methods. First, the method derives a mapping of values from the biased model to the reference data set. Remapping is then performed across the complete timeseries of the models, across the full distribution of values, and also for periods outside the reference period, such as for future projections.

Modelling hydrology

The effects of climate change on hydrology is modelled with the hydrological model S-HYPE (Lindström et al., 2010). S-HYPE calculates the river discharge and other hydrologically relevant parameters based on numerical representations of storage and flowprocesses in the ground, lakes and other water bodies.

The performance of S-HYPE can be explored in the SMHI web-service vattenwebb.

S-HYPE is delineated along the borders of about 40 000 catchments within Sweden, and in bordering upstream areas. The model encompasses also the influence of lakes and dams on the flow dynamics. Further, the model describes river regulations, with seasonally dependent settings of storage, production flow etc. S-HYPE also includes some artificial water transfers between catchments, which are mainly used for hydro power. The modelled regulation practices are stationary in time, i.e. they do not adjust to a changing climate, which needs to be considered when interpreting the results.

The service presents results for 262 river branches (BARO), which include complete or parts of water bodies, as well as coastal areas divided by municipalities (län). The larger areas have been constructed to fit the spatial detail suitable for the service and its users, and to account for the inherent spatial uncertainties in climate modelling. The presented results provide an overview of the expected mean changes of the different indicators across each BARO. Indicators based on river discharge are represented by a single outlet in the BARO. The outlet is sometimes upstream from the main outlet of the BARO so that the outlet point better describes the complete catchment area.

Scenarios are not forecasts

The results from calculations made with climate models presented here are called climate scenarios, they are not weather forecasts. Climate scenarios are based on assumptions of future emissions and amounts of greenhouse gases. They represent the statistical behaviour of the weather of a longer period of time - that is the climate.

Indicators

Hydrological climate indicators are calculated from standard hydrological parameters, such as river discharge.

Indicators
Name Definition Unit
Low flow days Number of days with river discharge < the mean low flow, which is defined as the mean of the day with the lowest river discharge of each year in 1971-2000. days
Low soil moisture days Number of days with soil moisture < the mean low soil moisture, which is defined as the mean of the day with the lowest soil moisture of each year in 1971-2000. days
Effektive precipitation (mean) The difference between precipitation and evapotranspiration from the ground and canopy over the day. Presented as the mean over the period, and the geographical area. mm/month
Air temperature (mean) The mean air temperature at 2m above the surface for the selected period, and the geographical area. °C
Runoff (mean) Mean river discharge for the selected period, and the geographical area. mm/month
Soil moisture (mean) Mean soil mositure for the selected period, and the geographical area. Unitless
Precipitation (mean) Mean precipitation for the selected period, and the geographical area. mm/day
River discharge (mean) Mean river discharge for the selected period, and the geographical area. m3/s
Flood recurrence (2,5,10, or 50-yr) The return level of river discharge for the selected return period, selected period, and geographical area. m3/s
Snow water (max) The mean of each year’s maximum snow water equivalent for the selected period and the geographical area. mm

Why are different reference periods used?

The present climate is described by the last finished 30-year normal period, 1991-2020. Previous climates are described by earlier normal periods, and for the study of climate change different normal periods are compared.

To study the ongoing climate change the reference normal period 1961-1990 is usually used. In some cases it might be relevant to use other normal periods, for example when pre-industrial time is used as a reference. Within climate science other periods are sometimes used. In the beginning of the 21st century the period 1961-1990 was not seen as representative for the current climate, and it was still too early to shift to 1991-2020.

Because of this IPCC has, for example, sometimes used the last 20-30 years as the reference to which the calculated climate change is compared with. Sometimes there are problems with choosing a period too close to the present. In many climate simulations the greenhouse gas emissions scenarios start already in 2006, which means that the future starts in 2006 and that the years from then on do not represent actual levels of greenhouse gases.

There may also be practical reasons not to use 1961-1990. Climate model simulations require large computer resources. Therefore, it sometimes happens, that simulations start after 1961 to save computer power. As a compromise the period 1971-2000 is then used as a compromise between computer power and historical reference.

Read more about reference periods (in Swedish)

Calculations for the end of the century

The end of the century period is defined as the years 2071-2100. A few climate models lack data for the last years of that period. Formally, they can therefore not be used to derive information for that period. However, the lack of a small part of the 30-year period makes very little difference to the end result. Therefore, the service uses a condition that at least 25 years must be available in the period to allow the model to be used in the calculations, and still be presented as valid for the period 2071-2100.

The hydrological modelling extends the data of these models to the year 2100 by repeating the previous years in the timeseries fed to the hydrological model. For example, in one case where the period December 2099 to December 2100 was lacking, the values for December 2098 was copied to December 2099, and then the complete year 2099 was copied to 2100. Because the hydrological calculations are affected by slow processes in the ground water, and storages, the resulting modelling will still differ between the two years. This approach is chosen because of the very few models that lack data, and the insignificant impact on the end result.

About ensembles

What is an ensemble? 

An ensemble is a collection of climate model scenarios, or ensemble members. The members in the ensemble can differ regarding the driving global climate model, the regional climate model, or the emission scenario. 

Why should one use model ensembles? 

An ensemble provides an overview of the spread in the future climate scenarios. The spread informs about uncertainty in the results, which is important for assessing the robustness and reliability of the information. The reliability of the result increases if several models produce a consistent result. 

The role of the climate model 

The climate models, both global and regional, will provide different responses to the forcing from a given emissions scenario. These differences arise due to the models’ approach to simulate different physical processes in the atmosphere, which relate to our incomplete understanding of certain processes and the ability to numerically simulate them.

The global climate models determine the main climate sensitivity, and the large-scale atmospheric circulation. They will therefore have the largest impact on the simulated large-scale climate, such as wintertime mean temperature. Other processes are determined at smaller scales, such as summertime precipitation which arise from local cloud processes, where the regional climate models to a large extent dictate the results (Kjellström et al., 2018; Sørland et al., 2018).

It is challenging to select a single best performing climate model, because the models differ in performance depending on the region, the variable and season studied. Including several models acts to increase the robustness and quality of the ensemble mean result, and the most common approach is to include as many models as possible in an ensemble.

Another kind of ensemble can be configured by making several simulations with a single global climate model, where the initial conditions at the start-up of the model leads to different results. This is an effect of the chaotic development of the atmosphere, which in itself leads to a spread in the results, which is related to the natural variability of the atmosphere.

Ensembles describe the reliability of the results 

The spread within the ensemble provides information on the reliability of the results. The set up of the ensemble can describe different uncertainties, such as model uncertainty or natural variability. 
This service presents two measures of robustness and model spread: the relative number of models that show and increase, and the standard deviation across the model ensemble.

The relative number of models that show and increase, presents the model agreement on the sign of the change. If most models agree on an increase, the result can be considered robust. Further, if most models disagree on an increase, it means that they agree on a decrease. If the models show little agreement on the sign of the change, the result is either not robust, or there is no significant change to be expected from the studied variable. The latter requires further investigations.
Standard deviation is a measure of the model spread. Even a robust signal can have significant spread when the magnitude of the change is uncertain. A large spread results in a large value of the standard deviation relative to the magnitude of the change.

The time period matters 

The number of models and scenarios to include in an ensemble partly depends on the time period studied. Generally, for the near future (a few decades into the future), the initial conditions influence the results and demands a large ensemble of different model and different initial states. At longer time horizons, the emission scenarios become more important than the initial conditions, but also the sensitivity of the climate model to the greenhouse gas forcing plays a large role.

Natural variability is important in the short term perspective

The climate system has it’s own natural variability in addition to anthropogenic effects. These natural fluctuations between years, or between decades complicates the analysis of calculated climate scenarios. Especially when studying climate change in the near term future.

Climate change until year 2100 is calculated to be big enough to be clearly visible even with large variations between years. The natural variations of the climate can not be exactly predicted (for example for a certain date). However, it is possible to study natural variability by creating an ensemble of several climate scenarios based on one emissions scenario, but with different initial conditions. By the end of the century the largest uncertainty comes from the choice of climate model and the choice of emissions scenario. 

Models used in this service

The climate model ensemble of the hydrological service differs somewhat from that of the meteorological service. This is mainly because of funding restraints which has limited the number of models to include. A sub-selection was made to include only those global and regional climate model combinations that include all three emission scenarios. This result in an ensemble of five different global models and seven regional models, and a complete set of 16 unique members. 

Skärmdump av tabell över klimatmodeller
Climate models included in the hydrological service. The notation 26, 45, 85 refer to the scenarios RCP2.6, RCP4.5 and RCP8.5, respectively. Further, notations r1, r3, r12 refer to different realizations of the same climate model, i.e. a model with the same setup but with different initial conditions.  Enlarge Image

Table of models (14 kB, xlsx)

Robustness of the results

Robustness is here a measure of how certain or uncertain the climate change signal is. Since the results are based on several different simulations with climate models it is possible to do statistical analyses of how robust (i.e. reliant) the results are. Here, two measures of robustness are used: the relative number of models that give an increase and the standard deviation.  

The relative number of models that give an increase is a measure of how large part of the ensemble that agrees on a positive change. The greater the share of models that point in the same direction, the more robust the result. If all models in the ensemble give a temperature increase that is, for example, a robust result. It is also a robust result if no models give an increase, i,e, they show a robust decrease. The result is not robust if, on the other hand, the models are relatively equally spread between increase and decrease. One example is shown in the left figure below. Dark green colours in northern Sweden mean that most models agree on an increase. However, in large parts of Southern Sweden, the models disagree on the sign of the change, except for the south-eastern parts where the models agree on a decrease.

Another measure of the robustness is the standard deviation. It shows how large the spread is between the models. The standard deviation is large when the spread is large, and small when the spread is small. The value of one standard deviation means that ca 68 % of the models are within that value. In the example displayed in the figure below, the spread is relatively small in most parts of northern Sweden, but substantial in the south. Several hydrological indicators are presented as relative changes. Standard deviations are presented despite the non-guaranteed normal distribution of such data. The results are most of the time not strongly affected by this, and the presented information is therefore deemed to be useful.

Två Sverikekartor där den vänstar visar robusthet och den högra standardavvikelse.
Caption: Robustness (left) is a measure of the relative number of models in the ensemble that give an increase, 100 % means that all models give an increase, 0 % means that no models give an increase. The standard deviation (right) shows the spread between members in the ensemble. Higher values mean larger spread.

Updates

Here, there is room for descriptions of updates made.

References

Johansson, B., Chen, D., 2003. The influence of wind and topography on precipitation distribution in Sweden: statistical analysis and modelling. Int. J. Climatol. 23 (12), 1523–1535.

Kjellström, E., Nikulin, G., Strandberg, G., Christensen, O. B., Jacob, D., Keuler, K., Lenderink, G., van Meijgaard, E., Schär, C., Somot, S., Sørland, S. L., Teichmann, C., and Vautard, R., 2018: European climate change at global mean temperature increases of 1.5 and 2 °C above pre-industrial conditions as simulated by the EURO-CORDEX regional climate models, Earth Syst. Dynam., 9, 459–478, 803 https://doi.org/10.5194/esd-9-459-2018.

Lindström, G., Pers, C.P., Rosberg, R., Strömqvist, J., and Arheimer, B. (2010): Development and test of the HYPE (Hydrological Predictions for the Environment) model – A water quality model for different spatial scales. Hydrology Research 41.3-4:295-319.

Moss, R. H. et al. 2010: The next generation of scenarios for climate change  research and assessment. Nature, Vol 463, 11 February 2010,  doi:10.1038/nature08823.

Rummukainen M., 2010: State-of-the-art with regional climate models. WIREs Clim Change 2010, 1:82–96, doi:10.1002/wcc.8.
Rummukainen, M. 2016: Added value in regional climate modeling, WIREs Clim Change 2016, 7:145–159. doi: 10.1002/wcc.378

Sørland, S. L., Schär, C., Lüthi, D. And Kjellström E., 2018: Bias patterns and climate change signals in GCM-RCM model chains, Environ. Res. Lett., 13, 074017, https://doi.org/10.1088/1748-9326/aacc7.

van Vuuren, D. P. et al., 2011: The representative concentration pathways: an overview, Climatic Change (2011) 109:5–31, DOI 10.1007/s10584-011-0148-z