Design and analysis of climate model experiments for the efficient estimation of anthropogenic signals

Presented herein is an experimental design that allows the effects of several radiative forcing factors on climate to be estimated as precisely as possible from a limited suite of atmosphere-only general circulation model (GCM) integrations. The forcings include the combined effect of observed changes in sea surface temperatures, sea ice extent, stratospheric (volcanic) aerosols, and solar output, plus the individual effects of several anthropogenic forcings. A single linear statistical model is used to estimate the forcing effects, each of which is represented by its global mean radiative forcing. The strong colinearity in time between the various anthropogenic forcings provides a technical problem that is overcome through the design of the experiment. This design uses every combination of anthropogenic forcing rather than having a few highly replicated ensembles, which is more commonly used in climate studies. Not only is this design highly efficient for a given number of integrations, but it also allows the estimation of (nonadditive) interactions between pairs of anthropogenic forcings. The simulated land surface air temperature changes since 1871 have been analyzed. The changes in natural and oceanic forcing, which itself contains some forcing from anthropogenic and natural influences, have the most influence. For the global mean, increasing greenhouse gases and the indirect aerosol effect had the largest anthropogenic effects. It was also found that an interaction between these two anthropogenic effects in the atmosphere-only GCM exists. This interaction is similar in magnitude to the individual effects of changing tropospheric and stratospheric ozone concentrations or to the direct (sulfate) aerosol effect. Various diagnostics are used to evaluate the fit of the statistical model. For the global mean, this shows that the land temperature response is proportional to the global mean radiative forcing, reinforcing the use of radiative forcing as a measure of climate change. The diagnostic tests also show that the linear model was suitable for analyses of land surface air temperature at each GCM grid point. Therefore, the linear model provides precise estimates of the space time signals for all forcing factors under consideration. For simulated 50-hPa temperatures, results show that tropospheric ozone increases have contributed to stratospheric cooling over the twentieth century almost as much as changes in well-mixed greenhouse gases.

EOQ, JIT and fixed costs in the Ready-mixed concrete industry

The successful implementation of just-in-time (JIT) purchasing policy in many industries has prompted many companies that still use the economic order quantity (EOQ) purchasing policy to ponder if they should switch to the JIT purchasing policy. Despite existing studies that directly compare the costs between the EOQ and JIT purchasing systems, this decision is, however, still difficult to be made, especially when price discount has to be considered. JIT purchasing may not always be successful even though plants that adopted JIT operations have experienced or can take advantage of physical space reduction. Hence, the objective of this study is to expand on a classical EOQ with a price discount model to derive the EOQ–JIT cost indifference point. The objective was tested and achieved through a survey and case study conducted in the ready-mixed concrete industry in Singapore.

In vitro fermentation of chitosan derivatives by mixed cultures of human faecal bacteria

Stirred, pH controlled batch cultures were carried out with faecal inocula and various chitosans to investigate the fermentation of chitosan derivatives by the human gut flora. Changes in bacterial levels and short chain fatty acids were measured over time. Low, medium and high molecular weight chitosan caused a decrease in bacteroides, bifidobacteria, clostridia and lactobacilli. A similar pattern was seen with chitosan oligosaccharide (COS). Butyrate levels also decreased. A three-stage fermentation model of the human colon was used for investigation of the metabolism of COS. In a region representing the proximal colon, clostridia decreased while lactobacilli increased. In the region representing the transverse colon, bacteroides and clostridia increased. Distally a small increase in bacteroides occurred. Butyrate levels increased. Under the highly competitive conditions of the human colon, many members of the microflora, are unable to compete for chitosans of low, medium or high molecular weight. COS were more easily utilised and when added to an in vitro colonic model led to increased production of butyrate, but some populations of potentially detrimental bacteria also increased. (c) 2005 Elsevier Ltd. All rights reserved.

A two-stage continuous culture system to study the effect of supplemental alpha-lactalbumin and glycomacropeptide on mixed cultures of human gut bacteria challenged with enteropathogenic Escherichia coli and Salmonella serotype Typhimurium

