Mechanisms of the negative shortwave cloud feedback in mid to high latitudes

Increases in cloud optical depth and liquid water path (LWP) are robust features of global warming model simulations in high latitudes, yielding a negative shortwave cloud feedback, but the mechanisms are still uncertain. We assess the importance of microphysical processes for the negative optical depth feedback by perturbing temperature in the microphysics schemes of two aquaplanet models, both of which have separate prognostic equations for liquid water and ice. We find that most of the LWP increase with warming is caused by a suppression of ice microphysical processes in mixed-phase clouds, resulting in reduced conversion efficiencies of liquid water to ice and precipitation. Perturbing the temperature-dependent phase partitioning of convective condensate also yields a small LWP increase. Together, the perturbations in large-scale microphysics and convective condensate partitioning explain more than two-thirds of the LWP response relative to a reference case with increased SSTs, and capture all of the vertical structure of the liquid water response. In support of these findings, we show the existence of a very robust positive relationship between monthly-mean LWP and temperature in CMIP5 models and observations in mixed-phase cloud regions only. In models, the historical LWP sensitivity to temperature is a good predictor of the forced global warming response poleward of about 45°, although models appear to overestimate the LWP response to warming compared to observations. We conclude that in climate models, the suppression of ice-phase microphysical processes that deplete cloud liquid water is a key driver of the LWP increase with warming and of the associated negative shortwave cloud feedback.

Regulation of early steps of GPVI signal transduction by phosphatases: a systems biology approach

We present a data-driven mathematical model of a key initiating step in platelet activation, a central process in the prevention of bleeding following Injury. In vascular disease, this process is activated inappropriately and causes thrombosis, heart attacks and stroke. The collagen receptor GPVI is the primary trigger for platelet activation at sites of injury. Understanding the complex molecular mechanisms initiated by this receptor is important for development of more effective antithrombotic medicines. In this work we developed a series of nonlinear ordinary differential equation models that are direct representations of biological hypotheses surrounding the initial steps in GPVI-stimulated signal transduction. At each stage model simulations were compared to our own quantitative, high-temporal experimental data that guides further experimental design, data collection and model refinement. Much is known about the linear forward reactions within platelet signalling pathways but knowledge of the roles of putative reverse reactions are poorly understood. An initial model, that includes a simple constitutively active phosphatase, was unable to explain experimental data. Model revisions, incorporating a complex pathway of interactions (and specifically the phosphatase TULA-2), provided a good description of the experimental data both based on observations of phosphorylation in samples from one donor and in those of a wider population. Our model was used to investigate the levels of proteins involved in regulating the pathway and the effect of low GPVI levels that have been associated with disease. Results indicate a clear separation in healthy and GPVI deficient states in respect of the signalling cascade dynamics associated with Syk tyrosine phosphorylation and activation. Our approach reveals the central importance of this negative feedback pathway that results in the temporal regulation of a specific class of protein tyrosine phosphatases in controlling the rate, and therefore extent, of GPVI-stimulated platelet activation.

Climate-driven expansion of blanket bogs in Britain during the Holocene

Blanket bog occupies approximately 6 % of the area of the UK today. The Holocene expansion of this hyperoceanic biome has previously been explained as a consequence of Neolithic forest clearance. However, the present distribution of blanket bog in Great Britain can be predicted accurately with a simple model (PeatStash) based on summer temperature and moisture index thresholds, and the same model correctly predicts the highly disjunct distribution of blanket bog worldwide. This finding suggests that climate, rather than land-use history, controls blanket-bog distribution in the UK and everywhere else. We set out to test this hypothesis for blanket bogs in the UK using bioclimate envelope modelling compared with a database of peat initiation age estimates. We used both pollen-based reconstructions and climate model simulations of climate changes between the mid-Holocene (6000 yr BP, 6 ka) and modern climate to drive PeatStash and predict areas of blanket bog. We compiled data on the timing of blanket-bog initiation, based on 228 age determinations at sites where peat directly overlies mineral soil. The model predicts large areas of northern Britain would have had blanket bog by 6000 yr BP, and the area suitable for peat growth extended to the south after this time. A similar pattern is shown by the basal peat ages and new blanket bog appeared over a larger area during the late Holocene, the greatest expansion being in Ireland, Wales and southwest England, as the model predicts. The expansion was driven by a summer cooling of about 2 °C, shown by both pollen-based reconstructions and climate models. The data show early Holocene (pre-Neolithic) blanket-bog initiation at over half of the sites in the core areas of Scotland, and northern England. The temporal patterns and concurrence of the bioclimate model predictions and initiation data suggest that climate change provides a parsimonious explanation for the early Holocene distribution and later expansion of blanket bogs in the UK, and it is not necessary to invoke anthropogenic activity as a driver of this major landscape change.

