Predicting microbial pollution concentrations in UK rivers in response to land use change

The Water Framework Directive has caused a paradigm shift towards the integrated management of recreational water quality through the development of drainage basin-wide programmes of measures. This has increased the need for a cost-effective diagnostic tool capable of accurately predicting riverine faecal indicator organism (FIO) concentrations. This paper outlines the application of models developed to fulfil this need, which represent the first transferrable generic FIO models to be developed for the UK to incorporate direct measures of key FIO sources (namely human and livestock population data) as predictor variables. We apply a recently developed transfer methodology, which enables the quantification of geometric mean presumptive faecal coliforms and presumptive intestinal enterococci concentrations for base- and high-flow during the summer bathing season in unmonitored UK watercourses, to predict FIO concentrations in the Humber river basin district. Because the FIO models incorporate explanatory variables which allow the effects of policy measures which influence livestock stocking rates to be assessed, we carry out empirical analysis of the differential effects of seven land use management and policy instruments (fiscal constraint, production constraint, cost intervention, area intervention, demand-side constraint, input constraint, and micro-level land use management) all of which can be used to reduce riverine FIO concentrations. This research provides insights into FIO source apportionment, explores a selection of pollution remediation strategies and the spatial differentiation of land use policies which could be implemented to deliver river quality improvements. All of the policy tools we model reduce FIO concentrations in rivers but our research suggests that the installation of streamside fencing in intensive milk producing areas may be the single most effective land management strategy to reduce riverine microbial pollution.

On the response of the Antarctic Circumpolar Current transport to climate change in coupled climate models

The CMIP3 (IPCC AR4) models show a consistent intensification and poleward shift of the westerly winds over the Southern Ocean during the 21st century. However, the responses of the Antarctic Circumpolar Currents (ACC) show great diversity in these models, with many even showing reductions in transport. To obtain some understanding of diverse responses in the ACC transport, we investigate both external atmospheric and internal oceanic processes that control the ACC transport responses in these models. While the strengthened westerlies act to increase the tilt of isopycnal surfaces and hence the ACC transport through Ekman pumping effects, the associated changes in buoyancy forcing generally tend to reduce the surface meridional density gradient. The steepening of isopycnal surfaces induced by increased wind forcing leads to enhanced (parameterized) eddy-induced transports that act to reduce the isopycnal slopes. There is also considerable narrowing of the ACC that tends to reduce the ACC transport, caused mainly by the poleward shifts of the subtropical gyres and to a lesser extent by the equatorward expansions of the subpolar gyres in some models. If the combined effect of these retarding processes is larger than that of enhanced Ekman pumping, the ACC transport will be reduced. In addition, the effect of Ekman pumping on the ACC is reduced in weakly stratified models. These findings give insight into the reliability of IPCC-class model predictions of the Southern Ocean circulation, and into the observed decadal-scale steady ACC transport.

A mechanism for the effect of tropospheric jet structure on the annular mode-like response to stratospheric forcing

For many climate forcings the dominant response of the extratropical circulation is a latitudinal shift of the tropospheric mid-latitude jets. The magnitude of this response appears to depend on climatological jet latitude in general circulation models (GCMs): lower latitude jets exhibit a larger shift. The reason for this latitude dependence is investigated for a particular forcing, heating of the equatorial stratosphere, which shifts the jet poleward. Spin-up ensembles with a simplified GCM are used to examine the evolution of the response for five different jet structures. These differ in the latitude of the eddy-driven jet, but have similar sub-tropical zonal winds. It is found that lower latitude jets exhibit a larger response due to stronger tropospheric eddy-mean flow feedbacks. A dominant feedback responsible for enhancing the poleward shift is an enhanced equatorward refraction of the eddies, resulting in an increased momentum flux, poleward of the low latitude critical line. The sensitivity of feedback strength to jet structure is associated with differences in the coherence of this behaviour across the spectrum of eddy phase speeds. In the configurations used, the higher latitude jets have a wider range of critical latitude locations. This reduces the coherence of the momentum flux anomalies associated with different phase speeds, with low phase speeds opposing the effect of high phase speeds. This suggests that, for a given sub-tropical zonal wind strength, the latitude of the eddy driven jet affects the feedback through its influence on the width of the region of westerly winds and the range of critical latitudes on the equatorward flank of the jet.

