Origin and impact of initialisation shocks in coupled atmosphere-ocean forecasts

Current methods for initialising coupled atmosphere-ocean forecasts often rely on the use of separate atmosphere and ocean analyses, the combination of which can leave the coupled system imbalanced at the beginning of the forecast, potentially accelerating the development of errors. Using a series of experiments with the European Centre for Medium-range Weather Forecasts coupled system, the magnitude and extent of these so-called initialisation shocks is quantified, and their impact on forecast skill measured. It is found that forecasts initialised by separate ocean and atmospheric analyses do exhibit initialisation shocks in lower atmospheric temperature, when compared to forecasts initialised using a coupled data assimilation method. These shocks result in as much as a doubling of root-mean-square error on the first day of the forecast in some regions, and in increases that are sustained for the duration of the 10-day forecasts performed here. However, the impacts of this choice of initialisation on forecast skill, assessed using independent datasets, were found to be negligible, at least over the limited period studied. Larger initialisation shocks are found to follow a change in either the atmospheric or ocean model component between the analysis and forecast phases: changes in the ocean component can lead to sea surface temperature shocks of more than 0.5K in some equatorial regions during the first day of the forecast. Implications for the development of coupled forecast systems, particularly with respect to coupled data assimilation methods, are discussed.

Metrics for linking emissions of gases and aerosols to global precipitation changes

Recent advances in understanding have made it possible to relate global precipitation changes directly to emissions of particular gases and aerosols that influence climate. Using these advances, new indices are developed here called the Global Precipitation-change Potential for pulse (GPP_P) and sustained (GPP_S) emissions, which measure the precipitation change per unit mass of emissions. The GPP can be used as a metric to compare the effects of different emissions. This is akin to the global warming potential (GWP) and the global temperature-change potential (GTP) which are used to place emissions on a common scale. Hence the GPP provides an additional perspective of the relative or absolute effects of emissions. It is however recognised that precipitation changes are predicted to be highly variable in size and sign between different regions and this limits the usefulness of a purely global metric. The GPP_P and GPP_S formulation consists of two terms, one dependent on the surface temperature change and the other dependent on the atmospheric component of the radiative forcing. For some forcing agents, and notably for CO2, these two terms oppose each other – as the forcing and temperature perturbations have different timescales, even the sign of the absolute GPP_P and GPP_S varies with time, and the opposing terms can make values sensitive to uncertainties in input parameters. This makes the choice of CO2 as a reference gas problematic, especially for the GPP_S at time horizons less than about 60 years. In addition, few studies have presented results for the surface/atmosphere partitioning of different forcings, leading to more uncertainty in quantifying the GPP than the GWP or GTP. Values of the GPP_P and GPP_S for five long- and short-lived forcing agents (CO2, CH4, N2O, sulphate and black carbon – BC) are presented, using illustrative values of required parameters. The resulting precipitation changes are given as the change at a specific time horizon (and hence they are end-point metrics) but it is noted that the GPPS can also be interpreted as the time-integrated effect of a pulse emission. Using CO2 as a references gas, the GPP_P and GPP_S for the non-CO2 species are larger than the corresponding GTP values. For BC emissions, the atmospheric forcing is sufficiently strong that the GPP_S is opposite in sign to the GTP_S. The sensitivity of these values to a number of input parameters is explored. The GPP can also be used to evaluate the contribution of different emissions to precipitation change during or after a period of emissions. As an illustration, the precipitation changes resulting from emissions in 2008 (using the GPP_P) and emissions sustained at 2008 levels (using the GPP_S) are presented. These indicate that for periods of 20 years (after the 2008 emissions) and 50 years (for sustained emissions at 2008 levels) methane is the dominant driver of positive precipitation changes due to those emissions. For sustained emissions, the sum of the effect of the five species included here does not become positive until after 50 years, by which time the global surface temperature increase exceeds 1 K.

Atmosphere-ocean coupled processes in the Madden-Julian oscillation

