Consistent and contrasting decadal Arctic sea ice thickness predictions from a highly optimized sea ice model

[1] Decadal hindcast simulations of Arctic Ocean sea ice thickness made by a modern dynamic-thermodynamic sea ice model and forced independently by both the ERA-40 and NCEP/NCAR reanalysis data sets are compared for the first time. Using comprehensive data sets of observations made between 1979 and 2001 of sea ice thickness, draft, extent, and speeds, we find that it is possible to tune model parameters to give satisfactory agreement with observed data, thereby highlighting the skill of modern sea ice models, though the parameter values chosen differ according to the model forcing used. We find a consistent decreasing trend in Arctic Ocean sea ice thickness since 1979, and a steady decline in the Eastern Arctic Ocean over the full 40-year period of comparison that accelerated after 1980, but the predictions of Western Arctic Ocean sea ice thickness between 1962 and 1980 differ substantially. The origins of differing thickness trends and variability were isolated not to parameter differences but to differences in the forcing fields applied, and in how they are applied. It is argued that uncertainty, differences and errors in sea ice model forcing sets complicate the use of models to determine the exact causes of the recently reported decline in Arctic sea ice thickness, but help in the determination of robust features if the models are tuned appropriately against observations.

Can polar lows be objectively identified and tracked in the ECMWF operational analysis and the ERA-Interim reanalysis?

Polar lows are maritime meso-cyclones associated with intense surface wind speeds and oceanic heat fluxes at high latitudes. The ability of the ERA-Interim (ERAI) reanalysis to represent polar lows in the North Atlantic is assessed by comparing ERAI and the ECMWF operational analysis for the period 2008-2011. First, the representation of a set of satellite observed polar lows over the Norwegian and Barents Seas in the operational analysis and ERAI is analysed. Then, the possibility of directly identifying and tracking the polar lows in the operational analysis and ERAI is explored using a tracking algorithm based on 850 hPa vorticity with objective identification criteria on cyclone dynamical intensity and atmospheric static stability. All but one of the satellite observed polar lows with a lifetime of at least 6 hours have an 850 hPa vorticity signature of a co-located mesocyclone in both the operational analysis and ERAI for most of their life cycles. However, the operational analysis has vorticity structures that better resemble the observed cloud patterns and stronger surface wind speed intensities compared to those in ERAI. By applying the objective identification criteria, about 55% of the satellite observed polar lows are identified and tracked in ERAI, while this fraction increases to about 70% in the operational analysis. Particularly in ERAI, the remaining observed polar lows are mainly not identified because they have too weak wind speed and vorticity intensity compared to the tested criteria. The implications of the tendency of ERAI to underestimate the polar low dynamical intensity for future studies of polar lows is discussed.

Evaluation of ERA-Interim reanalysis precipitation products using England and Wales observations

Precipitation forecast data from the ERA-Interim reanalysis (33 years) are evaluated using the daily England and Wales Precipitation (EWP) observations obtained from a rain gauge network. Observed and reanalysis daily precipitation data are both described well by Weibull distributions with indistinguishable shapes but different scale parameters, such that the reanalysis underestimates the observations by an average factor of 22%. The correlation between the observed and ERA-Interim time series of regional, daily precipitation is 0.91. ERA-Interim also captures the statistics of extreme precipitation including a slightly lower likelihood of the heaviest precipitation events (>15 mm day− 1 for the regional average) than indicated by the Weibull fit. ERA-Interim is also closer to EWP for the high precipitation events. Since these carry weight in longer accumulations, a smaller underestimation of 19% is found for monthly mean precipitation. The partition between convective and stratiform precipitation in the ERA-Interim forecast is also examined. In summer both components contribute equally to the total precipitation amount, while in winter the stratiform precipitation is approximately double convective. These results are expected to be relevant to other regions with low orography on the coast of a continent at the downstream end of mid-latitude stormtracks.

An aggregated wind power generation model based on MERRA reanalysis data: MATLAB model and example data for the April 2014 wind farm distribution of Great Britain.