Response of the North Atlantic storm track to climate change shaped by ocean–atmosphere coupling

A poleward shift of the mid-latitude storm tracks in response to anthropogenic greenhouse-gas forcing has been diagnosed in climate model simulations1, 2. Explanations of this effect have focused on atmospheric dynamics3, 4, 5, 6, 7. However, in contrast to storm tracks in other regions, the North Atlantic storm track responds by strengthening and extending farther east, in particular on its southern flank8. These adjustments are associated with an intensification and extension of the eddy-driven jet towards western Europe9 and are expected to have considerable societal impacts related to a rise in storminess in Europe10, 11, 12. Here, we apply a regression analysis to an ensemble of coupled climate model simulations to show that the coupling between ocean and atmosphere shapes the distinct storm-track response to greenhouse-gas forcing in the North Atlantic region. In the ensemble of simulations we analyse, at least half of the differences between the storm-track responses of different models are associated with uncertainties in ocean circulation changes. We compare the fully coupled simulations with both the associated slab model simulations and an ocean-forced experiment with one climate model to establish causality. We conclude that uncertainties in the response of the North Atlantic storm track to anthropogenic emissions could be reduced through tighter constraints on the future ocean circulation.

The Aqua-Planet Experiment (APE): Response to changed meridional SST profile

This paper explores the sensitivity of Atmospheric General Circulation Model (AGCM) simulations to changes in the meridional distribution of sea surface temperature (SST). The simulations are for an aqua-planet, a water covered Earth with no land, orography or sea-ice and with specified zonally symmetric SST. Simulations from 14 AGCMs developed for Numerical Weather Prediction and climate applications are compared. Four experiments are performed to study the sensitivity to the meridional SST profile. These profiles range from one in which the SST gradient continues to the equator to one which is flat approaching the equator, all with the same maximum SST at the equator. The zonal mean circulation of all models shows strong sensitivity to latitudinal distribution of SST. The Hadley circulation weakens and shifts poleward as the SST profile flattens in the tropics. One question of interest is the formation of a double versus a single ITCZ. There is a large variation between models of the strength of the ITCZ and where in the SST experiment sequence they transition from a single to double ITCZ. The SST profiles are defined such that as the equatorial SST gradient flattens, the maximum gradient increases and moves poleward. This leads to a weakening of the mid-latitude jet accompanied by a poleward shift of the jet core. Also considered are tropical wave activity and tropical precipitation frequency distributions. The details of each vary greatly between models, both with a given SST and in the response to the change in SST. One additional experiment is included to examine the sensitivity to an off-equatorial SST maximum. The upward branch of the Hadley circulation follows the SST maximum off the equator. The models that form a single precipitation maximum when the maximum SST is on the equator shift the precipitation maximum off equator and keep it centered over the SST maximum. Those that form a double with minimum on the equatorial maximum SST shift the double structure off the equator, keeping the minimum over the maximum SST. In both situations only modest changes appear in the shifted profile of zonal average precipitation. When the upward branch of the Hadley circulation moves into the hemisphere with SST maximum, the zonal average zonal, meridional and vertical winds all indicate that the Hadley cell in the other hemisphere dominates.

Assessment of surface winds over the Atlantic, Indian and Pacific Ocean sectors of the Southern Ocean in CMIP5 models: historical bias, forcing response, and state dependence

