In summer 2018, central and northern Europe were stricken by extreme drought and heat (DH2018). The DH2018 differed from previous events in being preceded by extreme spring warming and brightening, but moderate rainfall deficits, yet registering the fastest transition between wet winter conditions and extreme summer drought. Using 11 vegetation models, we show that spring conditions promoted increased vegetation growth, which, in turn, contributed to fast soil moisture depletion, amplifying the summer drought. We find regional asymmetries in summer ecosystem carbon fluxes: increased (reduced) sink in the northern (southern) areas affected by drought. These asymmetries can be explained by distinct legacy effects of spring growth and of water-use efficiency dynamics mediated by vegetation composition, rather than by distinct ecosystem responses to summer heat/drought. The asymmetries in carbon and water exchanges during spring and summer 2018 suggest that future land-management strategies could influence patterns of summer heat waves and droughts under long-term warming.

In Europe, three widespread extreme summer drought and heat (DH) events have occurred in 2003, 2010 and 2018. These events were comparable in magnitude but varied in their geographical distribution and biomes affected. In this study, we perform a comparative analysis of the impact of the DH events on ecosystem CO 2 fluxes over Europe based on an ensemble of 11 dynamic global vegetation models (DGVMs), and the observation-based FLUXCOM product. We find that all DH events were associated with decreases in net ecosystem productivity (NEP), but the gross summer flux anomalies differ between DGVMs and FLUXCOM. At the annual scale, FLUXCOM and DGVMs indicate close to neutral or above-average land CO 2 uptake in DH2003 and DH2018, due to increased productivity in spring and reduced respiration in autumn and winter compensating for less photosynthetic uptake in summer. Most DGVMs estimate lower gross primary production (GPP) sensitivity to soil moisture during extreme summers than FLUXCOM. Finally, we show that the different impacts of the DH events at continental-scale GPP are in part related to differences in vegetation composition of the regions affected and to regional compensating or offsetting effects from climate anomalies beyond the DH centres. This article is part of the theme issue ‘Impacts of the 2018 severe drought and heatwave in Europe: from site to continental scale’.

Abstract. In 2018 and 2019, central Europe was affected by two consecutive extreme dry and hot summers (DH18 and DH19). The DH18 event had severe impacts on ecosystems and likely affected vegetation activity in the subsequent year, for example through depletion of carbon reserves or damage from drought. Such legacies from drought and heat stress can further increase vegetation susceptibility to additional hazards. Temporally compound extremes such as DH18 and DH19 can, therefore, result in an amplification of impacts due to preconditioning effects of past disturbance legacies. Here, we evaluate how these two consecutive extreme summers impacted ecosystems in central Europe and how the vegetation responses to the first compound event (DH18) modulated the impacts of the second (DH19). To quantify changes in vegetation vulnerability to each compound event, we first train a set of statistical models for the period 2001–2017, which are then used to predict the impacts of DH18 and DH19 on enhanced vegetation index (EVI) anomalies from MODIS. These estimates correspond to expected EVI anomalies in DH18 and DH19 based on past sensitivity to climate. Large departures from the predicted values can indicate changes in vulnerability to dry and hot conditions and be used to identify modulating effects by vegetation activity and composition or other environmental factors on observed impacts. We find two regions in which the impacts of the two compound dry and hot (DH) events were significantly stronger than those expected based on previous climate–vegetation relationships. One region, largely dominated by grasslands and crops, showed much stronger impacts than expected in both DH events due to an amplification of their sensitivity to heat and drought, possibly linked to changing background CO2 and temperature conditions. A second region, dominated by forests and grasslands, showed browning from DH18 to DH19, even though dry and hot conditions were partly alleviated in 2019. This browning trajectory was mainly explained by the preconditioning role of DH18 on the impacts of DH19 due to interannual legacy effects and possibly by increased susceptibility to biotic disturbances, which are also promoted by warm conditions. Dry and hot summers are expected to become more frequent in the coming decades, posing a major threat to the stability of European forests. We show that state-of-the-art process-based models could not represent the decline in response to DH19 because they missed the interannual legacy effects from DH18 impacts. These gaps may result in an overestimation of the resilience and stability of temperate ecosystems in future model projections.

