Recent progress in understanding and projecting regional and global mean sea-level change

Considerable progress has been made in understanding the present and future regional and global sea level in the 2 years since the publication of the Fifth Assessment Report (AR5) of the Intergovernmental Panel on Climate Change. Here, we evaluate how the new results affect the AR5’s assessment of (i) historical sea level rise, including attribution of that rise and implications for the sea level budget, (ii) projections of the components and of total global mean sea level (GMSL), and (iii) projections of regional variability and emergence of the anthropogenic signal. In each of these cases, new work largely provides additional evidence in support of the AR5 assessment, providing greater confidence in those findings. Recent analyses confirm the twentieth century sea level rise, with some analyses showing a slightly smaller rate before 1990 and some a slightly larger value than reported in the AR5. There is now more evidence of an acceleration in the rate of rise. Ongoing ocean heat uptake and associated thermal expansion have continued since 2000, and are consistent with ocean thermal expansion reported in the AR5. A significant amount of heat is being stored deeper in the water column, with a larger rate of heat uptake since 2000 compared to the previous decades and with the largest storage in the Southern Ocean. The first formal detection studies for ocean thermal expansion and glacier mass loss since the AR5 have confirmed the AR5 finding of a significant anthropogenic contribution to sea level rise over the last 50 years. New projections of glacier loss from two regions suggest smaller contributions to GMSL rise from these regions than in studies assessed by the AR5; additional regional studies are required to further assess whether there are broader implications of these results. Mass loss from the Greenland Ice Sheet, primarily as a result of increased surface melting, and from the Antarctic Ice Sheet, primarily as a result of increased ice discharge, has accelerated. The largest estimates of acceleration in mass loss from the two ice sheets for 2003–2013 equal or exceed the acceleration of GMSL rise calculated from the satellite altimeter sea level record over the longer period of 1993–2014. However, when increased mass gain in land water storage and parts of East Antarctica, and decreased mass loss from glaciers in Alaska and some other regions are taken into account, the net acceleration in the ocean mass gain is consistent with the satellite altimeter record. New studies suggest that a marine ice sheet instability (MISI) may have been initiated in parts of the West Antarctic Ice Sheet (WAIS), but that it will affect only a limited number of ice streams in the twenty-first century. New projections of mass loss from the Greenland and Antarctic Ice Sheets by 2100, including a contribution from parts of WAIS undergoing unstable retreat, suggest a contribution that falls largely within the likely range (i.e., two thirds probability) of the AR5. These new results increase confidence in the AR5 likely range, indicating that there is a greater probability that sea level rise by 2100 will lie in this range with a corresponding decrease in the likelihood of an additional contribution of several tens of centimeters above the likely range. In view of the comparatively limited state of knowledge and understanding of rapid ice sheet dynamics, we continue to think that it is not yet possible to make reliable quantitative estimates of future GMSL rise outside the likely range. Projections of twenty-first century GMSL rise published since the AR5 depend on results from expert elicitation, but we have low confidence in conclusions based on these approaches. New work on regional projections and emergence of the anthropogenic signal suggests that the two commonly predicted features of future regional sea level change (the increasing tilt across the Antarctic Circumpolar Current and the dipole in the North Atlantic) are related to regional changes in wind stress and surface heat flux. Moreover, it is expected that sea level change in response to anthropogenic forcing, particularly in regions of relatively low unforced variability such as the low-latitude Atlantic, will be detectable over most of the ocean by 2040. The east-west contrast of sea level trends in the Pacific observed since the early 1990s cannot be satisfactorily accounted for by climate models, nor yet definitively attributed either to unforced variability or forced climate change.

Atmospheric and oceanic contributions to irreducible forecast uncertainty of Arctic surface climate

Uncertainty of Arctic seasonal to interannual predictions arising from model errors and initial state uncertainty has been widely discussed in the literature, whereas the irreducible forecast uncertainty (IFU) arising from the chaoticity of the climate system has received less attention. However, IFU provides important insights into the mechanisms through which predictability is lost, and hence can inform prioritization of model development and observations deployment. Here, we characterize how internal oceanic and surface atmospheric heat fluxes contribute to IFU of Arctic sea ice and upper ocean heat content in an Earth system model by analyzing a set of idealized ensemble prediction experiments. We find that atmospheric and oceanic heat flux are often equally important for driving unpredictable Arctic-wide changes in sea ice and surface water temperatures, and hence contribute equally to IFU. Atmospheric surface heat flux tends to dominate Arctic-wide changes for lead times of up to a year, whereas oceanic heat flux tends to dominate regionally and on interannual time scales. There is in general a strong negative covariance between surface heat flux and ocean vertical heat flux at depth, and anomalies of lateral ocean heat transport are wind-driven, which suggests that the unpredictable oceanic heat flux variability is mainly forced by the atmosphere. These results are qualitatively robust across different initial states, but substantial variations in the amplitude of IFU exist. We conclude that both atmospheric variability and the initial state of the upper ocean are key ingredients for predictions of Arctic surface climate on seasonal to interannual time scales.