An assessment of the fifth Coupled Models Intercomparison Project (CMIP5) models’ simulation of the near-surface westerly wind jet position and strength over the Atlantic, Indian and Pacific sectors of the Southern Ocean is presented. Compared with reanalysis climatologies there is an equatorward bias of 3.7° (inter-model standard deviation of ± 2.2°) in the ensemble mean position of the zonal mean jet. The ensemble mean strength is biased slightly too weak, with the largest biases over the Pacific sector (-1.6±1.1 m/s, 27 -22%). An analysis of atmosphere-only (AMIP) experiments indicates that 41% of the zonal mean position bias comes from coupling of the ocean/ice models to the atmosphere. The response to future emissions scenarios (RCP4.5 and RCP8.5) is characterized by two phases: (i) the period of most rapid ozone recovery (2000-2049) during which there is insignificant change in summer; and (ii) the period 2050-2098 during which RCP4.5 simulations show no significant change but RCP8.5 simulations show poleward shifts (0.30, 0.19 and 0.28°/decade over the Atlantic, Indian and Pacific sectors respectively), and increases in strength (0.06, 0.08 and 0.15 m/s/decade respectively). The models with larger equatorward position biases generally show larger poleward shifts (i.e. state dependence). This inter-model relationship is strongest over the Pacific sector (r=-0.89) and insignificant over the Atlantic sector (r=-0.50). However, an assessment of jet structure shows that over the Atlantic sector jet shift is significantly correlated with jet width whereas over the Pacific sector the distance between the sub-polar and sub-tropical westerly jets appears to be more important.

Dynamics of the lower stratospheric circulation response to ENSO

A robust feature of the observed response to El Nin˜o–Southern Oscillation (ENSO) is an altered circulation in the lower stratosphere. When sea surface temperatures (SSTs) in the tropical Pacific are warmer there is enhanced upwelling and cooling in the tropical lower stratosphere and downwelling and warming in the midlatitudes, while the opposite is true of cooler SSTs. The midlatitude lower stratospheric response to ENSO is larger in the Southern Hemisphere (SH) than in the Northern Hemisphere (NH). In this study the dynamical version of the Canadian Middle Atmosphere Model (CMAM) is used to simulate 25 realizations of the atmospheric response to the 1982/83 El Nin˜o and the 1973/74 La Nin˜ a. This version ofCMAMis a comprehensive high-top general circulation model that does not include interactive chemistry. The observed lower stratospheric response to ENSO is well reproduced by the simulations, allowing them to be used to investigate the mechanisms involved. Both the observed and simulated responses maximize in December–March and so this study focuses on understanding the mechanisms involved in that season. The response in tropical upwelling is predominantly driven by anomalous transient synoptic-scale wave drag in the SH subtropical lower stratosphere, which is also responsible for the compensating SH midlatitude response. This altered wave drag stems from an altered upward flux of wave activity from the troposphere into the lower stratosphere between 208 and 408S. The altered flux of wave activity can be divided into two distinct components. In the Pacific, the acceleration of the zonal wind in the subtropics from the warmer tropical SSTs results in a region between the midlatitude and subtropical jets where there is an enhanced source of low phase speed eddies. At other longitudes, an equatorward shift of the midlatitude jet from the extratropical tropospheric response to El Nin˜o results in an enhanced source of waves of higher phase speeds in the subtropics. The altered resolved wave drag is only apparent in the SH and the difference between the two hemispheres can be related to the difference in the climatological jet structures in this season and the projection of the wind anomalies associated with ENSO onto those structures.

Response of tropical cyclones to idealized climate change experiments in a global high resolution coupled general circulation model