Abstract. Quantifying the net carbon flux from land use and land cover changes (fLULCC) is critical for understanding the global carbon cycle and, hence, to support climate change mitigation. However, large-scale fLULCC is not directly measurable and has to be inferred from models instead, such as semi-empirical bookkeeping models and process-based dynamic global vegetation models (DGVMs). By definition, fLULCC estimates are not directly comparable between these two different model types. As an important example, DGVM-based fLULCC in the annual global carbon budgets is estimated under transient environmental forcing and includes the so-called loss of additional sink capacity (LASC). The LASC results from the impact of environmental changes on land carbon storage potential of managed land compared to potential vegetation and accumulates over time, which is not captured in bookkeeping models. The fLULCC from transient DGVM simulations, thus, strongly depends on the timing of land use and land cover changes mainly because LASC accumulation is cut off at the end of the simulated period. To estimate the LASC, the fLULCC from pre-industrial DGVM simulations, which is independent of changing environmental conditions, can be used. Additionally, DGVMs using constant present-day environmental forcing enable an approximation of bookkeeping estimates. Here, we analyse these three DGVM-derived fLULCC estimations (under transient, pre-industrial, and present-day forcing) for 12 models within 18 regions and quantify their differences as well as climate- and CO2-induced components and compare them to bookkeeping estimates. Averaged across the models, we find a global fLULCC (under transient conditions) of 2.0±0.6 PgC yr−1 for 2009–2018, of which ∼40 % are attributable to the LASC (0.8±0.3 PgC yr−1). From 1850 onward, the fLULCC accumulated to 189±56 PgC with 40±15 PgC from the LASC. Around 1960, the accumulating nature of the LASC causes global transient fLULCC estimates to exceed estimates under present-day conditions, despite generally increased carbon stocks in the latter. Regional hotspots of high cumulative and annual LASC values are found in the USA, China, Brazil, equatorial Africa, and Southeast Asia, mainly due to deforestation for cropland. Distinct negative LASC estimates in Europe (early reforestation) and from 2000 onward in the Ukraine (recultivation of post-Soviet abandoned agricultural land), indicate that fLULCC estimates in these regions are lower in transient DGVM compared to bookkeeping approaches. Our study unravels the strong dependence of fLULCC estimates on the time a certain land use and land cover change event happened to occur and on the chosen time period for the forcing of environmental conditions in the underlying simulations. We argue for an approach that provides an accounting of the fLULCC that is more robust against these choices, for example by estimating a mean DGVM ensemble fLULCC and LASC for a defined reference period and homogeneous environmental changes (CO2 only).

Some of the major drivers of extreme temperatures in Australia are related to large-scale variations in the climate. While many of the relationships between climate variability and heatwaves have been described, the preceding climate conditions and the circulation patterns responsible for interannual heatwave variations are still poorly understood. We have used rotated principal component analysis (PCA) to describe the variance of frequency, duration, magnitude and timing of Australian summer heatwaves. This provides a simple means of investigating both the spatial patterns of circulation and the preceding climate mode conditions of summer heatwaves. We find that spatial patterns of heatwave variability associated with El Niño-Southern Oscillation (ENSO) were focused in the north and northeast of Australia for heatwave occurrence. Patterns of heatwave amplitude in southeastern Australia were also weakly related to ENSO. The PCA did not find any direct, concurrent Indian Ocean Dipole (IOD) or Southern Annular Mode (SAM) signal. However, lead correlations indicate that winter IOD conditions can have a weak indirect effect on heatwave occurrence in northeast Australia and similarly, winter SAM conditions can affect heatwave timing in central eastern Australia. Regressions reveal that ENSO related heatwave duration and frequency in the northeast are related to the movement of the Walker circulation, while southeastern heatwave amplitude is associated with a north-south movement of the Tasman high pressure system. These results may lead to improvements in understanding the predictability and general preparedness for heatwaves in the summer season.