On the Damping Time Scale of EVP Sea Ice Dynamics

We propose to make the damping time scale, which governs the decay of pseudo-elastic waves in the Elastic Viscous Plastic (EVP) sea-ice solvers, independent of the external time step and large enough to warrant numerical stability for a moderate number of internal time steps. A necessary condition is that the forcing on sea ice varies slowly on the damping time scale, in which case an EVP solution may still approach a Viscous Plastic one, but on a time scale longer than a single external time step. In this case, the EVP method becomes very close to the recently proposed modified EVP (mEVP) method in terms of stability and simulated behavior. In a simple test case dealing with sea ice breaking under the forcing of a moving cyclone, the EVP method with an enlarged damping time scale can simulate linear kinematic features which are very similar to those from the traditional EVP implementation, although a much smaller number of internal time steps is used. There is more difference in sea-ice thickness and linear kinematic features simulated in a realistic Arctic configuration between using the traditional and our suggested choices of EVP damping time scales, but it is minor considering model uncertainties associated with choices of many other parameters in sea-ice models

Modified halocline water over the Laptev Sea continental margin : historical data analysis

Historical hydrographic data (1940s–2010) show a distinct cross-slope difference of the lower halocline water (LHW) over the Laptev Sea continental margins. Over the slope, the LHW is on average warmer and saltier by 0.2°C and 0.5 psu, respectively, relative to the off-slope LHW. The LHW temperature time series constructed from the on-slope historical records are related to the temperature of the Atlantic Water (AW) boundary current transporting warm water from the North Atlantic Ocean. In contrast, the on-slope LHW salinity is linked to the sea ice and wind forcing over the potential upstream source region in the Barents and northern Kara Seas, as also indicated by hydrodynamic model results. Over the Laptev Sea continental margin, saltier LHW favors weaker salinity stratification that, in turn, contributes to enhanced vertical mixing with underlying AW

Ocean Model Formulation Influences Transient Climate Response

The transient climate response (TCR) is 20% higher in the Alfred Wegener Institute Climate Model (AWI-CM) compared to the Max Planck Institute Earth System Model (MPI-ESM) whereas the equilibrium climate sensitivity (ECS) is by up to 10% higher in AWI-CM. These results are largely independent of the two considered model resolutions for each model. The two coupled CMIP6 models share the same atmosphere-land component ECHAM6.3 developed at the Max Planck Institute for Meteorology (MPI-M). However, ECHAM6.3 is coupled to two different ocean models, namely the MPIOM sea ice-ocean model developed at MPI-M and the FESOM sea ice-ocean model developed at the Alfred Wegener Institute, Helmholtz Centre for Polar and Marine Research (AWI). A reason for the different TCR is related to ocean heat uptake in response to greenhouse gas forcing. Specifically, AWI-CM simulations show stronger surface heating than MPI-ESM simulations while the latter accumulate more heat in the deeper ocean. The vertically integrated ocean heat content is increasing slower in AWI-CM model configurations compared to MPI-ESM model configurations in the high latitudes. Weaker vertical mixing in AWI-CM model configurations compared to MPI-ESM model configurations seems to be key for these differences. The strongest difference in vertical ocean mixing occurs inside the Weddell and Ross Gyres and the northern North Atlantic. Over the North Atlantic, these differences materialize in a lack of a warming hole in AWI-CM model configurations and the presence of a warming hole in MPI-ESM model configurations. All these differences occur largely independent of the considered model resolutions

Identifying climate change information needs for the himalayan region: Results from the GLACINDIA Stakeholder Workshop and Training Program

Here we present results of a workshop designed to bring together stakeholders from different states of the Indian side of the Himalayan arc and an international group of climate scientists in order to discuss how climate change research for this region can be tailored toward the needs of local communities. The stakeholder workshop was jointly organized by the Jawaharlal Nehru University (JNU), New Delhi, India, and the Climate Service Center 2.0, Hamburg, Germany, within the framework of the multidisciplinary international research project GLACINDIA. The project focuses on the water-related effects of changes in glacier mass balance and river runoff in western Himalayas. Given the research focus of the GLACINDIA project, the initial focus of the workshop was on glacier-related hydrological information. During stakeholder interactions the resulting discussion covered a much broader range of urgent climate change information needs for the Himalayan region.publishedVersio

On the impact of atmospheric vs oceanic resolutions on the representation of the sea surface temperature in the South Eastern Tropical Atlantic