The Madden-Julian oscillation (MJO) is a convectively coupled 30-70 day (intraseasonal) tropical atmospheric mode that drives variations in global weather, but which is poorly simulated in most atmospheric general circulation models. Over the past two decades, field campaigns and modeling experiments have suggested that tropical atmosphere-ocean interactions may sustain or amplify the pattern of enhanced and suppressed atmospheric convection that defines the MJO, and encourage its eastward propagation through the Indian and Pacific Oceans. New observations collected during the past decade have advanced our understand of the ocean response to atmospheric MJO forcing and the resulting intraseasonal sea surface temperature (SST) fluctuations. Numerous modeling studies have revealed a considerable impact of the mean state on MJO ocean-atmosphere coupled processes, as well as the importance of resolving the diurnal cycle of atmosphere--upper-ocean interactions. New diagnostic methods provide insight to atmospheric variability and physical processes associated with the MJO, but offer limited insight on the role of ocean feedbacks. Consequently, uncertainty remains concerning the role of the ocean in MJO theory. Our understanding of how atmosphere-ocean coupled processes affect the MJO can be improved by collecting observations in poorly sampled regions of MJO activity, assessing oceanic and atmospheric drivers of surface fluxes, improving the representation of upper-ocean mixing in coupled-model simulations, designing model experiments that minimize mean-state differences, and developing diagnostic tools to evaluate the nature and role of coupled ocean-atmosphere processes over the MJO cycle.

'Where water wells up from the earth’: excavations at the findspot of the Late Bronze Age hoard from Broadward, Shropshire

The paper begins by considering the importance of springs as a focus for votive deposits in Bronze Age Britain. This is not a new idea, but nowhere has this association been examined through the excavation of one of these features. The point is illustrated by excavation at the findspot of a famous group of Late Bronze Age weapons, the Broadward hoard, discovered in 1867. Little was known about the site, where it was found or the character of the original deposit, but a study of contemporary accounts of the hoard, combined with geophysical and topographical surveys, led to small-scale excavation in 2010, which showed that the deposit had most probably been buried in a pit on the edge of a spring. Other finds associated with the spring included an Early Bronze Age macehead, a Roman pot and various Saxon and medieval animal bones. The latest deposit, with a post-medieval carbon date, included a wooden knife or dagger. An adjacent palaeochannel provided an important environmental sequence for this part of the English–Welsh borderland and suggests that the Late Bronze Age hoard had been deposited not far from a settlement. A nearby earthwork enclosure was associated with a clay weight, which may be of similar date. Despite the limited scale of the fieldwork, it illustrates the potential for treating springs associated with artefact finds on the same terms as other archaeological deposits.

Time-lapse geophysical investigations over known archaeological features using electrical resistivity imaging and earth resistance

Electrical methods of geophysical survey are known to produce results that are hard to predict at different times of the year, and under differing weather conditions. This is a problem which can lead to misinterpretation of archaeological features under investigation. The dynamic relationship between a ‘natural’ soil matrix and an archaeological feature is a complex one, which greatly affects the success of the feature’s detection when using active electrical methods of geophysical survey. This study has monitored the gradual variation of measured resistivity over a selection of study areas. By targeting difficult to find, and often ‘missing’ electrical anomalies of known archaeological features, this study has increased the understanding of both the detection and interpretation capabilities of such geophysical surveys. A 16 month time-lapse study over 4 archaeological features has taken place to investigate the aforementioned detection problem across different soils and environments. In addition to the commonly used Twin-Probe earth resistance survey, electrical resistivity imaging (ERI) and quadrature electro-magnetic induction (EMI) were also utilised to explore the problem. Statistical analyses have provided a novel interpretation, which has yielded new insights into how the detection of archaeological features is influenced by the relationship between the target feature and the surrounding ‘natural’ soils. The study has highlighted both the complexity and previous misconceptions around the predictability of the electrical methods. The analysis has confirmed that each site provides an individual and nuanced situation, the variation clearly relating to the composition of the soils (particularly pore size) and the local weather history. The wide range of reasons behind survey success at each specific study site has been revealed. The outcomes have shown that a simplistic model of seasonality is not universally applicable to the electrical detection of archaeological features. This has led to the development of a method for quantifying survey success, enabling a deeper understanding of the unique way in which each site is affected by the interaction of local environmental and geological conditions.

Earth hummocks in north-east Okstindan, northern Norway: morphology, distribution and environmental constraints