Robust comparison of climate models with observations using blended land air and ocean sea surface temperatures

The level of agreement between climate model simulations and observed surface temperature change is a topic of scientific and policy concern. While the Earth system continues to accumulate energy due to anthropogenic and other radiative forcings, estimates of recent surface temperature evolution fall at the lower end of climate model projections. Global mean temperatures from climate model simulations are typically calculated using surface air temperatures, while the corresponding observations are based on a blend of air and sea surface temperatures. This work quantifies a systematic bias in model-observation comparisons arising from differential warming rates between sea surface temperatures and surface air temperatures over oceans. A further bias arises from the treatment of temperatures in regions where the sea ice boundary has changed. Applying the methodology of the HadCRUT4 record to climate model temperature fields accounts for 38% of the discrepancy in trend between models and observations over the period 1975–2014.

Irreducible uncertainty in near-term climate projections

Model simulations of the next few decades are widely used in assessments of climate change impacts and as guidance for adaptation. Their non-linear nature reveals a level of irreducible uncertainty which it is important to understand and quantify, especially for projections of near-term regional climate. Here we use large idealised initial condition ensembles of the FAMOUS global climate model with a 1 %/year compound increase in CO2 levels to quantify the range of future temperatures in model-based projections. These simulations explore the role of both atmospheric and oceanic initial conditions and are the largest such ensembles to date. Short-term simulated trends in global temperature are diverse, and cooling periods are more likely to be followed by larger warming rates. The spatial pattern of near-term temperature change varies considerably, but the proportion of the surface showing a warming is more consistent. In addition, ensemble spread in inter-annual temperature declines as the climate warms, especially in the North Atlantic. Over Europe, atmospheric initial condition uncertainty can, for certain ocean initial conditions, lead to 20 year trends in winter and summer in which every location can exhibit either strong cooling or rapid warming. However, the details of the distribution are highly sensitive to the ocean initial condition chosen and particularly the state of the Atlantic meridional overturning circulation. On longer timescales, the warming signal becomes more clear and consistent amongst different initial condition ensembles. An ensemble using a range of different oceanic initial conditions produces a larger spread in temperature trends than ensembles using a single ocean initial condition for all lead times. This highlights the potential benefits from initialising climate predictions from ocean states informed by observations. These results suggest that climate projections need to be performed with many more ensemble members than at present, using a range of ocean initial conditions, if the uncertainty in near-term regional climate is to be adequately quantified.

Simulating the effects of mid- to upper-tropospheric clouds on microwave emissions in EC-Earth using COSP

In this work, the Cloud Feedback Model Intercomparison (CFMIP) Observation Simulation Package (COSP) is expanded to include scattering and emission effects of clouds and precipitation at passive microwave frequencies. This represents an advancement over the official version of COSP (version 1.4.0) in which only clear-sky brightness temperatures are simulated. To highlight the potential utility of this new microwave simulator, COSP results generated using the climate model EC-Earth's version 3 atmosphere as input are compared with Microwave Humidity Sounder (MHS) channel (190.311 GHz) observations. Specifically, simulated seasonal brightness temperatures (TB) are contrasted with MHS observations for the period December 2005 to November 2006 to identify possible biases in EC-Earth's cloud and atmosphere fields. The EC-Earth's atmosphere closely reproduces the microwave signature of many of the major large-scale and regional scale features of the atmosphere and surface. Moreover, greater than 60 % of the simulated TB are within 3 K of the NOAA-18 observations. However, COSP is unable to simulate sufficiently low TB in areas of frequent deep convection. Within the Tropics, the model's atmosphere can yield an underestimation of TB by nearly 30 K for cloudy areas in the ITCZ. Possible reasons for this discrepancy include both incorrect amount of cloud ice water in the model simulations and incorrect ice particle scattering assumptions used in the COSP microwave forward model. These multiple sources of error highlight the non-unique nature of the simulated satellite measurements, a problem exacerbated by the fact that EC-Earth lacks detailed micro-physical parameters necessary for accurate forward model calculations. Such issues limit the robustness of our evaluation and suggest a general note of caution when making COSP-satellite observation evaluations.

Climate-driven expansion of blanket bogs in Britain during the Holocene