Processes driving intraseasonal displacements of the eastern edge of the warm pool: the contribution of westerly wind events

We investigate the processes responsible for the intraseasonal displacements of the eastern edge of the western Pacific warm pool (WPEE), which appear to play a role in the onset and development of El Niño events. We use 25 years of output from an ocean general circulation model experiment that is able to accurately capture the observed displacements of the WPEE, sea level anomalies, and upper ocean zonal currents at intraseasonal time scales in the western and central Pacific Ocean. Our results confirm that WPEE displacements driven by westerly wind events (WWEs) are largely controlled by zonal advection. This paper has also two novel findings: first, the zonal current anomalies responsible for the WPEE advection are driven primarily by local wind stress anomalies and not by intraseasonal wind-forced Kelvin waves as has been shown in most previous studies. Second, we find that intraseasonal WPEE fluctuations that are not related to WWEs are generally caused by intraseasonal variations in net heat flux, in contrast to interannual WPEE displacements that are largely driven by zonal advection. This study hence raises an interesting question: can surface heat flux-induced zonal WPEE motions contribute to El Niño–Southern Oscillation evolution, as WWEs have been shown to be able to do?

Reconstructing the subsurface ocean decadal variability using surface nudging in a perfect model framework

Initialising the ocean internal variability for decadal predictability studies is a new area of research and a variety of ad hoc methods are currently proposed. In this study, we explore how nudging with sea surface temperature (SST) and salinity (SSS) can reconstruct the three-dimensional variability of the ocean in a perfect model framework. This approach builds on the hypothesis that oceanic processes themselves will transport the surface information into the ocean interior as seen in ocean-only simulations. Five nudged simulations are designed to reconstruct a 150 years “target” simulation, defined as a portion of a long control simulation. The nudged simulations differ by the variables restored to, SST or SST + SSS, and by the area where the nudging is applied. The strength of the heat flux feedback is diagnosed from observations and the restoring coefficients for SSS use the same time-scale. We observed that this choice prevents spurious convection at high latitudes and near sea-ice border when nudging both SST and SSS. In the tropics, nudging the SST is enough to reconstruct the tropical atmosphere circulation and the associated dynamical and thermodynamical impacts on the underlying ocean. In the tropical Pacific Ocean, the profiles for temperature show a significant correlation from the surface down to 2,000 m, due to dynamical adjustment of the isopycnals. At mid-to-high latitudes, SSS nudging is required to reconstruct both the temperature and the salinity below the seasonal thermocline. This is particularly true in the North Atlantic where adding SSS nudging enables to reconstruct the deep convection regions of the target. By initiating a previously documented 20-year cycle of the model, the SST + SSS nudging is also able to reproduce most of the AMOC variations, a key source of decadal predictability. Reconstruction at depth does not significantly improve with amount of time spent nudging and the efficiency of the surface nudging rather depends on the period/events considered. The joint SST + SSS nudging applied everywhere is the most efficient approach. It ensures that the right water masses are formed at the right surface density, the subsequent circulation, subduction and deep convection further transporting them at depth. The results of this study underline the potential key role of SSS for decadal predictability and further make the case for sustained large-scale observations of this field.

ENSO representation in climate models: from CMIP3 to CMIP5

