Benefits of simulating precipitation characteristics over Africa with a regionally-coupled atmosphere–ocean model

High-quality climate information at appropriate spatial and temporal resolution is essential to develop and provide tailored climate services for Africa. A common method to produce regional climate change data is to dynamically downscale global climate projections by means of regional climate models (RCMs). Deficiencies in the representation of the sea surface temperatures (SSTs) in earth system models (ESMs) and missing atmosphere–ocean interactions in RCMs contribute to the precipitation bias. This study analyzes the influence of the regional atmosphere–ocean coupling on simulated precipitation and its characteristics over Africa, and identifies those regions providing an added value using the regionally coupled atmosphere–ocean model ROM. For the analysis, the MPI-ESM-LR historical simulation and emission scenario RCP8.5 were dynamically downscaled with ROM at a spatial resolution of 0.22° × 0.22° for the whole African continent, including the tropical Atlantic and the Southwest Indian Ocean. The results show that reduced SST warm biases in both oceans lead to more realistic simulated precipitation over most coastal regions of Sub-Saharan Africa and over southern Africa to varying degrees depending on the season. In particular, the annual precipitation cycles over the coastal regions of the Atlantic Ocean are closer to observations. Moreover, total precipitation and extreme precipitation indices in the coupled historical simulation are significantly lower and more realistic compared to observations over the majority of the analyzed sub-regions. Finally, atmosphere–ocean coupling can amplify or attenuate climate change signals from precipitation indices or even change their sign in a regional climate projection

East Asian summer precipitation in AWI‐CM3: Comparison with observations and CMIP6 models

Owing to the complicated spatial–temporal characteristics of East Asian precipitation (EAP), climate models have limited skills in simulating the modern Asian climate. This consequently leads to large uncertainties in simulations of the past EAP variation and future projections. Here, we explore the performance of the newly developed Alfred Wegener Institute Climate Model,version 3 (AWI-CM3) in simulating the climatological summer EAP. To test whether the model's skill depends on its atmosphere resolution, we design two AWI-CM3 simulations with different horizontal resolutions. The result shows that both simulations have acceptable performance in simulating the summer mean EAP, generally better than the majority of individual models participating in the Coupled Model Intercomparison Project (CMIP6). However, for the monthly EAP from June to August, AWI-CM3 exhibits a decayed skill, which is due to the subseasonal movement of the western Pacific subtropical high bias. The higher-resolution AWI-CM3 simulation shows an overall improvement relative to the one performed at a relatively lower resolution in all aspects taken into account regarding the EAP. We conclude that AWI-CM3 is a suitable tool for exploring the EAP for the observational period. Having verified the model's skill for modern climate, we suggest employing the AWI-CM3, especially with high atmosphere resolution, both for applications in paleoclimate studies and future projections

Multi-resolution climate modelling with the AWI Climate Model (AWI-CM)

The recently established AWI Climate Model (AWI-CM), a coupled configuration of the Finite Element Sea Ice-Ocean Model (FESOM) with the atmospheric model ECHAM6, uses a novel multi-resolution approach: Its ocean component builds on a finite element dynamical core supporting unstructured triangular surface grids, allowing to distribute the grid points in a flexible manner. This allows to concentrate resolution in dynamically important regions, with a continuous transition zone to the coarser resolution in other areas. The model is an ideal tool to study the influence of explicit resolution of smaller scales in dedicated experiments. The unique – spatially seamless – approach might also be of benefit when it comes to temporally seamless prediction, bridging the gap between numerical weather prediction and climate models. A first benchmark set-up of AWI-CM with moderate resolution in the atmosphere (T63) and 25km in key ocean areas, e.g. around the equator, achieved a similar overall simulation performance in a long control simulation compared to well-established CMIP5 models. In particular, the (isotropically) increased equatorial resolution considerably increased the realism of TIW activity and ENSO-related variability compared to standard resolutions. The potential of AWI-CM is further exploited within the EU project PRIMAVERA in the HighResMIP of CMIP6, where we plan to contribute simulations with eddy-resolving resolutions (1/12° or 9-10 km) in key areas of the global ocean, such as the Gulf Stream-North Atlantic Current region, the Agulhas retroflection zone, or the Arctic basin. First simulations show distinct improvements with respect to the development of deep temperature and salinity biases in the North Atlantic Ocean and an overall improvement of surface biases. At even higher resolutions of 4.5 km locally in the Arctic, linear kinematic features emerge in the simulated sea ice distribution with potentially strong impacts on air-sea fluxes in the coupled system. Although the tested set-ups are computationally very demanding (with numbers of grid points comparable to a regular 0.25° grid), the throughput is high at about 8 simulated years per day because of high scalability. In addition, we are about to finish the development of a finite volume version of the ocean model code (FESOM 2). It is already faster than the original FESOM version by a factor of two to three, which will further enlarge the set of computationally feasible applications

