On the climate response of the low-latitude Pacific Ocean to changes in the global freshwater cycle

Under global warming, the predicted intensification of the global freshwater cycle will modify the net freshwater flux at the ocean surface. Since the freshwater flux maintains ocean salinity structures, changes to the density-driven ocean circulation are likely. A modified ocean circulation could further alter the climate, potentially allowing rapid changes, as seen in the past. The relevant feedback mechanisms and timescales are poorly understood in detail, however, especially at low latitudes where the effects of salinity are relatively subtle. In an attempt to resolve some of these outstanding issues, we present an investigation of the climate response of the low-latitude Pacific region to changes in freshwater forcing. Initiated from the present-day thermohaline structure, a control run of a coupled ocean-atmosphere general circulation model is compared with a perturbation run in which the net freshwater flux is prescribed to be zero over the ocean. Such an extreme experiment helps to elucidate the general adjustment mechanisms and their timescales. The atmospheric greenhouse gas concentrations are held constant, and we restrict our attention to the adjustment of the upper 1,000 m of the Pacific Ocean between 40°N and 40°S, over 100 years. In the perturbation run, changes to the surface buoyancy, near-surface vertical mixing and mixed-layer depth are established within 1 year. Subsequently, relative to the control run, the surface of the low-latitude Pacific Ocean in the perturbation run warms by an average of 0.6°C, and the interior cools by up to 1.1°C, after a few decades. This vertical re-arrangement of the ocean heat content is shown to be achieved by a gradual shutdown of the heat flux due to isopycnal (i.e. along surfaces of constant density) mixing, the vertical component of which is downwards at low latitudes. This heat transfer depends crucially upon the existence of density-compensating temperature and salinity gradients on isopycnal surfaces. The timescale of the thermal changes in the perturbation run is therefore set by the timescale for the decay of isopycnal salinity gradients in response to the eliminated freshwater forcing, which we demonstrate to be around 10-20 years. Such isopycnal heat flux changes may play a role in the response of the low-latitude climate to a future accelerated freshwater cycle. Specifically, the mechanism appears to represent a weak negative sea surface temperature feedback, which we speculate might partially shield from view the anthropogenically-forced global warming signal at low latitudes. Furthermore, since the surface freshwater flux is shown to play a role in determining the ocean's thermal structure, it follows that evaporation and/or precipitation biases in general circulation models are likely to cause sea surface temperature biases.

The omega equation and potential vorticity

Two fundamental perspectives on the dynamics of midlatitude weather systems are provided by potential vorticity (PV) and the omega equation. The aim of this paper is to investigate the link between the two perspectives, which has so far received very little attention in the meteorological literature. It also aims to give a quantitative basis for discussion of quasi-geostrophic vertical motion in terms of components associated with system movement, maintaining a constant thermal structure, and with the development of that structure. The former links with the isentropic relative-flow analysis technique. Viewed in a moving frame of reference, the measured development of a system depends on the velocity of that frame of reference. The requirement that the development should be a minimum provides a quantitative method for determining the optimum system velocity. The component of vertical velocity associated with development is shown to satisfy an omega equation with forcing determined from the relative advection of interior PV and boundary temperature. The analysis carries through in the presence of diabatic heating provided the omega equation forcing is based on the interior PV and boundary thermal tendencies, including the heating effect. The analysis is shown to be possible also at the level of the semi-geostrophic approximation. The analysis technique is applied to a number of idealized problems that can be considered to be building blocks for midlatitude synoptic-scale dynamics. They focus on the influences of interior PV, boundary temperature, an interior boundary, baroclinic instability associated with two boundaries, and also diabatic heating. In each case, insights yielded by the new perspective are sought into the dynamical behaviour, especially that related to vertical motion. Copyright © 2003 Royal Meteorological Society

A model study of oceanic mechanisms affecting Equatorial Pacific sea surface temperature during the 1997-98 El Niño

