Rheology of discrete failure regimes of anisotropic sea ice

A rheological model of sea ice is presented that incorporates the orientational distribution of ice thickness in leads embedded in isotropic floe ice. Sea ice internal stress is determined by coulombic, ridging and tensile failure at orientations where corresponding failure criteria are satisfied at minimum stresses. Because sea ice traction increases in thinner leads and cohesion is finite, such failure line angles are determined by the orientational distribution of sea ice thickness relative to the imposed stresses. In contrast to the isotropic case, sea ice thickness anisotropy results in these failure lines becoming dependent on the stress magnitude. Although generally a given failure criteria type can be satisfied at many directions, only two at most are considered. The strain rate is determined by shearing along slip lines accompanied by dilatancy and closing or opening across orientations affected by ridging or tensile failure. The rheology is illustrated by a yield curve determined by combining coulombic and ridging failure for the case of two pairs of isotropically formed leads of different thicknesses rotated with regard to each other, which models two events of coulombic failure followed by dilatancy and refreezing. The yield curve consists of linear segments describing coulombic and ridging yield as failure switches from one lead to another as the stress grows. Because sliding along slip lines is accompanied by dilatancy, at typical Arctic sea ice deformation rates a one-day-long deformation event produces enough open water that these freshly formed slip lines are preferential places of ridging failure.

Effects of plant pathogens on population dynamics and community composition in grassland ecosystems: two case studies

Grassland ecosystems comprise a major portion of the earth’s terrestrial surface, ranging from high-input cultivated monocultures or simple species mixtures to relatively unmanaged but dynamic systems. Plant pathogens are a component of these systems with their impact dependent on many interacting factors, including grassland species population dynamics and community composition, the topics covered in this paper. Plant pathogens are affected by these interactions and also act reciprocally by modifying their nature. We review these features of disease in grasslands and then introduce the 150-year long-term Park Grass Experiment (PGE) at Rothamsted Research in the UK. We then consider in detail two plant-pathogen systems present in the PGE, Tragopogon pratensis-Puccinia hysterium and Holcus lanata-Puccinia coronata. These two systems have very different life history characteristics: the first, a biennial member of the Asteraceae infected by its host-specific, systemic rust; the second, a perennial grass infected by a host-non-specific rust. We illustrate how observational, experimental and modelling studies can contribute to a better understanding of population dynamics, competitive interactions and evolutionary outcomes. With Tragopogon pratensis-Puccinia hysterium, characterised as an “outbreak” species in the PGE, we show that pathogen-induced mortality is unlikely to be involved in host population regulation; and that the presence of even a short-lived seed-bank can affect the qualitative outcomes of the host-pathogen dynamics. With Holcus lanata-Puccinia coronata, we show how nutrient conditions can affect adaptation in terms of host defence mechanisms, and that co-existence of competing species affected by a common generalist pathogen is unlikely.

RuII(α-diimine) or RuIII(α-diimine·–)? Structural, spectroscopic, and theoretical evidence for the stabilization of a prominent metal-to-Ligand charge- transfer excited-state configuration in the ground state

The new compounds [Ru(R-DAB)(acac)2] (R-DAB = 1,4-diorganyl- 1,4-diazabuta-1,3-diene; R = tert-butyl, 4-methoxyphenyl, 2,6-dimethylphenyl; acac– = 2,4-pentanedionate) exhibit intrachelate ring bond lengths 1.297

The effect of phase-correlated returns and spatial smoothing on the accuracy of radar refractivity retrievals

Radar refractivity retrievals have the potential to accurately capture near-surface humidity fields from the phase change of ground clutter returns. In practice, phase changes are very noisy and the required smoothing will diminish large radial phase change gradients, leading to severe underestimates of large refractivity changes (ΔN). To mitigate this, the mean refractivity change over the field (ΔNfield) must be subtracted prior to smoothing. However, both observations and simulations indicate that highly correlated returns (e.g., when single targets straddle neighboring gates) result in underestimates of ΔNfield when pulse-pair processing is used. This may contribute to reported differences of up to 30 N units between surface observations and retrievals. This effect can be avoided if ΔNfield is estimated using a linear least squares fit to azimuthally averaged phase changes. Nevertheless, subsequent smoothing of the phase changes will still tend to diminish the all-important spatial perturbations in retrieved refractivity relative to ΔNfield; an iterative estimation approach may be required. The uncertainty in the target location within the range gate leads to additional phase noise proportional to ΔN, pulse length, and radar frequency. The use of short pulse lengths is recommended, not only to reduce this noise but to increase both the maximum detectable refractivity change and the number of suitable targets. Retrievals of refractivity fields must allow for large ΔN relative to an earlier reference field. This should be achievable for short pulses at S band, but phase noise due to target motion may prevent this at C band, while at X band even the retrieval of ΔN over shorter periods may at times be impossible.