Aims: Certain milk factors may promote the growth of a gastrointestinal microflora predominated by bifidobacteria and may aid in overcoming enteric infections. This may explain why breast-fed infants experience fewer intestinal infections than their formula-fed counterparts. The effect of formula supplementation with two such factors was investigated in this study. Methods and Results: Infant faecal specimens were used to ferment formulae supplemented with glycomacropeptide (GMP) and alpha-lactalbumin (alpha-la) in a two-stage compound continuous culture model. At steady state, all fermenter vessels were inoculated with 5 ml of 0.1 M phosphate-buffered saline (pH 7.2) containing 10(8) CFU ml(-1) of either enteropathogenic Escherichia coli 2348/69 (O127:H6) or Salmonella serotype Typhimurium (DSMZ 5569). Bacteriology was determined by independent fluorescence in situ hybridization. Vessels that contained breast milk (BM), as well as alpha-la and GMP supplemented formula had stable total counts of bifidobacteria while lactobacilli increased significantly only in vessels with breast milk. Bacteroides, clostridia and E. coli decreased significantly in all three groups prior to pathogen addition. Escherichia coli counts decreased in vessels containing BM and alpha-la while Salmonella decreased significantly in all vessels containing BM, alpha-la and GMP. Acetate was the predominant acid. Significance and Impact of the Study: Supplementation of infant formulae with appropriate milk proteins may be useful in mimicking the beneficial bacteriological effects of breast milk.

Crystal structures of [PhLi center dot(-)-sparteine](2), [PhOLi center dot(-)-sparteine](2), and the mixed aggregate [PhLi center dot PhOLi center dot 2(-)-sparteine]

Crystal structure determination of adducts of sparteine and PhLi, (-)-sparteine and PhOLi and of sparteine and PhLi/PhOLi reveal a four-membered ring with two lithium centers, each capped by a (-)-sparteine ligand, as central motif of all structure. Quantum-chemical calculations show that the mixed aggregate [PhLi center dot PhOLi center dot 2(-)-sparteine] is energetically more favorable than the model system {1/2[PhLi center dot(-)-sparteine](2) + 1/2[PhOLi center dot(-)-sparteine](2)}.

Model selection approaches for non-linear system identification: a review

The identification of non-linear systems using only observed finite datasets has become a mature research area over the last two decades. A class of linear-in-the-parameter models with universal approximation capabilities have been intensively studied and widely used due to the availability of many linear-learning algorithms and their inherent convergence conditions. This article presents a systematic overview of basic research on model selection approaches for linear-in-the-parameter models. One of the fundamental problems in non-linear system identification is to find the minimal model with the best model generalisation performance from observational data only. The important concepts in achieving good model generalisation used in various non-linear system-identification algorithms are first reviewed, including Bayesian parameter regularisation and models selective criteria based on the cross validation and experimental design. A significant advance in machine learning has been the development of the support vector machine as a means for identifying kernel models based on the structural risk minimisation principle. The developments on the convex optimisation-based model construction algorithms including the support vector regression algorithms are outlined. Input selection algorithms and on-line system identification algorithms are also included in this review. Finally, some industrial applications of non-linear models are discussed.

A refined LEED analysis of water on Ru{0001}: an experimental test of the partial dissociation model

Despite a number of earlier studies which seemed to confirm molecular adsorption of water on close-packed surfaces of late transition metals, new controversy has arisen over a recent theoretical work by Feibelman, according to which partial dissociation occurs on the Ru{0001} surface leading to a mixed (H2O + OH + H) superstructure. Here, we present a refined LEED-IV analysis of the (root3 x root3)R30degrees-D2O-Ru{0001} structure, testing explicitly this new model by Feibelman. Our results favour the model proposed earlier by Held and Menzel assuming intact water molecules with almost coplanar oxygen atoms and out-of-plane hydrogen atoms atop the slightly higher oxygen atoms. The partially dissociated model with an almost identical arrangement of oxygen atoms can, however, not unambiguously be excluded, especially when the single hydrogen atoms are not present in the surface unit cell. In contrast to the earlier LEED-IV analysis, we can, however, clearly exclude a buckled geometry of oxygen atoms.

Use of a reduced Schiff-Base ligand to prepare novel chloro-bridged chains of rare Cu(II) trinuclear complexes with mixed azido/oxo and chloro/oxo bridges

Two mixed bridged one-dimensional (1D) polynuclear complexes, [Cu3L2(mu(1,1)-N-3)(2)(mu-Cl)Cl](n) (1) and {[Cu3L2(mu-Cl)(3)Cl]center dot 0.46CH(3)OH}(n), (2), have been synthesized using the tridentate reduced Schiff-base ligand HL (2-[(2-dimethylamino-ethylamino)-methyl]-phenol). The complexes have been characterized by X-ray structural analyses and variable-temperature magnetic susceptibility measurements. In both complexes the basic trinuclear angular units are joined together by weak chloro bridges to form a 1D chain. The trinuclear structure of 1 is composed of two terminal square planar [Cu(L)(mu(1,1)-N-3)] units connected by a central Cu(II) atom through bridging nitrogen atoms of end-on azido ligands and the phenoxo oxygen atom of the tridentate ligand. These four coordinating atoms along with a chloride ion form a distorted trigonal bipyramidal geometry around the central Cu(II). The structure of 2 is similar; the only difference being a Cl bridge replacing the mu(1,1)-N-3 bridge in the trinuclear unit. The magnetic properties of both trinuclear complexes can be very well reproduced with a simple linear symmetrical trimer model (H = JS(i)S(i+1)) with only one intracluster exchange coupling (J) including a weak intertrimer interaction (.j) reproduced with the molecular field approximation. This model provides very satisfactory fits for both complexes in the whole temperature range with the following parameters: g = 2.136(3), J = 93.9(3) cm(-1) and zj= -0.90(3) cm(-1) (z = 2) for 1 and g = 2.073(7), J = -44.9(4) cm(-1) and zJ = -1.26(6) cm(-1) (z = 2) for 2.