Despite the efforts of the modelling community to improve the representation of the sea surface temperature (SST) over the South Eastern Tropical Atlantic, warm biases still persist. In this work we use four different configurations of the fully-coupled AWI Climate Model (AWI-CM) which allow us to gain physics-based insight into the role of the oceanic and atmospheric resolutions of the model in the regional distribution of the SST. Our results show that a sole refinement of the oceanic resolution reduces warm biases further than a single increase of the atmospheric component. An increased oceanic resolution is required (i) to simulate properly the Agulhas Current and its associated rings; (ii) to reinforce the northward-flowing Benguela Current and (iii) to intensify coastal upwelling. The best results are obtained when both resolutions are refined. However, even in that case, warm biases persist, reflecting that some processes and feedbacks are still not optimally resolved. Our results indicate that overheating is not due to insufficient upwelling, but rather due to upwelling of waters which are warmer than observations as a result of an erroneous representation of the vertical distribution of temperature. Errors in the representation of the vertical temperature profile are the consequence of a warm bias in the simulated climate state

On the role of horizontal resolution over the Tibetan Plateau in the REMO regional climate model

A number of studies have shown that added value is obtained by increasing the horizontal resolution of a regional climate model to capture additional fine-scale weather processes. However, the mechanisms leading to this added value are different over areas with complicated orographic features, such as the Tibetan Plateau (TP). To determine the role that horizontal resolution plays over the TP, a detailed comparison was made between the results from the REMO regional climate model at resolutions of 25 and 50 km for the period 1980–2007. The model was driven at the lateral boundaries by the European Centre for Medium-Range Weather Forecasts Interim Reanalysis data. The experiments differ only in representation of topography, all other land parameters (e.g., vegetation characteristics, soil texture) are the same. The results show that the high-resolution topography affects the regional air circulation near the ground surface around the edge of the TP, which leads to a redistribution of the transport of atmospheric water vapor, especially over the Brahmaputra and Irrawaddy valleys—the main water vapor paths for the southern TP—increasing the amount of atmospheric water vapor transported onto the TP by about 5. This, in turn, significantly decreases the temperature at 2 m by > 1.5 °C in winter in the high-resolution simulation of the southern TP. The impact of topography on the 2 m temperature over the TP is therefore by influencing the transport of atmospheric water vapor in the main water vapor paths. © 2018 Springer-Verlag GmbH Germany, part of Springer Natur

Sensitivity of deep ocean biases to horizontal resolution in prototype CMIP6 simulations with AWI-CM1.0

Models from phase 5 of the Coupled Model Inter-comparison Project (CMIP5) show substantial biases in the deep ocean that are larger than the level of natural variability and the response to enhanced greenhouse gas concentrations. Here, we analyze the influence of horizontal resolution in a hierarchy of five multi-resolution simulations with the AWI Climate Model (AWI-CM), the climate model used at the Al-fred Wegener Institute, Helmholtz Centre for Polar and Ma-rine Research, which employs a sea ice–ocean model com-ponent formulated on unstructured meshes. The ocean grid sizes considered range from a nominal resolution of ∼ 1◦ (CMIP5 type) up to locally eddy resolving. We show that increasing ocean resolution locally to resolve ocean eddies leads to reductions in deep ocean biases, although these im-provements are not strictly monotonic for the five different ocean grids. A detailed diagnosis of the simulations allows to identify the origins of the biases. We find that two key re-gions at the surface are responsible for the development of the deep bias in the Atlantic Ocean: the northeastern North Atlantic and the region adjacent to the Strait of Gibraltar. Furthermore, the Southern Ocean density structure is equally improved with locally explicitly resolved eddies compared to parameterized eddies. Part of the bias reduction can be traced back towards improved surface biases over outcrop-ping regions, which are in contact with deeper ocean layers along isopycnal surfaces. Our prototype simulations provide guidance for the optimal choice of ocean grids for AWI-CM to be used in the final runs for phase 6 of CMIP (CMIP6) and for the related flagship simulations in the High Resolution Model Intercomparison Project (HighResMIP). Quite remarkably, retaining resolution only in areas of high eddy activity along with excellent scalability characteristics of the unstructured-mesh sea ice–ocean model enables us to per-form the multi-centennial climate simulations needed in a CMIP context at (locally) eddy-resolving resolution with a throughput of 5–6 simulated years per day

Simulations for CMIP6 With the AWI Climate Model AWI‐CM‐1‐1

The Alfred Wegener Institute Climate Model (AWI‐CM) participates for the first time in the Coupled Model Intercomparison Project (CMIP), CMIP6. The sea ice‐ocean component, FESOM, runs on an unstructured mesh with horizontal resolutions ranging from 8 to 80 km. FESOM is coupled to the Max Planck Institute atmospheric model ECHAM 6.3 at a horizontal resolution of about 100 km. Using objective performance indices, it is shown that AWI‐CM performs better than the average of CMIP5 models. AWI‐CM shows an equilibrium climate sensitivity of 3.2°C, which is similar to the CMIP5 average, and a transient climate response of 2.1°C which is slightly higher than the CMIP5 average. The negative trend of Arctic sea‐ice extent in September over the past 30 years is 20–30% weaker in our simulations compared to observations. With the strongest emission scenario, the AMOC decreases by 25% until the end of the century which is less than the CMIP5 average of 40%. Patterns and even magnitude of simulated temperature and precipitation changes at the end of this century compared to present‐day climate under the strong emission scenario SSP585 are similar to the multi‐model CMIP5 mean. The simulations show a 11°C warming north of the Barents Sea and around 2°C to 3°C over most parts of the ocean as well as a wetting of the Arctic, subpolar, tropical, and Southern Ocean. Furthermore, in the northern middle latitudes in boreal summer and autumn as well as in the southern middle latitudes, a more zonal atmospheric flow is projected throughout the year