Quantifying errors due to frequency changes and target location uncertainty for radar refractivity retrievals

Radar refractivity retrievals can capture near-surface humidity changes, but noisy phase changes of the ground clutter returns limit the accuracy for both klystron- and magnetron-based systems. Observations with a C-band (5.6 cm) magnetron weather radar indicate that the correction for phase changes introduced by local oscillator frequency changes leads to refractivity errors no larger than 0.25 N units: equivalent to a relative humidity change of only 0.25% at 20°C. Requested stable local oscillator (STALO) frequency changes were accurate to 0.002 ppm based on laboratory measurements. More serious are the random phase change errors introduced when targets are not at the range-gate center and there are changes in the transmitter frequency (ΔfTx) or the refractivity (ΔN). Observations at C band with a 2-μs pulse show an additional 66° of phase change noise for a ΔfTx of 190 kHz (34 ppm); this allows the effect due to ΔN to be predicted. Even at S band with klystron transmitters, significant phase change noise should occur when a large ΔN develops relative to the reference period [e.g., ~55° when ΔN = 60 for the Next Generation Weather Radar (NEXRAD) radars]. At shorter wavelengths (e.g., C and X band) and with magnetron transmitters in particular, refractivity retrievals relative to an earlier reference period are even more difficult, and operational retrievals may be restricted to changes over shorter (e.g., hourly) periods of time. Target location errors can be reduced by using a shorter pulse or identified by a new technique making alternate measurements at two closely spaced frequencies, which could even be achieved with a dual–pulse repetition frequency (PRF) operation of a magnetron transmitter.

Impact of variable atmospheric and oceanic form drag on simulations of Arctic sea ice

Over Arctic sea ice, pressure ridges and floe andmelt pond edges all introduce discrete obstructions to the flow of air or water past the ice and are a source of form drag. In current climate models form drag is only accounted for by tuning the air–ice and ice–ocean drag coefficients, that is, by effectively altering the roughness length in a surface drag parameterization. The existing approach of the skin drag parameter tuning is poorly constrained by observations and fails to describe correctly the physics associated with the air–ice and ocean–ice drag. Here, the authors combine recent theoretical developments to deduce the total neutral form drag coefficients from properties of the ice cover such as ice concentration, vertical extent and area of the ridges, freeboard and floe draft, and the size of floes and melt ponds. The drag coefficients are incorporated into the Los Alamos Sea Ice Model (CICE) and show the influence of the new drag parameterization on the motion and state of the ice cover, with the most noticeable being a depletion of sea ice over the west boundary of the Arctic Ocean and over the Beaufort Sea. The new parameterization allows the drag coefficients to be coupled to the sea ice state and therefore to evolve spatially and temporally. It is found that the range of values predicted for the drag coefficients agree with the range of values measured in several regions of the Arctic. Finally, the implications of the new form drag formulation for the spinup or spindown of the Arctic Ocean are discussed.

Pan-Arctic and regional sea ice predictability: initialization month dependence

Seasonal-to-interannual predictions of Arctic sea ice may be important for Arctic communities and industries alike. Previous studies have suggested that Arctic sea ice is potentially predictable but that the skill of predictions of the September extent minimum, initialized in early summer, may be low. The authors demonstrate that a melt season “predictability barrier” and two predictability reemergence mechanisms, suggested by a previous study, are robust features of five global climate models. Analysis of idealized predictions with one of these models [Hadley Centre Global Environment Model, version 1.2 (HadGEM1.2)], initialized in January, May and July, demonstrates that this predictability barrier exists in initialized forecasts as well. As a result, the skill of sea ice extent and volume forecasts are strongly start date dependent and those that are initialized in May lose skill much faster than those initialized in January or July. Thus, in an operational setting, initializing predictions of extent and volume in July has strong advantages for the prediction of the September minimum when compared to predictions initialized in May. Furthermore, a regional analysis of sea ice predictability indicates that extent is predictable for longer in the seasonal ice zones of the North Atlantic and North Pacific than in the regions dominated by perennial ice in the central Arctic and marginal seas. In a number of the Eurasian shelf seas, which are important for Arctic shipping, only the forecasts initialized in July have continuous skill during the first summer. In contrast, predictability of ice volume persists for over 2 yr in the central Arctic but less in other regions.

The impacts of European and Asian anthropogenic sulphur dioxide emissions on Sahel rainfall