The MATLAB model is contained within the compressed folders (versions are available as .zip and .tgz). This model uses MERRA reanalysis data (>34 years available) to estimate the hourly aggregated wind power generation for a predefined (fixed) distribution of wind farms. A ready made example is included for the wind farm distribution of Great Britain, April 2014 ("CF.dat"). This consists of an hourly time series of GB-total capacity factor spanning the period 1980-2013 inclusive. Given the global nature of reanalysis data, the model can be applied to any specified distribution of wind farms in any region of the world. Users are, however, strongly advised to bear in mind the limitations of reanalysis data when using this model/data. This is discussed in our paper: Cannon, Brayshaw, Methven, Coker, Lenaghan. "Using reanalysis data to quantify extreme wind power generation statistics: a 33 year case study in Great Britain". Submitted to Renewable Energy in March, 2014. Additional information about the model is contained in the model code itself, in the accompanying ReadMe file, and on our website: http://www.met.reading.ac.uk/~energymet/data/Cannon2014/

Using reanalysis data to quantify extreme wind power generation statistics : a 33 year case study in Great Britain

With a rapidly increasing fraction of electricity generation being sourced from wind, extreme wind power generation events such as prolonged periods of low (or high) generation and ramps in generation, are a growing concern for the efficient and secure operation of national power systems. As extreme events occur infrequently, long and reliable meteorological records are required to accurately estimate their characteristics. Recent publications have begun to investigate the use of global meteorological “reanalysis” data sets for power system applications, many of which focus on long-term average statistics such as monthly-mean generation. Here we demonstrate that reanalysis data can also be used to estimate the frequency of relatively short-lived extreme events (including ramping on sub-daily time scales). Verification against 328 surface observation stations across the United Kingdom suggests that near-surface wind variability over spatiotemporal scales greater than around 300 km and 6 h can be faithfully reproduced using reanalysis, with no need for costly dynamical downscaling. A case study is presented in which a state-of-the-art, 33 year reanalysis data set (MERRA, from NASA-GMAO), is used to construct an hourly time series of nationally-aggregated wind power generation in Great Britain (GB), assuming a fixed, modern distribution of wind farms. The resultant generation estimates are highly correlated with recorded data from National Grid in the recent period, both for instantaneous hourly values and for variability over time intervals greater than around 6 h. This 33 year time series is then used to quantify the frequency with which different extreme GB-wide wind power generation events occur, as well as their seasonal and inter-annual variability. Several novel insights into the nature of extreme wind power generation events are described, including (i) that the number of prolonged low or high generation events is well approximated by a Poission-like random process, and (ii) whilst in general there is large seasonal variability, the magnitude of the most extreme ramps is similar in both summer and winter. An up-to-date version of the GB case study data as well as the underlying model are freely available for download from our website: http://www.met.reading.ac.uk/~energymet/data/Cannon2014/.

The WFDEI meteorological forcing data set: WATCH Forcing Data methodology applied to ERA-Interim reanalysis data

The WFDEI meteorological forcing data set has been generated using the same methodology as the widely used WATCH Forcing Data (WFD) by making use of the ERA-Interim reanalysis data. We discuss the specifics of how changes in the reanalysis and processing have led to improvement over the WFD. We attribute improvements in precipitation and wind speed to the latest reanalysis basis data and improved downward shortwave fluxes to the changes in the aerosol corrections. Covering 1979–2012, the WFDEI will allow more thorough comparisons of hydrological and Earth System model outputs with hydrologically and phenologically relevant satellite products than using the WFD.

ERA-Interim/Land: a global land surface reanalysis data set

ERA-Interim/Land is a global land surface reanalysis data set covering the period 1979–2010. It describes the evolution of soil moisture, soil temperature and snowpack. ERA-Interim/Land is the result of a single 32-year simulation with the latest ECMWF (European Centre for Medium-Range Weather Forecasts) land surface model driven by meteorological forcing from the ERA-Interim atmospheric reanalysis and precipitation adjustments based on monthly GPCP v2.1 (Global Precipitation Climatology Project). The horizontal resolution is about 80 km and the time frequency is 3-hourly. ERA-Interim/Land includes a number of parameterization improvements in the land surface scheme with respect to the original ERA-Interim data set, which makes it more suitable for climate studies involving land water resources. The quality of ERA-Interim/Land is assessed by comparing with ground-based and remote sensing observations. In particular, estimates of soil moisture, snow depth, surface albedo, turbulent latent and sensible fluxes, and river discharges are verified against a large number of site measurements. ERA-Interim/Land provides a global integrated and coherent estimate of soil moisture and snow water equivalent, which can also be used for the initialization of numerical weather prediction and climate models.

The effects of springtime mid-latitude storms on trace gas composition determined from the MACC reanalysis