The Relative Influence of Atmospheric and Oceanic Model Resolution on the Circulation of the North Atlantic Ocean in a Coupled Climate Model

It is often unclear how to optimally choose horizontal resolution for the oceanic and atmospheric components of coupled climate models, which has implications for their ability to make best use of available computational resources. Here we investigate the effect of using different combinations of horizontal resolutions in atmosphere and ocean on the simulated climate in a global coupled climate model (Alfred Wegener Institute Climate Model [AWI‐CM]). Particular attention is given to the Atlantic Meridional Overturning Circulation (AMOC). Four experiments with different combinations of relatively high and low resolutions in the ocean and atmosphere are conducted. We show that increases in atmospheric and oceanic resolution have clear impacts on the simulated AMOC, which are largely independent. Increased atmospheric resolution leads to a weaker AMOC. It also improves the simulated Gulf Stream separation; however, this is only the case if the ocean is locally eddy resolving and reacts to the improved atmosphere. We argue that our results can be explained by reduced mean winds caused by higher cyclone activity. Increased resolution of the ocean affects the AMOC in several ways, thereby locally increasing or reducing the AMOC. The finer topography (and reduced dissipation) in the vicinity of the Caribbean basin tends to locally increase the AMOC. However, there is a reduction in the AMOC around 45°N, which relates to the reduced mixed layer depth in the Labrador Sea in simulations with refined ocean and changes in the North Atlantic current pathway. Furthermore, the eddy‐induced changes in the Southern Ocean increase the strength of the deep cell

Regionally Coupled Atmosphere-Ocean-Marine Biogeochemistry Model ROM: 2. Studying the Climate Change Signal in the North Atlantic and Europe

Climate simulations for the North Atlantic and Europe for recent and future conditions simulated with the regionally coupled ROM model are analyzed and compared to the results from the MPI-ESM. The ROM simulations also include a biogeochemistry and ocean tides. For recent climate conditions, ROM generally improves the simulations compared to the driving model MPI-ESM. Reduced oceanic biases in the Northern Atlantic are found, as well as a better simulation of the atmospheric circulation, notably storm tracks and blocking. Regarding future climate projections for the 21st century following the RCP 4.5 and 8.5 scenarios, MPI-ESM and ROM largely agree qualitatively on the climate change signal over Europe. However, many important differences are identified. For example, ROM shows an SST cooling in the Subpolar Gyre, which is not present in MPI-ESM. Under the RCP8.5 scenario, ROM Arctic sea ice cover is thinner and reaches the seasonally ice-free state by 2055, well before MPI-ESM. This shows the decisive importance of higher ocean resolution and regional coupling for determining the regional responses to global warming trends. Regarding biogeochemistry, both ROM and MPI-ESM simulate a widespread decline in winter nutrient concentration in the North Atlantic of up to similar to 35%. On the other hand, the phytoplankton spring bloom in the Arctic and in the North-Western Atlantic starts earlier, and the yearly primary production is enhanced in the Arctic in the late 21st century. These results clearly demonstrate the added value of ROM to determine more detailed and more reliable climate projections at the regional scale. Plain Language Summary We downscale present climate and future climate change projections for the North Atlantic and Europe using a regionally coupled Earth System Model including atmosphere, ocean, river runoff, and ocean biogeochemistry components. This approach allows us to attain higher spatial resolution and to a more accurate representation of key physical processes, yielding a better simulation of present climate at regional and local scales when compared to the driving global climate model. Future climate change projections show more detail at regional and local scale, mostly related to the improvement in the representation of orography and bathymetry. These improvements along with a better representation of the interactive ocean-atmosphere coupling lead to other remarkable differences with the driving global climate model: (1) colder Sea Surface Temperature in the Subpolar Gyre region, indicating a local convection collapse and a Atlantic Meridional Overturning Circulation slowdown; (2) a seasonal free-ice Arctic is reached by 2055 under RCP8.5 scenario, well before projected by the driving global climate model; and (3) stronger reduction in nutrients in the North Atlantic by the end of the 21st century. These results clearly demonstrate the added value of the regionally coupled model system to determine more reliable climate projections at the regional scale