We analyse the ability of CMIP3 and CMIP5 coupled ocean–atmosphere general circulation models (CGCMs) to simulate the tropical Pacific mean state and El Niño-Southern Oscillation (ENSO). The CMIP5 multi-model ensemble displays an encouraging 30 % reduction of the pervasive cold bias in the western Pacific, but no quantum leap in ENSO performance compared to CMIP3. CMIP3 and CMIP5 can thus be considered as one large ensemble (CMIP3 + CMIP5) for multi-model ENSO analysis. The too large diversity in CMIP3 ENSO amplitude is however reduced by a factor of two in CMIP5 and the ENSO life cycle (location of surface temperature anomalies, seasonal phase locking) is modestly improved. Other fundamental ENSO characteristics such as central Pacific precipitation anomalies however remain poorly represented. The sea surface temperature (SST)-latent heat flux feedback is slightly improved in the CMIP5 ensemble but the wind-SST feedback is still underestimated by 20–50 % and the shortwave-SST feedbacks remain underestimated by a factor of two. The improvement in ENSO amplitudes might therefore result from error compensations. The ability of CMIP models to simulate the SST-shortwave feedback, a major source of erroneous ENSO in CGCMs, is further detailed. In observations, this feedback is strongly nonlinear because the real atmosphere switches from subsident (positive feedback) to convective (negative feedback) regimes under the effect of seasonal and interannual variations. Only one-third of CMIP3 + CMIP5 models reproduce this regime shift, with the other models remaining locked in one of the two regimes. The modelled shortwave feedback nonlinearity increases with ENSO amplitude and the amplitude of this feedback in the spring strongly relates with the models ability to simulate ENSO phase locking. In a final stage, a subset of metrics is proposed in order to synthesize the ability of each CMIP3 and CMIP5 models to simulate ENSO main characteristics and key atmospheric feedbacks.

Seasonal surface urban energy balance and wintertime stability simulated using three land-surface models in the high-latitude city Helsinki

The performance of three urban land surface models, run in offline mode, with their default external parameters, is evaluated for two distinctly different sites in Helsinki: Torni and Kumpula. The former is a dense city centre site with 22% vegetation, while the latter is a suburban site with over 50% vegetation. At both locations the models are compared against sensible and latent heat fluxes measured using the eddy covariance technique, along with snow depth observations. The cold climate experienced by the city causes strong seasonal variations that include snow cover and stable atmospheric conditions. Most of the time the three models are able to account for the differences between the study areas as well as the seasonal and diurnal variability of the energy balance components. However, the performances are not systematic across the modelled components, season and surface type. The net all-wave radiation is well simulated, with the greatest uncertainties related to snowmelt timing, when the fraction of snow cover has a key role, particularly in determining the surface albedo. For the turbulent fluxes, more variation between the models is seen which can partly be explained by the different methods in their calculation and partly by surface parameter values. For the sensible heat flux, simulation of wintertime values was the main problem, which also leads to issues in predicting near-surface stabilities particularly at the dense city centre site. All models have the most difficulties in simulating latent heat flux. This study particularly emphasizes that improvements are needed in the parameterization of anthropogenic heat flux and thermal parameters in winter, snow cover in spring and evapotranspiration in order to improve the surface energy balance modelling in cold climate cities.

Scale interactions between the MJO and the western maritime continent

State-of-the-art regional climate model simulations that are able to resolve key mesoscale circulations are used, for the first time, to understand the interaction between the large-scale convective environment of the MJO and processes governing the strong diurnal cycle over the islands of the Maritime Continent (MC). Convection is sustained in the late afternoon just inland of the coasts due to sea breeze convergence. Previous work has shown that the variability in MC rainfall associated with the MJO is manifested in changes to this diurnal cycle; land-based rainfall peaks before the active convective envelope of the MJO reaches the MC, whereas oceanic rainfall rates peak whilst the active envelope resides over the region. The model simulations show that the main controls on oceanic MC rainfall in the early active MJO phases are the large-scale environment and atmospheric stability, followed by high oceanic latent heat flux forced by high near-surface winds in the later active MJO phases. Over land, rainfall peaks before the main convective envelope arrives (in agreement with observations), even though the large-scale convective environment is only moderately favourable for convection. The causes of this early rainfall peak are convective triggers from land-sea breeze circulations that are strong due to high surface insolation and surface heating. During the peak MJO phases cloud cover increases and surface insolation decreases, which weakens the strength of the mesoscale circulations and reduces land-based rainfall, even though the large-scale environment remains favourable for convection at this time. Hence, scale interactions are an essential part of the MJO transition across the MC.

Impact of atmospheric forcing data on simulations of the Laptev Sea polynya dynamics using the sea-ice ocean model FESOM

