Research highlights 2023–24
These scientific articles, that reflect the Bolin Centre different Research Themes, were originally published in the Bolin Centre’s annual report 2023–24.

RT1 | ARTofMELT – A Stockholm University led expedition to investigate the role of Atmospheric rivers and the onset of Arctic melt

The collective global warming effects are larger in the Arctic than elsewhere on Earth, and for many aspects observed trends fall well outside expectations, considering current understanding and modelling skills. Much of these deficiencies are due to a lack of observations needed to better our understanding. The lengthening of the sea ice melt season is one aspect, especially the so-called shoulder seasons; the onsets of melt and freeze. The onset of the melt is especially difficult to observe. Most Arctic research expeditions occur in August and September when the ice is much easier to navigate than at the end of winter, because of the summer melt.
ARTofMELT was designed to alleviate this, targeting two aspects: first the onset of the sea ice melt itself but also so-called ‘atmospheric rivers’ (ARs). ARs are filaments of warm and moist air propagating into the Arctic and are associated with weather systems. The hypothesis was that the melt onset timing is governed by such AR events. To test this, ARTofMELT had to navigate the Arctic in late winter ice and also be able to position measuring systems and sampling devices to where ARs would appear. This navigation was facilitated by a virtual forecast centre set up at MISU, Stockholm, and at the European Centre for Medium Range Weather forecasts in Reading, UK, who provided daily forecast guidance.

Group photo of all participants in ARTofMELT. Photo: Paul Zieger
After intense preparations, the Swedish research icebreaker Oden left Longyearbyen, Svalbard, on May 8, 2023, heading north with 38 ARTofMELT scientists onboard, representing 10 nations from 19 universities or research institutes, also supported by onshore scientists from another 3 institutes. The Swedish Polar Research Secretariat supplied 13 support staff while Oden’s crew comprised another 23, totalling 74 individuals going into the Arctic to experience spring together. An underlying strategy was to cover a column starting in the upper ocean, up through ice and snow, to the top of the troposphere. As you need to be nimble to catch ARs, the core of the observations had to be made from Oden, but the plan also included mooring to ice floes to take measurements on the ice when ARs were encountered. The program was dominated by atmospheric science (meteorology, atmospheric chemistry, and physics) but also included work packages on oceanography, biogeochemistry, sea ice, and snow, and geophysical mapping was also conducted. In addition to the ship itself and the ice camps, a helicopter was used to transport scientists to nearby ice floes to sample snow and ice and to deploy the Helipad, an airborne instrument package with an extensive suite of instruments. The Helikite, a tethered and instrumented balloon, was deployed from Oden’s aft, and while researchers were at the ice camps, an underwater robot was surveying the ice from below.