We present an assessment of how tropical cyclone activity might change due to the influence of increased atmospheric carbon dioxide concentrations, using the UK’s High Resolution Global Environment Model (HiGEM) with N144 resolution (~90 km in the atmosphere and ~40 km in the ocean). Tropical cyclones are identified using a feature tracking algorithm applied to model output. Tropical cyclones from idealized 30-year 2×CO2 (2CO2) and 4×CO2 (4CO2) simulations are compared to those identified in a 150-year present-day simulation, which is separated into a 5-member ensemble of 30-year integrations. Tropical cyclones are shown to decrease in frequency globally by 9% in the 2CO2 and 26% in the 4CO2. Tropical cyclones only become more intese in the 4CO2, however uncoupled time slice experiments reveal an increase in intensity in the 2CO2. An investigation into the large-scale environmental conditions, known to influence tropical cyclone activity in the main development regions, is used to determine the response of tropical cyclone activity to increased atmospheric CO2. A weaker Walker circulation and a reduction in zonally averaged regions of updrafts lead to a shift in the location of tropical cyclones in the northern hemisphere. A decrease in mean ascent at 500 hPa contributes to the reduction of tropical cyclones in the 2CO2 in most basins. The larger reduction of tropical cyclones in the 4CO2 arises from further reduction of mean ascent at 500 hPa and a large enhancement of vertical wind shear, especially in the southern hemisphere, North Atlantic and North East Pacific.

Mid-Holocene climates of the Americas: a dynamical response to changed seasonality

Simulations of the climatic response to mid-Holocene (6 ka BP) orbital forcing with two coupled ocean–atmosphere models (FOAM and CSM) show enhancement of monsoonal precipitation in parts of the American Southwest, Central America and northernmost South America during Northern Hemisphere summer. The enhanced onshore flow that brings precipitation into Central America is caused by a northward displacement of the inter-tropical convergence zone, driven by cooling of the equatorial and warming of the northern subtropical and mid-latitude ocean. Ocean feedbacks also enhance precipitation over the American Southwest, although the increase in monsoon precipitation there is largely driven by increases in land-surface temperature. The northward shift in the equatorial precipitation band that causes enhanced precipitation in Central America and the American Southwest has a negative feedback effect on monsoonal precipitation in northern South America. The simulations demonstrate that mid-Holocene aridity in the mid-continent of North America is dynamically linked to the orbitally induced enhancement of the summer monsoon in the American Southwest, with a spatial structure (wet in the Southwest and dry in the mid-continent) similar to that found in strong monsoon years today. Changes in winter precipitation along the west coast of North America, in Central America and along the Gulf Coast, caused by southward-displacement of the westerly storm tracks, indicate that changes in the Northern Hemisphere winter monsoon also play a role in regional climate changes during the mid-Holocene. Although the simulations with FOAM and CSM differ in detail, the general mechanisms and patterns are common to both. The model results thus provide a coherent dynamical explanation for regional patterns of increased or decreased aridity shown by vegetation, lake status and aeolian data from the Americas

Intercomparison of simulated global vegetation distribution in response to 6 kyr BP orbital forcing

The response of ten atmospheric general circulation models to orbital forcing at 6 kyr BP has been investigated using the BIOME model, which predicts equilibrium vegetation distribution, as a diagnostic. Several common features emerge: (a) reduced tropical rain forest as a consequence of increased aridity in the equatorial zone, (b) expansion of moisture-demanding vegetation in the Old World subtropics as a consequence of the expansion of the Afro–Asian monsoon, (c) an increase in warm grass/shrub in the Northern Hemisphere continental interiors in response to warming and enhanced aridity, and (d) a northward shift in the tundra–forest boundary in response to a warmer growing season at high northern latitudes. These broadscale features are consistent from model to model, but there are differences in their expression at a regional scale. Vegetation changes associated with monsoon enhancement and high-latitude summer warming are consistent with palaeoenvironmental observations, but the simulated shifts in vegetation belts are too small in both cases. Vegetation changes due to warmer and more arid conditions in the midcontinents of the Northern Hemisphere are consistent with palaeoenvironmental data from North America, but data from Eurasia suggests conditions were wetter at 6 kyr BP than today. The models show quantitatively similar vegetation changes in the intertropical zone, and in the northern and southern extratropics. The small differences among models in the magnitude of the global vegetation response are not related to differences in global or zonal climate averages, but reflect differences in simulated regional features. Regional-scale analyses will therefore be necessary to identify the underlying causes of such differences among models.