The polynyas of the Laptev Sea are regions of particular interest due to the strong formation of Arctic sea-ice. In order to simulate the polynya dynamics and to quantify ice production, we apply the Finite Element Sea-Ice Ocean Model FESOM. In previous simulations FESOM has been forced with daily atmospheric NCEP (National Centers for Environmental Prediction) 1. For the periods 1 April to 9 May 2008 and 1 January to 8 February 2009 we examine the impact of different forcing data: daily and 6-hourly NCEP reanalyses 1 (1.875° x 1.875°), 6-hourly NCEP reanalyses 2 (1.875° x 1.875°), 6-hourly analyses from the GME (Global Model of the German Weather Service) (0.5° x 0.5°) and high-resolution hourly COSMO (Consortium for Small-Scale Modeling) data (5 km x 5 km). In all FESOM simulations, except for those with 6-hourly and daily NCEP 1 data, the openings and closings of polynyas are simulated in principle agreement with satellite products. Over the fast-ice area the wind fields of all atmospheric data are similar and close to in situ measurements. Over the polynya areas, however, there are strong differences between the forcing data with respect to air temperature and turbulent heat flux. These differences have a strong impact on sea-ice production rates. Depending on the forcing fields polynya ice production ranges from 1.4 km3 to 7.8 km3 during 1 April to 9 May 2011 and from 25.7 km3 to 66.2 km3 during 1 January to 8 February 2009. Therefore, atmospheric forcing data with high spatial and temporal resolution which account for the presence of the polynyas are needed to reduce the uncertainty in quantifying ice production in polynyas.

Impact of the Arctic Ocean Atlantic water layer on Siberian shelf hydrography

This paper examines the role of the Arctic Ocean Atlantic water (AW) in modifying the Laptev Sea shelf bottom hydrography on the basis of historical records from 1932 to 2008, field observations carried out in April–May 2008, and 2002–2009 cross‐slope measurements. A climatology of bottom hydrography demonstrates warming that extends offshore from the 30–50 m depth contour. Bottom layer temperature‐time series constructed from historical records links the Laptev Sea outer shelf to the AW boundary current transporting warm and saline water from the North Atlantic. The AW warming of the mid‐1990s and the mid‐2000s is consistent with outer shelf bottom temperature variability. For April–May 2008 we observed on‐shelf near‐bottom warm and saline water intrusions up to the 20 m isobath. These intrusions are typically about 0.2°C warmer and 1–1.5 practical salinity units saltier than ambient water. The 2002–2009 cross‐slope observations are suggestive for the continental slope upward heat flux from the AW to the overlying low‐halocline water (LHW). The lateral on‐shelf wind‐driven transport of the LHW then results in the bottom layer thermohaline anomalies recorded over the Laptev Sea shelf. We also found that polynya‐induced vertical mixing may act as a drainage of the bottom layer, permitting a relatively small portion of the AW heat to be directly released to the atmosphere. Finally, we see no significant warming (up until now) over the Laptev Sea shelf deeper than 10–15 m in the historical record. Future climate change, however, may bring more intrusions of Atlantic‐modified waters with potentially warmer temperature onto the shelf, which could have a critical impact on the stability of offshore submarine permafrost.

Processes controlling surface, bottom and lateral melt of Arctic sea ice in a state of the art sea ice model

We present a modelling study of processes controlling the summer melt of the Arctic sea ice cover. We perform a sensitivity study and focus our interest on the thermodynamics at the ice–atmosphere and ice–ocean interfaces. We use the Los Alamos community sea ice model CICE, and additionally implement and test three new parametrization schemes: (i) a prognostic mixed layer; (ii) a three equation boundary condition for the salt and heat flux at the ice–ocean interface; and (iii) a new lateral melt parametrization. Recent additions to the CICE model are also tested, including explicit melt ponds, a form drag parametrization and a halodynamic brine drainage scheme. The various sea ice parametrizations tested in this sensitivity study introduce a wide spread in the simulated sea ice characteristics. For each simulation, the total melt is decomposed into its surface, bottom and lateral melt components to assess the processes driving melt and how this varies regionally and temporally. Because this study quantifies the relative importance of several processes in driving the summer melt of sea ice, this work can serve as a guide for future research priorities.

Langmuir turbulence and surface heating in the ocean surface boundary layer

This study uses large-eddy simulation to investigate the structure of the ocean surface boundary layer (OSBL) in the presence of Langmuir turbulence and stabilizing surface heat fluxes. The OSBL consists of a weakly stratified layer, despite a surface heat flux, above a stratified thermocline. The weakly stratified (mixed) layer is maintained by a combination of a turbulent heat flux produced by the wave-driven Stokes drift and downgradient turbulent diffusion. The scaling of turbulence statistics, such as dissipation and vertical velocity variance, is only affected by the surface heat flux through changes in the mixed layer depth. Diagnostic models are proposed for the equilibrium boundary layer and mixed layer depths in the presence of surface heating. The models are a function of the initial mixed layer depth before heating is imposed and the Langmuir stability length. In the presence of radiative heating, the models are extended to account for the depth profile of the heating.