The relationship between springtime air pollution transport of ozone (O3) and carbon monoxide (CO) and mid-latitude cyclones is explored for the first time using the Monitoring Atmospheric Composition and Climate (MACC) reanalysis for the period 2003–2012. In this study, the most intense spring storms (95th percentile) are selected for two regions, the North Pacific (NP) and the North Atlantic (NA). These storms (∼60 storms over each region) often track over the major emission sources of East Asia and eastern North America. By compositing the storms, the distributions of O3 and CO within a "typical" intense storm are examined. We compare the storm-centered composite to background composites of "average conditions" created by sampling the reanalysis data of the previous year to the storm locations. Mid-latitude storms are found to redistribute concentrations of O3 and CO horizontally and vertically throughout the storm. This is clearly shown to occur through two main mechanisms: (1) vertical lifting of CO-rich and O3-poor air isentropically, from near the surface to the mid- to upper-troposphere in the region of the warm conveyor belt; and (2) descent of O3-rich and CO-poor air isentropically in the vicinity of the dry intrusion, from the stratosphere toward the mid-troposphere. This can be seen in the composite storm's life cycle as the storm intensifies, with area-averaged O3 (CO) increasing (decreasing) between 200 and 500 hPa. The influence of the storm dynamics compared to the background environment on the composition within an area around the storm center at the time of maximum intensity is as follows. Area-averaged O3 at 300 hPa is enhanced by 50 and 36% and by 11 and 7.6% at 500 hPa for the NP and NA regions, respectively. In contrast, area-averaged CO at 300 hPa decreases by 12% for NP and 5.5% for NA, and area-averaged CO at 500 hPa decreases by 2.4% for NP while there is little change over the NA region. From the mid-troposphere, O3-rich air is clearly seen to be transported toward the surface, but the downward transport of CO-poor air is not discernible due to the high levels of CO in the lower troposphere. Area-averaged O3 is slightly higher at 1000 hPa (3.5 and 1.8% for the NP and NA regions, respectively). There is an increase of CO at 1000 hPa for the NP region (3.3%) relative to the background composite and a~slight decrease in area-averaged CO for the NA region at 1000 hPa (-2.7%).

Combining satellite observations and reanalysis energy transports to estimate global net surface energy fluxes 1985-2012

Two methods are developed to estimate net surface energy fluxes based upon satellite-based reconstructions of radiative fluxes at the top of atmosphere and the atmospheric energy tendencies and transports from the ERA-Interim reanalysis. Method 1 applies the mass adjusted energy divergence from ERA-Interim while method 2 estimates energy divergence based upon the net energy difference at the top of atmosphere and the surface from ERA-Interim. To optimise the surface flux and its variability over ocean, the divergences over land are constrained to match the monthly area mean surface net energy flux variability derived from a simple relationship between the surface net energy flux and the surface temperature change. The energy divergences over the oceans are then adjusted to remove an unphysical residual global mean atmospheric energy divergence. The estimated net surface energy fluxes are compared with other data sets from reanalysis and atmospheric model simulations. The spatial correlation coefficients of multi-annual means between the estimations made here and other data sets are all around 0.9. There are good agreements in area mean anomaly variability over the global ocean, but discrepancies in the trend over the eastern Pacific are apparent.

Regionally coupled atmosphere-ocean-sea ice-marine biogeochemistry model ROM: 1. Description and validation

The general circulation models used to simulate global climate typically feature resolution too coarse to reproduce many smaller-scale processes, which are crucial to determining the regional responses to climate change. A novel approach to downscale climate change scenarios is presented which includes the interactions between the North Atlantic Ocean and the European shelves as well as their impact on the North Atlantic and European climate. The goal of this paper is to introduce the global ocean-regional atmosphere coupling concept and to show the potential benefits of this model system to simulate present-day climate. A global ocean-sea ice-marine biogeochemistry model (MPIOM/HAMOCC) with regionally high horizontal resolution is coupled to an atmospheric regional model (REMO) and global terrestrial hydrology model (HD) via the OASIS coupler. Moreover, results obtained with ROM using NCEP/NCAR reanalysis and ECHAM5/MPIOM CMIP3 historical simulations as boundary conditions are presented and discussed for the North Atlantic and North European region. The validation of all the model components, i.e., ocean, atmosphere, terrestrial hydrology, and ocean biogeochemistry is performed and discussed. The careful and detailed validation of ROM provides evidence that the proposed model system improves the simulation of many aspects of the regional climate, remarkably the ocean, even though some biases persist in other model components, thus leaving potential for future improvement. We conclude that ROM is a powerful tool to estimate possible impacts of climate change on the regional scale.

