Intercomparison of the Arctic sea ice cover in global ocean–sea ice reanalyses from the ORA-IP project

Ocean–sea ice reanalyses are crucial for assessing the variability and recent trends in the Arctic sea ice cover. This is especially true for sea ice volume, as long-term and large scale sea ice thickness observations are inexistent. Results from the Ocean ReAnalyses Intercomparison Project (ORA-IP) are presented, with a focus on Arctic sea ice fields reconstructed by state-of-the-art global ocean reanalyses. Differences between the various reanalyses are explored in terms of the effects of data assimilation, model physics and atmospheric forcing on properties of the sea ice cover, including concentration, thickness, velocity and snow. Amongst the 14 reanalyses studied here, 9 assimilate sea ice concentration, and none assimilate sea ice thickness data. The comparison reveals an overall agreement in the reconstructed concentration fields, mainly because of the constraints in surface temperature imposed by direct assimilation of ocean observations, prescribed or assimilated atmospheric forcing and assimilation of sea ice concentration. However, some spread still exists amongst the reanalyses, due to a variety of factors. In particular, a large spread in sea ice thickness is found within the ensemble of reanalyses, partially caused by the biases inherited from their sea ice model components. Biases are also affected by the assimilation of sea ice concentration and the treatment of sea ice thickness in the data assimilation process. An important outcome of this study is that the spatial distribution of ice volume varies widely between products, with no reanalysis standing out as clearly superior as compared to altimetry estimates. The ice thickness from systems without assimilation of sea ice concentration is not worse than that from systems constrained with sea ice observations. An evaluation of the sea ice velocity fields reveals that ice drifts too fast in most systems. As an ensemble, the ORA-IP reanalyses capture trends in Arctic sea ice area and extent relatively well. However, the ensemble can not be used to get a robust estimate of recent trends in the Arctic sea ice volume. Biases in the reanalyses certainly impact the simulated air–sea fluxes in the polar regions, and questions the suitability of current sea ice reanalyses to initialize seasonal forecasts.

Can climate models represent the precipitation associated with extratropical cyclones?

Extratropical cyclones produce the majority of precipitation in many regions of the extratropics. This study evaluates the ability of a climate model, HiGEM, to reproduce the precipitation associated with extratropical cyclones. The model is evaluated using the ERA-Interim reanalysis and GPCP dataset. The analysis employs a cyclone centred compositing technique, evaluates composites across a range of geographical areas and cyclone intensities and also investigates the ability of the model to reproduce the climatological distribution of cyclone associated precipitation across the Northern Hemisphere. Using this phenomena centred approach provides an ability to identify the processes which are responsible for climatological biases in the model. Composite precipitation intensities are found to be comparable when all cyclones across the Northern Hemisphere are included. When the cyclones are filtered by region or intensity, differences are found, in particular, HiGEM produces too much precipitation in its most intense cyclones relative to ERA-Interim and GPCP. Biases in the climatological distribution of cyclone associated precipitation are also found, with biases around the storm track regions associated with both the number of cyclones in HiGEM and also their average precipitation intensity. These results have implications for the reliability of future projections of extratropical precipitation from the model.

Spatio-temporal variability of polynya dynamics and ice production in the Laptev Sea between the winters of 1979/80 and 2007/08

Polynyas in the Laptev Sea are examined with respect to recurrence and interannual wintertime ice production.We use a polynya classification method based on passive microwave satellite data to derive daily polynya area from long-term sea-ice concentrations. This provides insight into the spatial and temporal variability of open-water and thin-ice regions on the Laptev Sea Shelf. Using thermal infrared satellite data to derive an empirical thin-ice distribution within the thickness range from 0 to 20 cm, we calculate daily average surface heat loss and the resulting wintertime ice formation within the Laptev Sea polynyas between 1979 and 2008 using reanalysis data supplied by the National Centers for Environmental Prediction, USA, as atmospheric forcing. Results indicate that previous studies significantly overestimate the contribution of polynyas to the ice production in the Laptev Sea. Average wintertime ice production in polynyas amounts to approximately 55 km39 27% and is mostly determined by the polynya area, wind speed and associated large-scale circulation patterns. No trend in ice production could be detected in the period from 1979/80 to 2007/08.

Evaluation of moisture sources for the Central European summer flood of May/June 2013 based on regional climate model simulations