In this study, the processes affecting sea surface temperature variability over the 1992–98 period, encompassing the very strong 1997–98 El Niño event, are analyzed. A tropical Pacific Ocean general circulation model, forced by a combination of weekly ERS1–2 and TAO wind stresses, and climatological heat and freshwater fluxes, is first validated against observations. The model reproduces the main features of the tropical Pacific mean state, despite a weaker than observed thermal stratification, a 0.1 m s−1 too strong (weak) South Equatorial Current (North Equatorial Countercurrent), and a slight underestimate of the Equatorial Undercurrent. Good agreement is found between the model dynamic height and TOPEX/Poseidon sea level variability, with correlation/rms differences of 0.80/4.7 cm on average in the 10°N–10°S band. The model sea surface temperature variability is a bit weak, but reproduces the main features of interannual variability during the 1992–98 period. The model compares well with the TAO current variability at the equator, with correlation/rms differences of 0.81/0.23 m s−1 for surface currents. The model therefore reproduces well the observed interannual variability, with wind stress as the only interannually varying forcing. This good agreement with observations provides confidence in the comprehensive three-dimensional circulation and thermal structure of the model. A close examination of mixed layer heat balance is thus undertaken, contrasting the mean seasonal cycle of the 1993–96 period and the 1997–98 El Niño. In the eastern Pacific, cooling by exchanges with the subsurface (vertical advection, mixing, and entrainment), the atmospheric forcing, and the eddies (mainly the tropical instability waves) are the three main contributors to the heat budget. In the central–western Pacific, the zonal advection by low-frequency currents becomes the main contributor. Westerly wind bursts (in December 1996 and March and June 1997) were found to play a decisive role in the onset of the 1997–98 El Niño. They contributed to the early warming in the eastern Pacific because the downwelling Kelvin waves that they excited diminished subsurface cooling there. But it is mainly through eastward advection of the warm pool that they generated temperature anomalies in the central Pacific. The end of El Niño can be linked to the large-scale easterly anomalies that developed in the western Pacific and spread eastward, from the end of 1997 onward. In the far-western Pacific, because of the shallower than normal thermocline, these easterlies cooled the SST by vertical processes. In the central Pacific, easterlies pushed the warm pool back to the west. In the east, they led to a shallower thermocline, which ultimately allowed subsurface cooling to resume and to quickly cool the surface layer.

How well do coupled models replicate ocean energetics relevant to ENSO?

Accurate replication of the processes associated with the energetics of the tropical ocean is necessary if coupled GCMs are to simulate the physics of ENSO correctly, including the transfer of energy from the winds to the ocean thermocline and energy dissipation during the ENSO cycle. Here, we analyze ocean energetics in coupled GCMs in terms of two integral parameters describing net energy loss in the system using the approach recently proposed by Brown and Fedorov (J Clim 23:1563–1580, 2010a) and Fedorov (J Clim 20:1108–1117, 2007). These parameters are (1) the efficiency c of the conversion of wind power into the buoyancy power that controls the rate of change of the available potential energy (APE) in the ocean and (2) the e-folding rate a that characterizes the damping of APE by turbulent diffusion and other processes. Estimating these two parameters for coupled models reveals potential deficiencies (and large differences) in how state-of-the-art coupled GCMs reproduce the ocean energetics as compared to ocean-only models and data assimilating models. The majority of the coupled models we analyzed show a lower efficiency (values of c in the range of 10–50% versus 50–60% for ocean-only simulations or reanalysis) and a relatively strong energy damping (values of a-1 in the range 0.4–1 years versus 0.9–1.2 years). These differences in the model energetics appear to reflect differences in the simulated thermal structure of the tropical ocean, the structure of ocean equatorial currents, and deficiencies in the way coupled models simulate ENSO.

The GODAE high resolution sea surface temperature pilot project

Sea surface temperature (SST) measurements are required by operational ocean and atmospheric forecasting systems to constrain modeled upper ocean circulation and thermal structure. The Global Ocean Data Assimilation Experiment (GODAE) High Resolution SST Pilot Project (GHRSST-PP) was initiated to address these needs by coordinating the provision of accurate, high-resolution, SST products for the global domain. The pilot project is now complete, but activities continue within the Group for High Resolution SST (GHRSST). The pilot project focused on harmonizing diverse satellite and in situ data streams that were indexed, processed, quality controlled, analyzed, and documented within a Regional/Global Task Sharing (R/GTS) framework implemented in an internationally distributed manner. Data with meaningful error estimates developed within GHRSST are provided by services within R/GTS. Currently, several terabytes of data are processed at international centers daily, creating more than 25 gigabytes of product. Ensemble SST analyses together with anomaly SST outputs are generated each day, providing confidence in SST analyses via diagnostic outputs. Diagnostic data sets are generated and Web interfaces are provided to monitor the quality of observation and analysis products. GHRSST research and development projects continue to tackle problems of instrument calibration, algorithm development, diurnal variability, skin temperature deviation, and validation/verification of GHRSST products. GHRSST also works closely with applications and users, providing a forum for discussion and feedback between SST users and producers on a regular basis. All data within the GHRSST R/GTS framework are freely available. This paper reviews the progress of GHRSST-PP, highlighting achievements that have been fundamental to the success of the pilot project.