The annual-cycle modulation of meridional asymmetry in ENSO’s atmospheric response and its dependence on ENSO zonal structure

Previous studies documented that a distinct southward shift of central-Pacific low-level wind anomalies occurring during the ENSO decaying phase, is caused by an interaction between the Western Pacific annual cycle and El Niño-Southern Oscillation (ENSO) variability. The present study finds that the meridional movement of the central-Pacific wind anomalies appears only during traditional Eastern-Pacific (or EP) El Niño events rather than in Central-Pacific (CP) El Niño events in which sea surface temperature (SST) anomalies are confined to the central Pacific. The zonal structure of ENSO-related SST anomalies therefore has an important effect on meridional asymmetry in the associated atmospheric response and its modulation by the annual cycle. In contrast to EP El Niño events, the SST anomalies of CP El Niño events extend further west towards to the warm pool region with its climatological warm SSTs. In the warm pool region, relatively small SST anomalies thus are able to excite convection anomalies on both sides of the equator, even with a meridionally asymmetric SST background state. Therefore, almost meridionally symmetric precipitation and wind anomalies are observed over the central Pacific during the decaying phase of CP El Niño events. The SST anomaly pattern of La Niña events is similar to CP El Niño events with a reversed sign. Accordingly, no distinct southward displacement of the atmospheric response occurs over the central Pacific during the La Niña decaying phase. These results have important implications for ENSO climate impacts over East Asia, since the anomalous low-level anticyclone over the western North Pacific is an integral part of the annual cycle-modulated ENSO response.

Resonance response of electrorheological fluids in vertical oscillation squeeze flow

The resonance effect of microcrystalline cellulose/castor oil electrorheological (ER) suspensions was studied in a compressed oscillatory squeeze flow under external electric fields. The resonance frequency first increases linearly with increasing external held, and then shift to high-field plateau. The amplitudes of resonance peak increase sharply with the applied fields in the range of 0.17-1.67kV/mm. The phase difference of the reduced displacement relative to the excitation force inverses in the case of resonance. A viscoelasticity model of the ER suspensions, which offers both the equivalent stiffness and the viscous damping, should be responsible for the appearance of resonance. The influence of the electric field on the resonance frequency and the resonance hump is consistent qualitatively with the interpretation of our proposed model. Storage modulus G' was presented for the purpose of investigating this influence.

Antarctic ocean and sea ice response to ozone depletion: a two timescale problem