Earth hummocks (also termed pounus or thúfur) are a common form of periglacial non-sorted patterned ground. The study objectives were to determine the morphology, distribution and development on slopes of earth hummocks in north-east Okstindan, Norway, an area with many hummocks but few documented accounts. The methodology involved detailed geomorphological mapping and precise measurement with a profileometer. The internal structure of the hummocks was investigated through excavations and sediment sample analyses. Fourteen sites with well-developed earth hummocks (accounting for over 650 individual hummock forms) were investigated. The sites have an average altitude of 750 m and occur on slopes with an average gradient of 7°. The hummock heights are in the range 0.11–0.52 m and their diameters 0.7–1.5 m, although coalescent forms are up to 5 m in length. The hummock morphology is characterised by a variable plan form, asymmetry with respect to upslope and downslope forms, downslope elongation, coalescence, and superimposed microtopography. The hummocks’ distribution appeared to have been controlled by the existence of a frost-susceptible ‘host’ sediment, but moisture availability and topographic position played a role. The authors conclude that differential frost heave and vegetation cover stability are critical for the hummocks’ longevity in the studied landscape.

Decadal predictability of regional-scale peak winds over Europe using the Earth System Model of the Max-Planck-Institute for Meteorology

The predictability of high impact weather events on multiple time scales is a crucial issue both in scientific and socio-economic terms. In this study, a statistical-dynamical downscaling (SDD) approach is applied to an ensemble of decadal hindcasts obtained with the Max-Planck-Institute Earth System Model (MPI-ESM) to estimate the decadal predictability of peak wind speeds (as a proxy for gusts) over Europe. Yearly initialized decadal ensemble simulations with ten members are investigated for the period 1979–2005. The SDD approach is trained with COSMO-CLM regional climate model simulations and ERA-Interim reanalysis data and applied to the MPI-ESM hindcasts. The simulations for the period 1990–1993, which was characterized by several windstorm clusters, are analyzed in detail. The anomalies of the 95 % peak wind quantile of the MPI-ESM hindcasts are in line with the positive anomalies in reanalysis data for this period. To evaluate both the skill of the decadal predictability system and the added value of the downscaling approach, quantile verification skill scores are calculated for both the MPI-ESM large-scale wind speeds and the SDD simulated regional peak winds. Skill scores are predominantly positive for the decadal predictability system, with the highest values for short lead times and for (peak) wind speeds equal or above the 75 % quantile. This provides evidence that the analyzed hindcasts and the downscaling technique are suitable for estimating wind and peak wind speeds over Central Europe on decadal time scales. The skill scores for SDD simulated peak winds are slightly lower than those for large-scale wind speeds. This behavior can be largely attributed to the fact that peak winds are a proxy for gusts, and thus have a higher variability than wind speeds. The introduced cost-efficient downscaling technique has the advantage of estimating not only wind speeds but also estimates peak winds (a proxy for gusts) and can be easily applied to large ensemble datasets like operational decadal prediction systems.

Current model capabilities for simulating black carbon and sulfate concentrations in the Arctic atmosphere: a multi-model evaluation using a comprehensive measurement data set

The concentrations of sulfate, black carbon (BC) and other aerosols in the Arctic are characterized by high values in late winter and spring (so-called Arctic Haze) and low values in summer. Models have long been struggling to capture this seasonality and especially the high concentrations associated with Arctic Haze. In this study, we evaluate sulfate and BC concentrations from eleven different models driven with the same emission inventory against a comprehensive pan-Arctic measurement data set over a time period of 2 years (2008–2009). The set of models consisted of one Lagrangian particle dispersion model, four chemistry transport models (CTMs), one atmospheric chemistry-weather forecast model and five chemistry climate models (CCMs), of which two were nudged to meteorological analyses and three were running freely. The measurement data set consisted of surface measurements of equivalent BC (eBC) from five stations (Alert, Barrow, Pallas, Tiksi and Zeppelin), elemental carbon (EC) from Station Nord and Alert and aircraft measurements of refractory BC (rBC) from six different campaigns. We find that the models generally captured the measured eBC or rBC and sulfate concentrations quite well, compared to previous comparisons. However, the aerosol seasonality at the surface is still too weak in most models. Concentrations of eBC and sulfate averaged over three surface sites are underestimated in winter/spring in all but one model (model means for January–March underestimated by 59 and 37 % for BC and sulfate, respectively), whereas concentrations in summer are overestimated in the model mean (by 88 and 44 % for July–September), but with overestimates as well as underestimates present in individual models. The most pronounced eBC underestimates, not included in the above multi-site average, are found for the station Tiksi in Siberia where the measured annual mean eBC concentration is 3 times higher than the average annual mean for all other stations. This suggests an underestimate of BC sources in Russia in the emission inventory used. Based on the campaign data, biomass burning was identified as another cause of the modeling problems. For sulfate, very large differences were found in the model ensemble, with an apparent anti-correlation between modeled surface concentrations and total atmospheric columns. There is a strong correlation between observed sulfate and eBC concentrations with consistent sulfate/eBC slopes found for all Arctic stations, indicating that the sources contributing to sulfate and BC are similar throughout the Arctic and that the aerosols are internally mixed and undergo similar removal. However, only three models reproduced this finding, whereas sulfate and BC are weakly correlated in the other models. Overall, no class of models (e.g., CTMs, CCMs) performed better than the others and differences are independent of model resolution.

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.