Climatology of Upper Tropospheric–Lower Stratospheric (UTLS) Jets and Tropopauses in MERRA

A global climatology (1979–2012) from the Modern-Era Retrospective Analysis for Research and Applications (MERRA) shows distributions and seasonal evolution of upper tropospheric jets and their relationships to the stratospheric subvortex and multiple tropopauses. The overall climatological patterns of upper tropospheric jets confirm those seen in previous studies, indicating accurate representation of jet stream dynamics in MERRA. The analysis shows a Northern Hemisphere (NH) upper tropospheric jet stretching nearly zonally from the mid-Atlantic across Africa and Asia. In winter–spring, this jet splits over the eastern Pacific, merges again over eastern North America, and then shifts poleward over the North Atlantic. The jets associated with tropical circulations are also captured, with upper tropospheric westerlies demarking cyclonic flow downstream from the Australian and Asian monsoon anticyclones and associated easterly jets. Multiple tropopauses associated with the thermal tropopause “break” commonly extend poleward from the subtropical upper tropospheric jet. In Southern Hemisphere (SH) summer, the tropopause break, along with a poleward-stretching secondary tropopause, often occurs across the tropical westerly jet downstream of the Australian monsoon region. SH high-latitude multiple tropopauses, nearly ubiquitous in June–July, are associated with the unique polar winter thermal structure. High-latitude multiple tropopauses in NH fall–winter are, however, sometimes associated with poleward-shifted upper tropospheric jets. The SH subvortex jet extends down near the level of the subtropical jet core in winter and spring. Most SH subvortex jets merge with an upper tropospheric jet between May and December; although much less persistent than in the SH, merged NH subvortex jets are common between November and April.

Further considerations of cosmic ray modulation of infra-red radiation in the atmosphere

Understanding effects of ionisation in the lower atmosphere is a new interdisciplinary area, crossing the traditionally distinct scientific boundaries between astro-particle and atmospheric physics and also requiring understanding of both heliospheric and magnetospheric influences on cosmic rays. Following the paper of Erlykin et al. (2014) we develop further the interpretation of our observed changes in long-wave (LW) radiation, Aplin and Lockwood (2013) by taking account of both cosmic ray ionisation yields and atmospheric radiative transfer. To demonstrate this, we show that the thermal structure of the whole atmosphere needs to be considered along with the vertical profile of ionisation. Allowing for, in particular, ionisation by all components of a cosmic ray shower and not just by the muons, reveals that the effect we have detected is certainly not inconsistent with laboratory observations of the LW absorption cross section. The analysis presented here, although very different from that of Erlykin et al., does come to the same conclusion that the events detected by AL were not caused by individual cosmic ray primaries – not because it is impossible on energetic grounds, but because events of the required energy are too infrequent for the 12 h_1 rate at which they were seen by the AL experiment. The present paper numerically models the effect of three different scenario changes to the primary GCR spectrum which all reproduce the required magnitude of the effect observed by AL. However, they cannot solely explain the observed delay in the peak effect which, if confirmed, would appear to open up a whole new and interesting area in the study of water oligomers and their effects on LW radiation. We argue that a technical artefact in the AL experiment is highly unlikely and that our initial observations merit both a wide-ranging follow-up experiment and more rigorous, self-consistent, three-dimensional radiative transfer modelling

The solsticial pause on Mars: 2 modelling and investigation of causes

The martian solsticial pause, presented in a companion paper (Lewis et al., this issue), was investigated further through a series of model runs using the UK version of the LMD/UK Mars Global Climate Model. It was found that the pause could not be adequately reproduced if radiatively active water ice clouds were omitted from the model. When clouds were used, along with a realistic time-dependent dust opacity distribution, a substantial minimum in near-surface transient eddy activity formed around solstice in both hemispheres. The net effect of the clouds in the model is, by altering the thermal structure of the atmosphere, to decrease the vertical shear of the westerly jet near the surface around solstice, and thus reduce baroclinic growth rates. A similar effect was seen under conditions of large dust loading, implying that northern midlatitude eddy activity will tend to become suppressed after a period of intense flushing storm formation around the northern cap edge. Suppression of baroclinic eddy generation by the barotropic component of the flow and via diabatic eddy dissipation were also investigated as possible mechanisms leading to the formation of the solsticial pause but were found not to make major contributions. Zonal variations in topography were found to be important, as their presence results in weakened transient eddies around winter solstice in both hemispheres, through modification of the near-surface flow. The zonal topographic asymmetry appears to be the primary reason for the weakness of eddy activity in the southern hemisphere relative to the northern hemisphere, and the ultimate cause of the solsticial pause in both hemispheres. The meridional topographic gradient was found to exert a much weaker influence on near-surface transient eddies.