In this study, the atmospheric component of a state-of-the-art climate model (HadGEM2-ES) has been used to investigate the impacts of regional anthropogenic sulphur dioxide emissions on boreal summer Sahel rainfall. The study focuses on the transient response of the West African monsoon (WAM) to a sudden change in regional anthropogenic sulphur dioxide emissions, including land surface feedbacks, but without sea surface temperature (SST) feedbacks. The response occurs in two distinct phases: 1) fast adjustment of the atmosphere on a time scale of days to weeks (up to 3 weeks) through aerosol-radiation and aerosol-cloud interactions with weak hydrological cycle changes and surface feedbacks. 2) adjustment of the atmosphere and land surface with significant local hydrological cycle changes and changes in atmospheric circulation (beyond 3 weeks). European emissions lead to an increase in shortwave (SW) scattering by increased sulphate burden, leading to a decrease in surface downward SW radiation which causes surface cooling over North Africa, a weakening of the Saharan heat low and WAM, and a decrease in Sahel precipitation. In contrast, Asian emissions lead to very little change in sulphate burden over North Africa, but they induce an adjustment of the Walker Circulation which leads again to a weakening of the WAM and a decrease in Sahel precipitation. The responses to European and Asian emissions during the second phase exhibit similar large scale patterns of anomalous atmospheric circulation and hydrological variables, suggesting a preferred response. The results support the idea that sulphate aerosol emissions contributed to the observed decline in Sahel precipitation in the second half of the twentieth century.

Cloud banding and winds in intense European cyclones: results from the DIAMET project

The DIAMET (DIAbatic influences on Mesoscale structures in ExTratropical storms) project aims to improve forecasts of high-impact weather in extratropical cyclones through field measurements, high-resolution numerical modeling, and improved design of ensemble forecasting and data assimilation systems. This article introduces DIAMET and presents some of the first results. Four field campaigns were conducted by the project, one of which, in late 2011, coincided with an exceptionally stormy period marked by an unusually strong, zonal North Atlantic jet stream and a succession of severe windstorms in northwest Europe. As a result, December 2011 had the highest monthly North Atlantic Oscillation index (2.52) of any December in the last 60 years. Detailed observations of several of these storms were gathered using the UK’s BAe146 research aircraft and extensive ground-based measurements. As an example of the results obtained during the campaign, observations are presented of cyclone Friedhelm on 8 December 2011, when surface winds with gusts exceeding 30 m s-1 crossed central Scotland, leading to widespread disruption to transportation and electricity supply. Friedhelm deepened 44 hPa in 24 hours and developed a pronounced bent-back front wrapping around the storm center. The strongest winds at 850 hPa and the surface occurred in the southern quadrant of the storm, and detailed measurements showed these to be most intense in clear air between bands of showers. High-resolution ensemble forecasts from the Met Office showed similar features, with the strongest winds aligned in linear swaths between the bands, suggesting that there is potential for improved skill in forecasts of damaging winds.

Compensation between resolved wave driving and parameterized orographic gravity wave driving of the Brewer–Dobson circulation and its response to climate change

Following recent findings, the interaction between resolved (Rossby) wave drag and parameterized orographic gravity wave drag (OGWD) is investigated, in terms of their driving of the Brewer–Dobson circulation (BDC), in a comprehensive climate model. To this end, the parameter that effectively determines the strength of OGWD in present-day and doubled CO2 simulations is varied. The authors focus on the Northern Hemisphere during winter when the largest response of the BDC to climate change is predicted to occur. It is found that increases in OGWD are to a remarkable degree compensated by a reduction in midlatitude resolved wave drag, thereby reducing the impact of changes in OGWD on the BDC. This compensation is also found for the response to climate change: changes in the OGWD contribution to the BDC response to climate change are compensated by opposite changes in the resolved wave drag contribution to the BDC response to climate change, thereby reducing the impact of changes in OGWD on the BDC response to climate change. By contrast, compensation does not occur at northern high latitudes, where resolved wave driving and the associated downwelling increase with increasing OGWD, both for the present-day climate and the response to climate change. These findings raise confidence in the credibility of climate model projections of the strengthened BDC.

The climate impact of past changes in halocarbons and CO2in the tropical UTLS region

A chemistry-climate model coupled to an ocean model is used to compare the climate impact of past (1960-2010) changes in concentrations of halocarbons with those of CO2 in the tropical upper troposphere and lower stratosphere. The halocarbon contribution to both upper troposphere warming and the associated increase in lower stratospheric upwelling is about 40% as large as that due to CO2. Trends in cold-point temperature and lower stratosphere water vapor are positive for both halocarbons and CO2, and are of about the same magnitude. Trends in lower stratosphere ozone are negative, due to the increased upwelling. These increases in water vapor and decreases in lower stratosphere ozone feed back on lower stratosphere temperature through radiative cooling. The radiative cooling from ozone is about a factor of two larger than that from water vapor in the vicinity of the cold-point tropopause, while water vapor dominates at heights above 50 hPa. For halocarbons this indirect radiative cooling more than offsets the direct radiative warming, and together with the adiabatic cooling accounts for the lack of a halocarbon-induced warming of the lower stratosphere. For CO2 the indirect cooling from increased water vapor and decreased ozone is of comparable magnitude to the direct warming from CO2 in the vicinity of the cold-point tropopause, and (together with the increased upwelling) lowers the height at which CO2 increases induce stratospheric cooling, thus explaining the relatively weak increase in cold-point temperature due to the CO2 increases.