Designing variable ocean model resolution based on the observed ocean variability

If unstructured meshes are refined to locally represent eddy dynamics in ocean circulation models, a practical question arises on how to vary the resolution and where to deploy the refinement. We propose to use the observed sea surface height variability as the refinement criterion. We explore the utility of this method (i) in a suite of idealized experiments simulating a wind-driven double gyre flow in a stratified circular basin and (ii) in simulations of global ocean circulation performed with FESOM. Two practical approaches of mesh refinement are compared. In the first approach the uniform refinement is confined within the areas where the observed variability exceeds a given threshold. In the second one the refinement varies linearly following the observed variability. The resolution is fixed in time. For the double gyre case it is shown that the variability obtained in a high-resolution reference run can be well captured on variable-resolution meshes if they are refined where the variability is high and additionally upstream the jet separation point. The second approach of mesh refinement proves to be more beneficial in terms of improvement downstream the midlatitude jet. Similarly, in global ocean simulations the mesh refinement based on the observed variability helps the model to simulate high variability at correct locations. The refinement also leads to a reduced bias in the upper-ocean temperatur

Delayed Antarctic sea-ice decline in high-resolution climate change simulations

Despite global warming and Arctic sea-ice loss, on average the Antarctic sea-ice extent has not declined since 1979 when satellite data became available. In contrast, climate model simulations tend to exhibit strong negative sea-ice trends for the same period. This Antarctic sea-ice paradox leads to low confidence in 21st-century sea-ice projections. Here we present multi-resolution climate change projections that account for Southern Ocean mesoscale eddies. The high-resolution configuration simulates stable September Antarctic sea-ice extent that is not projected to decline until the mid-21st century. We argue that one reason for this finding is a more realistic ocean circulation that increases the equatorward heat transport response to global warming. As a result, the ocean becomes more efficient at moderating the anthropogenic warming around Antarctica and hence at delaying sea-ice decline. Our study suggests that explicitly simulating Southern Ocean eddies is necessary for providing Antarctic sea-ice projections with higher confidence

Impact of ocean-atmosphere coupling on present and future Köppen-Geiger climate classification in Europe

The effect of air-sea coupling in the simulation of the European climate is assessed through a climate type classification that uses surface temperature and precipitation from a regional atmosphere-ocean coupled model and from its atmospheric component. The atmospheric setup in both models is the same, differing only in the representation of the oceanic fields. The simulations cover the present and future-time climate under the RCP8.5 CMIP5 scenario. Climate type distributions obtained from both coupled and uncoupled simulations are similar to those obtained from ERA5 for the 1976–2005 period. Both models simulate colder climate types for present-time in southern and northeastern regions compared to ERA5, possibly due to a weaker influence of the Atlantic circulation, and drier climate types in some western Mediterranean areas. Also, for present-time coupling leads to more humid winters (relatively drier summers) in some zones of north Spain and south France, and drier climates in some western Mediterranean spots. Based on simulations with these models under the RCP8.5 scenario, we find that by the end of the 21st century (2070–2099) the climate type distribution changes in more than 50% of the domain. While both models project the reduction of regions with cold climate types and the expansion of those with hot summers and hot arid climate types, these changes affect a larger area in the coupled simulation. These differences may be related to a drier signal in the coupled simulation, especially during summer, due to the influence of colder surface water in the North Atlantic Ocean and the Mediterranean Sea. Using a climate classification to evaluate the annual cycles of the simulated temperature and precipitation data provides a novel insight into the impact of air-ocean coupling on the representation of the climate, and consequently into the simulated impact on ecosystems and human activities in Europe