Heavy precipitation affected Central Europe in May/June 2013, triggering damaging floods both on the Danube and the Elbe rivers. Based on a modelling approach with COSMO-CLM, moisture fluxes, backward trajectories, cyclone tracks and precipitation fields are evaluated for the relevant time period 30 May–2 June 2013. We identify potential moisture sources and quantify their contribution to the flood event focusing on the Danube basin through sensitivity experiments: Control simulations are performed with undisturbed ERA-Interim boundary conditions, while multiple sensitivity experiments are driven with modified evaporation characteristics over selected marine and land areas. Two relevant cyclones are identified both in reanalysis and in our simulations, which moved counter-clockwise in a retrograde path from Southeastern Europe over Eastern Europe towards the northern slopes of the Alps. The control simulations represent the synoptic evolution of the event reasonably well. The evolution of the precipitation event in the control simulations shows some differences in terms of its spatial and temporal characteristics compared to observations. The main precipitation event can be separated into two phases concerning the moisture sources. Our modelling results provide evidence that the two main sources contributing to the event were the continental evapotranspiration (moisture recycling; both phases) and the North Atlantic Ocean (first phase only). The Mediterranean Sea played only a minor role as a moisture source. This study confirms the importance of continental moisture recycling for heavy precipitation events over Central Europe during the summer half year.

Understanding climate change projections for precipitation over western Europe with a weather typing approach

Precipitation over western Europe (WE) is projected to increase (decrease) roughly northward (equatorward) of 50°N during the 21st century. These changes are generally attributed to alterations in the regional large-scale circulation, e.g., jet stream, cyclone activity, and blocking frequencies. A novel weather typing within the sector (30°W–10°E, 25–70°N) is used for a more comprehensive dynamical interpretation of precipitation changes. A k-means clustering on daily mean sea level pressure was undertaken for ERA-Interim reanalysis (1979–2014). Eight weather types are identified: S1, S2, S3 (summertime types), W1, W2, W3 (wintertime types), B1, and B2 (blocking-like types). Their distinctive dynamical characteristics allow identifying the main large-scale precipitation-driving mechanisms. Simulations with 22 Coupled Model Intercomparison Project 5 models for recent climate conditions show biases in reproducing the observed seasonality of weather types. In particular, an overestimation of weather type frequencies associated with zonal airflow is identified. Considering projections following the (Representative Concentration Pathways) RCP8.5 scenario over 2071–2100, the frequencies of the three driest types (S1, B2, and W3) are projected to increase (mainly S1, +4%) in detriment of the rainiest types, particularly W1 (−3%). These changes explain most of the precipitation projections over WE. However, a weather type-independent background signal is identified (increase/decrease in precipitation over northern/southern WE), suggesting modifications in precipitation-generating processes and/or model inability to accurately simulate these processes. Despite these caveats in the precipitation scenarios for WE, which must be duly taken into account, our approach permits a better understanding of the projected trends for precipitation over WE.

Improvement and sensitivity analysis of thermal thin-ice thickness retrievals

Considering the sea ice decline in the Arctic during the last decades, polynyas are of high research interest since these features are core areas of new ice formation. The determination of ice formation requires accurate retrieval of polynya area and thin-ice thickness (TIT) distribution within the polynya.We use an established energy balance model to derive TITs with MODIS ice surface temperatures (Ts) and NCEP/DOE Reanalysis II in the Laptev Sea for two winter seasons. Improvements of the algorithm mainly concern the implementation of an iterative approach to calculate the atmospheric flux components taking the atmospheric stratification into account. Furthermore, a sensitivity study is performed to analyze the errors of the ice thickness. The results are the following: 1) 2-m air temperatures (Ta) and Ts have the highest impact on the retrieved ice thickness; 2) an overestimation of Ta yields smaller ice thickness errors as an underestimation of Ta; 3) NCEP Ta shows often a warm bias; and 4) the mean absolute error for ice thicknesses up to 20 cm is ±4.7 cm. Based on these results, we conclude that, despite the shortcomings of the NCEP data (coarse spatial resolution and no polynyas), this data set is appropriate in combination with MODIS Ts for the retrieval of TITs up to 20 cm in the Laptev Sea region. The TIT algorithm can be applied to other polynya regions and to past and future time periods. Our TIT product is a valuable data set for verification of other model and remote sensing ice thickness data.