Surface Urban Energy and Water Balance Scheme (SUEWS): development and evaluation at two UK sites

The Surface Urban Energy and Water Balance Scheme (SUEWS) is evaluated at two locations in the UK: a dense urban site in the centre of London and a residential suburban site in Swindon. Eddy covariance observations of the turbulent fluxes are used to assess model performance over a twoyear period (2011-2013). The distinct characteristics of the sites mean their surface energy exchanges differ considerably. The model suggests the largest differences can be attributed to surface cover (notably the proportion of vegetated versus impervious area) and the additional energy supplied by human activities. SUEWS performs better in summer than winter, and better at the suburban site than the dense urban site. One reason for this is the bias towards suburban summer field campaigns in observational data used to parameterise this (and other) model(s). The suitability of model parameters (such as albedo, energy use and water use) for the UK sites is considered and, where appropriate, alternative values are suggested. An alternative parameterisation for the surface conductance is implemented, which permits greater soil moisture deficits before evaporation is restricted at non-irrigated sites. Accounting for seasonal variation in the estimation of storage heat flux is necessary to obtain realistic wintertime fluxes.

Modeling study of the aspect ratio influence on urban canopy energy fluxes with a modified wall-canyon energy budget scheme

The influence of the aspect ratio (building height/street canyon width) and the mean building height of cities on local energy fluxes and temperatures is studied by means of an Urban Canopy Model (UCM) coupled with a one-dimensional second-order turbulence closure model. The UCM presented is similar to the Town Energy Balance (TEB) model in most of its features but differs in a few important aspects. In particular, the street canyon walls are treated separately which leads to a different budget of radiation within the street canyon walls. The UCM has been calibrated using observations of incoming global and diffuse solar radiation, incoming long-wave radiation and air temperature at a site in So Paulo, Brazil. Sensitivity studies with various aspect ratios have been performed to assess their impact on urban temperatures and energy fluxes at the top of the canopy layer. In these simulations, it is assumed that the anthropogenic heat flux and latent heat fluxes are negligible. Results show that the simulated net radiation and sensible heat fluxes at the top of the canopy decrease and the stored heat increases as the aspect ratio increases. The simulated air temperature follows the behavior of the sensible heat flux. (C) 2010 Elsevier Ltd. All rights reserved.

Cloud streets and land-water interactions in the Amazon

Cloud streets are common feature in the Amazon Basin. They form from the combination of the vertical trade wind stress and moist convection. Here, satellite imagery, data collected during the COBRA-PARA (Caxiuan Observations in the Biosphere, River and Atmosphere of Para) field campaign, and high resolution modeling are used to understand the streets` formation and behavior. The observations show that the streets have an aspect ratio of about 3.5 and they reach their maximum activity around 15:00 UTC when the wind shear is weaker, and the convective boundary layer reaches its maximum height. The simulations reveal that the cloud streets onset is caused by the local circulations and convection produced at the interfaces between forest and rivers of the Amazon. The satellite data and modeling show that the large rivers anchor the cloud streets producing a quasi-stationary horizontal pattern. The streets are associated with horizontal roll vortices parallel to the mean flow that organizes the turbulence causing advection of latent heat flux towards the upward branches. The streets have multiple warm plumes that promote a connection between the rolls. These spatial patterns allow fundamental insights on the interpretation of the Amazon exchanges between surface and atmosphere with important consequences for the climate change understanding.

The energy and water balance of a Eucalyptus plantation in southeast Brazil

The eddy covariance method was used to measure energy and water balance of a plantation of Eucalyptus (grandis x urophylla) hybrids over a 2 year period. The average daily evaporation rates were 5.4 (+/- 2.0) mm day(-1) in summer, but fell to 1.2 (+/- 0.3) mm day(-1) in winter. In contrast, the sensible heat flux was relatively low in summer but dominated the energy balance in winter. Evaporation accounted for 80% and 26% of the available energy, in summer and winter respectively. The annual evaporation was 82% (1124 mm) and 96% (1235 mm) of the annual rainfall recorded during the first and second year, respectively. Daily average canopy and aerodynamic conductance to water vapour were in the summer 51.9 (+/- 38.4) mm s(-1) 84.1 (+/- 25.6) mm s(-1), respectively; and in the winter 6.0 (+/- 10.5) mm s(-1) and 111.6 (+/- 24.6) mm s(-1), respectively. (C) 2010 Elsevier B.V. All rights reserved.