Atmospheric rivers in CMIP5 climate ensembles downscaled with a high-resolution regional climate model

Atmospheric rivers (ARs) are important drivers of hazardous precipitation levels and are often associated with intense floods. So far, the response of ARs to climate change in Europe has been investigated using global climate models within the CMIP5 framework. However, the spatial resolution of those models (1–3∘) is too coarse for an adequate assessment of local to regional precipitation patterns. Using a regional climate model with 0.22∘ resolution, we downscaled an ensemble consisting of 1 ERA-Interim (ERAI) reanalysis data hindcast simulation, 9 global historical, and 24 climate scenario simulations following greenhouse gas emission scenarios RCP2.6, RCP4.5, and RCP8.5. The performance of the climate model to simulate AR frequencies and AR-induced precipitation was tested against ERAI. Overall, we find a good agreement between the downscaled CMIP5 historical simulations and ERAI. However, the downscaled simulations better represented small-scale spatial characteristics. This was most evident over the terrain of the Iberian Peninsula, where the AR-induced precipitation pattern clearly reflected prominent east–west topographical elements, resulting in zonal bands of high and low AR impact. Over central Europe, the models simulated a smaller propagation distance of ARs toward eastern Europe than obtained using the ERAI data. Our models showed that ARs in a future warmer climate will be more frequent and more intense, especially in the higher-emission scenarios (RCP4.5, RCP8.5). However, assuming low emissions (RCP2.6), the related changes can be mostly mitigated. According to the high-emission scenario RCP8.5, AR-induced precipitation will increase by 20 %–40 % in western central Europe, whereas mean precipitation rates increase by a maximum of only 12 %. Over the Iberian Peninsula, AR-induced precipitation will slightly decrease (∼6 %) but the decrease in the mean rate will be larger (∼15 %). These changes will lead to an overall increased fractional contribution of ARs to heavy precipitation, with the greatest impact over the Iberian Peninsula (15 %–30 %) and western France (∼15 %). Likewise, the fractional share of yearly maximum precipitation attributable to ARs will increase over the Iberian Peninsula, the UK, and western France. Over Norway, average AR precipitation rates will decline by −5 % to −30 %, most likely due to dynamic changes, with ARs originating from latitudes > 60∘ N decreasing by up to 20 % and those originating south of 45∘ N increasing. This suggests that ARs over Norway will follow longer routes over the continent, such that additional moisture uptake will be impeded. By contrast, ARs from >60∘ N will take up moisture from the North Atlantic before making landfall over Norway. The found changes in the local AR pathway are probably driven by larger-scale circulation changes such as a change in dominating weather regimes and/or changes in the winter storm track over the North Atlantic

Impact of air–sea coupling on the simulation of Indian summer monsoon using a high-resolution Regional Earth System Model over CORDEX-SA