Abrupt transitions in the NAO control of explosive North Atlantic cyclone development

Explosive cyclones are intense extra-tropical low pressure systems featuring large deepening rates. In the Euro-Atlantic sector, they are a major source of life-threatening weather impacts due to their associated strong wind gusts, heavy precipitation and storm surges. The wintertime variability of the North Atlantic cyclonic activity is primarily modulated by the North Atlantic Oscillation (NAO). In this study, we investigate the interannual and multi-decadal variability of explosive North Atlantic cyclones using track density data from two reanalysis datasets (NCEP and ERA-40) and a control simulation of an atmosphere/ocean coupled General Circulation Model (GCM—ECHAM5/MPIOM1). The leading interannual and multi-decadal modes of variability of explosive cyclone track density are characterized by a strengthening/weakening pattern between Newfoundland and Iceland, which is mainly modulated by the NAO at both timescales. However, the NAO control of interannual cyclone variability is not stationary in time and abruptly fluctuates during periods of 20–25 years long both in NCEP and ECHAM5/MPIOM1. These transitions are accompanied by structural changes in the leading mode of explosive cyclone variability, and by decreased/enhanced baroclinicity over the sub-polar/sub-tropical North Atlantic. The influence of the ocean is apparently important for both the occurrence and persistence of such anomalous periods. In the GCM, the Atlantic Meridional Overturning Circulation appears to influence the large-scale baroclinicity and explosive cyclone development over the North Atlantic. These results permit a better understanding of explosive cyclogenesis variability at different climatic timescales and might help to improve predictions of these hazardous events.

Characterisation of synoptic conditions and cyclones associated with top ranking potential wind loss events over Iberia

Intense extra-tropical cyclones are often associated with strong winds, heavy precipitation and socio-economic impacts. Over southwestern Europe, such storms occur less often, but still cause high economic losses. We characterise the largescale atmospheric conditions and cyclone tracks during the top-100 potential losses over Iberia associated with wind events. Based on 65 years of reanalysis data,events are classified into four groups: (i) cyclone tracks crossing over Iberia on the event day (“Iberia”), (ii) cyclones crossing further north, typically southwest of the British Isles (“North”), (iii) cyclones crossing southwest to northeast near the northwest tip of Iberia (“West”), and (iv) so called “Hybrids”, characterised by a strong pressure gradient over Iberia due to the juxtaposition of low and high pressure centres. Generally, “Iberia” events are the most frequent (31% to 45% for top-100 vs.top-20), while “West” events are rare (10% to 12%). 70% of the events were primarily associated with a cyclone. Multi-decadal variability in the number of events is identified. While the peak in recent years is quite prominent, other comparably stormy periods occurred in the 1960s and 1980s. This study documents that damaging wind storms over Iberia are not rare events, and their frequency of occurrence undergoes strong multi-decadal variability.

Sensitivity of resolved and parameterized surface drag to changes in resolution and parameterization

The relative contribution of resolved and parameterized surface drag towards balancing the atmospheric angular momentum flux convergence (AMFC), and their sensitivity to horizontal resolution and parameterization, are investigated in an atmospheric model. This sensitivity can be difficult to elucidate in free-running climate models, in which the AMFC varies with changing climatologies and, as a result, the relative contributions of surface terms balancing the AMFC also vary. While the sensitivity question has previously been addressed using short-range forecasts, we demonstrate that a nudging framework is an effective method for constraining the AMFC. The Met Office Unified Model is integrated at three horizontal resolutions ranging from 130 km (N96) to 25 km (N512) while relaxing the model’s wind and temperature fields towards the ERAinterim reanalysis within the altitude regions of maximum AMFC. This method is validated against short range forecasts and good agreement is found. These experiments are then used to assess the fidelity of the exchange between parameterized and resolved orographic torques with changes in horizontal resolution. Although the parameterized orographic torque reduces substantially with increasing horizontal resolution, there is little change in resolved orographic torque over 20N to 50N. The tendencies produced by the nudging routine indicate that the additional drag at lower horizontal resolution is excessive. When parameterized orographic blocking is removed at the coarsest of these resolutions, there is a lack of compensation, and even compensation of the opposite sense, by the boundary layer and resolved torques which is particularly pronounced over 20N to 50N. This study demonstrates that there is strong sensitivity in the behaviour of the resolved and parameterized surface drag over this region.