Blanket bog occupies approximately 6% of the area of the UK today. The Holocene expansion of this hyperoceanic biome has previously been explained as a consequence of Neolithic forest clearance. However, the present distribution of blanket bog in Great Britain can be predicted accurately with a simple model (PeatStash) based on summer temperature and moisture index thresholds, and the same model correctly predicts the highly disjunct distribution of blanket bog worldwide. This finding suggests that climate, rather than land-use history, controls blanket-bog distribution in the UK and everywhere else. We set out to test this hypothesis for blanket bogs in the UK using bioclimate envelope modelling compared with a database of peat initiation age estimates. We used both pollen-based reconstructions and climate model simulations of climate changes between the mid-Holocene (6000 yr BP, 6 ka) and modern climate to drive PeatStash and predict areas of blanket bog. We compiled data on the timing of blanketbog initiation, based on 228 age determinations at sites where peat directly overlies mineral soil. The model predicts that large areas of northern Britain would have had blanket bog by 6000 yr BP, and the area suitable for peat growth extended to the south after this time. A similar pattern is shown by the basal peat ages and new blanket bog appeared over a larger area during the late Holocene, the greatest expansion being in Ireland,Wales, and southwest England, as the model predicts. The expansion was driven by a summer cooling of about 2 °C, shown by both pollen-based reconstructions and climate models. The data show early Holocene (pre- Neolithic) blanket-bog initiation at over half of the sites in the core areas of Scotland and northern England. The temporal patterns and concurrence of the bioclimate model predictions and initiation data suggest that climate change provides a parsimonious explanation for the early Holocene distribution and later expansion of blanket bogs in the UK, and it is not necessary to invoke anthropogenic activity as a driver of this major landscape change.

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.

Transatlantic flight times and climate change

Aircraft do not fly through a vacuum, but through an atmosphere whose meteorological characteristics are changing because of global warming. The impacts of aviation on climate change have long been recognised, but the impacts of climate change on aviation have only recently begun to emerge. These impacts include intensified turbulence and increased take-off weight restrictions. Here we investigate the influence of climate change on flight routes and journey times. We feed synthetic atmospheric wind fields generated from climate model simulations into a routing algorithm of the type used operationally by flight planners. We focus on transatlantic flights between London and New York, and how they change when the atmospheric concentration of carbon dioxide is doubled. We find that a strengthening of the prevailing jet-stream winds causes eastbound flights to significantly shorten and westbound flights to significantly lengthen in all seasons. Eastbound and westbound crossings in winter become approximately twice as likely to take under 5 h 20 min and over 7 h 00 min, respectively. For reasons that are explained using a conceptual model, the eastbound shortening and westbound lengthening do not cancel out, causing round-trip journey times to increase. Even assuming no future growth in aviation, the extrapolation of our results to all transatlantic traffic suggests that aircraft will collectively be airborne for an extra 2000 h each year, burning an extra 7.2 million gallons of jet fuel at a cost of US$ 22 million, and emitting an extra 70 million kg of carbon dioxide, which is equivalent to the annual emissions of 7100 average British homes. Our results provide further evidence of the two-way interaction between aviation and climate change.

Seasonality of submesoscale flows in the ocean surface boundary layer

A signature of submesoscale flows in the upper ocean is skewness in the distribution of relative vorticity. Expected to result for high Rossby-number flows, such skewness has implications for mixing, dissipation and stratification within the upper ocean. An array of moorings deployed in the Northeast Atlantic for one year as part of the OSMOSIS experiment reveals that relative vorticity is positively skewed during winter even though the scale of the Rossby number is less than 0.5. Furthermore, this skewness is reduced to zero during spring and autumn. There is also evidence of modest seasonal variations in the gradient Rossby number. The proposed mechanism by which relative vorticity is skewed is that the ratio of lateral to vertical buoyancy gradients, as summarized by the inverse gradient Richardson number, restricts its range during winter but less so at other times of the year. These results support recent observations and model simulations suggesting the upper ocean is host to a seasonal cycle in submesoscale turbulence.

The biogeophysical climatic impacts of anthropogenic land use change during the Holocene