Climate model calculations of the impact of aerosols from road transport and shipping

Road transport and shipping are copious sources of aerosols, which exert a 9 significant radiative forcing, compared to, for example, the CO2 emitted by these sectors. An 10 advanced atmospheric general circulation model, coupled to a mixed-layer ocean, is used to 11 calculate the climate response to the direct radiative forcing from such aerosols. The cases 12 considered include imposed distributions of black carbon and sulphate aerosols from road 13 transport, and sulphate aerosols from shipping; these are compared to the climate response 14 due to CO2 increases. The difficulties in calculating the climate response due to small 15 forcings are discussed, as the actual forcings have to be scaled by large amounts to enable a 16 climate response to be easily detected. Despite the much greater geographical inhomogeneity 17 in the sulphate forcing, the patterns of zonal and annual-mean surface temperature response 18 (although opposite in sign) closely resembles that resulting from homogeneous changes in 19 CO2. The surface temperature response to black carbon aerosols from road transport is shown 20 to be notably non-linear in scaling applied, probably due to the semi-direct response of clouds 21 to these aerosols. For the aerosol forcings considered here, the most widespread method of 22 calculating radiative forcing significantly overestimates their effect, relative to CO2, 23 compared to surface temperature changes calculated using the climate model.

Wind-driven mixing below the oceanic mixed layer

This study describes the turbulent processes in the upper ocean boundary layer forced by a constant surface stress in the absence of the Coriolis force using large-eddy simulation. The boundary layer that develops has a two-layer structure, a well-mixed layer above a stratified shear layer. The depth of the mixed layer is approximately constant, whereas the depth of the shear layer increases with time. The turbulent momentum flux varies approximately linearly from the surface to the base of the shear layer. There is a maximum in the production of turbulence through shear at the base of the mixed layer. The magnitude of the shear production increases with time. The increase is mainly a result of the increase in the turbulent momentum flux at the base of the mixed layer due to the increase in the depth of the boundary layer. The length scale for the shear turbulence is the boundary layer depth. A simple scaling is proposed for the magnitude of the shear production that depends on the surface forcing and the average mixed layer current. The scaling can be interpreted in terms of the divergence of a mean kinetic energy flux. A simple bulk model of the boundary layer is developed to obtain equations describing the variation of the mixed layer and boundary layer depths with time. The model shows that the rate at which the boundary layer deepens does not depend on the stratification of the thermocline. The bulk model shows that the variation in the mixed layer depth is small as long as the surface buoyancy flux is small.

Resolution of sharp fronts in the presence of model error in variational data assimilation

We show that the four-dimensional variational data assimilation method (4DVar) can be interpreted as a form of Tikhonov regularization, a very familiar method for solving ill-posed inverse problems. It is known from image restoration problems that L1-norm penalty regularization recovers sharp edges in the image more accurately than Tikhonov, or L2-norm, penalty regularization. We apply this idea from stationary inverse problems to 4DVar, a dynamical inverse problem, and give examples for an L1-norm penalty approach and a mixed total variation (TV) L1–L2-norm penalty approach. For problems with model error where sharp fronts are present and the background and observation error covariances are known, the mixed TV L1–L2-norm penalty performs better than either the L1-norm method or the strong constraint 4DVar (L2-norm)method. A strength of the mixed TV L1–L2-norm regularization is that in the case where a simplified form of the background error covariance matrix is used it produces a much more accurate analysis than 4DVar. The method thus has the potential in numerical weather prediction to overcome operational problems with poorly tuned background error covariance matrices.

The upper end of climate model temperature projections is inconsistent with past warming