An anatomy of the cooling of the North Atlantic Ocean in the 1960s and 1970s

In the 1960s and early 1970s sea surface temperatures in the North Atlantic Ocean cooled rapidly. There is still considerable uncertainty about the causes of this event, although various mechanisms have been proposed. In this observational study it is demonstrated that the cooling proceeded in several distinct stages. Cool anomalies initially appeared in the mid-1960s in the Nordic Seas and Gulf Stream Extension, before spreading to cover most of the Subpolar Gyre. Subsequently, cool anomalies spread into the tropical North Atlantic before retreating, in the late 1970s, back to the Subpolar Gyre. There is strong evidence that changes in atmospheric circulation, linked to a southward shift of the Atlantic ITCZ, played an important role in the event, particularly in the period 1972-76. Theories for the cooling event must account for its distinctive space-time evolution. Our analysis suggests that the most likely drivers were: 1) The “Great Salinity Anomaly” of the late 1960s; 2) An earlier warming of the subpolar North Atlantic, which may have led to a slow-down in the Atlantic Meridional Overturning Circulation; 3) An increase in anthropogenic sulphur dioxide emissions. Determining the relative importance of these factors is a key area for future work.

An investigation into the impact of using various techniques to estimate Arctic surface air temperature anomalies

Time series of global and regional mean Surface Air Temperature (SAT) anomalies are a common metric used to estimate recent climate change. Various techniques can be used to create these time series from meteorological station data. The degree of difference arising from using five different techniques, based on existing temperature anomaly dataset techniques, to estimate Arctic SAT anomalies over land and sea ice were investigated using reanalysis data as a testbed. Techniques which interpolated anomalies were found to result in smaller errors than non-interpolating techniques relative to the reanalysis reference. Kriging techniques provided the smallest errors in estimates of Arctic anomalies and Simple Kriging was often the best kriging method in this study, especially over sea ice. A linear interpolation technique had, on average, Root Mean Square Errors (RMSEs) up to 0.55 K larger than the two kriging techniques tested. Non-interpolating techniques provided the least representative anomaly estimates. Nonetheless, they serve as useful checks for confirming whether estimates from interpolating techniques are reasonable. The interaction of meteorological station coverage with estimation techniques between 1850 and 2011 was simulated using an ensemble dataset comprising repeated individual years (1979-2011). All techniques were found to have larger RMSEs for earlier station coverages. This supports calls for increased data sharing and data rescue, especially in sparsely observed regions such as the Arctic.

Simulation of the global ENSO–Tropical cyclone teleconnection by a high-resolution coupled general circulation model

This study assesses the influence of the El Niño–Southern Oscillation (ENSO) on global tropical cyclone activity using a 150-yr-long integration with a high-resolution coupled atmosphere–ocean general circulation model [High-Resolution Global Environmental Model (HiGEM); with N144 resolution: ~90 km in the atmosphere and ~40 km in the ocean]. Tropical cyclone activity is compared to an atmosphere-only simulation using the atmospheric component of HiGEM (HiGAM). Observations of tropical cyclones in the International Best Track Archive for Climate Stewardship (IBTrACS) and tropical cyclones identified in the Interim ECMWF Re-Analysis (ERA-Interim) are used to validate the models. Composite anomalies of tropical cyclone activity in El Niño and La Niña years are used. HiGEM is able to capture the shift in tropical cyclone locations to ENSO in the Pacific and Indian Oceans. However, HiGEM does not capture the expected ENSO–tropical cyclone teleconnection in the North Atlantic. HiGAM shows more skill in simulating the global ENSO–tropical cyclone teleconnection; however, variability in the Pacific is overpronounced. HiGAM is able to capture the ENSO–tropical cyclone teleconnection in the North Atlantic more accurately than HiGEM. An investigation into the large-scale environmental conditions, known to influence tropical cyclone activity, is used to further understand the response of tropical cyclone activity to ENSO in the North Atlantic and western North Pacific. The vertical wind shear response over the Caribbean is not captured in HiGEM compared to HiGAM and ERA-Interim. Biases in the mean ascent at 500 hPa in HiGEM remain in HiGAM over the western North Pacific; however, a more realistic low-level vorticity in HiGAM results in a more accurate ENSO–tropical cyclone teleconnection.

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.