The first agricultural societies were established around 10 ka BP and had spread across much of Europe and southern Asia by 5.5 ka BP with resultant anthropogenic deforestation for crop and pasture land. Various studies (e.g. Joos et al., 2004; Kaplan et al., 2011; Mitchell et al., 2013) have attempted to assess the biogeochemical implications for Holocene climate in terms of increased carbon dioxide and methane emissions. However, less work has been done to examine the biogeophysical impacts of this early land use change. In this study, global climate model simulations with Hadley Centre Coupled Model version 3 (HadCM3) were used to examine the biogeophysical effects of Holocene land cover change on climate, both globally and regionally, from the early Holocene (8 ka BP) to the early industrial era (1850 CE). Two experiments were performed with alternative descriptions of past vegetation: (i) one in which potential natural vegetation was simulated by Top-down Representation of Interactive Foliage and Flora Including Dynamics (TRIFFID) but without land use changes and (ii) one where the anthropogenic land use model Kaplan and Krumhardt 2010 (KK10; Kaplan et al., 2009, 2011) was used to set the HadCM3 crop regions. Snapshot simulations were run at 1000-year intervals to examine when the first signature of anthropogenic climate change can be detected both regionally, in the areas of land use change, and globally. Results from our model simulations indicate that in regions of early land disturbance such as Europe and south-east Asia detectable temperature changes, outside the normal range of variability, are encountered in the model as early as 7 ka BP in the June–July–August (JJA) season and throughout the entire annual cycle by 2–3 ka BP. Areas outside the regions of land disturbance are also affected, with virtually the whole globe experiencing significant temperature changes (predominantly cooling) by the early industrial period. The global annual mean temperature anomalies found in our single model simulations were −0.22 at 1850 CE, −0.11 at 2 ka BP, and −0.03 °C at 7 ka BP. Regionally, the largest temperature changes were in Europe with anomalies of −0.83 at 1850 CE, −0.58 at 2 ka BP, and −0.24 °C at 7 ka BP. Large-scale precipitation features such as the Indian monsoon, the Intertropical Convergence Zone (ITCZ), and the North Atlantic storm track are also impacted by local land use and remote teleconnections. We investigated how advection by surface winds, mean sea level pressure (MSLP) anomalies, and tropospheric stationary wave train disturbances in the mid- to high latitudes led to remote teleconnections.

Elevation-dependent warming in mountain regions of the world

There is growing evidence that the rate of warming is amplified with elevation, such that high-mountain environments experience more rapid changes in temperature than environments at lower elevations. Elevation-dependent warming (EDW) can accelerate the rate of change in mountain ecosystems, cryospheric systems, hydrological regimes and biodiversity. Here we review important mechanisms that contribute towards EDW: snow albedo and surface-based feedbacks; water vapour changes and latent heat release; surface water vapour and radiative flux changes; surface heat loss and temperature change; and aerosols. All lead to enhanced warming with elevation (or at a critical elevation), and it is believed that combinations of these mechanisms may account for contrasting regional patterns of EDW. We discuss future needs to increase knowledge of mountain temperature trends and their controlling mechanisms through improved observations, satellite-based remote sensing and model simulations.

Land surface anomaly simulations and predictions with a climate model: an El Nino Southern Oscillation case study

The ability of climate models to reproduce and predict land surface anomalies is an important but little-studied topic. In this study, an atmosphere and ocean assimilation scheme is used to determine whether HadCM3 can reproduce and predict snow water equivalent and soil moisture during the 1997–1998 El Nino Southern Oscillation event. Soil moisture is reproduced more successfully, though both snow and soil moisture show some predictability at 1- and 4-month lead times. This result suggests that land surface anomalies may be reasonably well initialized for climate model predictions and hydrological applications using atmospheric assimilation methods over a period of time.

Coupled simulations of the Indian monsoon intraseasonal oscillation using a fine-resolution mixed-layer ocean model

The intraseasonal variability (ISV) of the Indian summer monsoon is dominated by a 30–50 day oscillation between “active” and “break” events of enhanced and reduced rainfall over the subcontinent, respectively. These organized convective events form in the equatorial Indian Ocean and propagate north to India. Atmosphere–ocean coupled processes are thought to play a key role the intensity and propagation of these events. A high-resolution, coupled atmosphere–mixed-layer-oceanmodel is assembled: HadKPP. HadKPP comprises the Hadley Centre Atmospheric Model (HadAM3) and the K Profile Parameterization (KPP) mixed-layer ocean model. Following studies that upper-ocean vertical resolution and sub-diurnal coupling frequencies improve the simulation of ISV in SSTs, KPP is run at 1 m vertical resolution near the surface; the atmosphere and ocean are coupled every three hours. HadKPP accurately simulates the 30–50 day ISV in rainfall and SSTs over India and the Bay of Bengal, respectively, but suffers from low ISV on the equator. This is due to the HadAM3 convection scheme producing limited ISV in surface fluxes. HadKPP demonstrates little of the observed northward propagation of intraseasonal events, producing instead a standing oscillation. The lack of equatorial ISV in convection in HadAM3 constrains the ability of KPP to produce equatorial SST anomalies, which further weakens the ISV of convection. It is concluded that while atmosphere–ocean interactions are undoubtedly essential to an accurate simulation of ISV, they are not a panacea for model deficiencies. In regions where the atmospheric forcing is adequate, such as the Bay of Bengal, KPP produces SST anomalies that are comparable to the Tropical Rainfall Measuring Mission Microwave Imager (TMI) SST analyses in both their magnitude and their timing with respect to rainfall anomalies over India. HadKPP also displays a much-improved phase relationship between rainfall and SSTs over a HadAM3 ensemble forced by observed SSTs, when both are compared to observations. Coupling to mixed-layer models such as KPP has the potential to improve operational predictions of ISV, particularly when the persistence time of SST anomalies is shorter than the forecast lead time.