Introduction

Climate change will lead to widespread and diverse (positive and negative) impacts across human and environmental systems (IPCC 2014). The climate change experienced in a location will depend on global-scale (e.g. emissions of greenhouse gases—van Vuuren et al. 2011; van Vuuren and Carter 2014) to regional- and local-scale factors (such as interactions with orography and proximity to oceans). Consequently, the changes in climate projected across Europe vary between the oceanically influenced regions of northwest Europe to the continental climates of eastern and central Europe (Kovats et al. 2014; EEA 2018).

However, climate change will not occur in isolation and will be accompanied by socio-economic changes, such as in population, economies, societal and personal values, governance, and technology (O’Neill et al. 2017; Kok et al. 2018). The impacts of these highly uncertain socio-economic futures on key sectors such as water resources management and agriculture can be more significant than those due to climate change (Holman et al. 2017; Harrison et al. 2018).

The nations of Eastern and Central Europe have experienced major socio-economic changes over the past 30–40 years as they have transitioned to democratic nations and, in many cases, membership of the European Union. Nevertheless, economic convergence with western European countries is an ongoing process that will take many decades; so continuing socio-economic change is probable, meaning that understanding the combined impacts of future climate and socio-economic change is important for informing national adaptation policy.

We focus on the Czech Republic in central Europe, a region where future precipitation trends are uncertain (Kjellström et al. 2011; Kovats et al. 2014). Climate change in the Czech Republic is expected to lead to wetter winters and drier summers (Kyselý et al. 2011; Stepanek et al. 2016), more frequent extreme rainfall events (Kyselý and Beranová 2009; Kyselý et al. 2011; Rulfová et al. 2017), decreased spring and summer runoff (Hanel et al. 2012), increased flood damages (Rojas et al. 2013), increased pest pressures (Eitzinger et al. 2012), increased risk of land degradation (Trnka et al. 2016), changing crop suitability (Trnka et al. 2009; Machar et al. 2017) and crop yields (Trnka et al. 2004; Eitzinger et al. 2012), changing species suitability (e.g. Lučan et al. 2013; Kolář et al. 2017), and increased urban heatwave impacts (e.g. Lorencová et al. 2018).

Whilst there have been previous national climate change impact and vulnerability assessments in Europe (reviewed in EEA 2018), these have three key common limitations. Firstly, they tend to focus on the impacts of climate change alone, with quantitative scenarios of socio-economic change being used in less than half of the national assessments. Secondly, they have a sectoral or multi-sectoral (rather than inter-sectoral) perspective (EEA 2018) that neglects systemic interactions that can moderate or exacerbate impacts (Harrison et al. 2016). Finally, they tend to focus on the use of existing literature and expert workshop/interviews, with less than half of assessments taking a co-ordinated modelling approach. Therefore, the systematic and quantitative assessment of the cross-scale interactions or teleconnections that affect within-country impacts of climate change (e.g. due to agricultural trade, species migration, population (im)migration, transboundary river basin flows) is unlikely to be fully accounted for.

The aim of this paper therefore is to develop an improved understanding of the systemic impacts of climate change and socio-economic change on the Czech Republic, taking account of cross-sectoral interactions between agriculture, water, biodiversity, forestry, and land-use. Unlike previous studies, the trans-national interactions between the Czech Republic and the rest of Europe were explicitly considered. In doing so, we apply a spatially distributed European-scale integrated modelling tool (the IMPRESSIONS Integrated Assessment Platform2) but, for the first time, evaluate the impacts at a national, rather than continental, scale.

Methods

IMPRESSIONS Integrated Assessment Platform (IAP2)