A new high-resolution Regional Earth System Model, namely ROM, has been implemented over CORDEX-SA towards examining the impact of air–sea coupling on the Indian summer monsoon characteristics. ROM's simulated mean ISM rainfall and associated dynamical and thermodynamical processes, including the representation of northward and eastward propagating convention bands, are closer to observation than its standalone atmospheric model component (REMO), highlighting the advantage of air–sea coupling. However, the value addition of air–sea coupling varies spatially with more significant improvements over regions with large biases. Bay of Bengal and the eastern equatorial Indian Ocean are the most prominent region where the highest added value is observed with a significant reduction up to 50–500% precipitation bias. Most of the changes in precipitation over the ocean are associated with convective precipitation (CP) due to the suppression of convective activity caused by the negative feedback due to the inclusion of air–sea coupling. However, CP and large-scale precipitation (LP) improvements show east–west asymmetry over the Indian land region. The substantial LP bias reduction is noticed over the wet bias region of western central India due to its suppression, while enhanced CP over eastern central India contributed to the reduction of dry bias. An insignificant change is noticed over Tibetan Plateau, northern India, and Indo Gangetic plains. The weakening of moisture-laden low-level Somalia Jets causes the diminishing of moisture supply from the Arabian Sea (AS) towards Indian land regions resulting in suppressed precipitation, reducing wet bias, especially over western central India. The anomalous high kinetic energy over AS, wind shear, and tropospheric temperature gradient in REMO compared to observation is substantially reduced in the ROM, facilitating the favourable condition for suppressing moisture feeding and hence the wet bias over west-central India in ROM. The warmer midlatitude in ROM than REMO over eastern central India strengthens the convection, enhancing precipitation results in reducing the dry bias. Despite substantially improved ROM’performance, it still exhibits some systematic biases (wet/dry) partially associated with the persistent warm/cold SST bias and land–atmosphere interaction

Impact of air–sea coupling on the climate change signal over the Iberian Peninsula

In this work we use a regional atmosphere–ocean coupled model (RAOCM) and its stand-alone atmospheric component to gain insight into the impact of atmosphere–ocean coupling on the climate change signal over the Iberian Peninsula (IP). The IP climate is influenced by both the Atlantic Ocean and the Mediterranean sea. Complex interactions with the orography take place there and high-resolution models are required to realistically reproduce its current and future climate. We find that under the RCP8.5 scenario, the generalized 2-m air temperature (T2M) increase by the end of the twenty-first century (2070–2099) in the atmospheric-only simulation is tempered by the coupling. The impact of coupling is specially seen in summer, when the warming is stronger. Precipitation shows regionally-dependent changes in winter, whilst a drier climate is found in summer. The coupling generally reduces the magnitude of the changes. Differences in T2M and precipitation between the coupled and uncoupled simulations are caused by changes in the Atlantic large-scale circulation and in the Mediterranean Sea. Additionally, the differences in projected changes of T2M and precipitation with the RAOCM under the RCP8.5 and RCP4.5 scenarios are tackled. Results show that in winter and summer T2M increases less and precipitation changes are of a smaller magnitude with the RCP4.5. Whilst in summer changes present a similar regional distribution in both runs, in winter there are some differences in the NW of the IP due to differences in the North Atlantic circulation. The differences in the climate change signal from the RAOCM and the driving Global Coupled Model show that regionalization has an effect in terms of higher resolution over the land and ocean

Antarctic sea ice decline delayed well into the 21st century in a high-resolution climate projection

Despite ongoing global warming and strong sea ice decline in the Arctic, the sea ice extent around the Antarctic continent has not declined during the satellite era since 1979. This is in stark contrast to existing climate models that tend to show a strong negative sea ice trend for the same period; hence the confidence in projected Antarctic sea-ice changes is considered to be low. In the years since 2016, there has been significantly lower Antarctic sea ice extent, which some consider a sign of imminent change; however, others have argued that sea ice extent is expected to regress to the weak decadal trend in the near future. In this presentation, we show results from climate change projections with a new climate model that allows the simulation of mesoscale eddies in dynamically active ocean regions in a computationally efficient way. We find that the high-resolution configuration (HR) favours periods of stable Antarctic sea ice extent in September as observed over the satellite era. Sea ice is not projected to decline well into the 21st century in the HR simulations, which is similar to the delaying effect of, e.g., added glacial melt water in recent studies. The HR ocean configurations simulate an ocean heat transport that responds differently to global warming and is more efficient at moderating the anthropogenic warming of the Southern Ocean. As a consequence, decrease of Antarctic sea ice extent is significantly delayed, in contrast to what existing coarser-resolution climate models predict. Other explanations why current models simulate a non-observed decline of Antarctic sea-ice have been put forward, including the choice of included sea ice physics and underestimated simulated trends in westerly winds. Our results provide an alternative mechanism that might be strong enough to explain the gap between modeled and observed trends alone