The response of the Southern Ocean to a repeating seasonal cycle of ozone loss is studied in two coupled climate models and found to comprise both fast and slow processes. The fast response is similar to the inter-annual signature of the Southern Annular Mode (SAM) on Sea Surface Temperature (SST), on to which the ozone-hole forcing projects in the summer. It comprises enhanced northward Ekman drift inducing negative summertime SST anomalies around Antarctica, earlier sea ice freeze-up the following winter, and northward expansion of the sea ice edge year-round. The enhanced northward Ekman drift, however, results in upwelling of warm waters from below the mixed layer in the region of seasonal sea ice. With sustained bursts of westerly winds induced by ozone-hole depletion, this warming from below eventually dominates over the cooling from anomalous Ekman drift. The resulting slow-timescale response (years to decades) leads to warming of SSTs around Antarctica and ultimately a reduction in sea-ice cover year-round. This two-timescale behavior - rapid cooling followed by slow but persistent warming - is found in the two coupled models analysed, one with an idealized geometry, the other a complex global climate model with realistic geometry. Processes that control the timescale of the transition from cooling to warming, and their uncertainties are described. Finally we discuss the implications of our results for rationalizing previous studies of the effect of the ozone-hole on SST and sea-ice extent. %Interannual variability in the Southern Annular Mode (SAM) and sea ice covary such that an increase and southward shift in the surface westerlies (a positive phase of the SAM) coincides with a cooling of Sea Surface Temperature (SST) around 70-50$^\circ$S and an expansion of the sea ice cover, as seen in observations and models alike. Yet, in modeling studies, the Southern Ocean warms and sea ice extent decreases in response to sustained, multi-decadal positive SAM-like wind anomalies driven by 20th century ozone depletion. Why does the Southern Ocean appear to have disparate responses to SAM-like variability on interannual and multidecadal timescales? Here it is demonstrated that the response of the Southern Ocean to ozone depletion has a fast and a slow response. The fast response is similar to the interannual variability signature of the SAM. It is dominated by an enhanced northward Ekman drift, which transports heat northward and causes negative SST anomalies in summertime, earlier sea ice freeze-up the following winter, and northward expansion of the sea ice edge year round. The enhanced northward Ekman drift causes a region of Ekman divergence around 70-50$^\circ$S, which results in upwelling of warmer waters from below the mixed layer. With sustained westerly wind enhancement in that latitudinal band, the warming due to the anomalous upwelling of warm waters eventually dominates over the cooling from the anomalous Ekman drift. Hence, the slow response ultimately results in a positive SST anomaly and a reduction in the sea ice cover year round. We demonstrate this behavior in two models: one with an idealized geometry and another, more detailed, global climate model. However, the models disagree on the timescale of transition from the fast (cooling) to the slow (warming) response. Processes that controls this transition and their uncertainties are discussed.

Response of the North Atlantic wave climate to atmospheric modes of variability

This study investigates the relationship between the wind wave climate and the main climate modes of atmospheric variability in the North Atlantic Ocean. The modes considered are the North Atlantic Oscillation (NAO), the East Atlantic (EA) pattern, the East Atlantic Western Russian (EA/WR) pattern and the Scandinavian (SCAN) pattern. The wave dataset consists of buoys records, remote sensing altimetry observations and a numerical hindcast providing significant wave height (SWH), mean wave period (MWP) and mean wave direction (MWD) for the period 1989–2009. After evaluating the reliability of the hindcast, we focus on the impact of each mode on seasonal wave parameters and on the relative importance of wind-sea and swell components. Results demonstrate that the NAO and EA patterns are the most relevant, whereas EA/WR and SCAN patterns have a weaker impact on the North Atlantic wave climate variability. During their positive phases, both NAO and EA patterns are related to winter SWH at a rate that reaches 1 m per unit index along the Scottish coast (NAO) and Iberian coast (EA) patterns. In terms of winter MWD, the two modes induce a counterclockwise shift of up to 65° per negative NAO (positive EA) unit over west European coasts. They also increase the winter MWP in the North Sea and in the Bay of Biscay (up to 1 s per unit NAO) and along the western coasts of Europe and North Africa (1 s per unit EA). The impact of winter EA pattern on all wave parameters is mostly caused through the swell wave component.

The response of the Southern Hemispheric eddy-driven jet to future changes in shortwave radiation in CMIP5

A strong relationship is found between changes in the meridional gradient of absorbed shortwave radiation (ASR) and Southern Hemispheric jet shifts in 21st century climate simulations of CMIP5 (Coupled Model Intercomparison Project phase 5) coupled models. The relationship is such that models with increases in the meridional ASR gradient around the southern midlatitudes, and therefore increases in midlatitude baroclinicity, tend to produce a larger poleward jet shift. The ASR changes are shown to be dominated by changes in cloud properties, with sea ice declines playing a secondary role. We demonstrate that the ASR changes are the cause, and not the result, of the intermodel differences in jet response by comparing coupled simulations with experiments in which sea surface temperature increases are prescribed. Our results highlight the importance of reducing the uncertainty in cloud feedbacks in order to constrain future circulation changes.