The IMPRESSIONS Integrated Assessment Platform (IAP2), a further development of the widely published CLIMSAVE IAP (Harrison et al. 2013, 2015) to include the latest climate (dynamically downscaled CMIP5 climate model projections using Representative Concentration Pathways) and socio-economic scenarios (Shared Socio-economic Pathways), integrates a suite of sectoral models (water resources, crops, forestry, land-use allocation, species, floods) within a web-based platform. The IAP2 uses 10 arcmin × 10 arcmin (approximately 16 km × 16 km) grid cells across Europe (the EU28 plus Norway and Switzerland), although multiple land-uses are simulated within each cell. The models use average monthly to average annual input variables to simulate a wide range of impact and ecosystem service indicators. When simulating impacts, the models do not account for adaptation, with the exception of the land-use allocation model which includes autonomous adaptation to climate change due to changing crop yields. The linked models have been thoroughly validated and assessed through sensitivity (Kebede et al. 2015; Fronzek et al. 2018) and uncertainty analysis (Dunford et al. 2014; Brown et al. 2015). The model has been used in numerous climate change impacts and adaptation studies (Dunford et al. 2015; Harrison et al. 2016, 2018; Holman et al. 2017). A more detailed description of IAP2 is given in the Electronic Supplementary Material (ESM).

Exploratory scenarios of climate change and socio-economic change

The IAP2 includes three Representative Concentration Pathways (RCP 2.6, 4.5 and 8.5) and four European Shared Socio-economic Pathways (SSP 1, 3, 4, and 5) out to 2100. For each RCP, multiple CMIP5 global climate model (GCM) runs that had been dynamically downscaled using Regional Climate Models (RCMs) and are available within the Coordinated Regional Climate Downscaling experiment (CORDEX, Jacob et al. 2014) are used (Table S1 of ESM) to represent future climate model uncertainty in Europe (Madsen et al. 2016).