Research icebreaker Oden moored to an ice floe in the Fram Strait, to allow for work on the ice while receiving the first atmospheric river. Note the Helipad flying from the aft. Photo: Michael Tjernström
When the expedition ended in Longyearbyen on June 14, 2023 we had experienced only two AR, but the second triggered the onset of the melt, on June 10, 2023. AR episodes were fewer in number than expected but we also experienced two strong storms. The ice was much more difficult to navigate than expected and caused ARTofMELT to remain in the Fram Strait sea ice rather than situate itself farther north. Despite these setbacks, we covered relevant processes in all spheres from late winter through to the melt, collecting a wealth of unique observations and samples that will keep the science team occupied for years. The transition from winter to summer is more complex than first believed, and one main discovery was that it appeared under the ice before becoming obvious at the surface. We observed the onset of the spring algal bloom below the ice weeks before the surface reached melt; melt first occurred from below. While the surface melt onset was indeed triggered by an AR, the upper ocean and the ice was preconditioned by increasing solar radiation. A question for the future is how these timings are interconnected; what would have happened if the triggering AR had come earlier—or later?
ARTofMELT was co-lead by MISU and ACES.
RT1 | Code genealogy of climate models constructed
Almost ever since the advent of computers, climate models have been developed to simulate the Earth's climate. First, models were developed to scientifically understand the climate system and later to also understand climate change and predict future climate depending on greenhouse gas and aerosol emissions. Early climate models simulated only the atmosphere, later models added the ocean, and modern climate models, called ‘Earth system models’, also simulate the cryosphere; biosphere, and chemical processes in the atmosphere and ocean. Many countries started developing one or more state-funded climate models through research institutes. Over time, hundreds of such models have been developed, often borrowing computer code from one another. The immense complexity and size of the code often make starting a new model from scratch prohibitively expensive and time-consuming. This has led to a maze of model relationships, which until now have not been properly studied.
Many climate studies today rely on results coming from a set of models called a ‘model ensemble’; for example, by averaging the results of multiple models. This serves as protection from the deficiencies of any single model, but also allows scientists to quantify the uncertainty in their results and capture a greater amount of variability due to our uncertainty about climate processes and how to simulate them on a computer, with each model representing them differently. Including multiple models that share code in such an ensemble then leads to a certain amount of dependency, and certain groups of models sharing substantial amounts of code can be over-represented in the ensemble. Ideally, this should be known and quantified, and over-representing models should be avoided.
In our study, we constructed a code genealogy, or 'family tree', of 167 atmosphere-only, atmosphere–ocean, and Earth system models spanning several decades, including all models participating in the Coupled Model Intercomparison Project (CMIP) phases 3, 5, and 6. We have found that all of the studied models trace their history back to about 14 ancestral models, forming 14 model families. Within each family, models often tend to produce similar results, such as global temperature, so-called ‘climate feedbacks’, and sensitivity of temperature to greenhouse gas concentrations. We have found that only 3 model families comprise about 70% of all models in the latest CMIP phase.
We propose statistical weighting which takes into account the code relationships between models, so that, for example, more balanced averages can be calculated from model ensembles. This is in contrast to a simple 'model democracy', where every model in an ensemble is given the same statistical weight, as in ‘one model, one vote’, which would be only truly applicable if all models were completely independent.
When constructing the genealogy, our focus was on the component of models simulating the atmosphere, but a similar genealogy could be constructed with a focus on the ocean-simulating component. This was motivated by the fact that cloud processes are currently the leading source of uncertainty in future climate projections.
The projected sensitivity of temperature to greenhouse gases has increased on average in the latest generation of CMIP models. Our results contribute to explaining a part of this increase. Models in families with high sensitivity have proliferated more in the latest phase of CMIP compared to other model families and therefore have a greater weight in the overall result, even though they are often closely related. Our proposed weighting, accounting for dependencies, reduces the difference between the latest and the previous generation of CMIP models.
Future research could extend our analysis to the ocean component, quantify the strength of code relationships between models, or look at how the identified model families differ in their simulations of different aspects of the climate system. Overall, this could help explain why climate models do not always agree on their results, and thus reduce our uncertainty in future climate projections.
Our work was done at the Department of Meteorology and Bolin Centre for Climate Research at Stockholm University, and funded by the EU projects FORCeS and NextGEMS and the Swedish e-Science Research Centre. The results have been published in the Journal of Advances in Modelling Earth Systems (https://doi.org/10.1029/2022MS003588).

Figure from Kuma et al.(2023) illustrating the model code genealogy of models participating in the Coupled Model Intercomparison Project (CMIP) phase 3, 5, and 6, including their common ancestor models.
RT2 | Understanding the Impact of Remote Land-Use Changes on Wetland Hydrology
Wetlands are experiencing intense decline due to various factors, including climate change and human activities such as agricultural expansion, fragmentation by road infrastructure, water impoundment, and freshwater withdrawals for irrigation. Historically, research has focused on local and immediate threats to wetlands, such as pollution, direct land conversion, and water diversion. However, there is growing recognition that remote land-use changes, particularly those affecting upwind atmospheric moisture supply, can significantly influence wetland hydrologic regimes. Understanding these impacts is crucial, as they can alter precipitation patterns and moisture recycling, leading to changes in water availability and, consequently, wetland functionality.
A study by Fahrländer et al. (2024) evaluates the vulnerability of 40 wetlands protected by the Ramsar Convention and their hydrological basins to upwind precipitation shifts caused by land-use and hydroclimatic changes. The study investigates the precipitation source regions of these wetlands with a comprehensive approach that maps where precipitation predominantly comes from, its relationship to historical global land use change, and the hydroclimatic changes recently experienced by the wetlands. This approach allows us to assess how changes in upwind regions influence wetland moisture supply, which is a critical factor in maintaining their hydrologic balance and overall health Fahrländer et al. (2024) use a dataset of atmospheric moisture trajectories generated with the Lagrangian atmospheric moisture tracking model UTrack to trace the sources of atmospheric moisture contributing to precipitation in selected wetlands. They combine this model output with atmospheric reanalysis data, land use, and modelled evaporation and precipitation data to map the precipitation sources given their current land use and analyse changes in moisture supply through anthropogenic land use change. The results are then combined with trends in runoff, expressed as P–E, in a vulnerability indicator to assess the wetlands hydroclimatic vulnerability to upwind land use change.

Precipitation sheds (precipitation source regions) of the 40 Ramsar wetland basins that account for 70% of their total precipitation. Figure taken from Fahrländer et al. (2024) Figure 2.
They find that for 30 of the 40 wetlands, terrestrial sources contribute to at least half of the total annual precipitation. Notably, coastal wetland basins show the lowest terrestrial precipitation recycling ratios, primarily receiving ocean moisture. However, even these basins rely on land evaporation for 17%–23% of their annual precipitation. The findings reveal that historical land use changes have significantly reduced precipitation and moisture recycling in many regions, decreasing water availability in some wetlands. Regions such as Central Europe face substantial impacts to their hydrologic regimes due to high precipitation recycling ratios and significant land use changes.
The study also identifies ‘hotspot’ wetlands in South America and Central Asia which are particularly vulnerable. These regions exhibit a marked decrease in runoff and precipitation, driven by land-use-induced runoff impacts and climatic change. The declining surface water availability in these wetlands underscores the critical need to understand how climate change and land use changes impact ecosystems worldwide. The study highlights the significant influence of remote land use changes on wetland hydrology, emphasising the need to incorporate these effects into wetland restoration, management, and conservation strategies. Recognising the role of upwind land-use changes in shaping wetland hydrologic regimes is essential for developing effective protection strategies amidst evolving environmental conditions.
Moving forward, it is imperative to conduct more in-depth studies using high-resolution time-series data to better understand the dynamics of source-to-sink moisture flows. Such research will enhance our understanding of wetland moisture recycling and its broader implications for all water-dependent ecosystems.
Fahrländer, S. F., Wang-Erlandsson, L., Pranindita, A., & Jaramillo, F. (2024). Hydroclimatic vulnerability of wetlands to upwind land use changes. Earth's Future, 12, e2023EF003837. https://doi.org/10.1029/2023EF003837

Precipitation sheds (70% of P) of the eight identified “hotspot” wetland basins; ρterr : terrestrial precipitation recycling ratio; Figure taken from Fahrländer et al. (2024) – Supporting Information, Figure S5.
RT2 | Climate-active trace gas dynamics of CO2, CH4, and N2O in coastal ocean processes
Emissions of the trace gases carbon dioxide (CO2), methane (CH4), and nitrous oxide (N2O) from the coastal ocean, particularly from the shallow inshore areas, have generated strong interest regarding marine emissions. This because these coastal areas are highlighted as potential blue carbon sequestration areas, but are also among the most anthropogenically affected areas of the world thanks to human development and industrial use. RT2 has focused its research activities on a range of topics in the coastal environment by conducting local investigations in the Baltic Sea in collaboration with the Baltic Sea Centre and remote studies in the Arctic.
One of the highlights was the first-ever investigation of the shelf-wide assessment of N2O emissions from the Siberian shelf, looking at the relationship between N2O emissions and the land-derived discharge of nitrogen from large Russian rivers and coastal erosion in a rapidly warming region of the world, published by Wild et. al (2023). Our data shows that concentrations of dissolved nitrogen in the water and of total nitrogen in sediments increase toward large river mouths. On average, water-air N2O fluxes were low, but strong N2O sources and sinks were observed locally. An increase in water temperature could substantially reduce N2O solubility in ocean water, and add to biogeochemical and physical changes. This could alter N2O production in the shelf sediment, and convert the Siberian shelf to a N2O source in the future.
Coastal ecosystems can efficiently remove CO2 from the atmosphere and are therefore promoted for nature-based climate change mitigation. Natural CH4 emissions from these ecosystems may counterbalance atmospheric CO2 uptake. Still, knowledge of mechanisms sustaining such CH4 emissions and their contribution to net radiative forcing remains scarce for globally prevalent macroalgae, mixed vegetation, and surrounding depositional sediment habitats. A 2023 publication by Roth et. al (2023) showed that these habitats emit CH4 in the range of 0.1–2.9mg CH4m−2d−1 to the atmosphere. This was revealed by in situ CH4 emissions from macroalgae which were sustained by conventional methanogenic archaea, but also potentially associated with photosynthetic processes. The key finding was that over an annual cycle, CO2-equivalent CH4 emissions offset 28% and 35% of the carbon sink capacity attributed to atmospheric CO2 uptake in the macroalgae and mixed vegetation habitats offsetting the sink capacity of these potential blue carbon habitats. A multi-method, multi-instrument approach to tackle the complexity of trace gas emissions and carbon sequestration potential with an eye on coastal habitat diversity has risen to the forefront of RT2 investigations for the coming year.

From Roth et al., (2023). Seasonal net greenhouse gas balances. Daily mean net fluxes of CO2, CH4, and the net greenhouse gas balance (all expressed in CO2-eq. fluxes) (a), and the offset (in %) of the carbon sink capacity attributed to atmospheric CO2 uptake by CO2-eq CH4 emissions (b). Values in a are means ± standard error. Positive fluxes refer to an efflux from the water to the atmosphere (source), while negative fluxes depict an uptake of atmospheric GHGs (sink). CO2-equivalent CH4 fluxes were calculated using the sustained-flux global warming potential (SGWP) on a 100-year time horizon of 4521. The net greenhouse gas balance is calculated based on net CO2 and net CO2-eq. Original artwork by Elsa Wikander at Azote AB.
Ocean acidification, the other CO2 problem, continues to be a central research topic of Bolin Centre RT2 investigations, since increasing atmospheric CO2 drives ocean acidification globally. A recent review article by Gustafsson et al. (2023) highlights that in coastal seas, acidification trends can either be counteracted or enhanced depending on management strategies, and emphasises the importance of coastal monitoring studies to accompany model-based assessments. Ecosystem effects of acidification are so far small in the Baltic Sea, but changes should be anticipated unless CO2 emissions are curbed. Possible future acidification trends in the Baltic Sea, conditional on CO2 emissions, climate change, and changes in productivity can be assessed by means of model simulations. There are uncertainties regarding potential consequences for marine organisms, partly because of difficulties assigning critical thresholds, but also because of knowledge gaps regarding species’ capacity to adapt. Increased temporal and spatial monitoring of inorganic carbon system parameters would allow a better understanding of current acidification trends and improve capacity to predict possible future changes. An additional benefit is that such measurements also provide quantitative estimates of productivity. The technology required for precise measurements of the inorganic carbon system is readily available today. Regularly updated status evaluations of acidification, and of the inorganic carbon system in general, would support management when assessing climate change effects, eutrophication, or characteristics of the pelagic habitats. This would, however, have to be based on a spatially and temporally sufficient monitoring program.

Gustafsson et al. (2023). Simulated surface water pCO2 (upper panel) and pH (lower panel) in the Gotland Sea, using reconstructed forcing for the period 1851–1970 and climate change scenarios combined with nutrient load scenarios for the period 1971–2098. Thin pale lines indicate short-term variations while thick bright lines indicate annual means.
References
- Gustafsson, E., J. Carstensen, V. Fleming, B. G. Gustafsson, L. Hoikkala and G. Rehder (2023). Causes and consequences of acidification in the Baltic Sea: implications for monitoring and management. Scientific Reports 13(1): 16322.
- Roth, F., E. Broman, X. Sun, S. Bonaglia, F. Nascimento, J. Prytherch, V. Brüchert, M. Lundevall Zara, M. Brunberg, M. C. Geibel, C. Humborg and A. Norkko (2023). Methane emissions offset atmospheric carbon dioxide uptake in coastal macroalgae, mixed vegetation and sediment ecosystems. Nature Communications 14(1): 42.
- Wild, B., N. E. Ray, C. Lett, A. J. Davies, E. Kirillova, H. Holmstrand, E. Klevantceva, A. Osadchiev, I. Gangnus, E. Yakushev, D. Kosmach, O. Dudarev, Ö. Gustafsson, I. Semiletov and V. Brüchert (2023). Nitrous Oxide Dynamics in the Siberian Arctic Ocean and Vulnerability to Climate Change. Journal of Geophysical Research: Biogeosciences 128(5): e2022JG007326.
RT3 | Tracing the climatic and evolutionary history of Asian mammals: insights from dental traits and paleoclimate models
RT3 member Liping Liu and her team has explored how modern zoogeographic regions in Asia developed by examining the relationship between the dental traits of mammalian communities and climatic conditions. Using a method called redescription mining, Liu et. al (2023) have identified patterns of association between these traits and climate variables, providing insights into the history and distribution of mammal species in Asia.
Redescription mining methodology
Redescription mining is a key methodology used in this study. Redescriptions are pairs of logical statements about data variables that capture the interplay between the dental traits of mammals and their climatic conditions. Each logical statement, or query, selects a subset of localities where specific conditions are met, such as temperature or specific dental traits. The similarities between two subsets (one based on climate, the other on dental traits) is measured using the Jaccard coefficient, indicating the strength of their association.
Characterizing modern zoogeographic regions
The study identified nine significant redescriptions (patterns) that highlight the correlation between mammalian dental traits and climatic conditions across Asia, as seen in Figure 1. These redescriptions depict areas corresponding to prominent ecoregions and notable mammalian distribution patterns. This detailed characterisation helps us understand how specific regions are defined by both climate and the evolutionary traits of the mammals which inhabit them.
Emergence of modern zoogeographic patterns
Researchers traced the development of these patterns over five intervals of the Neogene period, which spans over the last 22 million years. They identified both north–south and southeast–northwest zoogeographic divisions in Asia. These patterns reflect changes in climate and tectonic activities, such as mountain ranges, which have significantly influenced the distribution and evolution of mammalian communities.
Impact of paleoclimate models
Various paleoclimate model simulations were used to understand the effect of different climate conditions on mammalian distribution. These models varied in terms of parameters like mountain elevation and climate data from different geological periods. By integrating these models, the researchers could analyse how past climate conditions impacted mammal distributions, and how these conditions have shaped current zoogeographic patterns.
Unique intergration of paleontology and climate science
One of the most interesting highlights of this study is its integration of paleontological data with climate science. The study combines data from fossil records with modern climate models, providing a comprehensive view of how mammalian communities have evolved in response to climatic changes. This interdisciplinary approach enhances our understanding of the historical context of modern zoogeographic regions.
Redescription mining as a tool
The adaptation of redescription mining, typically used in data science, to biological and environmental studies showcases its versatility in uncovering complex patterns in large datasets. This method allows researchers to identify and analyse associations between different variables, offering a powerful tool for exploring ecological and evolutionary dynamics.
Visual mapping of data
The study's results are visualised through maps, as seen in Figure 2, where different colours indicate whether the dental traits and climate conditions match the queries. This visual representation makes it easier to see the geographical distribution and overlap of these traits, providing a clear and intuitive understanding of the data.

Figure 1: Redescriptions rA-rl in the present-day dataset. Localities that satisfy both queries, only the dental traits query, only the climate query, and neither query, are drawn in purple, red, blue, and gray, respectively. For each redescription, we list the query over dental traits variables (qD), the query over bioclimatic variables (qC), the accuracy (J) as well as the size of its support as a percentage of the total number of localities (supp%)."
Insights into evolutionary biology
By linking dental traits, which evolve due to dietary needs, with climate conditions, the study provides insights into how environmental factors drive evolutionary changes in species. Understanding these links helps elucidate broader patterns of evolution and adaptation in response to changing climates over millions of years. This study sheds light on the historical development of mammalian communities in Asia through a novel approach that links climatic conditions with evolutionary traits. The findings highlight the dynamic nature of zoogeographic regions influenced by long-term climate changes and geological events. The use of redescription mining offers a powerful tool for exploring complex ecological and evolutionary patterns, contributing valuable knowledge to both palaeontology and climate science. This integrated approach provides a richer, more detailed picture of the evolution of biodiversity in response to environmental changes, offering important insights for future research in these fields.

Figure 2: Focus maps of redescriptions rB and rC (columns) evaluated on fossil localities from the different time intervals, considering the corresponding paleoclimate model simulation (rows).
This text is derived from the Liping Liu et. al (2023) publication The emergence of modern zoogeographic regions in Asia examined through climate–dental trait association patterns, published in Nature Communications, and can be found here: https://doi.org/10.1038/s41467-023-43807-w
RT3 | Understanding the effects of climate variability for historical grain harvest variations throughout Europe
Despite considerable differences between historical and modern agriculture, relationships between climate variability and harvest variations of past centuries can be informative when it comes to understanding general climatic risks to crop yields. However, such studies remain rare. A team of researchers representing several disciplines studied the climate signal in long harvest series from parts of Sweden, Switzerland, and Spain from the 16th to the 18th centuries. They found regionally consistent climate–harvest relationships of similar strength and character to modern ones for the same regions and crop types.
The understanding of associations between climate variability and grain harvest variations in early modern Europe (c. 1500–1800) has long been hampered by a lack of collaboration between historians and palaeoclimatologists and by a predominant focus on extreme climate events rather than long-term relationships. A group of historians and climate researchers, with experience in working on the climate–history interface collected a large dataset of already-published data covering grain harvests from the 16th and 18th centuries in central agricultural districts of east–central Sweden, northern Switzerland, and throughout Spain. The choice of regions was motivated by data availability as well as by a wish to obtain a north–south transect across Europe. State-of-the-art palaeoclimate reconstructions of seasonal temperature, precipitation, and soil moisture, based on tree-ring data and documentary sources, were used along with several statistical techniques to investigate the relationships between climate variability and grain harvest variations for barley, rye, oats and wheat. Different methods were used to investigate climate effects on inter-annual and longer timescales.
Grain harvests in east–central Sweden were found to have benefited mainly from wetter summers, in northern Switzerland from dry and warm winters, and in Spain from cooler springs and annual mean temperatures. The strengths of these relationships are comparable to modern climate–harvest relationships at comparable spatial scales. Spain exhibited particularly heterogeneous climate–harvest relationships due to its complex topography. For the summer season in Sweden and Switzerland, the relationship between temperature and harvest was found to be timescale dependent. Warmer summers reduced harvests (presumably through droughts, sub-optimal crop development rates, and plant pathogens) on inter-annual timescales, while this adverse effect of warm summers was not detected when considering longer timescales. Longer periods with warm summers instead increased harvests, at least in Sweden. Despite the northern location, the grain harvests from east–central Sweden (around Lake Mälaren and in Östergötland) were found to be more drought sensitive than those in Switzerland or even Spain. The reason for this is partly the low precipitation in spring and early summer in east–central Sweden, coinciding with the early part of the grain growing season, and the long hours of sunshine during this season. Grain in most of Spain, with its very hot and dry summers, has a key growing season from March to May, rendering the summer drought season largely irrelevant for crop harvests.
The study revealed that the relationship between climate variability and regional harvest variations during the 16th to 18th centuries were weaker than the relationship between climate variability and grain prices during the same period. This is because grain prices, through both trade and storage, represented an average of harvests from large parts of Europe for up to several years, leading to diffused impacts of local weather on harvests. Furthermore, climate variability only explains a portion of harvest variations. This is partly because poor harvests frequently resulted from single days or weeks of unfavourable weather conditions (and not seasonal climate averages). Harvest variations were partly due to factors related to seed quantity and quality, manure availability, access to draft animals, the availability of labour, and pest outbreaks.
While the strength and direction of the relationship between climate variability and grain harvest variations from the 16th to 18th centuries are surprisingly similar to those of today for the same regions and crops, the comparison between historical and modern agriculture is not without its limitations. The genetic variation of early modern grain types was much greater than it is today. Local varieties, suitable for local growth conditions, were used – particularly for barley and rye – complicating the comparisons of the relationship through time and between regions. More importantly, the longer time required to plant and harvest, and the limitations in technology used to dry the grain after harvest, posed other climatic hazards to historical grain production. Modern mechanisation, chemical fertilisers, pesticides, and higher-yielding seed varieties (which have resulted in multiplied harvest yields) have obviously changed the relationship between climate variability and grain harvest variations since pre-industrial times. Finally, while the relationships still remain similar to those of historical times, they are expected to change dramatically during the course of this century across Europe as a result of ongoing climatic warming and potential changes in precipitation patterns and their seasonality.

Swedish harvest correlations: Cross-correlation matrix for Sweden between climate series and grain tithe series at a county level for 10-year high-pass Gaussian filtered data (a) and line-arly detrended data (b), respectively. Correlations significant with a t test are marked with a dot (.).
The study Climatic signature in early modern European grain harvest yields was published in the scientific journal EGU Climate of the Past on 7 December 2023: https://doi.org/10.5194/cp-19-2463-2023
RT4 | The consequences of climate change for agroforestry in Arabica coffee's native range

Climate change is a major threat to agriculture, both directly by affecting crop growth and yield, but also indirectly by increasing pest and disease levels. Here, we present some key findings from a cross-disciplinary project involving meteorologists, pathologists, ecologists, and social scientists focused on coffee agroforestry in Arabica coffee's native range in southwestern Ethiopia, an area that also harbours the world’s only genetic reservoir of wild Arabica coffee genes.
Climate predictions
We used statistical downscaling to create predictions of the future climate of forest-covered areas in southwestern Ethiopia. To make these predictions easily accessible for stakeholders (e.g., policymakers, agricultural extension workers, and farmers), we created a visualisation tool on the Bolin Centre Database. Here, the user can select the climate model and emission scenarios, and the website shows dynamic maps, as well as the option to compare scenarios (https://bolin.su.se/data).
Climate-disease-yield relationships
We analysed the impact of past and present climates on disease dynamics, natural biocontrol, and yield, which allows us to learn from the past and present to predict the future. As one example, PhD student Biruk Ayalew focused on the fascinating interaction between the coffee leaf rust pathogen and its fungal hyperparasite. The fungal hyperparasite is widespread in the study area, and feeds on the coffee leaf rust pathogen. Biruk found that the coffee leaf rust pathogen and its hyperparasite differ in their climatic preferences, which suggests that there are opportunities to manage the microclimate in a way that increases natural biocontrol. A survey by a previous student at our department, Beyene Zewdie, indicated that the fungal hyperparasite controls disease levels under field conditions.
Perceptions of changes in climate and disease levels
As perceptions play an important role in adaptations to climate change, we interviewed smallholder farmers about their perceptions of changes in climate and disease levels. One interesting finding was that even though most scientific studies focus on increases in temperature, changes in precipitation and extreme weather events were at least equally important for farmers. Importantly, however, perceptions of climate change did not automatically translate into management adaptations.

An Arabica coffee leaf (green) with the coffee leaf rust pathogen (orange) and its fungal hyperparasite (white). Photo: Ayco Tack
Predicting changes at the landscape-level
In an ongoing part of the project, we are focusing on how future climatic changes and smallholder farmer perceptions will jointly reshape the landscape. At high elevations, we already observed how farmers are starting coffee production above the historic elevational range for growing coffee. We are currently examining how farmers at lower and middle elevations adapt to climatic changes.

A screenshot from the website where stakeholders can easily access the climate predictions, as hosted by the Bolin Climate Research Center database.
Taken together, climate-induced changes in coffee agroforestry will have major consequences for livelihoods, food security, forest cover, and biodiversity in southwestern Ethiopia.
References
- Ayalew, B., K. Hylander, G. Adugna, B. Zewdie and A. J. M. Tack. 2024. Impact of climate on a host-hyperparasite interaction on Arabica coffee in its native range. Journal of Applied Ecology 61: 538–550, https://doi.org/10.1111/1365-2664.14578
- Ayalew, B., K. Hylander, L. Börjeson, G. Adugna, D. Beche, F. Zignol and A. J. M. Tack. 2024. Using local knowledge to reconstruct climate-mediated changes in disease dynamics and yield – a case study on Arabica coffee in its native range. Plants, People, Planet. Accepted.
- Gomm, X., B. Ayalew, K. Hylander, F. Zignol, L. Börjeson and A. J. M. Tack. 2024. From climate perceptions to actions: case study on coffee farms in Ethiopia. Ambio, https://doi.org/10.1007/s13280-024-01990-0
RT4 | Soil drought weakens forest microclimate cooling
An investigation led by RT4 member Caroline Greiser, researcher at the Bolin Centre for Climate Research and the Physical Geography Department looked at the mechanisms which create cool microclimates beneath forest canopies during warm and dry summer days. The study reveals the importance of soil moisture in creating cooler forest microclimates compared to temperatures outside forests.

Caroline Greiser, leading author of the study, in the boreal forest in Sweden, where she also studies forest microclimate with the help of small temperature and soil moisture loggers. Photo: Marian Schönauer
Forests can buffer hot temperature extremes – a natural air conditioning effect. They do so by providing shade and by evaporating and transpiring water. One can imagine that forests sweat in the heat to keep their internal temperature low, and they need more water to be able to sweat. In other words, forests need soil water. Greiser and her team explored how soil moisture variability affects temperature offsets between outside and inside the forest on a daily basis, using data from temperate broadleaf forests in Central Europe spanning four subsequent summer seasons. They found that daily maximum temperatures inside forests are, on average, 2°C cooler than temperatures outside forests. Small tree seedlings, as well as much of the forest’s biodiversity, depend on these buffered forest microclimates.
Drier soils create a weaker cooling effect
The research team found that higher soil moisture levels improved the cooling effect in forests, emphasising the combined contributions of canopy shade, soil water evaporation, and plant transpiration to cooler microclimates.
“We used a network of small temperature and moisture loggers spread across different forest patches to link daily fluctuations of sub-canopy temperature to canopy cover and daily fluctuations of soil moisture at a given site,” says Caroline Greiser. She further adds: “Forest microclimate research often focuses on canopy cover as a major driver of understory cooling. Our study highlights the role of soil water in buffering understories from the rising heat.”
As climate change causes more disturbances to forest canopies and increases the risk of soil droughts, forests may lose their cooling function. The researchers therefore emphasise the significance of incorporating soil moisture into models predicting forest microclimate, biodiversity, and tree regeneration.

Tree canopies protect the forest understory from extremely hot temperatures. Photo: FelixMmittermeier @pexels.com
“The findings are alarming in the context of climate change as more frequent and more severe droughts may threaten the cooling functions of forests,” says lead author Caroline Greiser.
The study Higher soil moisture increases microclimate temperature buffering in temperate broadleaf forests was published in Agricultural and Forest Meteorology: doi.org/10.1016/j.agrformet.2023.109828
Last updated: 2025-11-10
Source: Bolin Centre for Climate Research