The solsticial pause on Mars: 1. A planetary wave reanalysis

Large-scale planetary waves are diagnosed from an analysis of profiles retrieved from the Thermal Emission Spectrometer aboard the Mars Global Surveyor spacecraft during its scientific mapping phase. The analysis is conducted by assimilating thermal profiles and total dust opacity retrievals into a Mars global circulation model. Transient waves are largest throughout the northern hemisphere autumn, winter and spring period and almost absent during the summer. The southern hemisphere exhibits generally weaker transient wave behaviour. A striking feature of the low-altitude transient waves in the analysis is that they show a broad subsidiary minimum in amplitude centred on the winter solstice, a period when the thermal contrast between the summer hemisphere and the winter pole is strongest and baroclinic wave activity might be expected to be strong. This behaviour, here called the ‘solsticial pause,’ is present in every year of the analysis. This strong pause is under-represented in many independent model experiments, which tend to produce relatively uniform baroclinic wave activity throughout the winter. This paper documents and diagnoses the transient wave solsticial pause found in the analysis; a companion paper investigates the origin of the phenomenon in a series of model experiments.

The Ocean Reanalysis Intercomparison project (ORA-IP)

Uncertainty in ocean analysis methods and deficiencies in the observing system are major obstacles for the reliable reconstruction of the past ocean climate. The variety of existing ocean reanalyses is exploited in a multi-reanalysis ensemble to improve the ocean state estimation and to gauge uncertainty levels. The ensemble-based analysis of signal-to-noise ratio allows the identification of ocean characteristics for which the estimation is robust (such as tropical mixed-layer-depth, upper ocean heat content), and where large uncertainty exists (deep ocean, Southern Ocean, sea ice thickness, salinity), providing guidance for future enhancement of the observing and data assimilation systems.

The role of local atmospheric forcing on the modulation of the ocean mixed layer depth in reanalysis and a coupled single column ocean model

The role of the local atmospheric forcing on the ocean mixed layer depth (MLD) over the global oceans is studied using ocean reanalysis data products and a single-column ocean model coupled to an atmospheric general circulation model. The focus of this study is on how the annual mean and the seasonal cycle of the MLD relate to various forcing characteristics in different parts of the world's ocean, and how anomalous variations in the monthly mean MLD relate to anomalous atmospheric forcings. By analysing both ocean reanalysis data and the single-column ocean model, regions with different dominant forcings and different mean and variability characteristics of the MLD can be identified. Many of the global oceans' MLD characteristics appear to be directly linked to different atmospheric forcing characteristics at different locations. Here, heating and wind-stress are identified as the main drivers; in some, mostly coastal, regions the atmospheric salinity forcing also contributes. The annual mean MLD is more closely related to the annual mean wind-stress and the MLD seasonality is more closely to the seasonality in heating. The single-column ocean model, however, also points out that the MLD characteristics over most global ocean regions, and in particular the tropics and subtropics, cannot be maintained by local atmospheric forcings only, but are also a result of ocean dynamics that are not simulated in a single-column ocean model. Thus, lateral ocean dynamics are essentially in correctly simulating observed MLD.

Climatology of polar lows over the Sea of Japan using the JRA-55 reanalysis

Polar lows are intense meso- a -scale cyclones that develop over the oceans poleward of the main baroclinic zone. A number of previous studies have reported polar low formation over the Sea of Japan within the East Asian winter monsoon. To understand the climatology of polar lows over the Sea of Japan, a tracking al- gorithm for polar lows is applied to the recent JRA-55 reanalysis. The polar low tracking is applied to 36 cold seasons (October–March) from October 1979 to March 2015. The polar lows over the Sea of Japan reach their maximum intensity on the southeastern side of the midline between the Japanese islands and the Asian continent. Consistent with previous case studies, composite analysis demonstrates that the polar low devel- opment is associated with the enhanced northerly flow on the western side of a synoptic-scale extratropical cyclone, with the cold trough in the midtroposphere and with increased heat fluxes from the sea surface. Furthermore, the present climatological study has revealed two dominant directions of motion of the polar lows: southward and eastward. Southward-moving polar lows are steered by a strong northerly flow in the lower troposphere, which is enhanced on the western side of synoptic-scale extratropical cyclones, while the eastward-moving polar lows occur within a planetary-scale westerly flow in the midlatitudes. Thus, the di- rection of polar low motion reflects the difference in planetary- and synoptic-scale conditions.