European versions (Kok et al. 2018) of four of the five global SSPs of O’Neill et al. (2017) capture the low and high challenges to both mitigation and adaptation. SSP2 (Middle of the Road) was omitted as a scenario lacking a distinct identity, as it is characterised by moderate changes on almost all its elements and intermediate challenges to both adaptation and mitigation (O’Neill et al. 2017; Kok et al. 2018). Quantified national values of population and GDP variables were provided from the global SSP database v1.0 (https://secure.iiasa.ac.at/web-apps/ene/SspDb/), whilst remaining scenario inputs were quantified using stakeholder estimates of change (See Table S2 and Table S3 of ESM and Pedde et al. (2018)). The European SSPs and their quantifications, which do not account for any adaptation to climate change, reflect contrasting but plausible futures to 2100:

  • SSP1—A high commitment to sustainable development through global cooperation results in reduced inequality and less resource intensive lifestyles. Investment in health and education and strong social cohesion enables communities to become more resilient to impacts and risks.

  • SSP3—Economic problems and increasing demand for resources lead to conflict between and within regions of Europe. Social break-down, declining living standards, and reducing long-term economic planning lead to a carbon intensive and highly unequal Europe.

  • SSP4—Global power increasingly belongs to a political and business elite who drive innovation and a high-tech green Europe alongside rising inequalities between the rich and poor.

  • SSP5—People place increasing faith in competitive markets, innovation, and participatory societies to produce rapid technological progress and development of human capital but lack environmental concern, as the development depends on fossil fuels.

Model runs and analyses

The IAP2 was run for:

  • Baseline (1 run): The baseline is represented by the CRU climatology for 1961–1990 (New et al. 2002) and socio-economic settings for 2000–2010.

  • Sensitivity analysis (110 runs): combinations of changes in annual average temperature (− 1 °C to + 11 °C in steps of 1 °C) and precipitation (− 60% to + 40% in steps of 10%) following (Fronzek et al. 2018).

  • Climate change scenarios alone (30 runs): RCP2.6, 4.5, and 8.5 for three to four GCM-RCMs and baseline socio-economic settings (Table S1) for the 2020s, 2050s, and 2080s time-slices.

  • Combined climate and socio-economic scenarios (99 runs): selected combinations of RCPs, GCM-RCMs, and SSPs (Table S1) for the 2020s, 2050s, and 2080s time-slices.

The simulated outputs for the Czech Republic for the baseline have been evaluated against CORINE land-cover (2000) and national land-cover statistics (2010), national crop areas (2010), and average annual crop yields for NUTS3 regions (for 2000–2017) from the Czech Statistical Office (2018). The sensitivity analysis runs were used to develop Impact Response Surfaces (IRS) for those grids within the Czech Republic, the larger ‘Continental’ bioclimatic classification region (Metzger et al. 2005) in which the Czech Republic lies, and the entire European region (Fig. S1). The IRS show the relative change in the selected IAP2 indicators relating to land-use, cropping, water, and biodiversity with respect to the 0 °C temperature change and 0% precipitation change simulation. These were evaluated to understand the relative sensitivity of the indicators within the Czech Republic to systematic changes in climate variables compared with the surrounding region and the continent. The IRS methodology allows us to conduct a “scenario-neutral” sensitivity analysis of the effects of temperature and precipitation changes. The methodology was not extended to socio-economic inputs, as these are examined in a scenario context in conjunction with the changing climate. Finally, we analysed simulated impacts under climate change–only and climate and socio-economic change to explore the changes within the selected output variables over time, by using multiple RCPs and climate models to account for the uncertainties in future emissions and climate model simulations.

Results

Evaluation of baseline performance in Czech Republic

The land-use classes simulated by the IAP2 were compared with the national land-cover data and the CORINE land-cover dataset for the baseline period (Fig. 1a). There is excellent agreement between the simulated and observed urban, agricultural, forest, and unmanaged land areas, although there are differences in the proportions of arable and grassland, with the IAP2 appearing to over-estimate grassland at the expense of arable areas. These differences can be attributed to the uncertainty in classifying grasslands of different types and management intensities. For example, the Czech national data only differentiate permanent grasslands (which may be extensively or intensively managed), with temporary or ley grasslands being included within their arable class (unlike the IAP2). Similarly, differences in class definition (Pérez-Hoyos et al. 2012), uncertainties inherent to satellite-based datasets (De Meij et al. 2018), and uncertainties relating to the derivation of ground-based and satellite-based products (Frolking et al. 1999) contribute to the differences with the CORINE estimates. The CORINE classes used for aggregation to IAP2 outputs can be found in Fig. S2 of ESM.

Fig. 1
figure 1

(a) Land-use classification for Integrated Assessment Platform 2 (IAP2) (baseline) and data from the Yearbook (2010) and Corine (2000) for the Czech Republic, (b,c) simulated and reported crop areas and yields for the four most extensive crops in the Czech Republic: (b) Observed (2010) and simulated (baseline) crop areas and (c) Boxplots of average annual crop yield [t/ha] within the NUTS3 areas. Each boxplot consists of 13 values (one per NUTS3 region), with the outliers of the sample shown as dots outside the boxplot. Observed values are the averages of the 2000–2017 period

Simulated crop areas and yields are compared with respective national data (Fig. 1b, c). We focus on four crops (wheat, rapeseed oil, barley, forage maize) that are the most widespread crops in the Czech Republic, and together cover about 50% of total agricultural area and about 70% of arable area. There is very good general agreement in the simulated and reported national crop areas (Fig. 1b). The reported median average annual yield for wheat and barley is within the inter-quartile range of the simulations among 13 NUTS regions (Fig. 1c). For oilseed rape and forage maize, although the above condition is not met, the model performance is still comparable with the reported yields, with ~ 15% underestimation and ~ 33% overestimation, respectively. Nevertheless, it is apparent that the range of reported average annual yields across the 13 NUTS regions in the Czech Republic is much narrower than for the baseline simulation, indicating greater spatial variability in the modelled crop yields (Fig. 1c). The values of simulated and observed crop yield for each of the 13 NUTS regions are given in Tables S4 and S5 of ESM respectively.

Sensitivity analysis

The results of the sensitivity analysis are shown as scenario-neutral IRS for a number of key outputs on crop yield, water use, biodiversity, and land-use for the Czech Republic, the Continental region, and Europe (Fig. 2). For the Czech Republic, simulated wheat yield decreases with temperature increases and all variations in precipitation. The IRS is similar to that of the Continental region, but shows generally larger yield decreases at high temperature increases than Europe, reflecting the reduced climatological range (and therefore lower variability in yield response to changes) at the smaller scales. The area of arable land generally increases with increasing temperature, which might seem counterintuitive given the decreasing yields and therefore gross margin (as costs remain constant). However, as yields decline in the Czech Republic and Europe, the area of arable land in Europe has to expand in order to meet the constant food demand. The larger arable area increases in the Czech Republic of up to 40% relative to the baseline compared with the Continental region and Europe as a whole indicate that the Czech Republic (despite yield decreases in crops such as wheat) has a relative competitive advantage under these climate conditions.

Fig. 2
figure 2

Impact response surface (IRS) analysis for (a) the Czech Republic, (b) the Continental region, and (c) the European domain. Precipitation changes are shown on the y axis and temperature changes on the x axis of each plot. Contour values refer to the percentage change in the indicator relative to the baseline

The changes in water availability and the Water Exploitation Index (WEI) (Fig. S3 of ESM) are highly sensitive to changes in precipitation with little sensitivity to temperature changes, due to the strong soil water control on actual evapotranspiration. Irrigation usage (Fig. 2) strongly increases for all temperature increases and with decreased precipitation in the Czech Republic reflecting the increased yields and gross margins of some irrigated crops compared with rainfed crops under conditions of increasing aridity combined with sufficient water resources to support increased irrigation. These increases are greater than in the Continental region and Europe as a whole reflecting the increased arable area and greater aridity. There is a gradual reduction in the number of species (Fig. 2) as temperature increases and precipitation decreases leading to major species loss for extreme levels of warming (> 5 °C) and drying (< − 50%) of more than half of simulated number of species. This arises from the combination of substantial changes in suitable climate and the dramatic land-use changes associated with agricultural expansion leading to reduced habitat space for non-agricultural (forest, montane, etc.) species. The proportional loss of species in the Czech Republic is higher than for the Continental region and Europe, reflecting the smaller area and reduced bioclimatic range.

Climate change impacts

Climate change impacts for the Czech Republic, without the influence of socio-economic changes, are shown in Fig. 3 for key IAP2 outputs (Table S6). Results are presented for each time period and RCP, and the uncertainty in the projections is shown by including all the GCM-RCM models of each ensemble (Table S1).

Fig. 3
figure 3

Crop, water, biodiversity, and land-use simulations from Integrated Assessment Platform 2 (IAP), for the baseline and future time periods, aggregated for the Czech Republic. Results are shown for each representative concentration pathway (RCP), with baseline socio-economics. Per time period, the dots represent the different models of the ensemble used in each RCP simulation

Generally, climate change impacts are, as expected, lower for the lower-end emission scenario (RCP2.6) and become more pronounced over time within the same emission scenario. Average rainfed wheat yields generally increase in the Czech Republic (Fig. 3) and Continental region (Fig. S4 of ESM), but decrease for Europe (Fig. S5 of ESM). However, there is uncertainty in the direction of change for many time periods in the Czech Republic due to climate model uncertainty. Average potato yields decrease for the projections of all RCPs and time periods to the extent that potatoes are no longer selected for some time period-RCP-climate model combinations (as shown by negligible or zero average yield). This reflects the loss of competitiveness of the Czech Republic for profitable potato cultivation, as potato yields remain largely constant in the Continental region (Fig. S4) and Europe (Fig. S5).

Irrigation usage increases as a consequence of both increased irrigation demand (caused by the warmer drier summer periods and increased arable areas) and the increased water availability which allows greater volumes to be allocated to agriculture (having allowed for the public, industrial, energy generation, and environmental water needs). The increased water demands are balanced by the increased water availability so that there are little changes in WEI over time in the Czech Republic.

The number of species with suitable climate and habitat space decreases slightly under RCP2.6, but further progressive decreases are projected over time with increasing emissions reflecting decreasing climate suitability for many species and future declines in some land-use classes within the Czech Republic, therefore restricting the habitat of some species. In particular, whilst there are increases in the area of intensive arable land and changes in the relative proportions of unmanaged and managed forests (increasing and decreasing, respectively), there are major losses of grassland—to the point of virtually disappearing in the 2080s under RCP8.5. Both intensive and extensive grassland decrease strongly with increasing climate change and time in the Czech Republic and Continental Region, but intensive grassland only slightly reduces at the European level. As the IAP2 aims to meet net European food demand under climate change by expanding/contracting production areas based on profitability levels, arable agriculture in the Czech Republic and Continental region generally becomes more profitable than dairy and livestock farming leading to land-use conversion, whereas dairying on intensive grassland maintains its profitability in other regions of Europe.

Climate and socio-economic change impacts

The response of key outputs to combined climate and socio-economic scenarios for the end of the twenty-first century time-slice (2080s), reported per RCP for the Czech Republic, is shown in Fig. 4. All the GCM-RCM models of each ensemble are shown to depict the range of uncertainty in the projections. It is apparent that the inclusion of the socio-economic scenarios leads to a greater diversity of outcomes for the 2080s than climate change alone.

Fig. 4
figure 4

Crop, water, biodiversity, and land-use simulations from Integrated Assessment Platform 2 (IAP), for the baseline (dashed line) and 2080s future time period (coloured dots), aggregated for the Czech Republic. For the 2080s, results are shown per representative concentration pathway (RCP), for each shared socio-economic pathway (SSP)xRCP scenario run. The dots of the same colour represent the different global climate model (GCM)-regional climate model (RCM) combinations of the ensemble used in each RCP simulation

First, looking at the agricultural indicators, wheat yields increase strongly within SSP4 and SSP5, which are both scenarios of strong technological and agronomic innovation, in contrast to the decreases in SSP3, where social and political break-down leads to productivity losses. Despite the innovation in SSP4 and SSP5, potato cultivation only remains profitable in the Czech Republic under SSP1, a scenario of increased European domestic production as the continent reduces food imports in order to reduce its external environmental burdens. This is the only SSP that leads to strong increases in the arable area in Europe in the 2080s (Fig. S8). Alternatively, in the Czech Republic, the arable area increases for at least one climate model for all RCP-SSP scenarios indicating that the country’s arable agriculture is a relative winner under most futures. As a consequence of this expansion in the arable area in the Czech Republic under the SSPs, there are large increases in irrigation usage compared with the projections forced with baseline socio-economics. The largest increases occur in SSP1 (under RCP2.6 and 4.5) due to the need to increase domestic food production and in SSP4 and 5 (under RCP8.5) where technological innovation drives down the cost of the irrigation, thereby making irrigated agriculture increasingly profitable compared with rainfed agriculture under the high-end climate change scenario.

Similarly to the climate change–only scenarios (Fig. 3), grassland agriculture declines under the SSPs except in SSP1 (Fig. 4) where the reduced dietary preference for red and white meat (Table S2) leads to increased demand for dairy products and consequently little change in the intensive grassland area. However, the expansion of arable areas and the increased demand for dairy products and therefore increased profitability of intensive grassland in SSP1 lead to a different spatial pattern of land-use changes compared with the other scenarios (Fig. S6).

In SSP1-RCP2.6, arable agriculture expands into the grassland areas and low elevation forest (Fig. S6). Consequently, the grassland areas displace much of the forest in the higher elevation areas near the northern, eastern, and south-western borders. In SSP3 and SSP5, the expansion of the arable area (Fig. 4) is similarly at the expense of the grassland areas and low elevation forest (Fig. 4 and Fig. S6), but the low profitability of Czech livestock farming in these scenarios compared with elsewhere in Europe means that the livestock farming is displaced from the country and forest remains. These large projected future changes in land-use in the Czech Republic due to the changing relative profitability of arable and grassland systems, relative to each other and to other parts of Europe, will have significant consequences for rural (agricultural and forestry) livelihoods and the country’s landscape.

WEI shows only marginal changes under all RCPs combined with baseline socio-economics but increases to levels indicating severe water scarcity (WEI > 0.4) for the Czech Republic under all socio-economic scenarios, although the underlying reasons for the increases differ between scenarios. For example, SSP5 is associated with the largest increase in WEI as a combination of large population and GDP increases, low behavioural interest in water saving, and low environmental regulation of abstraction leading to large water abstraction increases. In SSP3, the effects of the decreases in population and GDP on water demand are offset by low behavioural and technological water saving and reduced irrigation efficiency. SSP4 has the smallest impact on increasing WEI, as high levels of technological water savings offset much of the increased water demands. Similar changes in WEI are seen in the Continental region (Fig. S7 of ESM) and also Europe (Fig. S8 of ESM), albeit to a lower level due to the moderating influence of the wetter northern and western regions.

The number of simulated species present in the Czech Republic is projected to decrease in all three climate change scenarios, with the impacts and uncertainty getting more pronounced for the higher-end scenarios due to greater changes in a species’ climate and habitat space. The relative consistency across the SSPs for the same RCP reflects the common major changes in land-use across the SSPs (expansion of intensive arable and reduction/loss of grassland systems) despite SSP differences in the proportions of arable and forest land-uses (e.g. lower arable and higher forest areas in SSP4, compared with the other SSPs in Fig. 4).

Discussion

We have applied a trans-European integrated modelling platform (the IMPRESSIONS IAP2) to understand the consequences of climate and socio-economic change for the Czech Republic. First, we have shown that a range of outputs from the IAP2 can be acceptably validated at the national-scale against multiple indicators (land-use, crop yield, and crop area). Second, we have carried out a novel national-scale impact assessment that takes account of the consequences of cross-sectoral interactions at national to continental scales. Previous regional-scale applications of IAP2 (Clarke et al. 2017) did not focus on the inter-scale interactions of reported impacts. Moreover, previous national-scale climate change impact assessments for the Czech Republic have focused on single ‘sectors’ (e.g. hydrology by Dvorak et al. 1997; Dubrovsky et al. 2009; Hanel et al. 2012; agriculture by Trnka et al. 2004, 2009; Lorencová et al. 2013), without accounting for the impact flows between sectors and regions. Consequently, modelled impacts within a country can be over- or under-estimated if the systemic footprint of climate change, in terms of feedbacks and interactions with other sectors (Harrison et al. 2016) and cross-scale dependencies within the wider geographical area, is not taken into account. To the authors’ knowledge, this is the first time an integrated cross-sectoral assessment of climate change and socio-economic impacts for the Czech Republic, taking account of interactions with the rest of Europe, has been performed.

IRSs have been previously used to understand model response to simple system perturbations (e.g. Prudhomme et al. 2010; Remesan and Holman 2015; Fronzek et al. 2018). Here we use them to provide a first indication of the scale interdependencies of the modelled impacts under precipitation and temperature perturbations, by comparing national-, regional-, and continental-scale responses. Subtle but systematic scale-dependent differences were found between some indicator responses in the Czech Republic, Continental region, and Europe. These differences included slightly better response of average crop yields in the Czech Republic and Continental region compared with Europe and greater increases in the benefits of irrigated compared with rainfed cropping for farm profitability in the Czech Republic compared with both the Continental region and Europe (as seen in the much greater increase in irrigation usage). Such systematic scale-dependent differences meant that arable agriculture in the Czech Republic was projected to have a comparative profitability advantage over Europe as a whole across much of the IRS, although changing interannual variability which can be crucial for agricultural production is not addressed by the IAP2. Consequently, whilst the overall arable area in Europe was little changed across the IRS, the simulated arable area in the Czech Republic expanded by up to more than 20% at the expense of production in other parts of Europe (Fronzek et al. 2018). A partial consequence of these scale interdependencies in land-use changes is that the scenario-neutral climate changes lead to more pronounced negative impacts on biodiversity in the Czech Republic.

However, the scenario-neutral analysis does not show the full picture of scale interdependencies as these largely depend on socio-economic factors such as food imports and agricultural demands. Model projections based on combined climate and socio-economic scenarios lead to a far greater diversity of outcomes for the Czech Republic for the 2080s than those associated with climate change scenarios alone. Future changes in the socio-economic conditions of countries and continents will affect all aspects of human systems, e.g. levels of agronomic innovation, food demand (through changing population numbers, food imports, and dietary preferences), water demand (through technology uptake and personal behaviour), and regulation (through spatial planning and environmental protection). The impact of these factors on spatial patterns of European food productivity, food and water demands, and environmental constraints across Europe leads to complex spatially varying changes in the model outputs for the Czech Republic.

The projected changes that are largely consistent in their direction, if not magnitude, of change across most scenario combinations for the Czech Republic are an expansion of arable land, reduction or loss of grasslands, reduction in area of managed forests and increase in the area of unmanaged forest, water scarcity, and a reduction in the number of species. Some of these impacts would be expected from the direct impacts of climate and socio-economic change acting solely within the Czech Republic; however, the consideration of the country level interactions within the greater geographical area can reveal different patterns for some indicators. Koutroulis et al. (2018), in a pan-European assessment of water availability with no spatial interactions, also found increases in mean water availability for the Czech Republic under climate and socio-economic change. They also found decreases in freshwater vulnerability, indicating less water scarcity at the country level. However, we find that the country moves towards extreme water scarce conditions due largely to increases in non-agricultural water use (driven by socio-economic change), leading to greatly constrained increases in irrigation usage and a projected expansion of rainfed arable land in the Czech Republic. On the land-use aspect, Trnka et al. (2009) also report declines in grasslands under climate change; however, taking into account the growing importance of the Czech Republic as a crop-growing region for Europe, our results show far more dramatic declines for grassland areas.

Our findings can provide valuable insights to support adaptation actions and policy measures in the Czech Republic and more broadly. The significant changes in simulated land-use in response to future climate and socio-economic change (e.g. expansion of arable farming, reduction or loss of grassland systems, and the reduction in managed forestry) imply significant changes to the rural economy, landscape heterogeneity and character, water quality, and habitat availability and connectivity. In such circumstances, future policy will need to support measures for a sustainable transition (Blanco et al. 2017; Brown et al. 2018) that will also allow the protection of valuable ecosystem services, such as measures to prevent/minimise soil erosion, preserve key characteristics of the landscape that are valued for tourism, and to support biodiversity. Such measures include establishing protected areas for key habitats and biodiversity corridors to support migration of species within and into the Czech Republic (Keeley et al. 2018). However, to develop effective policies, more locally specific considerations—which are not covered by the existing modelling platform—should be taken into account. For example, the projected conversion of forest to arable land could, in reality, be limited by land ownership conflicts and/or local topography of the forest terrain.

Our analysis shows strong signals of change due to climate and socio-economic change, with the land-water-biodiversity system in the Czech Republic being particularly sensitive to socio-economic factors, reinforcing the findings of Harrison et al. (2018) and Fronzek et al. (2018). Socio-economic change therefore becomes the most important source of uncertainty in the future projections due to its dominant influence on all aspects of demand for, production of, trade-in and regulation of natural resources. It is imperative therefore that future impact assessments avoid either neglecting socio-future economic change and its uncertainties or assuming some form of Business-As-Usual trends. Immediate policy priorities should aim to transition the Czech and European societies and economies towards more sustainable futures (Hölscher et al. 2018; Tàbara et al. 2018) that reduce the systemic non-climate pressures from unsustainable behaviours, thereby increasing the effectiveness of adaptation actions to reduce climate impacts and support a positive change in the system (Dunford et al. 2015).

Conclusion

We used a novel approach to study the impacts of climate and socio-economic change at the national level of the Czech Republic, accounting for trans-national and cross-sectoral interactions. The major impacts, consistent across most scenarios, are the expansion of arable land, reduction or loss of grasslands and managed forests, increased water scarcity, and loss of species. These impacts are partially driven by the increased importance of the Czech Republic as a crop-growing region for Europe, due to an increasing competitive advantage within the continent.

Our analysis highlights the importance of considering wider geographical domain interactions in integrated national-scale assessments of climate and socio-economic change, as they may have considerable impacts on national-level results. Moreover, such an approach can reveal the changing role of a country within its wider geographical and economic area, which is important to inform relevant measures for adaptation. Finally, the importance of socio-economic drivers highlights the opportunity for adaptation actions to alleviate negative impacts and provide support for climate resilient and sustainable development.