Climate models predict a large range of possible future temperatures for a particular scenario of future emissions of greenhouse gases and other anthropogenic forcings of climate. Given that further warming in coming decades could threaten increasing risks of climatic disruption, it is important to determine whether model projections are consistent with temperature changes already observed. This can be achieved by quantifying the extent to which increases in well mixed greenhouse gases and changes in other anthropogenic and natural forcings have already altered temperature patterns around the globe. Here, for the first time, we combine multiple climate models into a single synthesized estimate of future warming rates consistent with past temperature changes. We show that the observed evolution of near-surface temperatures appears to indicate lower ranges (5–95%) for warming (0.35–0.82 K and 0.45–0.93 K by the 2020s (2020–9) relative to 1986–2005 under the RCP4.5 and 8.5 scenarios respectively) than the equivalent ranges projected by the CMIP5 climate models (0.48–1.00 K and 0.51–1.16 K respectively). Our results indicate that for each RCP the upper end of the range of CMIP5 climate model projections is inconsistent with past warming.

On the attribution of stratospheric ozone and temperature changes to changes in ozone-depleting substances and well-mixed greenhouse gases

The vertical profile of global-mean stratospheric temperature changes has traditionally represented an important diagnostic for the attribution of the cooling effects of stratospheric ozone depletion and CO2 increases. However, CO2-induced cooling alters ozone abundance by perturbing ozone chemistry, thereby coupling the stratospheric ozone and temperature responses to changes in CO2 and ozone-depleting substances (ODSs). Here we untangle the ozone-temperature coupling and show that the attribution of global-mean stratospheric temperature changes to CO2 and ODS changes (which are the true anthropogenic forcing agents) can be quite different from the traditional attribution to CO2 and ozone changes. The significance of these effects is quantified empirically using simulations from a three-dimensional chemistry-climate model. The results confirm the essential validity of the traditional approach in attributing changes during the past period of rapid ODS increases, although we find that about 10% of the upper stratospheric ozone decrease from ODS increases over the period 1975–1995 was offset by the increase in CO2, and the CO2-induced cooling in the upper stratosphere has been somewhat overestimated. When considering ozone recovery, however, the ozone-temperature coupling is a first-order effect; fully 2/5 of the upper stratospheric ozone increase projected to occur from 2010–2040 is attributable to CO2 increases. Thus, it has now become necessary to base attribution of global-mean stratospheric temperature changes on CO2 and ODS changes rather than on CO2 and ozone changes.

Heat and moisture budgets from airborne measurements and high-resolution model simulations

High-resolution simulations with a mesoscale model are performed to estimate heat and moisture budgets of a well-mixed boundary layer. The model budgets are validated against energy budgets obtained from airborne measurements over heterogeneous terrain in Western Germany. Time rate of change, vertical divergence, and horizontal advection for an atmospheric column of air are estimated. Results show that the time trend of specific humidity exhibits some deficiencies, while the potential temperature trend is matched accurately. Furthermore, the simulated turbulent surface fluxes of sensible and latent heat are comparable to the measured fluxes, leading to similar values of the vertical divergence. The analysis of different horizontal model resolutions exhibits improved surface fluxes with increased resolution, a fact attributed to a reduced aggregation effect. Scale-interaction effects could be identified: while time trends and advection are strongly influenced by mesoscale forcing, the turbulent surface fluxes are mainly controlled by microscale processes.

Climate determinism revisited: multiple equilibria in a complex climate model

Multiple equilibria in a coupled ocean–atmosphere–sea ice general circulation model (GCM) of an aquaplanet with many degrees of freedom are studied. Three different stable states are found for exactly the same set of parameters and external forcings: a cold state in which a polar sea ice cap extends into the midlatitudes; a warm state, which is ice free; and a completely sea ice–covered “snowball” state. Although low-order energy balance models of the climate are known to exhibit intransitivity (i.e., more than one climate state for a given set of governing equations), the results reported here are the first to demonstrate that this is a property of a complex coupled climate model with a consistent set of equations representing the 3D dynamics of the ocean and atmosphere. The coupled model notably includes atmospheric synoptic systems, large-scale circulation of the ocean, a fully active hydrological cycle, sea ice, and a seasonal cycle. There are no flux adjustments, with the system being solely forced by incoming solar radiation at the top of the atmosphere. It is demonstrated that the multiple equilibria owe their existence to the presence of meridional structure in ocean heat transport: namely, a large heat transport out of the tropics and a relatively weak high-latitude transport. The associated large midlatitude convergence of ocean heat transport leads to a preferred latitude at which the sea ice edge can rest. The mechanism operates in two very different ocean circulation regimes, suggesting that the stabilization of the large ice cap could be a robust feature of the climate system. Finally, the role of ocean heat convergence in permitting multiple equilibria is further explored in simpler models: an atmospheric GCM coupled to a slab mixed layer ocean and an energy balance model