Structure, and thermal, magnetic and chemical properties of copper(II)iodoacetate monohydrate

Reaction of iodoacetic acid with cupric carbonate in water in dimmed light yields green Cu(ICH2COO)(2 center dot)H2O (1). From X-ray crystallography, it is found to be a tetra-acetato bridged copper(II) dimer with the water molecules occupying the apical positions. In thermogravimetry, the coordinated water molecules are lost in the temperature range 50-100 degrees C. From magnetic susceptibility measurements in the temperature range 300-1.8 K, the exchange coupling constant J is found to be -142(1) cm(-1) and g = 2.18(2) with the spin Hamiltonian H = -2J{S-Cu1 center dot S-Cu2}. It reacts with 2,2'-bipyridine (bpy) to yield [Cu(bpy)(2)I]I. It oxidises thiophenol to Ph-S-S-Ph under dry N-2 atmosphere.

Synthesis, crystal structure, magnetic behavior and thermal property of three polynuclear complexes: [M(dca)(2)(H2O)2](n)center dot(hmt)(n) [M = Mn(II), Co(II)] and [Co(dca)(2)(bpds)](n) [dca, dicyanamide; hint, hexamethylenetetramine; bpds, 4,4 '-bipyridyl disulfide]

Three mu(1.5)-dicyanamide bridged Mn(II) and Co(II) complexes having molecular formula [Mn(dca)(2)(H2O)(2)](n)center dot(hmt)(n) (1), [Co(dca)(2) (H2O)(2)](n)center dot(hmt)(n) (2) and [Co(dca)(2)(bpds)](n) (3) [dca = dicyanamide; hmt = hexamethylenetetramine; bpds = 4,4'-bipyridyl disulfide] have been synthesized and characterized by single crystal X-ray diffraction study, low temperature (300-2 K) magnetic measurement and thermal behavior. The X-ray diffraction analysis of 1 and 2 reveals that they are isostructural, comprising of 1D coordination polymers [M(dca)(2)(H2O)(2)](n) [M = Mn(II), Co(II) for 1 and 2. respectively] with uncoordinated hmt molecules located among the chains. The [M(dca)(2)(H2O)(2)](n) chains and the lattice hint molecules are connected through H-bonds resulting in a 3D supramolecular architecture. The octahedral N4O2 chromophore surrounding the metal ion forms via two trans located water oxygens and four nitrogens from four nitrile dca. Complex 3 is a 1D chain formed by two mu(1.5)-dca and one bridging bpds. The octahedral N-6 coordination sphere surrounding the cobalt ions comprises four nitrogens from dca and two from bpds. Low temperature magnetic study indicates small antiferromagnetic coupling for all the complexes. Best fit parameters for 1: J = -0.17 cm(-1), g = -2.03 with R = 6.1 x 10(-4), for 2, J = -0.50 cm(-1), and for 3, J = -0.95 cm(-1). (c) 2005 Elsevier B.V. All rights reserved.

Synthesis, crystal structure and thermal studies of catena(L-glutamato)-aqua copper(II) monohydrate

The synthesis. crystal structure and thermal study of the blue catena-(L-glutamato)-aqua copper(II) monohydrate have been reported. The compound crystallizes in P2(1)2(1)2(1) space group and consists of a polymeric three-dimensional network of copper(II) which is coordinated with the amino nitrogen and the carboxylate oxygen Of L-glutamate, the side chain carboxylate oxygen of a neighbouring L-glutamate and the oxygen of a water molecule in the equatorial position. Weak coordination of two additional glutamate oxygen atoms to both the axial positions Completes a distorted octahedron. The crystal structure shows that the lattice water is stabilized by the formation of strong H-bonding network with the coordinated water molecule. Removal and reabsorption of the water molecule have been studied by thermal analysis.

A differential thermal expansion approach to crystal structure determination from powder diffraction data

Differential thermal expansion over the range 90-210 K has been applied successfully to determine the crystal structure of chlorothiazide from synchrotron powder diffraction data using direct methods. Key to the success of the approach is the use of a multi-data-set Pawley refinement to extract a set of reflection intensities that is more 'single-crystal-like' than those extracted from a single data set. The improvement in reflection intensity estimates is quantified by comparison with reference single-crystal intensities. (C) 2008 International Union of Crystallography Printed in Singapore - all rights reserved

Local and average structure in zinc cyanide: towards an understanding of the atomistic origin of negative thermal expansion

Neutron diffraction at 11.4 and 295 K and solid-state 67Zn NMR are used to determine both the local and average structures in the disordered, negative thermal expansion (NTE) material, Zn(CN)2. Solid-state NMR not only confirms that there is head-to-tail disorder of the C≡N groups present in the solid, but yields information about the relative abundances of the different Zn(CN)4-n(NC)n tetrahedral species, which do not follow a simple binomial distribution. The Zn(CN)4 and Zn(NC)4 species occur with much lower probabilities than are predicted by binomial theory, supporting the conclusion that they are of higher energy than the other local arrangements. The lowest energy arrangement is Zn(CN)2(NC)2. The use of total neutron diffraction at 11.4 K, with analysis of both the Bragg diffraction and the derived total correlation function, yields the first experimental determination of the individual Zn−N and Zn−C bond lengths as 1.969(2) and 2.030(2) Å, respectively. The very small difference in bond lengths, of ~0.06 Å, means that it is impossible to obtain these bond lengths using Bragg diffraction in isolation. Total neutron diffraction also provides information on both the average and local atomic displacements responsible for NTE in Zn(CN)2. The principal motions giving rise to NTE are shown to be those in which the carbon and nitrogen atoms within individual Zn−C≡N−Zn linkages are displaced to the same side of the Zn···Zn axis. Displacements of the carbon and nitrogen atoms to opposite sides of the Zn···Zn axis, suggested previously in X-ray studies as being responsible for NTE behavior, in fact make negligible contribution at temperatures up to 295 K.

Adsorption, orientation and thermal decomposition of copper(II) hexafluoroacetylacetonate on rutile TiO2(110)

We have investigated the adsorption and thermal decomposition of copper hexafluoroacetylacetonate (Cu-11(hfaC)(2)) on single crystal rutile TiO2(110). Low energy electron diffraction shows that room temperature saturation coverage of the Cu-II(hfac)(2) adsorbate forms an ordered (2 x 1) over-layer. X-ray and ultra-violet photoemission spectroscopy of the saturated surface were recorded as the sample was annealed in a sequential manner to reveal decomposition pathways. The results show that the molecule dissociatively adsorbs by detachment of one of the two ligands to form hfac and Cu-1(hfac) which chemisorb to the substrate at 298 K. These ligands only begin to decompose once the surface temperature exceeds 473 K where Cu core level shifts indicate metallisation. This reduction from Cu(I) to Cu(0) takes place in the absence of an external reducing agent and without disproportionation and is accompanied by the onset of decomposition of the hfac ligands. Finally, C K-edge near edge X-ray absorption fine structure experiments indicate that both the ligands adsorb aligned in the < 001 > direction and we propose a model in which the hfac ligands adsorb on the 5-fold coordinated Ti atoms and the Cu-1(hfac) moiety attaches to the bridging O atoms in a square planar geometry. The calculated tilt angle for these combined geometries is approximately 10 degrees to the surface normal.

Thermal expansion of deuterated hopeite, Zn-3(PO4)(2)center dot 4D(2)O

The lattice parameters extracted from Lebail analysis of neutron powder diffraction data collected between 2 and 300 K have been used to calculate the temperature evolution of the thermal expansion tensor for hopeite, Zn-3(PO4)(2)center dot 2H(2)O, Pnma,Z=4with a= 10.6065(4) angstrom, b = 18.2977(4) angstrom, c= 5.0257(2) A at 275 K. The a lattice parameter shows a negative thermal expansion, the b lattice parameter appears to saturate at 275 K while the c lattice parameter has a more typical positive thermal expansion. At 275 K, the magnitudes of the thermal expansion coefficients are alpha(a) = -1. 1(4) x 10(-5) K-1, alpha(b) = 2.4(9) x 10(-6) K-1 and alpha(c) = 3.6(2) x 10(-1) K-1. Under the conditions of these experiments, hopeite begins to dehydrate to the dihydrate between 300 and 325 K, and between 480 and 500 K the monohydrate is formed. The thermal expansion of the dihydrate has been calculated between 335 and 480 and at 480 K the magnitudes of the thermal expansion coefficients are alpha(a) = 1(2) x 10(-5) K-1, alpha(b) = 4(l) x 10(-6) K-1, alpha(c) = 4(2) x 10(-5) K-1, alpha(beta) = 1 (1) x 10(-1) K-1, and alpha(v) = 2(2) x 10(-1) K-1. The thermal expansion of hopeite is described in terms of its crystal structure and possible dehydration mechanisms for the alpha and beta modifications of hopeite are discussed.