The impact of oceanic heat transport on the atmospheric circulation

A general circulation model of intermediate complexity with an idealized Earth-like aquaplanet setup is used to study the impact of changes in the oceanic heat transport on the global atmospheric circulation. Focus is on the atmospheric mean meridional circulation and global thermodynamic properties. The atmosphere counterbalances to a large extent the imposed changes in the oceanic heat transport, but, nonetheless, significant modifications to the atmospheric general circulation are found. Increasing the strength of the oceanic heat transport up to 2.5 PW leads to an increase in the global mean near-surface temperature and to a decrease in its equator-to-pole gradient. For stronger transports, the gradient is reduced further, but the global mean remains approximately constant. This is linked to a cooling and a reversal of the temperature gradient in the tropics. Additionally, a stronger oceanic heat transport leads to a decline in the intensity and a poleward shift of the maxima of both the Hadley and Ferrel cells. Changes in zonal mean diabatic heating and friction impact the properties of the Hadley cell, while the behavior of the Ferrel cell is mostly controlled by friction. The efficiency of the climate machine, the intensity of the Lorenz energy cycle and the material entropy production of the system decline with increased oceanic heat transport. This suggests that the climate system becomes less efficient and turns into a state of reduced entropy production as the enhanced oceanic transport performs a stronger large-scale mixing between geophysical fluids with different temperatures, thus reducing the available energy in the climate system and bringing it closer to a state of thermal equilibrium.

Climate change projections using the IPSL-CM5 Earth System Model: from CMIP3 to CMIP5

We present the global general circulation model IPSL-CM5 developed to study the long-term response of the climate system to natural and anthropogenic forcings as part of the 5th Phase of the Coupled Model Intercomparison Project (CMIP5). This model includes an interactive carbon cycle, a representation of tropospheric and stratospheric chemistry, and a comprehensive representation of aerosols. As it represents the principal dynamical, physical, and bio-geochemical processes relevant to the climate system, it may be referred to as an Earth System Model. However, the IPSL-CM5 model may be used in a multitude of configurations associated with different boundary conditions and with a range of complexities in terms of processes and interactions. This paper presents an overview of the different model components and explains how they were coupled and used to simulate historical climate changes over the past 150 years and different scenarios of future climate change. A single version of the IPSL-CM5 model (IPSL-CM5A-LR) was used to provide climate projections associated with different socio-economic scenarios, including the different Representative Concentration Pathways considered by CMIP5 and several scenarios from the Special Report on Emission Scenarios considered by CMIP3. Results suggest that the magnitude of global warming projections primarily depends on the socio-economic scenario considered, that there is potential for an aggressive mitigation policy to limit global warming to about two degrees, and that the behavior of some components of the climate system such as the Arctic sea ice and the Atlantic Meridional Overturning Circulation may change drastically by the end of the twenty-first century in the case of a no climate policy scenario. Although the magnitude of regional temperature and precipitation changes depends fairly linearly on the magnitude of the projected global warming (and thus on the scenario considered), the geographical pattern of these changes is strikingly similar for the different scenarios. The representation of atmospheric physical processes in the model is shown to strongly influence the simulated climate variability and both the magnitude and pattern of the projected climate changes.

Ions in the atmosphere

In Earth’s atmosphere, an ion is a cluster of molecules carrying an overall charge, known as a molecular cluster ion. Such cluster ions, with dimensions of approximately one nanometre, have usually been referred to as small ions, and their motion in air constitutes a small electric current. Large ions (or Langevin ions), by comparison, are physically larger (tens to hundreds of nm) and consequently electrically less mobile. Usage of the term “ion” to represent these molecular clusters originates from the early history of atmospheric electricity, which spans the discovery of the electron and the elucidation of the structure of matter. The distinction between large and small ions originates from distinguishing ions that could be accelerated by atmospheric electric fields (and therefore directly contribute to the conductivity of air), and those (the large ions) which were insufficiently electrically mobile to contribute to electrical conduction in air.

Cosmic ray measurements in the atmosphere at several latitudes in October, 2014

Cosmic ray fluxes in the atmosphere were recorded during balloon flights in October 2014 in northern Murmansk region, Apatity (Russia; 67o33’N, 33o24’E), in Antarctica (observatory Mirny; 66o33’S, 93o00’E), in Moscow (Russia; 55o45’N, 37o37’E), in Reading (United King-dom; 51o27’N, 0o 58’W), in Mitzpe-Ramon (Israel; 30o36’N, 34o48’E) and in Zaragoza (Spain; 41o9’N, 0o54’W). Two type of cosmic ray detectors were used, namely, (1) the standard ra-diosonde and its modification constructed at the Lebedev Physical Institute (Moscow, Russia) and (2) the device manufactured at the Reading University (Reading, United Kingdom). We compare and analyze obtained data and focus on the estimation of the cosmic ray latitudinal effect in the atmosphere.

Selecting algorithms for Earth observation of climate within the European Space Agency Climate Change initiative: introduction to a special issue

This special issue is focused on the assessment of algorithms for the observation of Earth’s climate from environ- mental satellites. Climate data records derived by remote sensing are increasingly a key source of insight into the workings of and changes in Earth’s climate system. Producers of data sets must devote considerable effort and expertise to maximise the true climate signals in their products and minimise effects of data processing choices and changing sensors. A key choice is the selection of algorithm(s) for classification and/or retrieval of the climate variable. Within the European Space Agency Climate Change Initiative, science teams undertook systematic assessment of algorithms for a range of essential climate variables. The papers in the special issue report some of these exercises (for ocean colour, aerosol, ozone, greenhouse gases, clouds, soil moisture, sea surface temper- ature and glaciers). The contributions show that assessment exercises must be designed with care, considering issues such as the relative importance of different aspects of data quality (accuracy, precision, stability, sensitivity, coverage, etc.), the availability and degree of independence of validation data and the limitations of validation in characterising some important aspects of data (such as long-term stability or spatial coherence). As well as re- quiring a significant investment of expertise and effort, systematic comparisons are found to be highly valuable. They reveal the relative strengths and weaknesses of different algorithmic approaches under different observa- tional contexts, and help ensure that scientific conclusions drawn from climate data records are not influenced by observational artifacts, but are robust.

The EarthCARE satellite: the next step forward in global measurements of clouds, aerosols, precipitation, and radiation

The collective representation within global models of aerosol, cloud, precipitation, and their radiative properties remains unsatisfactory. They constitute the largest source of uncertainty in predictions of climatic change and hamper the ability of numerical weather prediction models to forecast high-impact weather events. The joint European Space Agency (ESA)–Japan Aerospace Exploration Agency (JAXA) Earth Clouds, Aerosol and Radiation Explorer (EarthCARE) satellite mission, scheduled for launch in 2018, will help to resolve these weaknesses by providing global profiles of cloud, aerosol, precipitation, and associated radiative properties inferred from a combination of measurements made by its collocated active and passive sensors. EarthCARE will improve our understanding of cloud and aerosol processes by extending the invaluable dataset acquired by the A-Train satellites CloudSat, Cloud–Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO), and Aqua. Specifically, EarthCARE’s cloud profiling radar, with 7 dB more sensitivity than CloudSat, will detect more thin clouds and its Doppler capability will provide novel information on convection, precipitating ice particle, and raindrop fall speeds. EarthCARE’s 355-nm high-spectral-resolution lidar will measure directly and accurately cloud and aerosol extinction and optical depth. Combining this with backscatter and polarization information should lead to an unprecedented ability to identify aerosol type. The multispectral imager will provide a context for, and the ability to construct, the cloud and aerosol distribution in 3D domains around the narrow 2D retrieved cross section. The consistency of the retrievals will be assessed to within a target of ±10 W m–2 on the (10 km)2 scale by comparing the multiview broadband radiometer observations to the top-of-atmosphere fluxes estimated by 3D radiative transfer models acting on retrieved 3D domains.